## IDENTIFIKASI DUGAAN SITUS PURBAKALA MENGGUNAKAN METODE GEOLISTRIK RESISTIVITAS KONFIGURASI WENNER DI SITUS ALASSUMUR, KABUPATEN BONDOWOSO
Juan Pandu Gya Nur Rochman*, Amien Widodo, Thufeil Amr Adausy dan Hanif Dhiyaz Ulhaq F
Departemen Teknik Geofisika, Fakultas Teknik Sipil, Perencanaan, dan Kebumian, Institut Teknologi Sepuluh Nopember, Surabaya, Indonesia
* juanpandu@geofisika.its.ac.id
Received 13-08-2021, Revised 05-03-2022, Accepted 16-03-2022 Available Online 17-03-2022, Published Regularly April 2022

## ABSTRACT
The discovery of ancient brick structures in Alas Sumur Village, Pujer District, Bondowoso Regency is interesting to study, because so far Bondowoso Regency is generally in the form of relics of the large stone age, namely Megalithic Sites. The purpose of this study was to identify the brick structure at the site. The method used is the Wenner configuration geoelectric method which can describe the subsurface conditions and the estimated structure of the bricks based on differences in rock resistivity values. The number of measurement tracks is 4 tracks, with a length of 1.3, and 4 tracks that are 48 meters long with a space of 1 meter, while track 2 is 96 meters long with a space of 2 meters. Data processing results in the form of cross-sectional data resistivity values. Based on the processing results, obtained lithology of breccia rocks, andesite rocks, and sandy clay. The distribution of the resistivity values of breccia has a range of values (34 – 101 meters), andesite rocks have a range of (102 – 141 meters), sandy clay rocks have a range (26 – 33 meters), anomalies are detected with a range of values (19 – 25Ωmeters) as bricks. . These anomalies are found on paths 1, 3, and 4. The discovery of these anomalies is expected to help related parties to carry out further development processes.
Keywords : alassumu; bricks; resistivity; wenner.

## ABSTRAK
Penemuan  struktur  batubata  kuno  di  Desa  Alas  Sumur,  Kecamatan  Pujer  Kabupaten  Bondowoso menjadi menarik untuk diteliti, karena selama ini Kabupaten Bondowoso umumnya berupa peninggalan zaman batu besar yaitu Situs Megalithikum. Tujuan penelitian ini adalah untuk identifikasi struktur batu bata di situs tersebut. Metode yang digunakan adalah metode geolistrik konfigurasi wenner yang dapat menggambarkan kondisi bawah permukaan  dan dugaan stuktur batu bata  berdasarkan perbedaan nilai resistivitas batuan . Jumlah Lintasan pengukuran sebanyak 4 lintasan, dengan panjang untuk lintasan 1,3, dan 4 sepanjang 48 meter dengan  spasi 1 meter sedangkan lintasan 2 sepanjang 96 meter dengan spasi 2 meter. Data hasil pengolahan berupa data penampang nilai resistivitas. Berdasarkan hasil pengolahan, didapatkan  litologi  batuan  breksi,  batuan  andesit,  dan  lempung  pasiran.  Persebaran  nilai  resistivitas batuan  breksi  memiliki  range  nilai  (34  –  101  Ωmeter),  batuan  andesit  memiliki  range  (102  –  141 Ωmeter), batuan lempung pasiran memiliki range ( 26 – 33 Ωmeter), anomali terdeteksi dengan range nilai (19 - 25Ωmeter) sebagai batubata. Anomali tersebut terdapat pada lintasan 1, 3, dan 4. Penemuan anomali  tersebut  diharapkan  dapat  membantu  pihak  terkait  untuk  melakukan  proses  pengembangan selanjutnya.
Kata kunci : alassumur; batubata; resistivitas; wenner.

