## PEMODELAN PERUBAHAN PENGGUNAAN LAHAN KABUPATEN HUMBANG HASUNDUTAN MENGGUNAKAN CLUE-S

## Modeling Land-Use Change in Humbang Hasundutan District Using
Clue-S
Toga Pandapotan Sinurat 1)* , Khursatul Munibah 2) , dan Dwi Putro Tejo Baskoro 2)
1) Dinas Kehutanan dan Lingkungan Hidup, Jl. Raya Bonandolok KM. 3,5 Kompleks Perkantoran Purba Dolok- Doloksanggul Kab. Humbang Hasundutan Prov. Sumatera Utara 22457
2) Departemen Ilmu Tanah dan Sumberdaya Lahan, Fakultas Pertanian, Jl. Meranti Kampus IPB Dramaga, Bogor 16680

## ABSTRACT
The growth of a region has changed land-use pattern, not only in urban areas but also rural areas such as Humbang Hasundutan regency, North Sumatra Province. As one district that newly developed since its establishment in 2003, land-use change  has  affected  the  environment  such as  the  emergence  of  degraded  land due  to  deforestation  and other  environmental vulnerability issues in some sub districts.  CLUE-S (the Conversion of Land Use and its Effects at Small region extents) model with a combination of dynamic systems and quantification of empirical relationships between land-use and its driving factors was  chosen  to  predict  land-use  change  that  improves  an  understanding  of  land  allocation  planning.  The  aims  of  this  study were to analyze land-use change in Humbang Hasundutan district and predict the changes in year 2033. The results showed Humbang  Hasundutan  district  has  5,362  hectares  of  forest  cover  loss  in  period  of  years  2003  to  2013.    Land  use  change simulation  by  using  scenario with  restricted  area  policy  and  forest  land  rehabilitation  can  maintain  forest  cover  over  30  % and contribute to augment forest vegetation on degraded lands by 47.87 % in year 2033.
Keywords: CLUE-S, degraded land, forest cover, land-use change

## ABSTRAK
Pertumbuhan suatu wilayah merubah pola penggunaan lahan tidak hanya di kawasan perkotaan (urban) akan tetapi juga terjadi  pada  wilayah  dengan  karakteristik  perdesaan  (rural)  seperti  Kabupaten  Humbang  Hasundutan.  Sebagai  salah  satu kabupaten  yang  baru  berkembang  sejak  terbentuk  tahun  2003,  perubahan  penggunaan  lahan  di  Kabupaten  Humbang Hasundutan  mengindikasikan  terjadinya  dampak  negatif  terhadap  lingkungan  seperti  munculnya  lahan  terdegradasi  akibat deforestasi  dan  isu  kerentanan  lingkungan  di  beberapa  kecamatan  yaitu  erosi  dan  bahaya  longsor.  Model  CLUE-S  ( the Conversion  of  Land  Use  and  its  Effects  at  Small  region  extents )  adalah  suatu  pemodelan  spasial  untuk  mensimulasikan perubahan penggunaan lahan berdasarkan kuantifikasi hubungan secara empirik berbagai faktor pendorong untuk memprediksi perubahan  penggunaan  lahan  di  masa  yang  akan  datang.  Hasil  analisis  perubahan  penggunaan  lahan  menunjukkan  bahwa  di Kabupaten Humbang Hasundutan terjadi deforestasi seluas 5,362 hektar selama periode tahun 2003-2013. Proyeksi perubahan penggunaan lahan berdasarkan laju perubahan alami menunjukkan luas penggunaan lahan hutan pada akhir tahun 2033 hanya tersisa  sebesar  27.5%  dari  luas  wilayah  Kabupaten  Humbang  Hasundutan.  Kondisi  ini  menjadi  tidak  sesuai  dengan  proporsi luas  hutan  yang  dialokasikan  menurut  pola  ruang  dalam    penyusunan  Rencana  Tata  Ruang  Wilayah  (RTRW)  yaitu  sebesar 30%  dari  luas  wilayah.  Pemodelan  menggunakan  skenario  kawasan  terbatas  dan  alokasi  rehabilitasi  lahan  hutan  dapat mempertahankan  hutan  lebih  dari  30%  dan  berkontribusi  mengembalikan  penutupan  vegetasi  hutan  pada  lahan  terdegradasi sebesar 47.87%.
Kata kunci: CLUE-S, lahan terdegradasi, hutan, perubahan penggunaan lahan

