## ANALISIS RESPON HIDROLOGI TERHADAP PENERAPAN TEKNIK KONSERVASI TANAH DI SUB DAS LENGKONG MENGGUNAKAN MODEL SWAT

## Hydrological Response Analysis for the Application of Soil Conservation Techniques in Sub-watershed Lengkong Using SWAT Model
Gunadi Firdaus 1)* , Oteng Haridjaja 2) , dan Suria Darma Tarigan 2)
1) Alumni Program Studi Ilmu Pengelolaan Daerah Aliran Sungai, Sekolah Pascasarjana IPB, Jl. Raya Darmaga,
Gedung Andi Hakim Nasoetion Kampus IPB Darmaga Bogor 16680
2) Departemen Ilmu Tanah dan Sumberdaya Lahan Fakultas Pertanian IPB, Jln. Meranti Kampus IPB Dramaga Bogor 16680

## ABSTRACT
One  cause  of  flooding  and  erosion  is  the  biophysical  conditions  of  the  upstream  watershed  that  can  not  support  the hydrological function of the watershed.  It is necessary to implement appropriate  measures of soil conservation techniques to improve the hydrological functions of the watershed. This  research aims to:  (1) analyze the hydrological response based on biophysical  condition  of  the  watershed  on  a  meso-scale  level,  and  (2)  analyze  the  hydrological  response  based  on  soil conservation  practices  scenario  application.  The  study  area  located  in  Lengkong  watershed  in  the  upper  Cisadane  basin covering  an  area  of  1,788  ha.  This  research  was  carried  out  using  Soil  and  Water  Assessment  Tool  (SWAT)  hydrological modeling. Hydrological response analysis resulted from soil conservation scenario application showed that the ratio of Qmax and Qmin value of 149.71 for scenario 1 (strip cropping) and 2 (contour cropping), 149.80 for scenario 3 (terracing), 150.25 for scenario 4 (combine scenario 1 and 2), and 149.31 for scenario 5 (combine scenario 1, 2 and 3) respectively. Based on the potential erosion ratio with tolerable erosion, an erosion hazard index of 2.63 (moderate) for scenario 1; 2.57  (moderate) for scenario  2;  0.60  (low)  for  scenario  3;  2.45  (moderate)  for  scenario  4;  and  0.44  (low)  for  scenario  5.  Application  of  soil conservation  techniques  that  combine  the  contour  cropping  and  strip  cropping  with  terracing  collectively,  was  the  best  soil conservation techniques to improve  hydrological response of the watershed. This techniques are recommended to be applied especially in Lengkong watershed and generally in the upper of Cisadane basin.
Keywords: Flow-out, sediment-yield, soil conservation, SWAT, watershed

## ABSTRAK
Salah  satu  penyebab  banjir  dan  erosi  adalah  kondisi  biofisik  di  bagian  hulu  Daerah  Aliran  Sungai  (DAS)  yang  sudah tidak  dapat  mendukung  fungsi  hidrologis  DAS.  Perlu  dilakukan  upaya-upaya  penerapan  teknik  konservasi  tanah  yang  tepat untuk  memperbaiki  fungsi  hidrologis  DAS  tersebut.  Penelitian  ini  bertujuan  untuk:  (1)  menganalisis  respon  hidrologi berdasarkan  kondisi  biofisik  DAS  pada  tingkat  analisis  skala  meso,  dan  (2)  menganalisis  respon  hidrologi  berdasarkan penerapan  skenario  teknik  konservasi  tanah.  Wilayah  kajian  untuk  penelitian  adalah  di  sub  DAS  Lengkong  yang  terletak  di bagian  hulu  DAS  Cisadane  seluas  1,788  ha.  Penelitian  ini  dilakukan  dengan  menggunakan  pemodelan  hidrologi Soil  and Water Assessment Tool (SWAT). Hasil analisis respon hidrologi terhadap penerapan skenario konservasi tanah  menunjukkan Koefisien  Regim  Sungai  sebesar  149.71  untuk  skenario  1  (penanaman  tanaman  strip)  dan  2  (penanaman  pohon  mengikuti kontur), sebesar 149.80 untuk skenario 3 (pembuatan teras), sebesar 150.25 untuk skenario 4 (gabungan skenario 1 dan 2), dan sebesar 149.31 untuk skenario 5 (gabungan skenario 1, 2 dan 3). Berdasarkan nisbah erosi potensial dengan erosi yang dapat ditoleransi, diperoleh indeks bahaya erosi untuk skenario 1 sebesar 2.63 (sedang), skenario 2 sebesar 2.57 (sedang), skenario 3 sebesar  0.60  (rendah),  skenario  4  sebesar  2.45  (sedang),  dan  skenario  5  sebesar  0.44  (rendah).  Penerapan  teknik  konservasi tanah  yang  mengkombinasikan  penanaman  pohon  mengikuti  kontur  dan  penanaman  tanaman  strip  dengan  pembuatan  teras secara bersamaan, merupakan teknik konservasi tanah yang terbaik untuk memperbaiki respon hidrologi DAS, sehingga dapat direkomendasikan untuk diterapkan di sub DAS Lengkong pada khususnya, dan di DAS Cisadane bagian hulu pada umumnya.
Kata kunci: Debit, sedimen,  konservasi tanah, SWAT, DAS