## PENDAHULUAN
Kabupaten  Bondowoso  merupakan  daerah  yang  terkenal  akan  peninggalan  masa megalitikum [1] .  Hal  ini  dibuktikan  dengan  ditemukannya  pecahan-pecahan  benda peninggalan  berupa  sarkofagus,  kubur  bilik,  menhir,  dolmen  maupun  batu  dalkon  yang tersebar  di  beberapa  kecamatan [2] .  Baru-baru  ini  ditemukan  struktur  yang  diduga  situs peninggalan purbakala  di  Desa Alas Sumur, Kecamatan Pujer. Penemuan bermula setelah seorang  warga  menggali  sumur  untuk  mencari  air  bersih,  namun  sesuatu  janggal  dengan penemuan  batu  bata  kuno  dengan  arah  barat  laut  -  tenggara.Hal  ini  menjadi  penemuan menarik karena baru saat ini ditemukan jejak kesinambungan peradaban di Bondowoso yang tidak  hanya  berhenti  pada  zaman  megalitikum.  Oleh  karena  itu,  penelitian  lebih  lanjut mengenai  singkapan  struktur  batu  bata  ini  perlu  dilakukan,  salah  satunya  dengan menggunakan  metode  dalam  geofisika.  Metode  yang  umum  digunakan  adalah  metode geolistrik. Geolistrik adalah salah satu metode geofisika yang mempelajari sifat aliran listrik di dalam bumi dengan cara mendeteksinya di bawah permukaan bumi.
Prinsip kerja geolistrik yaitu dengan menginjeksikan arus listrik ke permukaan tanah melalui elektroda dan mengukur beda potensial dengan elektroda lain [3] . Metode ini menghasilkan nilai resistivitas bawah permukaan melalui pengukuran beda potensial yang ditimbulkan dari injeksi  arus  listrik  ke  dalam  bumi.  Berdasarkan  nilai  resistivitas  hasil  pengukuran,  dapat diketahui kondisi bawah permukaan bumi. Metode geolistrik telah terbukti efektif dan relatif murah digunakan untuk memetakan situs – situs arkeologi [4] , diantaranya Identifikasi Situs Maelang Banyuwangi Jawa Timur [5] , Situs Kadipaten Terung Sidoarjo [6,7] , dan Identifikasi Batuan Pondasi Candi Badut [8] . Penggunaan metode geolistrik ini diharapkan mampu untuk mendeteksi  struktur  bata  sehingga  mempermudah  dalam  perencaan  pelestarian  situs arkeologi ini.

## GEOLOGI REGIONAL
Lokasi  penelitian  termasuk  ke  dalam  formasi  batuan  gunungapi  raung  muda  (Qhvr) [9] . Litologi penyusunnya merupakan batuan breksi dan tuff. Batuan gunungapi tersusun atas tufa dan breksi gunung api batuapung yang berlapis tebal. Tufa berwarna coklat atau kemerahan, sebagian besar lapuk. Breksi  gunungapi  terdiri atas  komponen basal,  andesit dan obsidian pada matriks batupasir tufaan kasar, fragmen berukuran angular sub dan sortasi buruk. Pada Gambar 1 merupakan peta geologi daerah penelitian.
Gambar 1. Geologi Daerah Penelitian [9] .

## GEOLISTRIK RESISTIVITAS
Metode geolistrik pada penelitian ini digunakan untuk mengetahui susunan, kedalaman, dan penyebaran  lapisan  arkeologi  di  bawah  permukaan  berdasarkan  harga  tahanan  jenis  yang diperoleh.
Gambar 2. Konsep Geolistrik Resistivitas
Prinsip metode geolistrik ini dengan cara mengalirkan arus listrik searah (DC) ke dalam bumi melalui sepasang elektrode arus (C 1 dan C 2 ), yang kemudian diterima oleh sepasang elektrode potensial (P 1 dan P 2 ) [10] , seperti pada Gambar 2. Dari data pengukuran di lapangan, data yang dihasilkan merupakan data resistivitas semu dengan persamaan dibawah :
𝜌 𝛼 = 2𝜋 [( 1 𝑟1 − 1 𝑟2 )−( 1 𝑟3 − 1 𝑟4 )] ∆𝑉 𝐼 (1)
𝜌 𝛼 = 𝐾 ∆𝑉 𝐼 (2)
Dimana : 𝜌 𝛼 : Tahanan jenis semu (Ωmeter)
𝑟 1 , 𝑟 2, 𝑟 3, 𝑟 4 : Jarak antar elektroda arus dan potensial (meter) ∆𝑉 : Beda Potensial (mV) I : Kuat arus yang dialirkan (milliampere) K : Faktor geometri
Salah  satu  konfigurasi  yang  baik  untuk  studi  kasus  arkeologi  yang  dangkah  adalah konfigurasi Wenner karena  memiliki  resolusi  bagus  secara  lateral [3]  [11] . Pada  konfigurasi wenner jarak elektroda antara potensial dan arus (ABMN) memiliki jarak yang sama adalah a untuk layer 1, untuk layer 2 jaraknya 2a dan seterusnya. Untuk menghitung faktor geometri konfigurasi wenner R 1 = a rincian (A dan M), R 2 = 2a dengan rincian (AN), R 3 = dengan rincian (BM), dan R 4 = (BN) seperti pada Gambar 3.
Gambar 2. Konfigurasi Wenner
Faktor  geometri  (K)  konfigurasi Wenner berdasarkan  persamaan  1  dan  2  yaitu  sebagai berikut :
𝐾 = 2𝜋𝑎 (3)

