Abstrak —  Indonesia  merupakan  negara  yang  mempunyai kondisi  geologis  yang  unik  karena  berada  pada  pertemuan  3 lempeng  tektonik  besar,  yaitu  Lempeng  Indo-Australia, Lempeng  Eurasia  dan  Lempeng  Pasifik. Indonesia juga berada pada zona “Cincin Api Pasifik”, yaitu daerah barisan gunung api vulkanik  yang  aktif.  Wilayah Indonesia yang berada di daerah cincin api pasifik tersebut membawa manfaat yang sangat besar. S alah  satu manfaatnya adalah potensi energi panasbumi. Area Lahendong merupakan lapangan panas bumi pertama dan satu- satunya  di  S ulawesi.  Cadangan    terbukti reservoir Lahendong adalah  sebesar  80  MW  dengan  potensi  pengembangan sebesar 150  MW.  Area  Lahendong  mempunyai  karakteristik  reservoir low  permeability namun  memiliki  temperatur  sangat  tinggi. S eperti energi-energi lainnya, energi panasbumi juga memiliki kelebihan dan kekurangan yang jamak,  karena tidak ada energi yang  benar-benar  sempurna  dampak  kebermanfaatnya.  S alah satu  dampak  negatif  yang  dapat  terjadi    adalah  terjadinya deformasi yang disebabkan oleh ekspl oitasi fluida panasbumi itu sendiri. Penelitian  ini bertujuan untuk memantau aktivitas dari eksploitasi  panasbumi  di  PLTP  Lahendong,  S ulawesi  Utara dengan    metode Differential  Interferometric    Synthetic  Aperture Radar (DInSAR) . S edangkan untuk pengolahannya menggunakan  perangkat  lunak  bebas  terbuka  ( open  source software ) GMTSAR untuk mengolah data satelit ALOS-PALSAR daerah  eksploitasi  panasbumi  Lahendong.  Metode  yang digunakan adalah  dengan  melakukan  proses  pemfokusan  data SAR ,    transformasi    koordinat    ke sistem radar menggunakan infromasi  orbit  yang  teliti, image  alignment , interferome dan phase    unwrapping menggunakan   algoritma SNAPHU . Hasil metode DInSAR mengindikasikan  bahwa  telah  terjadi  proses deformasi  di  sekitar  unit  1  dan  2  area  eksploitasi  panasbumi Lahendong  berupa  penurunan  muka  tanah sebesar 3 sampai 4 cm..
Kata Kunci — ALOS PALSAR , DInSAR, Eksploitasi Panasbumi, GMTSAR, Line of Sight

## I. PENDAHULUAN
AMPAI saat ini, Indonesia diketahui memiliki ketersediaan  potensi  energi  geothermal  yang  sangat melimpah,  yaitu  sekitar  40%  dari  cadangan  energi geothermal di  seluruh  dunia [1] .  Terdapat  265  titik  lokasi prospek geothermal yang tersebar di sepanjang jalur vulkanik yang  membentang  dari  Pulau  Sumatera,  Jawa,  Bali,  Nusa Tenggara,  Sulawesi,  Maluku,  dan  Papua. PLTP  Lahendong salah  satunya,  merupakan  lapangan  panas   bumi  pertama  dan satu-satunya  di  Sulawesi.  Berdasarkan  data  tahun  2010,  total
pembangkitan  PLTP Lahendong sebesar 60 MW  (3 x 20MW) dapat menyumbang hingga 40%   kebutuhan listrik di Sulawesi Utara [2] . Seiring  dengan  perkembangan  teknologi  yang sangat  pesat,  saat  ini  telah  memunculkan  berbagai  metode untuk  melakukan  pemantauan  terhadap  aktifitas  deformasi. InSAR ( Interferometric  Synthetic Aperture Radar ) merupakan salah satu teknik yang mampu digunakan untuk mengidentifikasi  terjadinya  deformasi  di  permukaan  tanah dengan memanfaatkan data dari satelit penginderaan jauh aktif yaitu RADAR ( Radio  Detectionand  Ranging ). Pada teknik ini dilakukan  pembentukan  interferogram  yang  didapat  dari  dua buah  atau  lebih  citra SAR ( master dan slave )  yang  direkam pada  objek  yang  sama di permukaan bumi namun pada waktu yang  berbeda.  Jika  suatu  titik  di  tanah  begerak,  maka  jarak antara sensor dengan titik akan berubah, dan berpengaruh pada nilai phase yang  direkam  oleh  sensor SAR .  Nilai  beda phase inilah yang menunjukkan adanya ground  movement . Interferogram  yang  terbentuk  memiliki  efek  topografi  dan atmosferik,  oleh  karena  itu  perlu  dilakukan  pengurangan dengan  proses differential  inferometry  SAR ( DInSAR ), sehingga  dihasilkan  informasi  deformasi serta  dapat diketahui besar pergerakan tanah ( subsidance atau uplift ) yang terjadi. Data yang digunakan dalam penelitian ini adalah data PALSAR ( Phased  Array  type  L  band  Synthetic  Aperture Radar ). PALSAR merupakan  salah  satu  sensor  aktif  dengan sistem SAR yang  dibawa  oleh  satelit ALOS ( Advanced  Land Observing Satellite ). Sensor yang membawa saluran L dengan panjang  gelombang  23,6  cm  ini  dapat  digunakan  pada  siang maupun malam  hari [3] . Dengan  menerapkan  teknik DInSAR dan  data ALOS PALSAR sebagai  data  pendukung,  maka  pada  penelitian  ini dilakukan  pengestimasian  besar  deformasi  di  area  eksploitasi geothermal akibat  eksploitasi fluida geothermal dengan studi kasus area eksploitasi geothermal Lahendong, Sulawesi Utara.