## PENDAHULUAN
DAS  Cisadane  termasuk  dalam  DAS  Prioritas Nasional (Kemenhut, 2009). Hasil review lahan kritis yang dilakukan  oleh  Kementerian  Kehutanan  melalui  Balai Pengelolaan  Daerah  Aliran  Sungai  (BPDAS)  Citarum-
Ciliwung pada tahun 2009 menunjukkan bahwa luas lahan kritis  di  DAS  Cisadane  adalah  55,685  ha  atau  mencapai 36%  dari  luas  total  DAS  Cisadane  154,654  ha.  Sampai dengan akhir tahun 2012, penanganan lahan kritis di DAS Cisadane  tidak  banyak  dilakukan,  sehingga  luasnya  tidak berkurang,  malah  cenderung  bertambah  seiring  dengan
*) Penulis Korespondensi: Telp. +628128356338; Email. gunadifirdaus@yahoo.co.id
berlangsungnya  gangguan  dan  kerusakan  yang  dilakukan oleh  berbagai  aktivitas  manusia  terutama  di  bagian  hulu DAS. Tercatat jumlah erosi di Cisadane bagian hulu seluas 42,837  ha  pada  tahun  2011  adalah  sebesar  6.2  juta  ton tahun -1 dengan laju erosi mencapai 144.78 ton ha -1 tahun -1 (BPDAS Citarum-Ciliwung, 2012).
Dalam fokus wilayah yang lebih sempit, salah satu Sub DAS Cisadane yang berada di bagian hulu DAS yaitu Sub DAS Lengkong dengan luas 1,778 ha atau sekitar 4% dari  luas  DAS  Cisadane  bagian  hulu,  mempunyai    laju erosi  sebesar  165.38  ton  ha -1 tahun -1 atau  20.6  ton  ha -1 tahun -1 lebih  tinggi  dari  laju  erosi  DAS  Cisadane  bagian hulu.  Laju  erosi  yang  cukup  tinggi  di  DAS  Cisadane bagian hulu pada umumnya dan Sub DAS Lengkong pada khususnya,  disebabkan  terutama  oleh  sistem  budidaya pertanian  lahan  kering  seperti  ladang  dan  pertanian campuran  yang  tidak  dibarengi  dengan  penerapan  teknik
konservasi  tanah  yang  baik  dan  benar  (BPDAS  Citarum- Ciliwung, 2012). Penelitian ini bertujuan untuk menganalisis respon hidrologi berdasarkan kondisi biofisik DAS  pada  tingkat  analisis  skala  meso,  dan  menganalisis respon  hidrologi  berdasarkan  penerapan  skenario  teknik konservasi tanah.

## BAHAN DAN METODE
Penelitian dilakukan di Sub DAS Lengkong yang terletak  dibagian  hulu  DAS  Cisadane  yang  secara administrasi masuk ke dalam Desa Pasirbuncir Kecamatan Caringin  Kabupaten  Bogor  Propinsi  Jawa  Barat  (Gambar 1) seluas 1,788 ha dengan lokasi outlet yang menjadi fokus penelitian  terletak  pada  koordinat  106°  49’  45.4”  Bujur Timur, dan 6° 45’ 28.3” Lintang Selatan.
Gambar 1. Peta situasi wilayah penelitian