## RESISTIVITAS BATUAN
Nilai  resistivitas  batuan  cukup  bervariatif.  Perbedaan  nilai  resistivitas  diakibatkan materialnya  yang  terkandung  di  dalam  batuan  tersebut,  antara  lain  :  densitas,  porositas, ukuran  dan  bentuk  pori-pori  batuan,  serta  seberapa  tahan  terhadap  kandungan  air.  Nilai resistivitas beberapa litogi batuan dapat dilihat pada Tabel 1.
Tabel 1. Nilai Resistivitas Batuan [12] No. Litologi Nilai Resistivitas (Ωmeter) 1. Lempung 10 – 200 2. Pasir 100 – 600 3. Alluvium 10 - 800 4. Breksi 100 - 500 5. Andesit 100 - 1000

## METODE
MULAI Studi Literatur dan Survei Lapangan Akuisisi Data Data Geolistrik Pengaturan Data Sesuai Format Res2dinv Input Data ke dalam Software Res2dinv Pengaturan Parameter Inversi RMS Error < 5% ? Penampang Resistivitas 2D Perubahan Parameter Inversi Interpretasi Data Penampang 2D Analisis SELESAI Ya Tidak Data Geologi Regional
Gambar 3. Diagram alir penelitian
Diagram  alir  penelitian  dapat  dilihat  pada  Gambar  4.  Penelitian  di  awali  dengan  mencari studi literature mengenai kondisi geologi dan bebrapa referensi terkait. Kemudian dilakukan survei pendahuluan guna menentukan penentuan lintasan yang sesuai. Selanjutnya dilakukan akuisisi  data  sesuai  berdasarkan  lintasan  yang  ditentukan . Data  yang  diperoleh,  akan dilakukan pengolahan dengan menggunakan inversi Smoothness-Constrained Least-Square. Setelah dilakukan pengolahan, langkah selanjutnya merupakan intepretasi data berdasarkan data pendukung yang ada.
Gambar 4. Desain akuisisi penelitian
Lokasi  penelitian  berada  di  desa  Alas  Sumur,  Kecamatan  Pujer,  Kabupaten  Bondowoso. Pengukuran dilaksanan pada tanggal 20 – 22 November 2020 selama 3 hari, dengan  total 4 lintasan. Lokasi penelitian berada di pemukiman warga yang padat dan ladang warga yang cukup lebat. Sehingga penentuan lokasi lintasan, panjang lintasan dan spasi masing – masing lintasan disesuaikan dengan kondisi real di lapangan (Gambar 5). Pada lintasan 1, 3 , dan 4 digunakan spasi pengukuran 1 meter, sedangkan pada lintasan 2 digunakan spasi 2 m (Tabel 2).  Spasi yang rapat ini dibutuhkan untuk mendapatkan resolusi hasil yang detail.
Tabel 2. Data lintasan akuisisi Lintasan Spasi Elektroda (m) Latitude Longitude Panjang Lintasan (m) Arah Lintasan 1 1 7°59'53.63"S  113°52'28.09"E 48 Timur - Barat 7°59'53.53"S  113°52'26.45"E Lintasan 2 2 7°59'53.60"S  113°52'30.14"E 96 Selatan - Utara 7°59'56.73"S  113°52'30.17"E Lintasan 3 1 7°59'55.19"S  113°52'26.28"E 48 Selatan - Utara 7°59'55.02"S  113°52'24.74"E Lintasan 4 1 7°59'52.31"S  113°52'26.73"E 48 Timur - Barat 7°59'53.85"S  113°52'26.66"E Sumur
Copyright © 2022 Universitas Sebelas Maret
Gambar 5. Alat ukur Geolistrik Multichannel MAE X612-EM
Alat  yang  digunakan  dalam  penelitian  ini  adalah  geolistrik multichannel MAE  X612-EM (Gambar 6). Dari pengukuran data yang diperoleh  antara lain beda pontensial (V), besar arus listrik (I). Hasil perhitungan  yang diperoleh dari akuisisi data  didapatkan resistivitas semu (ρ),  dikarenakan  bumi  (bawah  permukaan)  tidak  bersifat  homogen. Pengolahan  data dilakukan  dengan  menghitung  faktor  geometri  dari  konfigurasi  Wenner  untuk menghilangkan  pengaruh  letak  elektroda  potensial  terhadap  letak  kedua  elektroda  arus. Setelah  diperoleh  hasil  faktor  geometri  dari  konfigurasi  Wenner  kemudian  menghitung resistivitas  semu  dan  melakukan  proses  inversi  dengan  sotfware  Res2Dinv  untuk memperoleh penampang 2D.
Metode  inversi  yang  digunakan  adalah  metode Smoothness-Constrained  Least-Square. Metode ini dapat meminimalkan perbedaan antara data lapangan dan model prediksi yang dapat dilihat dari penampang 2D dan menghasilkan model dengan variasi nilai resistivitas yang smooth sehingga  cocok  untuk  studi  ketebalan  sedimen  atau  litologi  bawah permukaan [13] . Proses  iterasi  digunakan untuk  mendapatkan  error hingga  di  bawah 5% [14] . Dari  pengolahan  data,  Langkah  selanjutnya  merupakan  interpretasi  data.  Interpretasi  data dilakukan  dengan  menggunakan  pendekatan  parameter  nilai-nilai  resistivitas  material (benda)  dari  penampang  2D  hasil  Res2Dinv.  Penampang  2D  yang  didapatkan,  lalu diinterpretasikan dengan data pendukung berupa peta geologi regional lembar Besuki, peta geologi  Gunungapi  Raung,  stratigrafi  terukur  yang  dilakukan  di  dalam  sumur  gali,  dan penelitian yang terkait.