## PENDAHULUAN
Pengalokasian lahan berdasarkan keterkaitan fungsi dan  kemampuan  lahan  diharapkan  dapat  mengendalikan perubahan penggunaan lahan serta meminimalkan dampak negatif terhadap lingkungan. Penggunaan lahan yang tidak terkendali  sebagaimana  dikemukakan  Rossiter  (1996) dapat  mengarahkan  pada  inefisiensi  pemanfaatan  sumber
daya  alam,  rusaknya  sumber  daya  lahan,  kemiskinan  dan masalah sosial lainnya.
Pertumbuhan wilayah Kabupaten Humbang Hasundutan  terbilang  relatif  pesat  dibandingkan  dengan kabupaten  di  sekitarnya.  Hal  ini  ditunjukkan  dengan adanya  peningkatan  positif  nilai  PDRB  dari  5.32%  tahun 2009  naik  menjadi  6.03%    tahun  2013  (BPS,  2013). Perkembangan  ini  pada  dasarnya  telah  merubah  pola penggunaan  lahan  di  Kabupaten  Humbang  Hasundutan
*) Penulis Korespondensi: Telp. +6285214138588;  Email. togapanda@gmail.com
ditandai  dengan  munculnya  pusat-pusat  kegiatan  di  setiap kecamatan  dan  terbentuknya  permukiman-permukiman baru, akan tetapi perubahan tersebut juga mengindikasikan adanya  dampak  negatif  terhadap  lingkungan  yaitu timbulnya  lahan  terdegradasi.  Menurut  data  spasial  lahan kritis  BPDAS  WU  (2013),  di  Kabupaten  Humbang Hasundutan  terdapat  lahan  terdegradasi  seluas  27,870 hektar  (11.2%)  dengan  kategori  sangat  kritis;  68,274 hektar (27.5%) kategori kritis; 88,890 hektar (35.9%) agak kritis;  40,453  hektar  (16.3%)  potensial  kritis;  dan  22,365 hektar  (9%) termasuk tidak kritis.
Munculnya lahan terdegradasi menyebabkan turunnya  kualitas  suatu  lahan  dan  mengarah  pada rendahnya produktivitas dari lahan tersebut. Hal ini sesuai dengan  data  BPS  (2013)  tahun  2007  s.d  2011,  bahwa produksi  tanaman  pangan  di  Kabupaten  Humbang Hasundutan  mengalami  penurunan  khususnya  produksi padi  (gabah  kering  giling)  yang  semula  111,181  ton menjadi  103,410  ton.  Indikasi  lainnya  terkait  isu kerentanan  lingkungan  di  beberapa  lokasi  pada  daerah perbukitan seperti Kecamatan Pakkat, Parlilitan,
Tarabintang, Onanganjang dan Baktiraja pernah dinyatakan  waspada  longsor  oleh  BPBD  Humbang Hasundutan.
Kim et  al. (2011)  menyebutkan  suatu  prediksi perubahan  penggunaan  lahan  dapat  menyediakan  peluang untuk menghadapi permasalahan lingkungan dan perubahan  regional  di  masa  depan.    Belum  diketahuinya perubahan  penggunaan  lahan  di  Kabupaten  Humbang Hasundutan  dan  faktor-faktor  yang  mempengaruhi perubahan  tersebut  mendasari  pentingnya  penelitian dilakukan sekaligus memproyeksikan perubahan penggunaan lahan di masa yang akan datang dalam rangka membangun  kebijakan  terkait  perencanaan  tata  ruang  dan alokasi penggunaan lahan.
Kerangka awal model Conversion of Land Use and Its  Effects (CLUE)  yang  dikembangkan  Veldkamp  dan Fresco  (1996),  Verburg et  al. (1999),  dan  Verburg et  al . (2002) dibuat untuk mensimulasikan perubahan penggunaan  lahan  untuk  wilayah  dengan  cakupan  yang sangat  luas  (skala  global).  Oleh  karena  perbedaan representasi  data  dan  fitur  lain  untuk  diaplikasikan  secara regional,  model  CLUE  tidak  dapat  diterapkan,  sehingga dikembangkan  pendekatan  pemodelan  CLUE-S  ( the Conversion  of  Land  Uses  and  its  Effects  at  Small  region extent )  dengan  menggunakan  resolusi  spasial  tidak  lebih dari  1  x  1  km.  Simulasi  perubahan  penggunaan  lahan menggunakan  CLUE-S  dapat  dijalankan  dengan  berbagai kondisi  skenario.  Penerapan  menggunakan  skenario  yang berbeda-beda  dapat  mendukung  dalam  mengeksplorasi perubahan  penggunaan  lahan  di  masa  yang  akan  datang (Verburg et al., 2002).
Penelitian  ini  bertujuan  untuk  memproyeksikan perubahan  penggunaan  lahan  Kabupaten  Humbang Hasundutan pada Tahun 2033 dengan menerapkan kondisi skenario  yang  dapat  mempertahankan  proporsi  hutan sesuai ketentuan pola ruang berdasarkan UU  No 26 tahun 2007 tentang Penataan Ruang.