## II. METODOLOGI PENELITIAN

## A. Lok asi Penelitian
Lokasi  penelitian  ini  mengambil  studi  kasus  wilayah  kerja pengusahaan geothermal Lahendong, Provinsi Sulawesi Utara. Posisi  Geografisnya  terletak  pada  1 0 13’12”  –  1 0 18’36”  Lintang Utara dan 124 0 48’0” – 124 0 53’24” Bujur Timur dengan luas wilayah tercatat 12 km 2 [4] .

## Pemanfaatan  Metode Differential Intermerometry Synthetic Aperture Radar ( DInSAR ) untuk Pemantauan  Deformasi  Akibat Aktivitas  Eksploitasi Panasbumi
(Studi Kasus:  PLTP Lahendong  Unit 1 dan 2)
Roni Kurniawan dan Ira Mutiara Anjasmara Jurusan Teknik Geomatika, Fakultas Teknik Sipil  dan Perencanaan, Institut Teknologi  Sepuluh  Nopember (ITS)
Jl. Arief Rahman Hakim, Surabaya 60111 Indonesia
e-mail : ira@geodesy.its.ac.id
S
Gambar 1. Peta T opografi Area Lahendong [5

## B. Data yang Digunak an
Penelitian ini menggunakan data citra satelit ALOS PALSAR
Level  1.0  dengan Fine  Beam  Dual  Plorazitaion (FBD)  pada tanggal  20  Oktober  2007,  22  Oktober  2008,  25  Oktober  2009, dan  29  Oktober  2010.  Setelah  itu  dilakukan  proses pembentukan  citra interferogram dengan menggabungkan dua buah  citra  yang  dimana  salah  satu  citra  disebut  citra master dan  yang  citra  lainnya  disebut  citra slave. Penelitian  ini menggunakan  data  citra  tahun  2010  sebagai  citra master dan ketiga citra lainnya sebagai citra slave . Tabe l 1. Spesifikasi Citra ALOS PALSAR yang Digunakan
No Nam a Citra Perekam an Orbit Fram e Arah 1 ALPSRP 092560010 20 Okt 2007 9256 10 Ascending 2 ALPSRP 146240010 22 Okt 2008 14624 10 Ascending 3 ALPSRP 199920010 25 Okt 2009 19992 10 Ascending 4 ALPSRP 253600010 28 Okt 2010 25360 10 Ascending
Selain  itu,  pada  metode DInSAR ini  juga  diperlukan  data DEM yang  digunakan  untuk  proses differential  2-pass interferometry.  DEM yang  digunakan  dalam  penelitian  ini adalah DEM SRTM  1’ yang memiliki  ketelitian  horizontal 0,00027 o ,  sebanding  dengan  ~30  m  dan  ketelitian vertical dari DEM  SRTM  1’ ini  tidak  melebihi  16  m. DEM ini diunduh dari website http://earthexplorer.usgs.gov/ . DEM  SRTM  1’ digunakan  untuk  mengurangi phase topografi  yang  terdapat  pada phase interferogram  sehingga didapatkan phase deformasi [6] .  Dikarenakan  resolusi  citra SAR dan DEM berbeda, maka phase topografi yang ada tidak seluruhnya  hilang. DEM  SRTM  1 disimulasikan  dengan menggunakan  data parameter dari citra master yang terbentuk dari  tahap preprocess pengolahan  data SAR .  Tampilan DEM yang digunakan dapat dilihat pada Gambar  3

## C. Proses Pengolahan Citra dengan GMTSAR
Gambar 3. Diagram Alir Pengolahan Citra dengan GMT SAR
Preprocessing Alignment Citra Slave (SLC) Interfere Filtering Geocode Data DEM SRTM 1’ Dem2topophase File DEM grid: Topo_ra.grd ALOS PALSAR (Master) ALOS PALSAR (Slave)
Data geocode: amp_ll.grd, phase_ll.grd, corr_ll.grd, dan unwrap_ll.grd
File Data *.PRM Dan *.raw File Data *.PRM Dan *.raw Focussing Citra Slave (SLC) Citra Master (SLC) Data Interferogram: real.grd dan imag.grd
Gambar 2.T ampilan DEM SRTM 1’ Sebagian Provinsi Sulawesi  Utara
Penjelasan tahap pengolahan data SAR dengan GMTSAR : 1. Tahap preprocess
Tahap  ini  dilakukan  untuk  mendapatkan  informasi  orbit dan  mengetahui  parameter ascii  file yang berisi parameter untuk InSAR  processing dengan  mengekstrak orbital satellite dan  kecepatannya  dari leader  file  *.LED , image alignment ,  dan  mengestimasi doppler  centroid dari  data raw yang terdapat pada masing-masing data SAR . 2. Tahap focussing dan alignment
Tahap focus dilakukan  untuk  membuat  citra Single  Look Complex (SLC)  dari  masing-masing data SAR . Sedangkan tahap align dilakukan  untuk  mencari  area  tampalan  dan parameter yang didapatkan dari tahap preprocess oleh citra master terhadap citra slave . 3. Tahap Transformasi Koordinat (DEM2topophase) DEM yang  digunakan  adalah DEM  SRTM  1’ yang memiliki  sistem  koordinat  geografis.  Pemrosesan  yang dilakukan  mayoritas  dalam  sistem  koordinat  radar,  maka harus  dilakukan  transformasi  dari  sistem  koordinat geografis ke sistem koordinat radar.
Hasil dari  tahapan ini ada tiga jenis data, yaitu:
a. Data topo_ra.grd ,  merupakan    koordinat range dan azimuth dalam  sistem radar.
b. Postscript image DEM.ps dan topo_ra.ps.
c. Data trans.dat ,  merupakan  koordinat  sekutu  antara radar  dan DEM yakni  koordinat range , azimuth dan koordinat geografis 4. Tahap Interfere
Tahap  ini  dilakukan  untuk  membuat  interferogram  dan secara  bersamaan  mengurangi  efek  topografi  dengan menggunakan DEM yang ada. Tahap interfere dilakukan melalui komen intf.csh . Tahapan tersebut akan menghasilkan  dua file yaitu imag.grd dan real.grd yang kemudian  dihapus setelah tidak diperlukan lagi.
5. Tahap Filtering Tahap  ini  dilakukan  untuk  menaikkan  nilai signal  noise ratio (SNR) yang  terdapat  pada  interferogram.  Jenis filtering yang  digunakan  oleh GMTSAR adalah  metode Gaussian .  Terdapat    beberapa    jenis  filter yang bervariasi dari  100  m  hingga  700  m.  Dari  tahap  pemfilteran  ini, diperoleh file interferogram  sebelum dan setelah filtering . 6. Tahap Snaphu (Unwrap Processing) Citra  yang  dihasilkan  dari  proses interfere dan filtering masih  dalam  satuan radian (satuan phase sudut)  dalam rentang  -π  sampai  dengan  π.  Tahap snaphu dilakukan melalui  komen snaphu.csh . 7. Tahap Geocode
Tahap geocode dilakukan  melalui  komen geocode.csh . Dengan  kemampuan  software GMT ,  hasil  tersebut  dapat ditampilkan  dalam  bentuk phostscript dan k ml  Google Earth .  Walaupun  telah  bergeoreferensi,  grid,  skala,  dan informasi lainnya belum ada. Untuk itu, skill menggunakan GMT perlu  dilakukan  sehingga  diperoleh hasil yang  artistik  akan ditampilkan  di bagian hasil.

## III. HASIL DAN ANALISA

## A. Interferogram Phase dan Amplitude
Citra Single  Look   Complex ( SLC )  diolah  menjadi interferogram phase dan amplitude dengan  tahapan  mencari area  tampalan  ( offset )  dan  komputasi.  Akan  tetapi,  dengan menggunakan  perangkat  lunak GMTSAR5 ,  citra SLC yang dihasilkan tidak dapat ditampilkan.
Hasil  interferogram  yang  telah  dikurangkan  dengan  efek topografi yang berupa citra phase dan amplitude dalam sistem koordinat radar dapat dilihat pada Gambar  4, 5, dan 6.
Gambar 4. Interferogram Phase T ahun 2007-2010 (a), Interferogram Am plitude T ahun 2007-2010 (b)
Gambar 5.  Interferogram Phase T ahun 2008-2010 (c), Interferogram Am plitude T ahun 2008-2010 (d)
Gambar 6.  Interferogram Phase T ahun 2009-2010 (a), Interferogram Am plitude T ahun 2009-2010 (b)

## B. Koherensi Citra
Untuk  melihat  tingkat koherensi citra secara kualitatif dapat dilihat dari citra  koherensi pada Gambar  7.
Objek  daratan  pada  citra  memiliki  nilai  koherensi  rata-rata kurang  dari  0,5.  Rendahnya  nilai  koherensi  ini  karena sebagian  besar  area  studi  kasus  memiliki  kerapatan  vegetasi yang  tinggi.  Kerapatan  vegetasi  ini  akan  mempengaruhi  nilai koherensi  citra  karena  pergerakan  dan  perubahan  vegetasi akan  berpengaruh  pada  pantulan back scatter -nya.  Apabila
(a) (b) (c) (d) (e) (f)
citra  koherensi  memiliki  tingkat  koherensi  yang  tinggi,  maka citra akan menunjukan warna yang terang.
Gambar 7. Citra Koherensi T ahun 2007-2010 (a), Citra Koherensi T ahun 2008-2010 (b),  Citra Koherensi T ahun 2009-2010 (c)
C. Citra Interferogram Setelah Filtering
Filtering interferogram ini menggunakan gauss_alos_200m yang  merupakan  default  dari software  GMTSAR .  Filter  ini dipilih  karena  mampu  memberikan  hasil  yang  baik  dengan waktu relatif  cepat.
Pada  Gambar  8,  9,  dan  10  adalah  hasil  interferogram sebelum  dan  sesudah  melalui  proses filtering. Pada  hasil interferogram  sebelum  melalui  proses filtering ,    menunjukan adanya noise yang  relatif  besar  pada  area  yang  memiliki tingkat  koherensi  yang  tinggi.  Sedangkan  interferogram setelah filtering -nya  menunjukan  peningkatan  ketajaman area yang masih memiliki  koherensi yang tinggi.
Gambar 10. Interferogram Phase Sebelum Filtering T ahun 2009-2010 (e), Interferogram Phase Setelah Filtering T ahun 2009-2010 (f)
D. Citra Interferogram Setelah Unwrapping
Setelah  dilakukan unwrapping pada  interferogram,  maka seperti  yang  terlihat  pada  Gambar  11,  pola  deformasi  area sudah dapat diketahui meski masih dalam satuan phase . Untuk mengestimasi  besar  deformasi  dalam  satuan  metrik,  menurut Gabriel (1989) [7] harus  dilakukan  kalkulasi  dengan menggunakan  rumus displacement  of  earth’s  surface sepanjang sensor Line of Sight ( LOS ), yaitu:
∆∅𝑑𝑒𝑓𝑜 = 4𝜋𝛿R/λ (1)
Dimana λ adalah panjang gelombang SAR dan ∆∅𝑑𝑒𝑓𝑜 adalah nilai  perbedaan phase .

## E. Citra Interferogram Setelah Geocode
Tahap  pengolahan  terakhir  adalah  melakukan georeference terhadap citra . Hasilnya dapat dilihat pada Gambar  12.
(a) (b) (c) (a) (b)
Gambar 8. Interferogram Phase Sebelum Filtering T ahun 2007-2010 (a), Interferogram Phase Setelah Filtering T ahun 2007-2010 (b)
(c) (d)
Gambar 9. Interferogram Phase Sebelum Filtering T ahun 2008-2010 (c), Interferogram Phase Setelah Filtering T ahun 2008-2010 (d)
(e) (f) (a) (b) (c)
Gambar 11. Interferogram Phase Citra Setelah Unwrapping T ahun 2007-2010 (a), T ahun 2008-2010 (b), dan T ahun 2009-2010 (c)

## F. Keak urasian Estimasi Panjang Baseline
Panjang baseline merupakan  faktor  terpenting  dalam mengukur  besar surface deformation dengan InSAR . Menurut Ferreti  (2007) [8] ,  untuk  mengetahui  penyebaran  titik  orbit dari data pengamatan SAR ,  dapat  memilih baseline perpendicular pendek  (≤  150  m)  atau baseline temporal yang pendek.  Hal  ini  juga  merupakan  faktor  didapatkannya  nilai koherensi yang tinggi antara citra master dan slave . Berdasarkan  perhitungan  pada GMTSAR ,  estimasi  panjang baseline  perpendicular data SAR area  Lahendong  adalah 178,475270  m  untuk  pasangan  citra  tahun  2007  dan  2010; 183,664538  m  untuk  pasangan  citra  tahun  2008  dan  2010; serta 174,048545  m  untuk pasangan citra 2009  dan 2010.

## G. Nilai Koherensi
Telah  dijelaskan  pada pembahasan sebelumnya bahwa nilai koherensi suatu citra berada antara 0  - 1. Proses interferomety dapat  dikatakan  baik  dan  akurat  apabila  nilai  koherensi citranya  berada  diantara  0,5  –  1,0.  Apabila  nilainya  dibawah 0,5,  maka  hasil  citra  dari  proses interferometry masih memiliki  informasi  yang  berarti,  namun  citra  dengan  nilai koherensi  tersebut  juga  menampilkan  peningkatan level noise yang sebanding dengan semakin kecilnya nilai koherensi citra. Pada  Gambar  13  merupakan  grafik  nilai  koherensi  yang menunjukan bahwa besar koherensi dari citra koherensi sangat rendah (< 0,5  ).
Berdasarkan  gambar  grafik  koherensi  di  atas,  dapat disimpulkan  hasil  plotting  besar  koherensi  citra  rata-rata adalah  antara  0.05-0.15.  Rendahnya  nilai  koherensi  ini menyebabkan  banyaknya  piksel  yang  kosong  setelah  citra melalui  proses unwrap .

## H. Estimasi Besar Deformasi
Pada  Gambar  14  menunjukkan  hasil   akhir pengolahan citra dengan  metode DInSAR yang  dimana  telah  di- overlay -kan dengan nilai estimasi deformasi  pada setiap pikselnya.
(a) (b) (c) Gambar 12. Interferogram Phase Citra Setelah Geocoding T ahun 2007-2010 (a), T ahun 2008-2010 (b), dan T ahun 2009-2010 (c)
(a) (b) (c)
Gambar 13. Grafik Koherensi Citra T ahun 2007-2010 (a), T ahun 2008-2010 (b), dan T ahun 2009-2010 (c)
(a)
Kecilnya  nilai  koherensi  menjadi  salah  satu  penyebab banyaknya nilai-nilai  piksel  yang kosong.

## IV. KESIMPULAN DAN SARAN
Berdasarkan  hasil  analisa  deformasi  permukaan  tanah akibat  eksploitasi  panasbumi  di  PLTP  Lahendong,  maka didapatkan beberapa kesimpulan,  yaitu:
1. Dari  hasil  pengolahan  dengan  menggunakan  metode DInSAR didapatkan estimasi besar deformasi berdasarkan Line of Sight (LOS) dengan nilai 3 sampai 4 cm.
2. Berdasarkan  hasil unwrap pada  pengolahan  data, menunjukan  bahwa  proses  deformasi  sudah  terjadi  di area  eksploitasi  panasbumi  Lahendong  dan  itu  berupa penurunan muka tanah (land subsidence) .
Kemudian  berdasarkan  hasil  pengolahan  data  dan kesimpulan  yang  diperoleh,  maka  penulis  memberikan
beberapa  saran  yang  dapat  diberikan  untuk  memaksimalkan penelitian selanjutnya adalah sebagai berikut:
1. Untuk  mencegah  rendahnya  nilai  koherensi  pada  citra akibat  dekorelasi  temporal,  maka  diperlukan  data  yang cukup  banyak  dengan  jarak  temporal dan baseline yang berdekatan.  Untuk  jarak  temporal  dapat  dilakukan akuisisi  data  setiap  6  bulan  dengan jarak baseline ≤ 150
m.
2. Selain  tenik DInSAR ini  masih  memiliki  kekurangan berupa  dekorelasi  temporal,  atmosferik,  dan noise ,  maka untuk  menambah  keakuratan  dari  hasil pengolahan data, dapat  digunakan  teknik  terbaru  dalam  pengolahan  data SAR yaitu  teknik PS-InSAR  (Persistent  Scattere Interferometry  SAR) ,  seperti yang pernah dilakukan oleh
Eneva  (2012) [9] yang  meneliti  deformasi  permukaan  di Imperial  Valley
3. Untuk  penelitian  selanjutnya  dapat  diberikan  tambahan data  pengukuran  langsung  di  lapangan  panasbumi Lahendong,  berupa  data  GPS  serta  melakukan  uji validasi  untuk  mengetahui  keakuratan  nilai  deformasi hasil  pengolahan SAR dengan  menggunakan  teknik DInSAR .

## UCAPAN TERIMA KASIH
Penulis mengucapkan terima kasih kepada Ibu Ira Mutiara Anjasmara  yang  sudah  berkenan  memberikan  bimbingannya dalam  menyelesaikan  penelitian  ini  dan  kepada  rekan -rekan
Laboratorium
Geodinamika dan Lingkungan Teknik
Geomatika  ITS  yang  sudah  memberikan  masukan -masukan yang membangun terhadap penelitian ini.

## DAFTAR PUSTAKA
[1] Nenny.  2009. Sekilas  Tentang  Panas  Bum i .  Bandung:  Institut T eknologi Bandung.
[2] Prabowo,  T .  2011. Evaluasi  Reservoir  Lapangan  Geotherm al Lahendong-Sulawesi Utara Setelah Berproduksi Selam a 10 Tahun . Bandar  Lampung,  Indonesia:  Proceeding  of  Annual  Indonesian Geothermal Association Meeting & Conference
[3] Japan Space Systems. 2012. PALSAR User's Guide. Japan: Japan Space Systems (J-space systems).
[4] Sardjito. 2013. Pengembangan Panas Bumi Lahendong.
<http://www.dunia-energi.com/pge-lanjutkan-pengembangan-panas- bumi-lahendong/> dikunjungi tanggal 20 Desember 2015 pukul 11.34
[5] Brehme, M. 2015. The Role of Fault Zone on Structure, Operation and Prospect  of  Geotherm al  Reservoirs  (Case  Study:  Lahendong, Indonesia). Göttingen: Georg-August University School of Science. [6] Hanssen, R. (2001). Radar Interferom etry: Data Interpretation and Error Analysis. London: Kluwer Academic Publisher. The Netherlands: Delft University of T echnology.
[7] Gabriel,  A.  G.,  Goldstein,  R. M. & Zebker, H. A. (1989). Mapping sm all  elevation  changes  over  large  area s:  Differential  radar interferom etry . J. Geophys. Res. 94, 9183–9191.
[8] Ferretti, A., Andrea Monti-Guarnieri, Claudio Prati, dan Fabio Rocca. (2007a). InSAR  Principles  :  Guidelines  For  SAR  Interferom etry Processing And Interpretation . T he  Netherlands: ESA Publications. [8] Ferretti, A., Andrea Monti-Guarnieri, Claudio Prati, dan Fabio Rocca. (2007b). InSAR Processing : A Practical Approach . T he Netherlands: ESA Publications.
[9] Eneva, M.. (2012). Surface Deformation in Imperial Valley, CA, from Satellite  Radar  Interferom etry. Geothermal  Resource  Council T ransaction, Vol. 36.
Gambar 14. Estimasi Besar Deformasi Area PLTP Lahendong Tahun 2007-2010 (a), T ahun 2008-2010 (b), dan T ahun 2009-2010 (c)
(b) (c)