## HASIL DAN PEMBAHASAN
Penelitian ini menggunakan data stratigrafi terukur (sumur gali) sebagai data pendukung dan sebagai data acuan. Sumur terletak pada  koordinat 113°87’42.19”S, 7.99’81.45”E. Data ini diambil  dengan  cara  turun  ke  sumur  dengan  tali  seadanya  dan  dilakukan  dengan  metode meraba  dengan  mengambil  sampel  yang  bisa  diambil  dan  melihat  warna    tiap  perbedaan litologi. Pengamatan stratigrafi ini membuahkan  hasil bahwa dari kedalaman 0 meter dari
permukaan hingga 8 meter tersusun atas batuan breksi vulkanik dan lempung pasiran. Pada Gambar 7 ditunjukkan susunan stratigrafi terukur hasil pengamatan.
Gambar 7. (a). Hasil data stratigrafi terukur, (b). Sumur pengambilan data stratigrafi
Hasil dari stratigrafi terukur merupakan acuan untuk melakukan intepretasi lanjut dari data pengukuran geolistrik. Dengan adanya data tersebut diharapkan bisa sebagai data tambahan selain  berdasarkan  data  geologi  regional  lembar  Besuki  maupun  peta  geologi  gunungapi Raung. Selain itu terdapat data pendukung yaitu penelitian di daerah Grujugan, Bondowoso dengan persebaran nilai 10,1 Ωm – 32 Ωm sebagai Alluvium, lempung dan pasir. Rentang nilai  101  Ωm  –  321  Ωm  sebagai  andesit (2) .  Dari  data  pendukung  tersebut  dilakukan interpretasi  sebagai  referensi  nilai  resitivitas  dan  litologi  hasil  penampang.  Pada  saat intepretasi  data,  dari  tiap  penampang range resistivitas  disamakan  guna  mempermudah proses intepretasi. Proses penyamaan range resistivitas melihat persebaran tiap penampang. Berikut merupakan pembahasan dari penampang tiap lintasan.
Gambar 8. Penampang lintasan 1
Penampang lintasan 1 yang berarah timur-barat (Gambar 8) dengan nilai rms error sebesar 3%. Persebaran nilai resistivitas yang ditunjukan pada penampang lintasan 1 memiliki range 12.5 Ωm – 143 Ωm. Nilai resistivitas 34 Ωm – 101 Ωm diinterpretasi sebagai batuan breksi. Daerah  penelitian  termasuk  daerah  batuan  gunung  api  raung  muda  (Qhvr)  yang  terdapat
(a) (b)
batuan  breksi  dan  divalidasi  dengan  adanya  data  pengukuran  stratigrafi  terukur.  Nilai resistivitas  26  Ωm  –  33  Ωm  diinterpretasi  merupakan  batuan  lempung  pasiran.  Daerah penelitian merupakan lahar pujer berdasarkan pada data pendukung peta geologi gunungapi Raung  dengan  terdapat  adanya  pengendapan  yang  terbawa  oleh  aliran  lahar  gunung  api Raung. Anomali terdapat pada penampang dengan range nilai resistivitas 19 Ωm – 25 Ωm. Anomali  tersebut  dihimpit  oleh  litologi  lempung  pasiran,  dan  secara  desain  akuisisi memotong tepat pada meter ke 22 hingga 24 dan jelas terdapat anomali dengan kedalaman tersebut. Dari hasil interpretasi lintasan 1, maka dijadikan sebagai acuan untuk interpretasi lintasan  selanjutnya.  Nilai  resistivitas  batu  bata  ini  yang  menjadi  kunci  dalam  identifikasi lintasan  lainnya.  Pada  lingkaran  kuning  pada  hasil  penampang  merupakan  bukan  anomali dikarenakan pada elektroda 32 dan 40 terdapat paving sehingga menggunakan elektroda plat dan hasilnya nilai resistivitas sangat rendah, sehingga didefinisikan bukan zona interest yang dicari.