## BAHAN DAN METODE
Kabupaten Humbang Hasundutan  merupakan  salah satu kabupaten di wilayah Provinsi Sumatera Utara dengan
luas wilayah 247,848 hektar terdiri atas 10 kecamatan, 153 desa  dan  1  kelurahan,  dengan  jumlah  penduduk  sebanyak 176,429  jiwa  (BPS,  2013).  Wilayah  secara  geografis terletak  pada  posisi  98°  15.08’  -  98°  57.88’  Bujur  Timur dan  1°  59.91’  -  2°  30.84’  Lintang  Utara  (Gambar  1) dengan  karakteristik  wilayah  dataran  tinggi  pegunungan pada ketinggian 65 s.d 2,025 m dpl.
Bahan  yang  digunakan  dalam  penelitian  ini  adalah citra Landsat 7 ETM+ Tahun 2003 dan citra Landsat 8 OLI Tahun 2013. Peta administrasi, jaringan jalan, curah hujan diperoleh  dari  data  tematik  RTRW  Kabupaten  Humbang Hasundutan  Tahun  2011-2013  skala  1  :  100,000.  Data kependudukan  dan  jumlah  rumah  tangga  pertanian diperoleh  dari  BPS  Humbang  Hasundutan  Tahun  2013. Data elevasi dan lereng diperoleh dari data SRTM ( Shuttle Radar  Topography  Mission )  resolusi  90  meter.  Peta penetapan  moratorium  izin  baru  pada  hutan  primer  dan lahan  gambut  skala  1  :  250,000  bersumber  dari Kementrian  Lingkungan  Hidup  dan  Kehutanan  Tahun 2015.  Peta  Satuan  Lahan  dan  Tanah  Skala  1  :  250,000 Lembar  0618  Sidikalang,  Sumatera  bersumber  dari Puslittan  Tahun  1989.  Peta  Sungai  dan  Peta  Lahan  Kritis Skala  1  :  100,000  bersumber  dari  BPDAS  Wampu  Sei- Ular, Medan. Data sasaran kegiatan Rehabilitasi Hutan dan Lahan  (RHL)  menurut  Buku  Rencana  Pengelolaan Rehabilitasi  Hutan  dan  Lahan  (RP-RHL)  Kabupaten Humbang  Hasundutan  Tahun  2013  –  2017.    Alat  yang digunakan  dalam  penelitian  ini  adalah  perangkat  lunak Dyna CLUE-S v2.0, Notepad , pengolah data penginderaan jauh, pengolah peta dan pengolah data statistik.
Klasifikasi  penggunaan  lahan  diperoleh  melalui interpretasi  visual  terhadap  citra  satelit  dalam  dua  titik tahun  yaitu  Tahun  2003  dan  Tahun  2013  dan  dievaluasi berdasarkan  nilai overall  accuracy dan kappa  coefficient melalui  konfirmasi  100  titik  uji  ( groundtruth )  terhadap aktual di lapangan dan citra resolusi tinggi Google Earth . Analisis faktor-faktor yang mempengaruhi perubahan  penggunaan  lahan  dilakukan  menggunakan regresi  logistik  biner  dijalankan  dengan  metode  forward stepwise yang bertujuan untuk mengeliminasi faktor-faktor pengaruh  yang  tidak  signifikan  dalam  model.  Regresi logistik  biner  adalah  model  regresi  dimana  variabel dependen  bersifat  biner  atau  dikotomi.  Fungsi  regresi merupakan  suatu  respon  monotonic  garis  kurva  dibatasi antara 0 dan 1 dengan fungsi logistik pada persamaan,
( ) ( ) ( )
dimana  pi  adalah  nilai  probabilitas  peningkatan  sel penggunaan  lahan  sedangkan  E(Y)  perolehan  nilai  dari variabel  dependen biner Y,  β0 adalah nilai  konstan  yang diestimasi,  βn  adalah  koefisien  yang  diestimasi  dari masing-masing  variabel  independen  Xi  (Pontius  dan Schneider,  2001).  Variabel  dependen  yang  terdiri  dari setiap  kelas  penggunaan  lahan  ditransformasikan  dalam peta raster biner (0 dan 1) dengan resolusi piksel 100 x 100 meter.  Nilai  1  menunjukkan  satu  kelas  penggunaan  lahan sedangkan  nilai  0  tidak  terdapat  penggunaan  lahan. Variabel  independen  berupa  faktor-faktor  biofisik  dan sosial ekonomi diklasifikasikan berdasarkan kelas kategori masing-masing tipe data dari nilai terkecil sampai terbesar dengan  metode  peng-kelasan  natural  breaks  dalam  format
raster  dengan  resolusi  piksel  yang  sama  dengan  kelas penggunaan  lahan.  Untuk  memperoleh  nilai  numerik  tiap sel  raster  dari  masing-masing  variabel,  data  raster seluruhnya  dikonversi  ke  dalam  format  ASCII  text  dan kombinasi nilai sel antar variabel ditabulasikan berdasarkan  formulasi  regresi  logistik  serta  selanjutnya diolah menggunakan aplikasi statistik.
Proyeksi  perubahan  penggunaan  lahan  didasarkan pada pola perubahan pada minimal dua titik waktu (t0 dan t1) untuk menghasilkan pola perubahan penggunaan lahan pada  waktu (t2). Simulasi perubahan menggunakan model CLUE-S  diproyeksi  sampai  Tahun  2033  ( baseline 20 tahun  berdasarkan  jangka  waktu  Rencana  Tata  Ruang Wilayah).  Pengaturan  dilakuan  terhadap  paramater  yang menjadi  input  dalam  aplikasi  CLUE-S  terdiri  atas  kelas penggunaan  lahan,  matriks  konversi,  elastisitas  konversi, kebutuhan  lahan,  karakteristik  lokasi  penggunaan  lahan dan  kondisi  skenario  (Verburg,  2010).  Seluruh  rangkaian proses analisis diilustrasikan dengan Gambar 2.

## HASIL DAN PEMBAHASAN

