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

J. Tanah Lingk., 16 (1) April 2014: 16-23

ISSN 1410-7333

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 Firdaus1)*, Oteng Haridjaja2), dan Suria Darma Tarigan2) 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

1)

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

16Penulis Korespondensi: Telp. +628128356338; Email. [email protected] *)

Analisis Respon Hidrologi Terhadap Penerapan Teknik Konservasi Tanah (Firdaus, G., O. Haridjaja, dan S.D. Tarigan)

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 CitarumCiliwung, 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 DEM Tutupan lahan Klasifikasi tanah

Skala 1:25,000 1:25,000 1:100,000

Sumber Citra SRTM 30 m Interpretasi citra wordview Puslitanak, Bogor

Tahun 2013 2011 1992

17


ISSN 1410-7333

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 (R2) dan efisiensi Nash-Sutcliffe (NS) digunakan perangkat lunak SWATPlot. Untuk perhitungan R2 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).

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):

∑ ∑

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 1060 82’ 87.3’’ bujur timur dan 60 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

( (

) )

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. 1 2 3 4

Parameter CN2 ALPHA_BF GW_DELAY GWQMN

5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28

GW_REVAP ESCO CH_N2 CH_K2 ALPHA_BNK SOL_AWC SOL_K SOL_BD RCHRG_DP SURLAG USLE_P REVAPMN GWHT SHALLST DEEPST GW_SPYLD SOL_ZMX SOL_CRK SLSUBBSN HRU_SLP OV_N LAT_TTIME LAT_SED SLSOIL

Definisi Curve Number in moisture condition II Baseflow alpha factor. Groundwater delay. Threshold depth of water in the shallow aquifer required for return flow to occur. Groundwater "revap" coefficient. Soil evaporation compensation factor. Manning's "n" value for the main channel. Effective hydraulic conductivity in main channel alluvium. Baseflow alpha factor for bank storage. Available water capacity of the soil layer. Saturated hydraulic conductivity. Moist bulk density. Deep aquifer percolation fraction. Surface runoff lag time. USLE support practice faktor Threshold depth of water in the shallow aquifer for "revap" to occur. Initial groundwater height. Initial depth of water in the shallow aquifer. Initial depth of water in the deep aquifer . Specific yield of the shallow aquifer. Maximum rooting depth of soil profile. Crack volume potential of soil. Average slope length. Average slope steepness Manning's "n" value for overland flow. Lateral flow travel time. Sediment concentration in lateral flow and groundwater flow . Slope length for lateral subsurface flow.

Selain itu, untuk melihat keakuratan pola hasil keluaran model dengan hasil observasi lapangan digunakan koefisien deterministik atau persamaan linier:

(

̅) (

18

( ̅)

)

Unit hari hari mm mm jam-1 mm mm-1 mm jam-1 hari mm m mm mm m3 m-3 mm m3 m-3 mm m m-1 hari mg l-1 m

dimana X adalah besarnya variable pengamatan, X adalah variable rata-rata pengamatan dan Y adalah variable perhitungan model. Hasil perhitungan R2 menunjukkan evaluasi kelayakan model tersebut, apabila R2 mendekati 1 maka terdapat hubungan yang erat antara hasil prediksi model dengan hasil observasi lapangan. Apabila hasil


pengujian R2 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 R2 ≥ 0.75 maka model sudah dapat digunakan untuk menjalankan skenario penerapan konservasi tanah.

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:

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 (m3 detik-1) dan Qmin adalah debit rata-rata harian minimum (m3 detik-1).

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

19


ISSN 1410-7333

Tabel 3. Kelompok hidrologi, morfologi dan sifat-sifat tanah Deskripsi Jumlah Horison tanah

Satuan -

Andic Humitropepts 3

Typic Hapludands 3

Kelompok hidrologi tanah

-

B

A

Kedalaman tanah

mm

1200

1400

Lapisan Tanah keKedalaman lapisan tanah

mm

Bobot isi tanah (Bulk Density)

g cm-3 -1

Kapasitas air tanah tersedia

mm mm

Kandungan bahan organik

1 200

2 400

3 600

1 300

2 500

3 600

0.8

1.0

1.1

0.7

0.9

1.0

0.1

0.1

0.1

0.17

0.17

0.17

%

4.5

3.2

3.0

2.7

2.1

2.0

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

Fraksi

0.21

0.22

0.23

0.21

0.22

0.23

-

0.32

0.28

0.27

0.28

0.26

0.24

Konduktivitas hidrolik

Albedo tanah Nilai erodibilitas tanah (USLE-K)

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 Ladang Pertanian campuran Hutan Semak belukar Padang rumput Sawah Pemukiman

Kode SWAT Luas (ha) AGRL 87 AGRC 107 FRSE 1050 FRST 469 PAST 3 RICE 61 URML 11 Jumlah: 1788 *) pada kelompok hidrologi tanah B, Sumber: Neitsch et al. (2005)