## Gambar 9. Penampang lintasan 2
Pada hasil pengolahan lintasan 2 (Gambar 9) memiliki rms yang paling rendah dibandingkan dengan lintasan yang lain yaitu sebesar 2.5%. Lintasan ini memiliki orientasi arah selatan – utara, dengan rentang persebaran nilai resistivitas 34 Ωm – 143 Ωm. Intepretasi pada lintasan ini memiliki 2 jenis litologi yaitu breksi dengan nilai resistivitas 34 Ωm – 101 Ωm dan andesit dengan nilai resistivitas 102 Ωm  – 141 Ωm. Pada lintasan kedua tidak ditemukan adanya anomali dugaan batu bata, dikarenakan persebaran nilai resistivitas cukup tinggi. Diduga kuat karena  tumbuhan  yang  cukup  lebat  pada  daerah  penelitian  dan  juga  dijustifikasi  bahwa lintasan  ini  tidak  terdapat  adanya  pengendapan  yang  terjadi  pada  stratigrafi  endapan sekunder. Sehingga, litologi  yang ada pada penampang nilai  resistivitas  tinggi  merupakan batuan yang memiliki matriks fragmen kompak.
Gambar 10. Penampang lintasan 3
Intepretasi pada lintasan ketiga ini terdapat 2 litologi yang mendominasi yaitu breksi dengan nilai resistivitas 34 Ωm – 101 Ωm dan lempung pasiran dengan nilai resistivitas 26 Ωm – 33 Ωm. Hal ini diduga pada saat proses pengendapan, tidak mengalami proses yang merata hal ini  dibuktikan  dengan  penampang  pada  Gambar  10.  Anomali  terdapat  pada  penampang dengan range nilai  resistivitas 19 Ωm  – 25 Ωm.  Lingkaran  warna  kuning bukan anomali, dikarenakan  fakta  di  lapangan  menunjukkan  adanya  selokan  yang  tergenang  air  pada elektroda  meter  30  –  32,  nilai  resistivitasnya  sangat  rendah,  sehingga  didefinisikan  bukan interest yang  dicari.  Pada  pengolahan  lintasan  3  memiliki rms  error sebesar  2.8%.  Pada lintasan ini memiliki orientasi lintasannya  berarah timur – barat.
Gambar 11. Penampang lintasan 4
Persebaran  nilai  resistivitas  yang  ditunjukan  pada  penampang  lintasan  4  (Gambar  11) memiliki range 12.5 Ωm – 143 Ωm. Interpretasi nilai resistivitas 34 Ωm – 101 Ωm sebagai breksi, nilai resistivitas 26 Ωm  – 33 Ωm sebagai  lempung pasiran. Anomali terdapat pada penampang dengan terdapat pada range nilai resistivitas 19 Ωm – 25 Ωm. Berdasarkan hasil inversi anomali ditemukan pada tengah lintasan yang sesuai pada desain akuisisi memotong tepat pada dugaan arah kemenerusan tenggara-barat laut. Hal ini sudah sesuai dengan lintasan pertama  yang  merupakan  validator.  Faktor  cuaca  juga  mempengaruhi  hasil  inversi  yang dilakukan, dengan hasil data akhir membuat litologi berubah menjadi nilai resistivitas rendah pada  dari  elektroda  32  hingga  48.  Hal  itu  membuat  faktor  pembacaan  inversi  menjadikan penampang memiliki nilai resistivitas yang rendah. Pada lintasan 4 memiliki arah lintasan selatan – utara.
Dari  hasil  pengolahan  dan  interpretasi  pada  keempat  lintasan,  anomali  ditemukan  pada litologi  lempung  pasiran.  Anomali  terkubur,  diduga  berasal  dari  adanya  pengendapan  di daerah  penelitian  akibat  aliran  lahar  gunung  api  Raung.  Berdasarkan  fasies  gunung  api, daerah  penelitian  termasuk  ke  dalam  fasies  distal.  Dibuktikan  fasies  distal  sebagai  daerah pengendapan terjauh dari sumber dan didominasi oleh endapan rombakan gunung api seperti halnya  breksi  lahar,  breksi  fluvial,  konglomerat,  batupasir,  dan  batulanau [15] .  Hal  tersebut diperkuat dengan adanya kali Sampean yang dekat dengan daerah penelitian merupakan jalur dari aliran lahar pujer.