## Perubahan Penggunaan Lahan Kabupaten Humbang Hasundutan
Interpretasi penggunaan lahan tahun 2003 dan 2013 menghasilkan  sembilan  kelas  penggunaan  lahan  yaitu belukar,  rawa/  gambut,  hutan,  hutan  tanaman,  lahan terbuka,  permukiman,  pertanian  lahan  kering,  sawah  dan tubuh  air  (Gambar  3).  Masing-masing  nilai overall accuracy dan kappa  coefficient klasifikasi  penggunaan lahan  terkonfirmasi  untuk  tahun  2003  sebesar  0.95%  dan 0.91, sedangkan untuk tahun 2013 sebesar 0.93% dan 0.90. Hal  ini  menunjukkan  hasil  interpretasi  sangat  baik  dan dapat dilanjutkan untuk analisis selanjutnya.
Perubahan  penggunaan  lahan  di  Kabupaten Humbang  Hasundutan  periode  2003  s.d 2013  terjadi  pada areal  seluas  13,510  hektar  yang  terdiri  dari  penggunaan lahan  yang  mengalami  penurunan  luas  terjadi  hutan sebesar 6.4%, rawa/gambut sebesar 44% dan lahan terbuka sebesar  5%  sedangkan  peningkatan  luas  terjadi  pertanian lahan  kering  sebesar  5.8%,  belukar  sebesar  19%,  lahan terbangun sebesar 39.4% dan sawah sebesar 3%.
Gambar 1.  Wilayah studi
Gambar 2.  Metode analisis perubahan pengunaan dengan Model CLUE-S
" J ! ! ! ! ! ! ! ! ! D A N A U   T O B A KEC. PARLILITAN KEC. PAKKAT KEC. POLLUNG KEC. TARA BINTANG KEC. ONAN GANJANG KEC. DOLOK SANGGUL KEC. SIJAMA POLANG KEC. LINTONG NIHUTA KEC. PARANGINAN KEC. BAKTI RAJA KABUPATEN HUMBANG HASUNDUTAN KABUPATEN TAPANULI UTARA KABUPATEN TAPANULI TENGAH KABUPATEN PAK-PAK BHARAT KABUPATEN SAMOSIR KABUPATEN TOBA SAMOSIR Pakkat Pollung Bakt ir aja Parlilitan Paran gin an Tar abintan g Sijamap ola ng On an Ganjang Doloksa nggul Lin tong N ih uta 99°0'E 99°0'E 98°45'E 98°45'E 98°30'E 98°30'E 98°15'E 98°15'E 2° 30 'N 2° 30 'N 2° 15 'N 2° 15 'N 2° 0' N 2° 0' N S A M U D E R A IN D O N E S IA RIAU NAD SUMATERA UTARA 10 0 10 20 30 5 KM Perairan Jalan Raya " J Ibukota Kabupaten ! Ibukota Kecamatan KETERANGAN : µ
Pemodelan Perubahan Penggunaan Lahan dengan Menggunakan Clue-S (Sinurat, T.P., K. Munibah, dan D.P.T. Baskoro)
(a) Tahun 2003
Perubahan  penggunaan  lahan  pada  periode  ini menujukkan  terjadinya  deforestasi  hutan  seluas  5,362 hektar dengan laju penyusutan 536.2 hektar per tahunnya. Pola  perubahan  penggunaan  lahan  hutan  adalah  berubah menjadi pertanian lahan kering, belukar, lahan terbuka dan sawah.  Perbandingan perubahan luas masing-masing kelas penggunaan lahan disajikan dalam Tabel 1.

## Faktor-faktor yang Mempengaruhi Perubahan Penggunaan Lahan
Faktor-faktor yang mempengaruhi perubahan penggunaan  lahan  adalah  faktor  yang  diseleksi  dari  aspek sosial  ekonomi  yaitu  jumlah  penduduk  (jiwa)  dan  rumah tangga  pertanian  (RT),  sedangkan  aspek  biofisik  terdiri atas,  curah  hujan  (mm  thn -1 ),  elevasi  (m  dpl), slope (%) dan jenis tanah serta aspek aksesibilitas jarak dari ibukota kabupaten,  jarak  dari  kota  kecamatan,  jarak  dari  jalan utama,  jarak  dari  sungai  besar  yang  dihasilkan menggunakan metode euclidean distance .
Hubungan  berbagai  faktor  yang  mempengaruhi peningkatan  sel  tiap  jenis  penggunaan  lahan  dapat diketahui  dari  nilai  koefisien  beta  (β)  pada  Tabel  2a  dan nilai  peningkatan  peluang  pada  Exp  (β)  pada  Tabel  2b. Ketepatan  hasil  regresi  dievaluasi  berdasarkan  nilai  ROC ( Relative Operating Characteristic ) masing-masing penggunaan  lahan.  Nilai  ROC  >  0.5  menunjukkan pengaruh  faktor-faktor penduga dalam  model tidak terjadi secara kebetulan (Pontius dan Schneider, 2001).
(b) Tahun 2013
Nilai  ROC  dalam  model  regresi  diperoleh  untuk rawa/gambut  (0.995),  lahan  terbangun  (0.949),  hutan (0.905),  lahan  terbuka  (0.896),  sawah  (0.877),  pertanian lahan  kering  (0.861),  dan  belukar  (0.834).  Hutan  tanaman dan  tubuh  air  tidak  diregresikan  karena  bersifat  statis (tidak  mengalami  perubahan  luas)  dan  nilai  konstanta ditetapkan sebesar 0.7 (Verburg, 2010).