## Bahan
Data yang digunakan dibagi menjadi dua jenis yaitu data  spasial  dan  data  numerik.  Data  spasial  yang digunakan  adalah:  (1)  peta  tutupan  lahan  skala  1  :  25,000 berdasarkan  interpretasi  citra  satelit wordview ,  (2)  peta klasifikasi  tanah  skala  1  :  100,000  dari  Pusat  Penelitian Tanah dan Agroklimat (Puslitanak) Bogor, dan (3) Digital Elevation Model (DEM) dari Citra SRTM (Tabel 1).  Data spasial  digunakan  untuk  keperluan  pembentukan  jaringan sungai,  pembentukan outlet ,  batas  DAS  dan  HRU.  Data numerik yang digunakan yaitu: (1) data curah hujan harian dan suhu maksimum dan minimum harian selama 3 tahun (2009-2011)  didapatkan  dari  BPDAS  Citarum-Ciliwung dan  Balai  Pengelolaan  Sumber  Daya  Air  (BPSDA) Ciliwung-Cisadane, (2) data iklim didapatkan dari Stasiun Badan  Meteorologi  Klimatologi  dan  Geofisika  (BMKG)
Dramaga  yang  merupakan  hasil  pengolahan  data  cuaca bulanan  (curah  hujan,  suhu  udara,  kelembaban  udara, radiasi  surya  dan  kecepatan  angin)  selama  10  tahun  yaitu tahun  2001  sampai  2010,  dan  (3)  data  debit  dan  sedimen rata-rata  harian  selama  3  tahun  (2009-2011)  dari  BPDAS Citarum-Ciliwung,  yang  akan  digunakan  untuk  kalibrasi dan validasi model.
Tabel  1.  Kualifikasi  data  spasial  yang  digunakan  di  Sub  DAS Lengkong
Data Spasial Skala Sumber Tahun DEM Tutupan lahan Klasifikasi tanah 1:25,000 1:25,000 1:100,000 Citra SRTM 30 m Interpretasi citra wordview Puslitanak, Bogor 2013 2011 1992

## Alat
Peralatan  yang  digunakan  dalam  penelitian  ini adalah  seperangkat  komputer  dengan  perangkat  lunak ArcGIS  9.3,  ArcSWAT  versi  2009.93.5  dan  Microsoft Office  2010.  Untuk  perhitungan  koefisien  deterministik (R 2 ) dan efisiensi Nash-Sutcliffe (NS) digunakan perangkat lunak  SWATPlot.  Untuk  perhitungan  R 2 dengan  grafik sebaran  seri  data  secara  linier  digunakan  program MS.excel.  Untuk  proses  kalibrasi  dan  validasi  model digunakan perangkat lunak SWAT-CUP dengan algoritma Sequential Uncertainty Fitting Version 2 (SUFI2).

## Analisis Data

## Penetapan HRU dan Menjalankan Model SWAT
Penetapan Hidrology  Respons  Unit (HRU)  diawali oleh deliniasi batas Sub DAS Lengkong berdasarkan input peta  DEM,  dengan  titik outlet pada  koordinat  106 0 82’ 87.3 ’’ bujur  timur  dan  6 0 75’  78 ’’ lintang  selatan  dengan menggunakan model SWAT. HRU dibentuk menggunakan model SWAT dengan cara menggabungkan ( overlay ) peta tanah, peta liputan lahan, dan peta lereng. HRU yang telah
terbentuk  digabungkan  dengan  data-data  iklim  meliputi data  iklim  global,  data  curah  hujan  harian  rata-rata  serta data  suhu  maksimum  dan  minimum  harian  rata-rata, kemudian  SWAT  dijalankan.  Model  SWAT  yang  telah dijalankan  akan  menghasilkan output model  diantaranya debit  harian  rata-rata  (FLOW_OUT)  dan  hasil  sedimen (SYLD).