## KESIMPULAN
Berdasarkan hasil penelitian yang telah dilakukan, maka dapat disimpulkan bahwa litologi penyusun di daerah penelitian berupa batuan Breksi, Lempung pasiran, dan Andesit. Masing masing memiliki range yang berbeda dengan nilai 26 – 33 Ωm sebagai lempung pasiran, 34 – 101 Ωm sebagai breksi, 102 – 143 Ωm sebagai andesit. Anomali yang berasosiasi dengan batu bata mempunyai rentang nilai resistivitas 19 – 25 Ωm , yaitu terindikasi pada lintasan 1 dengan  kedalaman  3  –  5  meter,  lintasan  3  dengan  kedalaman  3  –  5  meter,  dan  lintasan  4 terdapat  2  posisi  masing-masing  dengan  kedalaman  2,5  –  4  meter.  Hasil  tersebut menunjukkan  kesesuaian  dengan  stratigrafi  daerah  penelitian  yang  divalidasi  melalui pengukuran  stratigrafi  secara  langsung.  Penemuan  struktur  anomali  tersebut  dapat dimanfaatkan  oleh  dinas  terkait  untuk  dikaji  lebih  lanjut  atau  dilakukan  pengembangan proses selanjutnya.

## UCAPAN TERIMA KASIH
Terima  kasih  kepada  DPRD  Kabupaten  Bondowoso,  Dinas  Pendidikan  Kabupaten Bondowoso  dan  Balai  Pelestarian  Cagar  Budaya  (BPCB)  Jawa  Timur  atas  kerjasama  nya dalam penelitian di Situs Alas Sumur Kabupatern Bondowoso ini.