## Model Perubahan Penggunaan Lahan
Penggunaan  lahan  tahun  2003  dan  tahun  2013 ditranformasikan  menjadi  ASCII  text  dengan  penamaan file  (cov_all.0).  Matriks  konversi  adalah  nilai  perubahan yang  diperkenankan  terjadi  pada  setiap  jenis  penggunaan lahan,  dihasilkan  dari  tabulasi  silang  penggunaan  lahan pada dua titik tahun dan disusun dalam format text dengan nama  (allow.txt).  Nilai  0  menandakan  jenis  penggunaan lahan  tidak  boleh  berubah  sedangkan  nilai  1  dapat berubah.  Matriks  konversi  penggunaan  lahan  disajikan pada  Tabel  3.  Elastistitas  konversi  berhubungan  dengan mudah  tidaknya  suatu  penggunaan  lahan  terkonversi  dan mengalami  perubahan  kembali  ke  sifat  semula.  Nilai elastisitas  berada  pada  rentang  0  dan  1,  mendekati  1 diartikan  jenis  penggunaan  lahan  tersebut  tidak  mudah berubah  sebaliknya  nilai  elastisitas  0  menandakan  jenis tersebut  mudah  terkonversi  (Verburg,  2010).  Nilai elastisitas masing-masing penggunaaan lahan pada Tabel 4 diatur dalam file (alloc1.reg).
JENIS PENGGUNAAN LAHAN Belukar Rawa/ Gambut Hutan Hutan Tanaman Lahan Terbuka Lahan Terbangun Pertanian Lhn Kering Sawah Tubuh Air ¯ ¯
5 0 5 10 15 20 2.5 KM 5 0 5 10 15 20 2.5 KM KETERANGAN Red:    Layer_6 Green: Layer_5 Blue:   Layer_4
JENIS PENGGUNAAN LAHAN Belukar Rawa/ Gambut Hutan Hutan Tanaman Lahan Terbuka Lahan Terbangun Pertanian Lhn Kering Sawah Tubuh Air ¯ ¯ 5 0 5 10 15 20 2.5 KM 5 0 5 10 15 20 2.5 KM
No. Jenis Penggunaan Lahan 2003 2013 Luas Perubahan (ha) % (ha) % (ha) % 1 Belukar 8,074 3.3 9,606 3.9 1,532 19.0 2 Rawa/ Gambut 2,330 0.9 1,304 0.5 -1,026 -44.0 3 Hutan 83,768 33.8 78,406 31.6 -5,362 -6.4 4 Hutan Tanaman 51,613 20.8 51,613 20.8 0 0.0 5 Terbuka 7,351 3.0 6,984 2.8 -367 -5.0 6 Terbangun 340 0.1 474 0.2 134 39.4 7 Pertanian Lahan Kering 80,733 32.6 85,410 34.5 4,677 5.8 8 Sawah 13,520 5.5 13,932 5.6 412 3.0 9 Air 119 0.0 119 0.0 0 0.0 Total  Luas  (ha) 247,848 100 247,848 100 13,510 5.5 Tabel 1.  Perubahan luas penggunaan lahan periode tahun 2003 s.d 2013