## Uji Akurasi Hasil Model, Kalibrasi, dan Validasi
Output model  diuji  keakuratannya  dengan menggunakan  persamaan  efisiensi  model Nash-Sutcliffe (NS):
∑ (     ) ∑ (     )
dimana  y  adalah variable rata-rata  yang  terukur  (mm),  ŷ adalah variable hasil simulasi (mm), dan y adalah variable terukur.  Efisiensi  model  NS  dikelompokkan  menjadi  tiga kelas  (Nash et  al., 1970)  yaitu  baik  (NS  ≥  0.75), memuaskan  (0.75  >  NS  >  0.36)  dan  kurang  memuaskan (NS < 0.36).
Tabel 2. Parameter input model untuk kalibrasi
No. Parameter Definisi Unit 1 CN2 Curve Number in moisture condition II - 2 ALPHA_BF Baseflow alpha factor. hari 3 GW_DELAY Groundwater delay. hari 4 GWQMN Threshold depth of water in the shallow aquifer required for return flow to occur. mm 5 GW_REVAP Groundwater "revap" coefficient. - 6 ESCO Soil evaporation compensation factor. - 7 CH_N2 Manning's "n" value for the main channel. - 8 CH_K2 Effective hydraulic conductivity in main channel alluvium. mm jam -1 9 ALPHA_BNK Baseflow alpha factor for bank storage. - 10 SOL_AWC Available water capacity of the soil layer. mm mm -1 11 SOL_K Saturated hydraulic conductivity. mm jam -1 12 SOL_BD Moist bulk density. - 13 RCHRG_DP Deep aquifer percolation fraction. - 14 SURLAG Surface runoff lag time. hari 15 USLE_P USLE support practice faktor - 16 REVAPMN Threshold depth of water in the shallow aquifer for "revap" to occur. mm 17 GWHT Initial groundwater height. m 18 SHALLST Initial depth of water in the shallow aquifer. mm 19 DEEPST Initial depth of water in the deep aquifer . mm 20 GW_SPYLD Specific yield of the shallow aquifer. m 3 m -3 21 SOL_ZMX Maximum rooting depth of soil profile. mm 22 SOL_CRK Crack volume potential of soil. m 3 m -3 23 SLSUBBSN Average slope length. mm 24 HRU_SLP Average slope steepness m m -1 25 OV_N Manning's "n" value for overland flow. - 26 LAT_TTIME Lateral flow travel time. hari 27 LAT_SED Sediment concentration in lateral flow and groundwater flow . mg l -1 28 SLSOIL Slope length for lateral subsurface flow. m
Selain  itu,  untuk  melihat  keakuratan  pola  hasil keluaran  model  dengan  hasil  observasi  lapangan digunakan koefisien deterministik atau persamaan linier:
(     ̅) (     )
(     ̅)
dimana X adalah  besarnya variable pengamatan, X adalah variable rata-rata  pengamatan  dan Y adalah variable perhitungan  model.  Hasil  perhitungan  R 2 menunjukkan evaluasi kelayakan model tersebut, apabila R 2 mendekati 1 maka  terdapat  hubungan  yang  erat  antara  hasil  prediksi model  dengan  hasil  observasi  lapangan.  Apabila  hasil
pengujian  R 2 dan  NS  <  0.75  maka  dilakuan  kalibrasi menggunakan  SWAT-CUP  terhadap  28  parameter  (Tabel 2)  yang  dipilih  berdasarkan  hasil  penelitian  Mulyana (2012), Latifah (2013) dan Yustika (2013 ).
Nilai  terbaik  untuk  masing-masing  parameter digunakan  untuk  menggantikan  nilai  parameter  model sebelumnya,  kemudian  dilakukan  validasi  dengan  cara menjalankan  model  yang  telah  terkalibrasi  tersebut terhadap  tahun  yang  berbeda.  Apabila  uji  akurasi  model hasil  validasi  mencapai  nilai  NS  dan  R 2 ≥  0.75  maka model sudah dapat digunakan untuk menjalankan skenario penerapan konservasi tanah.

## Skenario Penerapan Konservasi Tanah
Berdasarkan
pada
Rencana
Pengelolaan
Rehabilitasi  Hutan  dan  Lahan  (RP-RHL)  Kabupaten Bogor  Tahun  2011–2015,  teknik  konservasi  tanah  yang dijadikan  skenario  dalam  penelitian  ini  adalah:  (1) strip cropping atau  penanaman  tanaman  strip,  (2) contour cropping atau  penerapan  teknik  penanaman  pohon mengikuti  kontur,  (3) terracing atau  pembuatan  teras,  (4) kombinasi skenario 1 dan 2, dan (5) kombinasi skenario 1, 2  dan  3.    Semua  skenario  teknik  konservasi  tanah  itu dilakukan  di  semua  HRU  pada  kelerengan  <  45%  dengan pengecualian  tutupan  lahan  hutan  dan  pemukiman.  Dari lima  skenario  tersebut,  dapat  dilihat  masing-masing hasilnya  dan  perbandingan  respon  hidrologinya  berupa Koefisien  Regim  Sungai  (KRS),  dan  Indek  Bahaya  Erosi (IBE).
Untuk perhitungan KRS digunakan persamaan:
dimana KRS adalah  koefisien  regim  sungai, Qmak adalah debit  rata-rata  harian  maksimum  (m 3 detik -1 )  dan Qmin adalah  debit  rata-rata  harian  minimum  (m 3 detik -1 ).
Semakin  kecil  KRS  semakin  menunjukkan  kondisi  DAS dalam keadaan baik (Asdak, 2007).
Untuk  perhitungan  indeks  bahaya  erosi  digunakan persamaan:
dimana IBE adalah  indeks  bahaya  erosi, EP adalah  erosi potensial (ton ha -1 tahun -1 ) dan TSL atau Etol adalah erosi yang  dapat  ditoleransi  (ton  ha -1 tahun -1 ).  Adapun  untuk memperoleh  nilai  erosi  yang  masih  dapat  ditoleransi digunakan persamaan:
dimana DE adalah  kedalaman  solum  tanah  (mm), Fd adalah  nilai  faktor  kedalaman  tanah,  dan T adalah  umur guna sumber daya tanah (tahun). Untuk mendapatkan nilai Etol dalam  satuan  ton  ha -1 tahun -1 maka  nilai Etol dikalikan  dengan  bobot  isi  tanah  (g  cm -3 )  dikalikan  10 (Arsyad, 2010).