## DAFTAR PUSTAKA
1. Amalia,  R.,  Megalitik,  S.,  Kecamatan,  P.,  Amalia,  R.,  &  Swastika,  K.  2013. Situs Megalitik di Desa Penanggungan Kecamatan Maesan Kabupaten Bondowoso sebagai Sumber  Pembelajaran  Sejarah  (  Megalitic  Site  at  Penanggungan  Village  Maesan Subdistrict Bondowoso Regency as History Learning Source ) .
2. Karisma,  U. 2013. Pola  Distribusi  Resistivitas  Bawah  Permukaan  Situs  Megalitikum Dengan  Metode  Geolistrik  RES3D  di  Kecamatan  Grujugan  Kabupaten  Bondowoso. Skripsi , FMIPA Universitas Jember, Jember.
3. Loke, M. H. 1999. A Practical Guide to 2D and 3D Surveys . Electrical Imaging Surveys for Environmental and Engineering Studies, 8-10.
4. Isdarmadi,  K.,  Susilo,  A.,  Sunaryo,  S.  2013.  Estimation  of  Structures  Under  Ancient Archaeological  Sites  Site  of  Jabung  Temple,  Probolinggo,  Indonesia  Using  Geolistrik Resistivity Method. Nat. B J. Health Environ. Sci. , 2, 50–57.
5. Fauzan,  M.,  &  Purwanto,  M.  S.  2017.  Penerapan  Metode  Resistivitas  2D  untuk Identifikasi  Bawah  Permukaan  Situs  Maelang  Banyuwangi  Jawa  Timur. Jurnal  Sains Dan Seni ITS , 6 (2).
6. Widodo, A., Syaifuddin, F., Vinca, V., Warnana, D. D., Rochman, J. P. G. N., Ariyanti, N., & Lestari, W. 2019. Data Acquisition of 2D Geophysical Resistivity Methods with Dipole-Dipole  Configuration  for  Identification  the  Subsurface  Brick  Stone  Sites  of Kadipaten Terung Sidoarjo. IOP Conference Series: Materials Science and Engineering , 546 (2).
7. Widodo,  A.,  Syaifuddin,  F.,  Mudhofar,  A.,  Warnana,  D.  D.,  Rochman,  J.  P.  G.  N., Ariyanti, N., & Lestari, W. 2019. Identification the subsurface structures of kadipaten terung  site  using  surface  3d  resistivity  methods. IOP  Conference  Series:  Materials Science and Engineering , 546 (3).
8. Luthfin, A., Cahyadi, H., & Jufri. 2020. Identifikasi Batuan Pondasi Candi ( Andesit ) di Bawah  Permukaan  Sekitar  Candi  Badut  dengan  Metode  Geolistrik  Resistivitas. Indonesian Journal of Applied Physics , 10 (2), 106–115.
9. Pendowo, B., Samodra, H. 1997 . Peta Geologi Lembar Besuki Jawa. Pusat Penelitian dan Pengembangan Geologi . Indonesia
10. Reynolds,  J.M.  1997. An  Intruduction  to  Applied  and  Enviromental  Geophysics . Chichester: Jhon Wiley & Sons Ltd.
11. Rochman,  J.  P.  G.  N.,  Widodo,  A.,  Bahri,  A.  S.,  Syamsuddin,  F.,  &  Wien,  L.  2017. Aplikasi  Metode  Geolistrik  Tahanan  Jenis  untuk  Mengetahui  Bawah  Permukaan. Geosaintek , 1, 93–98.
12. Telford.  W.  M,  Geldart.  L.  P,  &  Sheriff.  R.  E.  1990. Applied  Geophysics.  Applied Geophysics. Second Edi . New York: Cambridge and Hall.
13. Mia Azhari, Akmam, M. 2016. Analisis Jenis Batuan Menggunakan Metode Geolistrik Tahanan  Jenis  Konfigurasi  Wenner  di  Bukit  Apit  Puhun  Kecamatan  Guguk  Panjang Kota Bukittinggi. Pillar Of Physics , 34(1), 1–4.
14. Loke, M. H. 2014. RES2DINV Manual: Rapid 2-D Resistivity & IP inversion using the least-squares  method:  Geoelectrical  Imaging  2-D  &  3D  GEOTOMO  SOFTWARE . www.geoelectrical.co.
15. Bronto,  S.  (2006).  Fasies  gunung  api  dan  aplikasinya. Indonesian  Journal  on Geoscience , 1(2), 59–71.