## Gambar 3.  Klasifikasi penggunaan lahan
Pengaruh Faktor Jenis Penggunaan Lahan BLK 1 GBT 2 HUT 3 HTI 4 TBK 5 PM 6 PLK 7 SW 8 A 9  Jumlah Penduduk (jiwa) -1.033 3.097 -0.817 - 0.292 - - 0.496 -  Jumlah Rumah Tangga Pertanian (RT) 0.466 -1.704 0.871 - - 0.422 0.163 -0.462 -  Curah Hujan (mm/ thn) 1.053 -4.765 1.291 - -1.389 - 0.113 -0.733 -  Elevasi (m dpl) -0.482 - 0.993 - -0.575 - -1.352 -1.232 -  Slope (%) -0.032 -2.057 0.426 - 0.891 -0.430 -0.276 -0.506 -  Jarak dari Ibukota (m) -0.666 1.714 0.298 - -0.207 -0.109 -0.064 0.100 -  Jarak dari Kota Kecamatan (m) -0.428 -1.725 0.383 - -1.129 -0.654 -0.144 -0.230 -  Jarak dari Jalan Utama (m) 0.213 -1.116 0.074 - 0.556 -2.153 -0.483 -0.612 -  Jarak dari Sungai (m) -0.068 -0.110 0.150 - -0.161 - -0.030 0.236 -  Asosiasi tanah  dystropepts, distrandepts, tropudults - -2.850 -0.180 - - - 0.487 -0.627 -  Asosiasi tanah dystropepts, eutrandepts, dystrandepts - - -0.109 - - - 0.380 -0.593 -  Asosiasi tanah  dystropepts, tropudults, troportents - - 1.511 - - - -1.069 -1.179 - Constant 0.605 -0.859 -12.860 0.7 1.291 -1.967 4.795 4.513 0.7 ROC 0.834 0.995 0.905 - 0.896 0.949 0.861 0.877 - 1 BLK  (Belukar), 2 GBT  (Rawa/  gambut), 3 HUT  (Hutan), 4 HTI  (Hutan  Tanaman), 5 TBK  (Lahan  Terbuka), 6 PM (Lahan  Terbangun/  Permukiman), 7 PLK (Pertanian Lahan Kering), 8 SW (Sawah), 9 A (Tubuh Air)
Tabel 2a.  Nilai koefisien beta (β) pengaruh berbagai faktor terhadap kelas penggunaan lahan
Pengaruh Faktor Jenis Penggunaan Lahan
BLK 1 GBT 2 HUT 3
HTI 4
TBK 5
PM 6
PLK 7
SW 8 A 9  Jumlah Penduduk (jiwa) 0.356 22.14 0.441 - 1.338 - - 1.642 -  Jumlah Rumah Tangga Pertanian (RT) 1.593 0.182 2.388 - - 1.525 1.177 0.63 -  Curah Hujan (mm/ thn) 2.866 0.009 3.636 - 0.249 - 1.12 0.48 -  Elevasi (m dpl) 0.618 2.698 - 0.563 - 0.259 0.292 -  Slope (%) 0.969 0.128 1.531 - 2.436 0.651 0.759 0.603 -  Jarak dari Ibukota (m) 0.514 5.552 1.346 - 0.813 0.897 0.938 1.105 -  Jarak dari Kota Kecamatan (m) 0.652 0.178 1.467 - 0.323 0.52 0.866 0.795 -  Jarak dari Jalan Utama (m) 1.237 0.328 1.077 - 1.744 0.116 0.617 0.542 -  Jarak dari Sungai (m) 0.935 0.895 1.161 - 0.852 - 0.971 1.266 -  Asosiasi tanah dystropepts, distrandepts, tropudults - 0.058 0.835 - - - 1.628 0.534 -  Asosiasi tanah dystropepts, eutrandepts, dystrandepts - - 0.896 - - - 1.463 0.553 -  Asosiasi  tanah dystropepts, tropudults, troportents - - 4.531 - - - 0.343 0.307 -
1 BLK  (Belukar), 2 GBT  (Rawa/  gambut), 3 HUT  (Hutan), 4 HTI  (Hutan  Tanaman), 5 TBK  (Lahan  Terbuka), 6 PM  (Lahan  Terbangun/  Permukiman), 7 PLK (Pertanian Lahan Kering), 8 SW (Sawah), 9 A (Tubuh Air)
Tabel 2b.  Nilai peluang Exp (β) kelas penggunaan lahan
Penggunaan Lahan BLK 1
GBT 2 HUT 3 HTI 4 TBK 5 PM 6 PLK 7 SW 8 A 9 BLK 1 1 0 0 0 1
0 0 0 0 GBT 2 0 1 0 0 0
1 1 1 0 HUT 3 1 0 1 0 1 0 1 1 0 HTI 4 0 0 0 1
0 0 0 0 0 TBK 5 0 0 0 0 1
1 1 0 0 PM 6 0 0 0 0 0
1 0 0 0 PLK 7 1 0 1 0 1
1 1 1 0 SW 8 0 0 0 0 0
1 1 1 0 A 9 0 0 0 0 0
0 0 0 1 1 BLK (Belukar), 2 GBT (Rawa/ gambut), 3 HUT (Hutan), 4 HTI (Hutan Tanaman), 5 TBK (Lahan Terbuka), 6 PM (Lahan Terbangun/ Permukiman),
7 PLK (Pertanian Lahan Kering), 8 SW (Sawah), 9 A (Tubuh Air)
Tabel 3.  Matriks konversi penggunaan lahan
Kelas Penggunaan Lahan BLK 1
GBT 2 HUT 3 HTI 4 TBK 5 PM 6 PLK 7 SW 8 A 9 Elastisitas Konversi 0.5 0.4 0.7 1 0.3 1 0.5 0.8 1 1 BLK (Belukar), 2 GBT (Rawa/ gambut), 3 HUT (Hutan), 4 HTI (Hutan Tanaman), 5 TBK (Lahan Terbuka), 6 PM (Lahan Terbangun/ Permukiman), 7 PLK (Pertanian Lahan Kering), 8 SW (Sawah), 9 A (Tubuh Air)
Tabel 4.  Nilai elastisitas jenis penggunaan lahan
Moratorium region_park1
-9,998 0 Keterangan : locspec 0
1
Kebutuhan  penggunaan  lahan  adalah  tabulasi  rata- rata  luas  perubahan  penggunaan  lahan  per  tahun  prediksi yang  disusun  dalam  format  text  dengan  penamaan (demand.in)  dan  dimodifikasi  sesuai  dengan  kebutuhan skenario.

## Validasi Model
Validasi model CLUE-S menggunakan penggunaan lahan  tahun  2003  memprediksi  tahun  2013  dan dibandingkan  dengan  penggunaan  lahan  aktual  tahun 2013. Simulasi ini menghasilkan nilai Kappa sebesar 0.96 dari  tabulasi  silang  masing-masing  kelas  penggunaan lahan.  Hal  ini  menunjukkan  model  sangat  baik  dan  dapat dilanjutkan untuk simulasi proyeksi perubahan penggunaan lahan. Penggunaan lahan hasil simulasi untuk validasi model ditunjukan pada Gambar 4.