## HASIL DAN PEMBAHASAN

## Identifikasi Kondisi Biofisik

## Jenis Tanah
Hasil  klasifikasi  tanah  (Gambar  2),  berdasarkan peta  tanah  dari  Pusat  Penelitian  Tanah  dan  Agroklimat (Puslitanak)  Bogor  skala  1  :  100,000,  jenis  tanah  di  sub DAS  Lengkong  diklasifikasikan  menjadi  dua  yaitu  Andic Humitropepts  seluas  113  ha  (6%)  yang  berada  di  bagian hilir  sub  DAS  Lengkong  dan  Typic  Hapludands  seluas  1 675  (94%)  berada  di  bagian  hulu,  tengah  dan  sedikit  di bagian  hilir  Sub  DAS  Lengkong. Database peta  tanah spasial diberikan tambahan informasi  kelompok hidrologi, morfologi,  dan  sifat-sifat  tanah,  hasil  uji  laboratorium  2 sample profil tanah oleh Puslitanak (Tabel 3).
Gambar 2. Peta klasifikasi tanah berdasarkan model SWAT

## Tabel 3. Kelompok hidrologi, morfologi dan sifat-sifat tanah
Deskripsi Satuan Andic Humitropepts Typic Hapludands Jumlah Horison  tanah - 3 3 Kelompok hidrologi tanah - B A Kedalaman tanah mm 1200 1400 Lapisan Tanah ke- 1 2 3 1 2 3 Kedalaman lapisan tanah mm 200 400 600 300 500 600 Bobot isi tanah ( Bulk  Density ) g cm -3 0.8 1.0 1.1 0.7 0.9 1.0 Kapasitas air tanah tersedia mm mm -1 0.1 0.1 0.1 0.17 0.17 0.17 Kandungan bahan organik % 4.5 3.2 3.0 2.7 2.1 2.0 Konduktivitas hidrolik mm jam -1 26.2 26.2 26.2 29.7 29.7 29.7 Kandungan liat % 16 16 16 16 14 13 Kandungan debu % 21 21 21 21 18 17 Kandungan pasir % 53 53 53 58 56 55 Kandungan batu % 10 10 10 5 12 15 Albedo tanah Fraksi 0.21 0.22 0.23 0.21 0.22 0.23 Nilai erodibilitas tanah (USLE-K) - 0.32 0.28 0.27 0.28 0.26 0.24

## Tutupan Lahan
Berdasarkan  interpretasi  citra  satelit,  tutupan  lahan di  Sub  DAS  Lengkong  digolongkan  dalam  tujuh  kelas pada  model  SWAT  yaitu:  (1)  ladang,  (2)  petanian  lahan kering campuran, (3) hutan, (4) semak belukar, (5) padang rumput,  (6)  sawah,  dan  (7)  pemukiman  dengan  nilai  CN
tertinggi  92  untuk  pemukiman  dan  C-USLE  tertinggi  0.2 untuk ladang (Tabel 4).
Tutupan  lahan  hutan  berada  di  bagian  paling  hulu dengan luas mencapai 58.7% dari luas Sub DAS Lengkong (Gambar  3).  Di  bagian  hilir  Sub  DAS  Lengkong didominasi oleh tutupan lahan ladang sedangkan di bagian tengah DAS didominasi oleh tutupan lahan semak belukar 468 ha dan pertanian campuran 107 ha.
Tabel 4. Hasil identifikasi tutupan lahan berdasarkan model SWAT
Tutupan lahan Kode SWAT Luas (ha) Luas (%) Nilai CN*) C-USLE Ladang Pertanian campuran Hutan Semak belukar Padang rumput Sawah Pemukiman AGRL AGRC FRSE FRST PAST RICE URML 87 107 1050 469 3 61 11 4.85 6.00 58.70 26.25 0.15 3.42 0.64 77 73 55 60 69 73 92 0.200 0.030 0.001 0.010 0.030 0.030 - Jumlah: - 1788 100.00 - - *) pada kelompok hidrologi tanah B, Sumber: Neitsch et al. (2005)