Luas (%) 4.85 6.00 58.70 26.25 0.15 3.42 0.64 100.00

Gambar 3. Peta klasifikasi tutupan lahan berdasarkan model SWAT

20

Nilai CN*) 77 73 55 60 69 73 92 -

C-USLE 0.200 0.030 0.001 0.010 0.030 0.030 -


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 (%) 0-8 8-15 15-25 25-40 >40 Jumlah:

Definisi Landai Bergelombang Berbukit Curam Sangat Curam -

Luas (ha) 17 55 185 422 1,110 1,788

Luas (%) 0.93 3.06 10.36 23.60 62.50 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 R2 sebesar 0.46 dan NS sebesar 0.41, sedangkan pengujian terhadap hasil sedimen model awal sebelum kalibrasi diperoleh R2 dan NS masing-masing 0.42 dan 0.39.

Tabel 6. Nilai kalibrasi menggunakan model SWAT-CUP No

Parameter *)

Nilai input Nilai hasil kalibrasi awal Fitted_value Min_value 1 r__CN2.mgt 1.00 -0.12 -0.20 2 v__ALPHA_BF.gw 0.05 0.20 0.10 3 v__GW_DELAY.gw 31.00 39.35 0.00 4 v__GWQMN.gw 2,500.00 1,122.50 0.00 5 v__GW_REVAP.gw 0.02 0.01 0.00 6 v__ESCO.hru 0.95 0.96 0.01 7 v__CH_N2.rte 0.014 0.21 0.00 8 v__CH_K2.rte 0.00 109.13 0.00 9 v__ALPHA_BNK.rte 0.00 0.02 0.00 10 r__SOL_AWC(1).sol 1.00 0.12 -0.20 11 r__SOL_K(1).sol 1.00 -0.30 -0.80 12 r__SOL_BD(1).sol 1.00 0.33 -0.50 13 v__RCHRG_DP.gw 0.00 0.36 0.00 14 v__SURLAG.bsn 4.00 19.96 1.00 15 r__USLE_P.mgt 1.00 -0.03 -0.10 16 v__REVAPMN.gw 10.00 128.75 0.00 17 v__GWHT.gw 1.00 15.51 0.00 18 v__SHALLST.gw 400.00 556.50 0.00 19 v__DEEPST.gw 1,000.00 2,197.50 0.00 20 v__GW_SPYLD.gw 0.03.00 0.12 0.00 21 r__SOL_ZMX.sol 1.00 -0.15 -0.20 22 r__SOL_CRK.sol 1.00 0.14 -0.20 23 r__SLSUBBSN.hru 1.00 -0.09 -0.20 24 r__HRU_SLP.hru 1.00 0.37 -0.20 25 r__OV_N.hru 1.00 -0.08 -0.20 26 r__LAT_TTIME.hru 1.00 -0.07 -0.20 27 r__LAT_SED.hru 1.00 -0.10 -0.20 28 r__SLSOIL.hru 1.00 0.22 -0.20 *) v = nilai awal diganti fitted_value; r = 1 + fitted_value dikalikan nilai awal (Abbaspour, 2011)

Max_value 0.20 1.00 45.00 5,000.00 0.20 1.00 0.30 150.00 1.00 0.40 0.80 0.60 1.00 24.00 0.10 500.00 25.00 1,000.00 3,000.00 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40

21


ISSN 1410-7333

Kalibrasi dilakukan dengan menggunakan SWATCUP 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 R2 sebesar 0.86 dan NS sebesar 0.85. Hasil uji keakuratan sedimen hasil validasi model berdasarkan SWATPlot, R2 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 masingmasing 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 Qmak (m3 detik-1)

Qmin (m3 detik-1)

KRS (Qmak/Qmin)

0. tanpa skenario 1. strip cropping 2. contour cropping 3. terracing 4. (1 + 2)

15.42 15.24

0.10 0.10

152.82 149.71

15.24

0.10

149.71

15.25 15.28

0.10 0.10

149.80 150.25

5. (1 + 2 + 3)

15.23

0.10

149.31

Skenario

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 0. tanpa skenario 1. strip cropping 2. contour cropping 3. terracing 4. (1 + 2) 5. (1 + 2 + 3)

Erosi potensial (ton ha-1 tahun-1) 129.07 84.78 82.83 19.32 79.01 14.29

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 22

Erosi yang dapat ditoleransi (ton ha-1 tahun-1) 32.2 32.2 32.2 32.2 32.2 32.2

Indeks Bahaya Erosi 4.01 2.63 2.57 0.60 2.45 0.44

Kategori Tinggi Sedang Sedang Rendah Sedang Rendah

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.

23

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

Recommend Documents