## Skenario dalam Simulasi
Model  simulasi  dijalankan  dengan  dua  kondisi skenario.  Skenario  0  merupakan  kondisi  laju  perubahan alami  tanpa  adanya  larangan  konversi  untuk  masing- masing  penggunaan  lahan.  Skenario  1  disusun  dengan menetapkan  areal  yang  tidak  diperkenankan  terjadi perubahan  berdasarkan  peta  moratorium  izin  baru  pada hutan primer dan kawasan bergambut serta mengalokasikan rehabilitasi hutan pada lahan sangat kritis.
Kabupaten  Humbang  Hasundutan  ditetapkan moratorium  izin  baru  (Gambar  5)  seluas  2,202.6  hektar untuk  kawasan  bergambut  dan  seluas  60,085.6  hektar berupa  kawasan  hutan  primer  (Kemenlinghut,  2015). Rehabilitasi  hutan  dialokasikan  seluas  500  hektar  per tahun yang ditetapkan berdasarkan data sasaran rehabilitasi hutan  menurut  Buku  Rencana  Pengelolaan  Rehabilitasi Hutan dan Lahan Tahun 2013 – 2017 Kabupaten Humbang Hasundutan (Dishutlindup, 2013). Gambar 6 menunjukkan lokasi  spesifik  rehabilitasi  hutan  pada  lahan  sangat  kritis. Perbandingan  kebutuhan  luas  dan  alokasi  lahan  setiap skenario  diilustrasikan  dalam  Gambar  7  menurut  hasil tabulasi  rata-rata  perubahan  penggunaan  lahan  per tahunnya sesuai jangka waktu proyeksi.

## Proyeksi Perubahan Penggunaan Lahan Tahun 2033
Hasil  simulasi  perubahan  penggunaan  lahan  untuk proyeksi  tahun  2033  disajikan  pada  Gambar  8.  Tabel  5 menunjukkan bahwa simulasi dengan laju perubahan alami (skenario  0)  mengalami  perubahan  areal  terluas  yaitu sebesar  22,400  hektar  sedangkan  perubahan  lahan  pada skenario 1 hanya seluas 2,036 hektar.
¯ ¯ 5 0 5 10 15 20 2.5 KM 5 0 5 10 15 20 2.5 KM ¯ ¯ 5 0 5 10 15 20 2.5 KM 5 0 5 10 15 20 2.5 KM KETERANGAN Belukar Rawa/ Gambut Hutan Hutan Tanaman Lahan Terbuka Lahan Terbangun Pertanian Lahan Kering Sawah Tubuh Air
Gambar 4.  Penggunaan lahan hasil validasi model
Gambar 5.  Kawasan  moratorium  izin  pada  hutan  primer  dan kawasan bergambut Gambar 6.  Lahan terdegradasi dengan kategori sangat kritis Moratorium region_park1 -9,998 0 Keterangan : locspec 0 1
demand.in0 (Skenario 0)

## (a) Skenario 0
Melalui  pemodelan spasial menggunakan CLUE-S, proyeksi  perubahan  penggunaan  lahan  pada  tahun  2033 berdasarkan  skenario  0  menunjukkan  bahwa  deforestasi meningkat  seluas  10,274  hektar  ke  wilayah  barat Kabupaten  Humbang  Hasundutan  sehingga  luas  areal berhutan  tersisa  sebesar  27.5  %.  Namun,  skenario  1 menunjukkan deforestasi dapat ditekan menjadi seluas 286 hektar dan luas areal berhutan dapat dipertahankan sampai 31.5 %.
demand.in1(Skenario 1) (b) Skenario 1
Simulasi  CLUE-S  dengan  kebijakan  kawasan terbatas  dan  rehabilitasi  hutan  (skenario  1)  dapat mengembalikan  penutupan  vegetasi  hutan  pada  lahan sangat  sangat  kritis  sebesar  42.87  %.    Menurut  Peraturan Menteri Kehutanan Nomor P.32/ Menhut-II/ 2009 tentang Tata Cara Penyusunan Rencana Teknik Rehabilitasi Hutan dan Lahan Daerah Aliran Sungai (RTKRHL-DAS) bahwa penilaian kriteria lahan kritis berdasarkan tutupan vegetasi mencapai 41 – 60% termasuk dalam kelas kategori sedang. Kontribusi  tutupan  vegetasi  masing-masing  skenario disajikan dalam Tabel 6.

## KETERANGAN :
Belukar Rawa/ Gambut Hutan Hutan Tanaman Lahan Terbuka Lahan Terbangun
Pertanian Lahan Kering Sawah Tubuh Air 4.5 0 4.5 9 13.5 18 2.25 KM ¯ KETERANGAN : Belukar Rawa/ Gambut Hutan Hutan Tanaman Lahan Terbuka Lahan Terbangun Pertanian Lahan Kering Sawah Tubuh Air 4.5 0 4.5 9 13.5 18 2.25 KM ¯
Gambar 7.  Kebutuhan penggunaan lahan proyeksi tahun 2033
Gambar 8.  Proyeksi penggunaan lahan tahun 2033
No PL 1
Skenario 0 Prb 2 (ha) Skenario 1 Prb 2 (ha) (ha) % (ha) % 1 BLK 3 11,914 4.8 2,308 9,791 4.0 185 2 GBT 4 1,304 0.5 0 1,285 0.5 -19 3 HUT 5 68,132 27.5 -10,274 78,120 31.5 -286 4 HTI 6 51,613 20.8 0 51,613 20.8 0 5 TBK 7 6,058 2.4 -926 6,271 2.5 -713 6 PM 8 690 0.3 216 688 0.3 214 7 PLK 9 93,251 37.6 7,841 85,561 34.5 151 8 SW 10 14,767 6.0 835 14,400 5.8 468 9 A 11 119 0.0 0 119 0.0 0 Total (ha) 247,848 100 22,400 247,848 100 2,036 1 PL  (Penggunaan  Lahan), 2 Prb  (Luas  Perubahan), 3 BLK  (Belukar), 4 GBT  (Rawa/  gambut), 5 HUT  (Hutan), 6 HTI  (Hutan  Tanaman), 7 TBK  (Lahan Terbuka), 8 PM (Lahan Terbangun/ Permukiman), 9 PLK (Pertanian Lahan Kering), 10 SW (Sawah), 11 A (Tubuh Air)
Tabel 5.  Perubahan luas penggunaan lahan pada masing-masing kondisi skenario