## Kelas Kelerengan
Berdasarkan  hasil  model  SWAT,  kelas  kelerengan Sub DAS  Lengkong didominasi oleh lereng sangat  curam (>40%) seluas 1,110 ha (63%) khas daerah hulu DAS dan
kelas kelerengan tersempit adalah landai (0-8%) seluas 17 ha (Tabel 5). Lereng yang sangat curam dan curam berada di bagian hulu DAS (Gambar 4).
Tabel 5. Hasil identifikasi kelas kelerengan model SWAT Slope (%) Definisi Luas (ha) Luas (%) 0-8 8-15 15-25 25-40 >40 Landai Bergelombang Berbukit Curam Sangat Curam 17 55 185 422 1,110 0.93 3.06 10.36 23.60 62.50 Jumlah: - 1,788 100.00
Gambar 4. Peta klasifikasi lereng berdasarkan model SWAT

## Uji Akurasi Model, Kalibrasi dan Validasi
Hasil  pengujian  akurasi  debit  model  sebelum kalibrasi terhadap debit hasil observasi lapangan selama 3 tahun  (2009  –  2011)  dengan  menggunakan  SWATPlot
diperoleh R 2 sebesar 0.46 dan NS sebesar 0.41, sedangkan pengujian  terhadap  hasil  sedimen  model  awal  sebelum kalibrasi  diperoleh  R 2 dan  NS  masing-masing  0.42  dan 0.39.

## Tabel 6. Nilai kalibrasi menggunakan model SWAT-CUP
No Parameter *) Nilai input awal Nilai hasil kalibrasi Fitted_value Min_value Max_value 1 r__CN2.mgt 1.00 -0.12 -0.20 0.20 2 v__ALPHA_BF.gw 0.05 0.20 0.10 1.00 3 v__GW_DELAY.gw 31.00 39.35 0.00 45.00 4 v__GWQMN.gw 2,500.00 1,122.50 0.00 5,000.00 5 v__GW_REVAP.gw 0.02 0.01 0.00 0.20 6 v__ESCO.hru 0.95 0.96 0.01 1.00 7 v__CH_N2.rte 0.014 0.21 0.00 0.30 8 v__CH_K2.rte 0.00 109.13 0.00 150.00 9 v__ALPHA_BNK.rte 0.00 0.02 0.00 1.00 10 r__SOL_AWC(1).sol 1.00 0.12 -0.20 0.40 11 r__SOL_K(1).sol 1.00 -0.30 -0.80 0.80 12 r__SOL_BD(1).sol 1.00 0.33 -0.50 0.60 13 v__RCHRG_DP.gw 0.00 0.36 0.00 1.00 14 v__SURLAG.bsn 4.00 19.96 1.00 24.00 15 r__USLE_P.mgt 1.00 -0.03 -0.10 0.10 16 v__REVAPMN.gw 10.00 128.75 0.00 500.00 17 v__GWHT.gw 1.00 15.51 0.00 25.00 18 v__SHALLST.gw 400.00 556.50 0.00 1,000.00 19 v__DEEPST.gw 1,000.00 2,197.50 0.00 3,000.00 20 v__GW_SPYLD.gw 0.03.00 0.12 0.00 0.40 21 r__SOL_ZMX.sol 1.00 -0.15 -0.20 0.40 22 r__SOL_CRK.sol 1.00 0.14 -0.20 0.40 23 r__SLSUBBSN.hru 1.00 -0.09 -0.20 0.40 24 r__HRU_SLP.hru 1.00 0.37 -0.20 0.40 25 r__OV_N.hru 1.00 -0.08 -0.20 0.40 26 r__LAT_TTIME.hru 1.00 -0.07 -0.20 0.40 27 r__LAT_SED.hru 1.00 -0.10 -0.20 0.40 28 r__SLSOIL.hru 1.00 0.22 -0.20 0.40 *) v = nilai awal diganti fitted_value ; r  = 1 + fitted_value dikalikan nilai awal (Abbaspour, 2011)
Kalibrasi  dilakukan  dengan  menggunakan  SWAT- CUP  terhadap  debit  dan  sedimen  model  selama  periode  3 tahun  (2009  –  2011)  sebanyak  1095  data  debit  dan sedimen (n=1095) dengan iterasi yang dilakukan sebanyak 1000  simulasi.  Hasil  kalibrasi  diperoleh  nilai  terbaik ( fitted_value )  untuk  masing-masing  parameter input yang digunakan untuk menggantikan nilai awal (Tabel 6). Hasil uji keakuratan debit hasil validasi  model tahun 2011 pada SWATPlot R 2 sebesar 0.86 dan NS sebesar 0.85. Hasil uji keakuratan  sedimen  hasil  validasi  model  berdasarkan SWATPlot, R 2 sebesar 0.85 dan NS sebesar 0.76. Dengan demikian  model  sudah  dapat  digunakan  untuk  melakukan penerapan skenario teknik konservasi tanah.

## Respon Hidrologi terhadap Kondisi Biofisik
Hasil  analisis  pada  model  SWAT  menunjukkan bahwa  pemukiman  memiliki  tebal  aliran  permukaan sebesar  2,314.32  mm  pada  tahun  2011,  diikuti  dengan tutupan lahan sawah sebesar 1,530.01 mm, ladang sebesar 1,034.10  mm,  pertanian  campuran  sebesar  686.40  mm, semak belukar sebesar 597.92 mm padang rumput sebesar 344.23  mm,  dan  hutan  sebesar  176.35  mm.  Angka  ini menunjukkan bahwa infiltrasi yang terjadi pada pemukiman  dan  sawah  nilainya  kecil  karena  pemukiman didominasi oleh areal  yang  kedap air dan sawah  memiliki mulsa  air  yang  juga  kedap  air,  sehingga  pada  kedua  tipe tutupan  lahan  ini  aliran  permukaan  menjadi  tinggi. Berdasarkan  tebal  curah  hujan  tahun  2011  sebesar  3,447 mm,  maka  koefisien  aliran  untuk  masing-masing penutupan  lahan  berturut-turut  mulai  dari  nilai  yang tertinggi  yaitu  pemukiman  sebesar  0.67;  sawah  sebesar 0.44;  ladang  sebesar  0.30;  pertanian  campuran  sebesar 0.20;  semak  belukar  sebesar  0.17;  padang  rumput  sebesar 0.10 dan hutan sebesar 0.05.
Untuk  nilai  erosi  yang  ditimbulkan  pada  masing- masing  tutupan  lahan,  hasil  analisis  model  SWAT menunjukkan  bahwa  ladang  menghasilkan  erosi  sebesar 263.91 ton ha -1 pada tahun 2011, diikuti dengan pertanian campuran  sebesar  165.23  ton  ha -1 ,  pemukiman  sebesar 89.42 ton ha -1 , sawah sebesar 35.9 ton ha -1 , semak belukar sebesar  20.77  ton  ha -1 ,  hutan  sebesar  12.46  ton  ha -1 dan
padang rumput sebesar 7.08 ton ha -1 . Ladang dan pertanian campuran  mempunyai  nilai  erosi  yang  tinggi  karena adanya  pengolahan  tanah  pada  musim  tanam  yang dilakukan pada musim hujan.

## Respon Hidrologi terhadap Skenario Konservasi Tanah

## Fluktuasi Debit Aliran
Hasil  perhitungan  KRS  untuk  5  skenario  teknik konservasi  tanah  yang  diterapkan  beruturut-turut  dari mulai  yang  nilainya  kecil  adalah  skenario  5  sebesar 149.31, skenario 1 dan skenario 2 sebesar 149.71, skenario 3 sebesar 149.80 dan skenario 4 sebesar 150.25 (Tabel 7).
Tabel 7. KRS pada beberapa skenario teknik konservasi tanah
Skenario Qmak (m 3 detik -1 ) Qmin (m 3 detik -1 ) KRS (Qmak/Qmin) 0. tanpa skenario 15.42 0.10 152.82 1. strip cropping 15.24 0.10 149.71 2. contour cropping 15.24 0.10 149.71 3. terracing 15.25 0.10 149.80 4. (1 + 2) 15.28 0.10 150.25 5. (1 + 2 + 3) 15.23 0.10 149.31
Keterangan : Qmak = debit rata-rata harian maksimum dan Qmin = debit rata-rata harian minimum.

## Indeks Bahaya Erosi
Berdasarkan  perbandingan  antara  erosi  potensial dengan  erosi  yang  dapat  ditoleransi  (TSL),  diperoleh indeks bahaya erosi (Tabel 8) untuk skenario1, 2, 3, 4 dan 5  berturut-turut  adalah  2.63  (sedang),  2.57  (sedang),  0.60 (rendah),  2.45  (sedang)  dan  0.44  (rendah).  Untuk perhitungan indeks bahaya erosi ini, nilai TSL adalah 32.2 ton  ha  tahun -1 berdasarkan  kedalaman  solum  tanah  1,400 mm, nilai faktor kedalaman tanah 0.8 dan umur guna tanah 400 tahun (Arsyad, 2010) serta bobot isi tanah 1.15 g cm -3 .
Tabel 8. Indeks bahaya erosi  pada lima skenario konservasi tanah
Skenario Erosi potensial (ton ha -1 tahun -1 ) Erosi yang dapat ditoleransi (ton ha -1 tahun -1 ) Indeks Bahaya Erosi Kategori 0. tanpa skenario 129.07 32.2 4.01 Tinggi 1. strip cropping 84.78 32.2 2.63 Sedang 2. contour cropping 82.83 32.2 2.57 Sedang 3. terracing 19.32 32.2 0.60 Rendah 4. (1 + 2) 79.01 32.2 2.45 Sedang 5. (1 + 2 + 3) 14.29 32.2 0.44 Rendah