## SIMPULAN
Model  perubahan  penggunaan  lahan  menggunakan CLUE-S  untuk  mempertahankan  proporsi  luas  hutan sebesar  30  %  sampai  akhir  proyeksi  tahun  2033 sebagaimana  disyaratkan  dalam  penyusunan  pola  ruang berdasarkan  UU  No  26  Tahun  2007  tentang  Penataan Ruang  adalah  dengan  menerapkan  kebijakan  larangan konversi lahan pada hutan primer dan  lahan gambut serta pengalokasian rehabilitasi hutan pada lahan sangat kritis.

## UCAPAN TERIMA KASIH
Ucapan terima kasih kepada Pusdiklatren BAPPENAS  RI  yang  telah  membiayai  penelitian  ini melalui  anggaran  Beasiswa  BAPPENAS  tahun  2014, kepada Institut  for  Environmental  Studies ,  Universitas Amsterdam  atas  pengembangan  perangkat  lunak  Dyna CLUE-S  v2.0  yang  dapat  digunakan  secara  terbuka  dan Bupati  Humbang  Hasundutan  beserta  jajarannya  atas  izin penelitian dan bantuan data-data dalam penelitian.

## DAFTAR PUSTAKA

## [BPDAS] Balai Pengelolaan Daerah Aliran Sungai.  2013.
Up  Dating/  Review  Lahan  Kritis  Wilayah  Kerja BPDAS Wampu Sei Ular . Balai Pengelolaan Daerah Aliran Sungai Wampu Sei Ular, Medan.
[BPS]  Badan  Pusat  Statistik.  2013. Kabupaten  Humbang Hasundutan  dalam  Angka. Humbang  Hasundutan, Doloksanggul.
[Dishutlindup]  Dinas  Kehutanan  dan  Lingkungan  Hidup. 2013. Rencana Pengelolaan Rehabilitasi Hutan dan Lahan  Tahun  2013-2017 .  Dinas  Kehutanan  dan
Lingkungan Hidup Kabupaten Humbang Hasundutan,   Doloksanggul.
[Kemenlinghut]  Kementerian  Lingkungan  Hidup  dan Kehutanan.  2015. Keputusan  Menteri  Kehutanan dan Lingkungan Hidup Nomor : SK. 2312/ Menhut- VII/ IPSDH/ 2015 tentang Penetapan Peta Indikatif Penundaan Izin Baru Pemanfaatan Hutan, Penggunaan  Kawasan  dan  Perubahan  Peruntukan Kawasan  Hutan  dan  Areal  Penggunaan  Lain (Revisi  VIII) .  Kementerian  Lingkungan  Hidup  dan Kehutanan RI, Jakarta.
Kim,  I.,  G.Y.  Jeong,  S.  Park,  and    J.  Tenhunen.  2011. Predicted  Land  Use  Change  in  the  Soyang  River Basin . 2011 Terreco Science Conference. Karlsruhe Institute  of  Technology,  Garmisch-Partenkirchen, Germany. p. 17-24.
Pontius,  R.G.,  and    L.C.  Schneider.  2001.  Land-cover change model validation by an ROC method for the Ipswich watershed, Massachusetts, USA. Agriculture,  Ecosystems  &  Environment ,  85:  239- 248.
Rossiter,  D.G.  1996.  A  theoretical  framework  for  land evaluation. Geoderma , 72: 165-190.
Veldkamp,  A.,  and  L.O.  Fresco.  1996.  CLUE-CR:  An integrated  multi-scale  model  to  simulate  land  use change  scenarios  in  Costa  Rica. Ecological Modeling , 91:231-248.
Veldkamp,  A.,  and  E. Lambin. 2001. Predicting land-use change. Agriculture,  Ecosystems  &  Environment , 85: 1-6.
Verburg,  P.H.,  W.  Soepboer,  A.  Veldkamp,  R.  Limpiada, V. Espaldon, and S.S.A. Mastura. 2002. Modelling the  spatial  dynamics  of  regional  land  uses:  The CLUE-S  model. Environmental  Management , 30:391-405.
Verburg,  P.H.  2010. The  CLUE  Modelling  Framework. Course Material . University of Asterdam,
Amsterdam.
No. Kondisi Simulasi Lahan Terdegradasi (2013) Kontribusi Penutupan Vegetasi Kategori Lahan Kritis (P.32/Menhut-II/09) Skor Luas (ha) % 1 Skenario 0 27,877 7,627 27.36 Buruk 100 2 Skenario 1 27,877 11,952 42.87 Sedang 150 Tabel 6.  Kontribusi terhadap penutupan vegetasi hutan pada lahan terdegradasi