## SIMPULAN
Hasil  analisis  respon  hidrologi  terhadap  kondisi biofisik  DAS  menunjukkan  bahwa  tutupan  lahan pemukiman dan sawah mempunyai tebal aliran permukaan lebih  tinggi  dibandingkan  dengan  tutupan  lahan  lainnya. Hal ini disebabkan pemukiman mempunyai areal kedap air yang  menghambat  terjadinya  proses  infiltrasi.  Adapun pada  sawah,  tebal  aliran  permukaan  yang  tinggi
disebabkan  oleh  adanya  lapisan  tapak  bajak  yang menghambat  proses  infiltrasi.  Erosi  yang  ditimbulkan untuk  masing-masing  tutupan  lahan,  menunjukkan  bahwa ladang  dan  pertanian  campuran  mempunyai  erosi  lebih tinggi  dibandingkan  tutupan  lahan  lainnya.  Hal  ini disebabkan adanya kegiatan pengolahan tanah pada musim tanam yang dilakukan sewaktu musim hujan.
Hasil  analisis  respon  hidrologi  terhadap  penerapan skenario  teknik  konservasi  tanah  menunjukkan  bahwa
skenario teknik konservasi tanah yang mengkombinasikan penanaman  tanaman  strip,  penanaman  pohon  mengikuti kontur  dan  pembuatan  teras  secara  bersamaan  merupakan konservasi  tanah  terbaik  dengan  KRS  dan  indek  bahaya erosi  yang  lebih  rendah  dibandingkan  dengan  skenario lainnya.  Penerapan  teknik  konservasi  tanah  yang mengkombinasikan  penanaman  pohon  mengikuti  kontur, dan  penanaman  tanaman  strip  dengan  pembuatan  teras secara  bersamaan,  merupakan  teknik  konservasi  tanah yang  paling  baik  untuk  memperbaiki  respon  hidrologi DAS  yaitu  kondisi  fluktuasi  debit  yang  lebih  stabil  dan sedimen yang lebih menurun, sehingga dapat direkomendasikan  untuk  dapat  diterapkan  di  Sub  DAS Lengkong  pada  khususnya  dan  di  DAS  Cisadane  bagian hulu pada umumnya.

## DAFTAR PUSTAKA
Abbaspour, K.C. 2011. SWAT-CUP4: SWAT Calibration And  Uncertainty  Programs  –  User  Manual. Eawag: Swiss  Federal  Institute  of  Aquatic  Science  and Technology. Swiss.
Arsyad,  S.  2010. Konservasi  Tanah  dan  Air  Edisi  Kedua . IPB Press. Bogor.
Asdak,  C.  2007. Hidrologi  dan  Pengelolaan  Daerah Aliran  Sungai. Gajah  Mada  University  Press.
Yogyakarta.
[BPDASCTW]  Balai  Pengelolaan  Daerah  Aliran  Sungai Citarum-Ciliwung. 2012.
Rencana Tindak
Pengelolaan DAS Cisadane . Bogor.
[Kemenhut]  Kementerian  Kehutanan  Republik  Indonesia,
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-1014 . Jakarta.
Latifah, I. 2013. Analisis ketersediaan air, sedimentasi, dan karbon  organik  dengan  model  SWAT  di  Hulu  DAS Jeneberang, Sulawesi Selatan [Tesis]. IPB. Bogor.
Mulyana,  N.  2012. Analisis  luas  tutupan  hutan  terhadap ketersediaan green water dan blue water di Sub DAS Gumbasa  dan  Sub  DAS  Cisadane  Hulu  dengan aplikasi model SWAT [Disertasi]. IPB. Bogor.
Nash, J.E. and J.V. Sutcliffe. 1970. River flow forecasting through  conceptual  models  part  I  –  discussion  of principles. Journal of Hydrology , 10: 282-190.
Neitsch, S.L., J.G. Arnold, J.R. Kiniry, J.R. Williams, and K.W.  King.  2005. Soil  and Water  Assessment  Tool: Theoritical Documentation Version 2005 . Texas.
Yustika,  R.D.  2013.  Pengelolaan  lahan  terbaik  hasil simulasi  model  SWAT  untuk  mengurangi  aliran permukaan di Sub DAS Ciliwung Hulu [Tesis]. IPB. Bogor.