Pemetaan Zona Kerentanan Air Tanah di Pulau Bakalan, Kabupaten Banggai Kepulauan, Indonesia .............................. (Febriarta, et al)

## PEMETAAN ZONA KERENTANAN AIR TANAH DI PULAU BAKALAN, KABUPATEN BANGGAI KEPULAUAN, INDONESIA
(Mapping of Groundwater Vulnerability Zones in Bakalan Island, Banggai Kepulauan, Indonesia)
Erik Febriarta 1 , Ahmad Cahyadi 2 , dan Septian Vienastra 3
1 Magister Pengelolaan Pesisir dan Daerah Aliran Sungai, Fakultas Geografi, Universitas Gadjah Mada
2 Departemen Geografi Lingkungan, Fakultas Geografi, Universitas Gadjah Mada
3 Teknik Geologi, Fakultas Teknologi Mineral, Institut Sains dan Teknologi AKPRIND Yogyakarta Jl. Kaliurang Sekip Utara, Bulaksumur, Sleman, Daerah Istimewa Yogyakarta E-mail: e.febriarta@gmail.com
Diterima: 22 Februari 2021; Direvisi:5 April 2021; Disetujui untuk Dipublikasikan: 19 April 2021

## ABSTRAK
Keterdapatan  sumber  daya  airtanah  di  kepulauan  sangat  tergantung  oleh  faktor  meteorologi  (curah hujan)  sebagai  sumber  air  dan  faktor  geologi  yaitu  formasi  penyusun  batuan  tempat  cadangan  airtanah tersimpan  (akuifer).  Pulau  Bakalan  berada  di  Kepulauan  Banggai,  yang  secara  kondisi  geologi  regional berada  di  Formasi  Batugamping  koral  kuarter  (Ql).  Karakteristik  airtanah  di  daerah  batugamping  sangat dinamis,  terbatas,  dipengaruhi  oleh  perkembangan  dari  aliran  celah  dan  pelarutan.  Sehingga  dapat dikatagorikan  daerah  airtanah  langka,  dan  bersifat  setempat  atau  lokal.  Keterbatasan  airtanah  juga dipengaruhi  oleh  aktifitas  pengunaannya.  Semakin  besar  penggunaan  dan  semakin  intensif  pemakaian airtanah,  maka  dapat  mengakibatkan  penurunan  kuantitas  dan  penurunan  kualitas  airtanah.  Penurunan kualitas  airtanah  dapat  menimbulkan  permasalahan  kesehatan.  Penurunan  kualitas  air  dapat  disebabkan tercemarnya sistem akuifer oleh sumber pencemar melalui pori-pori batuan pada formasi batugamping yang mudah  terlarut.  Untuk  mengetahui  distribusi  kerentanan  airtanah,  maka  tujuan  dari  penelitian  ini  adalah menentukan  sebaran  zona  kerentanan  airtanah  pada  batugamping  terumbu.  Untuk  mencapai  tujuan tersebut dilakukan analisis spasial multi-kriteria dengan pendekatan DRASTIC. Pendekatan DRASTIC dapat memberikan  infromasi  zona  kerentanan  berdasarkan  sifat  hidrogeologi  pada  batugamping.  Analisis pemberian  indekss  pembobotan  parameter  dengan  nilai  tinggi,  mengambarkan  informasi  faktor  yang dominan  atau  yang  paling  mempengaruhi  zona  kerentanan  airtanah.  Hasil  skor  perhitungan  setiap parameter  diperoleh  dari  perkalian  linier.  Indekss  kerentanan  air  tersebut  dikelaskan  ke  dalam  lima  tingkat kerentanan dengan pendekatan equal interval. Hasil analisis menunjukkan kerentanan sangat rendah seluas 0,003%,  kerentanan  rendah  seluas  18,79%,  kerentanan  sedang  seluas  18,48%,  kerentanan  tinggi  seluas 13,53% dan kerentanan sangat tinggi seluas 49,16% dari luas pulau 40,7 km 2 .
Kata kunci : airtanah, Banggai kepulauan, batugamping koral, DRASTIC, kerentanan

## ABSTRACT
Groundwater availability in a group of islands is highly dependent on meteorological factors (rainfall) as water  sources  and  geological  factors  (rock  formations  where  groundwater  reserves  are  stored,  termed aquifers). Bakalan Island, part of the Banggai Kepulauan, geologically lies on the Quaternary coral limestone formation (Ql). Groundwater in limestone areas is very dynamic, limited in quantity, and susceptible to fissure flow  development and dissolution. For these reasons, the  island is categorized as an area  with scarce and local groundwater. The prevailing activities are also influenced because more extensive and intensive use of groundwater can decrease its quantity and quality. Groundwater  quality  decline  is attributable to an  aquifer system polluted by contaminants that flow through rock pores in highly dissoluble limestone formation and, thus,  can  cause  health  problems.  This  research  set  out  to  determine  the  distribution  of  groundwater vulnerability zones in coral limestones using a multi-criteria spatial analysis with the DRASTIC approach. The DRASTIC approach produces information on vulnerability zones based on the hydrogeological properties of limestone. A parameter with a high weight index  reflects a dominant or  the most influencing factor in these zones and its score was obtained from linear multiplication. The groundwater vulnerability index values were differentiated into five levels using the equal interval approach. The results showed that the island area, 40.7 km 2 , was 0.003% very low vulnerability, 18.79% low vulnerability, 18.48% medium vulnerability, 13.53% high vulnerability, and 49.16% very high vulnerability.
Keywords: Banggai kepulauan, coral limestone, DRASTIC, groundwater, vulnerability

## PENDAHULUAN
Sumber  daya  air  di  pulau-pulau  kecil  relatif terbatas,  hal  tersebut  dipengaruhi  oleh  siklus hidrologi  yang  pendek  dan  daerah  tangkapan  air yang  relatif  kecil,  kemudian  mengalir  menuju kembali  ke  laut.  Faktor  meteorologis  atau  kondisi intensitas dan curah hujan mempengaruhi imbuhan  air  yang  masuk  kedalam  tanah  menjadi cadangan  airtanah  (Febriarta  et  al.,  2018).  Faktor geologis  dipengaruhi  oleh  kondisi  penysusun batuan  dari  sifat  keterdapatan  airtanah,  yaitu  sifat batuan dalam menahan (menyimpan) dan melolosakan  airtanah  (akuifer)  (Singhal  &  Gupta, 2010; Todd & Mays, 2005). Kemampuan air dalam menyimpan  air  sebagai  cadangan  airtanah  dapat diketahui  dengan  nilai  potensi  porositas  efektif ( specific  yield ).  Nilai specific  yield (sy)  atau porositas  efektif,  merupakan  perbandingan  dalam persen  (%)  air  yang  dapat  diambil  dari  tanah  atau batuan  yang  jenuh  air  dibandingkan  dengan volume  total  batuan  atau  tanah  (Todd  &  Mays, 2005).  Semakin  besar  nilai  porositas  efektif  (sy) maka  semakin  besar  volume  cadangan  air  yang dapat tersimpan dan dapat diloloskan oleh formasi batuan,  dan  sebaliknya  bila  semakin  kecil,  maka nilai porositas efektif (sy)  maka akan memiliki nilai rendah hingga langka (Singhal & Gupta, 2010).
Karakteristik  airtanah  dengan  parameter porositas  efektif  (sy)  tinggi  memiliki  korelasi dengan  kondisi  keterusan  airtanah  dan  potensi debit  air  yang  tinggi  juga  (Febriarta  &  Larasati, 2020). Potensi airtanah pada kawasan batugamping  memiliki  sifat  keterusan  yang  kecil hingga sedang, dengan material penyusun batuannya  berupa  celah  atau  rekah,  dengan  sifat keterusan  (melolosakan  airtanah)  dengan  potensi sedang  hingga  tinggi  bila  sudah  berkembang  dan berada pada media aliran  berupa saluran maupun pelarutan  hingga  terbentuknya  saluran-saluran (sungai  bawah  tanah)  (Singhal  &  Gupta,  2010). Pada  perkembangan  tersebut  sangat  dipengaruhi oleh  daerah  imbuhan  airtanah  yang  luas  dan sumber air (curah hujan) dengan intensitas kondisi curah  hujan  yang  tinggi  (Febriarta  &  Oktama, 2020; Purnama et al., 2019).

## Geologi Regional
Secara  geologi  regional  Pulau  Bakalan berada pada Formasi  Terumbu Koral Kuarter (Ql). Penyusun  batuan  formasi  tersebut  antara  lain batugamping  terumbu  dan  batugamping  koral bersisipan napal (KESDM, 2012). Secara stratigrafi  pengendapan  umur  batuan  tersebut Holosen. Kenampakan morfologi perbukitan batugamping dengan litologi batugamping terumbu pada  bagian  utara  dari  Pulau  Bakalan  secara karakteristik keterdapatan airtanah tidak terdistrubusi merata. Dilihat dari stratigrafi
pengendapan,  di  daerah  tersebut  merupakan bagian  dari  endapan  permukaan  dengan  potensi pelapukan  yang  tinggi.  Batuan  penysun  pada bagian selatan dengan morfologi dataran tersusun atas  batugamping bioklastik  dan  rombakan batugamping  koral  ( Gambar  2 ).  Morfolgi  pada batuan  barat  terdapat  perbukitan  karst  dengan penyusun batuan batugamping terumbu perselingan napal dan batugamping koral fragmental  hingga  tengah  pulau,  sedangkan bagian  timur  pulau  memiliki  morfologi  datar, dengan  kemiringan  lereng  landau  hingga  datar, dengan  penyusun  batuan  batugamping  koral ( Gambar  3 ).  Batugamping  bioklastik  terbentuk pada  lingkungan  pengendapan  lingkungan  laut yang  dangkal  dan  dekat  dengan  pantai  dengan partikel  yang  terombak  oleh  proses  abrasi. Endapan  tersebut  dapat  terbentuk  dari  akibat gelombang yang terabrasi, sedangkan perkembangan  terumbu  merupakan  hasil  dari pelongsoran  koral.  (BouDagher-Fadel  &  Marcelle, 2013).  Batugamping  bioklastik  (koral  bioklastik, fragmental)  yang  terdapat  pada  bagian  selatan pulau ( Gambar 1 ), tersusun atas material fargmen (material  rombakan)  kerangka  organisme  dengan sifat  lepas-lepas  akibat  proses  tertrasportasi. Material  lepas-lepas  memiliki  media  lubang  antar butir  yang  dengan  persentase  tingi  sehingga memiliki  nilai  percepatan  aliran  airtanah  yang cepat  (Todd  &  Mays,  2005;  Wicaksono  et  al., 2020).  Karakteristik  penyusun  material  tersebut memiliki  nilai  percepatan  pengaliran  yang  tinggi dan  dapat  mempercepat  pencemaran  terhadap sumber  pencemar  (polutan)  ke  dalam  sistem akuifer (Aller et al., 1987a; Gunawan et al, 2012).
Berdasarkan  peta  Kawasan  Bentang  Alam Karst (KBAK), Pulau Bakalan mempunyai kawasan karst  seluas  42,85%  ( Gambar  1 , Gambar  2 ,  dan Gambar  3 )  dari  luas  wilayah  Pulau  Bakalan  (40,7 km 2 )  yang  merupakan  bagian  dari  KBAK  Banggai Kepulaun  (RI,  2010).  Definisi  kawasan  karst menurut  Permen  ESDM  nomor  17  (2012),  adalah bentang  alam  yang  terbentuk  akibat  pelarutan  air pada  batugamping  dan/atau  dolomit  dengan perkembangan  ciri-ciri  eksokarst  dan  endokarst tertentu  (RI,  2010;  Singhal  &  Gupta,  2010). Kawasan  karst  tersebut  berada  di  bagian  utara dengan  ketinggian  tempat  0-78  mdpal  ( Gambar  2 dan Gambar  3 ).  Keterdapatan  airtanah  pada formasi  batugamping  sangat  dinamis  yang  sangat dipengaruhi  oleh  pelarutan  batuan,  dampak  dari sifat  dinamis  tersebut  antara  lain;  sifat  aliran  alir tanah  bersifat  lokal,  keterdapatan  airtanah  dapat berpindah  diakibatkan  oleh  pelarutan  batuan  yang berada  jauh  semakin  kedalam  permukaan, hingga potensi  pencemaran  yang  relatif  tinggi  karena langsung  terhubung  dengan  cadangan  sumber  air (Febriarta  et  al,  2020;  Fetter,  2014;  Kalhor  et  al., 2019).
Gambar 1 . Geologi regional batugamping terumbu dan kondisi hidrogeologi Pulau Bakalan.
Gambar 2 . Topografi utara-selatan dengan morfologi bergelombang.

## Hidrogeologi Regional
Pulau  Bakalan  merupakan  bagian  dari gugusan kepulauan di kabupaten Banggai Kepulauan,  Provinsi  Sulawesi  Tengah.  Secara pembentukannya  Pulau  Bakalan  terbentuk  oleh proses  teknogenesis  pengangkatan  yang  dicirikan dengan  struktur  utama  yaitu  sesar  dan  lipatan (KESDM,  2015b;  Surono  &  Sukido,  1993). Sehingga  potensi  airtanah  sangat  dikontrol  oleh sifat  batuan  dan  morfologi  lereng,  dengan  potensi keterdapatan  sumber  daya  air  sedang  (terbatas) pada  pegunungan  karst  batugamping  terumbu (BPLH,  2014).  Karakteristik  hidrogeologi  pada stratigrafi  kuarter  dengan  penyususn  batuan Formasi  Batugamping  Terumbu  (Ql)  berupa akuifer dengan aliran celah, lubang hasil pelarutan bukan  termasuk  akuifer  dan  tidak  termasuk  dalam kawasan  Cekungan  Airtanah  (CAT)  (BPLH,  2013, 2014; KESDM, 2018).
Berdasakan  sifat  hidrogeologi  maka,  Pulau Bakalan  dapat  dikategorikan  mempunyai  sumber daya air yang sangat terbatas hingga langka (tidak terdapat airtanah). Berdasarkan kondisi hidrogeologi  terdapat  tiga  potensi  akuifer  yaitu, zona dengan produktivitas sedang, dengan akiufer media ruang rekah atau saluran dengan keterusan sedang,  penyusun  batuan  berupa  batugamping terumbu  (Febriarta  et  al,  2020;  KESDM,  2015b). Zona  produktivitas  rendah  hingga  sedang  pada akuifer dengan media aliran melaui ruang, rekahan dan  saluran,  zona  tersebut  berada  di  bagian tengah,  dan  sebagian  berada  di  kawasan  Karst Banggai  Kepulauan.  Zona  produktivitas  kecil dengan  airtanah  langka  pada  penyusun  batuan batugamping koral, zona tersebut berada di bagian utara  Pulau  Bakalan  dan  masuk  kawasan  karst, yaitu  Karst  Banggai  Kepulauan.  Distribusi  kondisi hidrogeologi  di  Pulau  Bakalan  disajikan  pada Gambar 1 .

## Kerentanan Airtanah
Dinamika airtanah di pesisir sangat dipengruhi oleh kondisi formasi batuan dengan sifat parameter  akuifer  dan  pengaruh  air  laut  di  daerah sepandan  pantai  atau  dataran  pantai  (Cahyadi, 2019;  Febriarta  &  Widyastuti,  2020).  Pulau  koral memiliki karakteristik morfologi  yang landai hingga datar,  sehingga  imbuhan  air  yang  meresap kedalam  tanah  relatf  sedikit  (Purnama  &  Cahyadi, 2019;  Trinanda,  2017).  Keterdapatan  airtanah berkembang  pada  daerah  batugamping  terumbu atau  akuifer  karst.  Akuifer  karst  adalah  akuifer yang  memiliki  batuan  induk  yang  signifikan kelarutan  dalam  air  (White  et  al.,  2017).  Potensi perkembangan rongga, celah mapun saluran pada batuan  yang  terbentuk  oleh  pelarutan  dapat mengakibatkan sumber polutan (sumber pencemar)  masuk  kedalam  sistem  airtanah (Ferreira  et  al.,  2004;  White  et  al.,  2017).  Sifat akuifer  karst  sangat  permeabel  dan  rentan
terhadap  kontaminasi  (Linggasari  et  al.,  2019; Thapa et al., 2018; Vrba & Zaporozec, 1994).
Suatu  potensi  atau  kondisi  dimana  sistem airtanah dapat bertahan terhadap sumber pencemar dari permukaan hingga kedalam akuifer dapat  didefiniskan  sebagai  kerentanan  airtanah (Harter  &  Walker,  2001;  Putranto  et  al.,  2016). Percepatan pencemaran airtanah dapat dipengaruhi  oleh  kondisi  hidrogeologi  dari  formasi batuan  dalam  sifat kelulusan  airtanah  (Febriarta  & Oktama  et  al,  2020;  Voutchkova  et  al.,  2021). Sumber pencemar secara dinamis dapat dipengaruhi  oleh  sumber  pencemar  serta  kondisi infiltrasi hingga mencapai zona jenuh air, sehingga dapat  mengakibatkan  risiko  pencemaran  rendah pada  daerah  dengan  sumber  pencemar  (polutan) yang tinggi (DPU, 1996; Putranto et al., 2020; Vrba & Zaporozec, 1994).
Lathamani  et al.,  (2015)  menyebutkan  bahwa penilaian kerentanan airtanah dapat diperoleh dari kombinasi  parameter  sifat  pembawa  air  (kondisi geologi). Putranto et al., (2020) dan Voutchkova et al., (2021) menyatakan bahwa parameter hidrogeologi merupakan parameter yang mempengaruhi percepatan penyebaran pencemaran.  Oleh  karena  hal  tersebut  penilaian kerentanan  airtanah  terhadap  pencemaran  dapat diperoleh  dari  perhitungan  multi  kriteria  dari kecepatan  kelulusan  airtanah,  dengan  metode DRASTIC  (Aller  et  al.,  1987a;  Sugianti  et  al, 2016).Berkaitan dengan peningkatan akses informasi  dalam  pembangunan  berkelanjutan, maka  diperlukan  pemetaan  tentang  sumber  daya air  secara  keruangan  (spasial)  mengenai  potensi, lokasi, sebaran airtanah yang dangkal, dan kondisi kerentanan yang dipengaruhi oleh faktor penysusn batuan  yaitu  dominasi  batugamping  bioklastik fragmental  dan  terumbu,  yang  rentan  terhadap sumber  pencemar  (polutan)  di  Pulau  Bakalan. Maka  tujuan  dari  penelitian  ini  adalah  mengetahui potensi  tingkat  kerentanan  airtanah  terhadap pencemaran  pada  formasi  batugamping  koral kuarter (Ql), dengan pendekatan metode DRASTIC di Pulau Bakalan di Kabupaten Banggai Kepulauan.

## METODE
Dalam perhitungan pendekatan zona kerentanan  airtanah  dengan  DRASTIC,  terdapat beberapa  tahapan  yang    perlu    dilakukan    seperti yang    tertera    pada  diagram  alir  penelitian ( Gambar  4 ).  Pendekatan  kerentanan  dengan metode  DRASTIC  merupakan  sebuah  evaluasi untuk  airtanah  terhadap  potensi  pencemaran dengan  memperhatikan  sifat  hidrogeologi  (Aller  et al.,  1987a;  Ferreira  et  al.,  2004;  Harter  &  Walker, 2001;  White  et  al.,  2017).  Sifat  hidrogeologi merupakan  parameter  akuifer  dalam  menyimpan dan  meloloskan  (mengalirkan)  yang  berhubungan dengan waktu dan kecepatan aliran dan ketebalan terhadap  zona  tak  jenuh    (Fetter,  2014;  Singhal  & Gupta, 2010).
Gambar 4 . Langkah kerja kajian kerentanan airtanah terhadap pencemaran
Parameter yang digunakan dalam pendekatan kerentanan  airtanah  terhadap  pencemaran  antara lain: depth to water (D)/kedalaman muka airtanah, recharge (R)/imbuhan  air ,  aquifer (A)/jenis  akufer , soil media (S)/tekstur tanah , topography (T)/kemiringan  lereng ,  impact  of  vadose  zone media (I)/jenis  zona  tak  jenuh,  dan conductivity hydraulic (C)/konduktivitas hidraulis akuifer ( Tabel  1 ). Depth  to  water (D)/kedalaman  muka airtanah  merupakan  parameter  kedalaman  muka airtanah  dari  permukaan.  Semakin  dangkal  muka airtanah,  memiliki  potensi  pencemaran  yang  lebih tinggi  atau  mudah  mengalirkan  pencemaran.  Nilai Recharge (R) / imbuhan air adalah kondisi airtanah yang  masuk  kedalam  tanah  (Purnama  et  al., 2019).  Untuk  mendapatkan  nilai  imbuhan  airtanah diperoleh  dari  pendekatan  rumus  keseimbangan air ( water balance ) berdasarkan metode Thornthwaite-Matter (BSN,  2002),  seperti  pada Persamaan 1 .
Recharge = Curah hujan (presipitasi)- aliran permukaan (runoff) -  evapotranspirasi........... ......... (1)
Aquifer  (A)  /  jenis  akuifer  merupakan pendekatan  yang  dipengaruhi  oleh  sifat  pembawa airtanah  dengan  rongga  antar  butir  atau  media aliran  air.  Akuifer  dengan  aliran  media  antar  butir yang semakin rapat atau semakin kecil mempunyai nilai kerentanan yang rendah, sedangkan  kondisi  akuifer  dengan  media  rongga antar  butir  semakin  besar  hingga  rongga  dan percelahan  memiliki  nilai  yang  semakin  tinggi  dan semakin  rentan  terhadap  kerentanan  airtanah (Vienastra  &  Febriarta,  2018). Soil  media (S)/tekstur tanah merupakan parameter kerentanan dari sifat peresapan  air (imbuhan) dari permukaan  tanah  hingga  masuk  ke  sistem  akuifer (Lathamani  et  al.,  2015).  Kondisi  tekstur  semakin lepas-lepas mempengaruhi nilai kerentanan, selain
tekstur,  ketebalan  tanah  mempengaruhi  proses filtrasi (Putranto et al, 2016).
Parameter topography (T)/kemiringan  lereng merupakan  faktor  yang  mempengaruhi  dalam meresapkan  imbuhan  airtanah.  Daerah  dengan kemiringan  lereng  yang  semakin  datar  hingga landa akan mempengaruhi proses filtrasi (Febriarta,  et  al.,  2020). Impact  of  vadose  zone media (I)/jenis zona tak jenuh, dipengaruhi oleh sifat  aliran  dari  setiap  batuan.  Parameter  terakhir adalah Conductivity  hydraulic (C)/konduktivitas hidraulis  akuifer,  merupakan  kemampuan  batuan dalam  melolosakan  airtanah  yang  dinyatakan dalam satuan m/hari (BSN, 2016). Nilai konduktivitas  hidraulis  akuifer  yang semakin cepat akan  mempunyai  nilai  kerentanan  yang  tinggi terhadap kerentanan pencemaran airtanah (Febriarta  &  Widyastuti,  2020).  Nilai  potensi konduktivitas  hidraulis,  mengacu  persamaan  nilai terhadap  batuan  (Todd  &  Mays,  2005)  ( Tabel  2 ). Penilaian  dalam  setiap  kelas  dinilai  dengan  nlai terendah 1 dan terbesar atau yang mempengaruhi faktor  kerentanan  dengan  nilai  10.  Kelas  dan  nilai setiap  parameter  kerentanan  airtanah  terhadap pencemaran  dengan  metode  DRASTIC  disajikan pada Tabel 1 .
Setiap  parameter  kerentanan  memiliki  bobot dominan terhadap kerentanan yang paling mempunyai dampak terhadap pencemaran. Penilaian  setiap  kelas  dari  komponen  parameter kerentanan dikalikan dengan nilai bobot ( Tabel  3 )  untuk  mendapatkan  informasi  skala pengaruh  kerentanan  penecemaran.  Nilai  skala pembobotan adalah 1-5, dengan angka pembobotan  bernilai  besar  (nilai  5)  menunjukkan faktor  yang  penting  dan  paling  berpengaruh, sedangkan  niali  kecil  (nilai  1)  menunjukkan  faktor yang paling rendah dalam mempengaruhi kerentanan.
Tabel 1 . Kelas dan nilai (R) parameter kerentanan air tanah terhadap pencemaran dengan metode DRASTIC (D) Kedalaman Air tanah (m) (R) Imbuhan Air (mm) (A) Jenis Akuifer (S) Tekstur Tanah (T) Kemiringan Lereng (%) (I) Jenis Zona Tak Jenuh (C) Konduktivitas Hidraulis (m/hari) Kelas Nilai Kelas Nilai Kelas Nilai Kelas Nilai Kelas Nilai Kelas Nilai Kelas Nilai 0-3  10 200-300  9 Batu gamping, karst 10 kerikil 10 0-2 10 lapisan batas 1 >82 10 3-4,5  9 150-200  8 basal 9 pasir 9 2-6 9 lempung/ lanau 3 41-82 8 4,5-9  7 120-150  6 Batu gamping masif, pasir, kerikil 8 gambut 8 6-12 5 serpih 3 28,8-41 6 9-15  5 100-120  5 Batu pasir masif 6 agregat lempung 7 12-18 3 Batu gamping 3 12,3-28,7 4 15-22  3 80-100 3 Batu pasir, batu gamping, serpih 6 lempung pasiran 6 >18 1 Batu pasir 6 4,1-12,3 2 22-30  2 60-80 2 glacial till / sedimen oleh pencairan es 5 lempung 5 pasir kerikil dan lanau 6 0,04-4,1 1 >30  1 lapukan metamorf/ beku 4 lempung lanau 4 pasir dan kerikil 8 metamorf/ beku 3 lempung liat 3 basal 9 Serpih masif 2 Campuran 2 Batu gamping, karst 10 lempung non agregat 1
Sumber: Aller et al, (1987) ; White et al, (2017); Putranto et al, (2019); Voutchkovaa et al, (2021)
Tabel 2. Konduktivitas hidraulis pada batuan. Material batuan Konduktivitas Hidrolik (m/hari) Material batuan Konduktivitas Hidrolik (m/hari) Kerikil, kasar 150 Gumuk Pasir 20 Kerikil, halus 250 Gambut 5,7 Pasir, kasar 45 Skis / batuan metamorf 0,2 Pasir, sedang 12 Slate /batuan metamorf halus 0,00008 Pasir, halus 2,5 Till sedimen, pasiran 0,49 Lanau 0,08 Till sedimen, berkerikil 30 Tanah liat 0,0002 Tuff 0,2 Batu pasir, berbutir halus 0,2 Basal 0,01 Batu pasir, berbutir sedang 3,2 Gabbro, lapuk 0,2 Batugamping 0,94 Granit, lapuk 1,4 Dolomit 0,001 Sumber: Todd & Mays, (2005); Singhal & Gupta, (2010); Fetter, (2014) Tabel 3. Bobot (W) parameter DRASTIC. Parameter Nilai Bobot (D) Depth to water / kedalaman muka airtanah 5 (R) Recharge / Imbuhan Air 4 (A) Aquifer / jenis akufer 3 (S) Soil media / tekstur tanah 2 (T) Topography / kemiringan lereng, 1 (I) Impact of vadose zone media / jenis zona tak jenuh,
5 (C) Conductivity hydraulic / konduktivitas hidraulis akuifer 3
Sumber: Aller et al, (1987); Sugianti et al, (2016); Devianto et al, (2019); Putranto et al, (2019)
Perhitungan  nilai  dan  bobot  dari  setiap parameter kerentanan DRASTIC, dilakukan secara keruangan (spasial) untuk mengambarkan informasi  distribusi  kondisi  sertiap  parameter kerentanan  DRASTIC.  Fungsi  matematika  yang digunakan  untuk  menghitung  seluruh  parameter kerentanan (Tabel 1) dengan pembobotan sekenario  (Tabel  3)  digunakan  fungsi  perhitungan linier  setiap  variabel.  Hasil  dari  perhitungan menghasilkan  indeks  kerentanan.  Pendekatan rumus yang digunakan untuk menghitung nilai total dari indeks kerentanan (Aller et al., 1987b), seperti pada Persamaan 2 .
Kerentanan  =  D R D W +  R R R W +  A R A W +  S R S W + T R T W + I R I W + C R C W ........................... .................. (2)
dimana:
( R ) : nilai untuk masing-masing parameter ( W ) : bobot untuk masing-masing parameter (D) : Depth to water / kedalaman muka airtanah (R) : Recharge /Imbuhan aAir (A) : Aquifer /jenis akufer
(S) : Soil media /tekstur tanah (T) : Topography /kemiringan lereng,
(I)  : Impact  of  vadose  zone  media /jenis  zona  tak jenuh,
(C): Conductivity  hydraulic /konduktivitas  hidraulis akuifer
Tabel 4. Kelas kerentanan airtanah DRASTIC. Kelas kerentanan Nilai indeks DRASTIC Sangat Rendah 86-100 Rendah 101-120 Sedang 121-140 Tinggi 141-148 Sangat Tinggi >140
Sumber: Pacheco et al, (2015); Thapa et al, (2018)
Penentuan  kelas  kerentanan  pencemaran airtanah  dengan  metode  DRASTIC,  mengacu pada  pembagian  kelas  menurut  Pacheco  et  al., (2015)  dan Thapa  et al., (2018).  Pembagian kelas kerentanan  airtanah  tersebut,  dibagi  menjadi  lima (5)  kelas,  yaitu  sangat  rendah,  rendah,  sedang, tinggi  dan  sangat  tinggi  dengan  rentang  nilai seperti pada Tabel 4 .

## HASIL DAN PEMBAHASAN

## Kedalaman Muka Airtanah (D)
Kedalaman  muka  airtanah  merupakan  selisih nilai  kedalaman  dari  permukaan  tanah  hingga  ke tinggi muka airtanah dalam satuan meter (m). Nilai kedalaman  muka  airtanah  merupakan  data sekunder  yang  diperoleh  dari  hasil  kajian  BPLH, (2013),  dan  BPLH,  (2014).  Kondisi  kedalaman airtanah  pada  perbukitan  karst  (bagian  utara) memiliki  kedalaman  9-11  m.  Terdapat  beberapa lokasi  yang  tergenang  pada  cekungan  dan  rawa dengan  kondisi  kedalaman  airtanah  <2  m.
Sedangkan  kondisi  kedalaman  airtanah  pada dataran  secara  umum  kedalaman  muka  airtanah relatif  dangkal  dengan  nilai  rata-rata  0-4  meter. Daerah pesisir sepadan pantai memiliki kedalaman muka  airtanah  0  m,  hal  tersebut  dipengaruhi  oleh tinggi  muka  laut  (0  mdpal).  Penilaian  nilai  dan bobot, distribusi kedalaman muka airtanah disajikan pada Tabel 5 .
Tabel 5. Perhitungan nilai kedalaman muka airtanah (D). Kedalaman muka airtanah (m) Kelas (m) Nilai  Bobot Skor 2 0-3 10 5 50 3 3-4,5 9 5 45 5 4,5-9 7 5 35 10 9-15 5 5 25 15 15-22 3 5 15

## Imbuhan Airtanah (R)
Kondisi  imbuhan  airtanah  diperoleh  dari perhitungan  neraca  air  meteorologis,  dengan mengunakan  data  curah  hujan  dan  temperatur rata-rata  selama  10  tahun.  Di  Kabupaten  Banggai Kepulauan  tidak  terdapat  stasiun  hujan  sebagai sumber  data  pengolahan,  sehingga  sumber  data sekunder  menggunakan  data  curah  hujan  dari Alfiandy  et  al,  (2020),  BPS  (2020),  BPLH,  (2014) BMKG, (2020b), BMKG (2020a), dan Paski, (2018) yang  dianalisis  dengan  metode Thornthwaite- Matter , dan pendekatan pada Persamaan 1 untuk mengetahui imbuhan airtanah. Berdasarkan pendekatan  tersebut,  diketahui  hasil  perhitungan nilai  potensi  air  imbuhan  di  Pulau  Bakalan  merata dengn  nilai  sebesar  235  mm/tahun.  Menurut Pacheco et al., (2015) kondisi curah hujan tersebut kategori  lebat  atau  tinggi.  Nilai  tersebut  termasuk kelas  nilai  yang  tinggi  dari  kelas  kerentanan ( Tabel 6 ) dalam mempengaruhi potensi pencemaran airtanah.
Tabel 6. Perhitungan nilai imbuhan airtanah (R). Imbuhan air (mm/tahun) Kelas (m) Nilai  Bobot Skor 235 200-300 9 4 36
Jenis Akuifer (A )
Parameter  jenis  akuifer  merupakan  jenis batuan  yang  mempunyai  sifat  akuifer  atau  dapat menyimpan dan meloloskan airtanah. Jenis akuifer di  Pulau  Bakalan  secara  umum  didominasi  oleh batugamping yang merupakan bagian dari Formasi Batugamping  koral  (Ql).  Menurut  KESDM  (2015a) dan  KESDM  (2015b)  jenis  batuan  (litologi)  akuifer terdiri atas batugamping terumbu di bagian tengah pulau  dan  berada  di  kawasan  karst,  dan batugamping koral.
Batugamping  terumbu  memiliki  faktor  yang lebih  mempengaruhi  dalam  kerentanan  airtanah dengan  nilai  10  poin.  Secara  umum  jenis  akuifer tersebut, memiliki nilai yang tinggi dipengaruhi oleh sifat  rongga  antar  butir  atau  media  aliran  air dengan  sifat  rekah,  celah  maupun  saluran pelarutan,  yang  mempunyai  respon  yang  cepat kedalam  sistem  akuifer  airtanah.  Jenis  akuifer pada  bagian  selatan,  tersusun  atas  batuan fragmental  dari  rombakan  batugamping  bioklastik, nilai  tersebut  memiliki  nilai  yang  tinggi  (8  poin) dalam  parameter  kerentanan.  Jenis  litologi  akuifer berupa batugamping dengan media celah memiliki potensi  mengalirkan  air  dengan  cepat  (Harter  & Walker,  2001).  Seperti  pada  hasil  penelitian Vienastra  &  Febriarta  (2020)  dan  Febriarta  et  al., (2020)  penilaian  kerentanan  pada  batuagamping memiliki  kerentanan  sangat  tinggi.  Nilai  jenis akuifer  di  Pulau  Bakalan  disajikan  pada Tabel  7 dan  distribusi  spasial  (keruangan)  disajikan  pada Gambar 5 .
Tabel 7.
Perhitungan nilai jenis akuifer (A). Jenis akuifer Kelas    Nilai  Bobot Skor Batugamping terumbu 10 10 3 30 Batugamping bioklastik, fragmental 8 8 3 24

## Tektur Tanah (S)
Tekstur  tanah  dalam  parameter  kerentanan airtanah,  dipengaruhi  oleh  sifat  filtrasi  airtanah masuk kedalam tanah.  BPLH, (2014) dan Sisultan (2017)  menyebutkan  bahwa,  jenis  perkembangan tanah di Pulau Bakalan adalah berasosiasi dengan Eutric Regosols . Perkembangan tanah ini memiliki tebal  0-30  cm  dengan  dominiasi  tekstur  lempung sedang.  Tanah  tersebut  meniliki  kondisi  drainase sedang.  Perkembangan  tanah  pada  daerah  pada cekungan  didominiasi  menjadi  daerah  yang tergenang  air  ( Gambar  5 )  yang  memiliki  tingkat kerentanan  rendah,  sedangkan  pada  tektur  tanah dengan  tekstur  campuran  pasiran  memiliki  nilai kerentanan  lebih  tinggi  dengan  kondisi  drainase baik.  Nilai  kerentanan  airtanah  dari  parameter tekstur  tanah  disajikan  pada Tabel  8 dan  secara keruangan (spasial) disajikan pada Gambar 5 .
Tabel 8.
Perhitungan nilai tektur tanah (S). Tekstur tanah Kelas    Nilai  Bobot Skor Pasir 9 9 2 18 Gambut/rawa (genangan) 8 8 2 16 Lempung pasiran 6 6 2 12 Lempung 2 2 2 4

## Kemiringan Lereng (T)
Kemiringan  lereng  diperoleh  dari  analisis secara keruangan (spasial) dengan sistem informasi  geografis  mengubah  nilai  ketinggian
menjadi  nilai  kemiringan  menggunan  data digital elevation  model nasional  (DEMNAS)  dari  BIG, (2018).  Secara  morofologi,  di  Pulau  Bangkalan bagian  utara  merupakan  perbukitan  karst  dengan lokasi  tertinggi  76  mdpal.  Morfologi  bergelombang tersebut  memiliki  rentang  kermiringan  lereng  12- 18% atau 8-14°. Sedangkan bagian tengah hingga selatan, mempunyai morfologi relatif datar  dengan ketinggian  rata-rata  25  mdapl.  Kondisi  kemiringan lereng secara umum berkisar 2-6% atau 2-4°, yang termasuk  landai.  Bagian  selatan  Pulau  Bakalan, banyak  terdapat  cekungan  pada  ketinggian  ±  5 mdpal,  daerah  cekungan  tersebut  berkembang menjadi  daerah  rawa/gambut  (genangan),  dan pada cekungan yang relatif besar terbentuk danau.
Terdapat  tiga  (3)  danau  besar  yaitu  Danau Petinding,  Danau  Bakalan,  dan  Danau  Kobit ( Gambar  1 ).  Topografi  disekitar  danau,  relatif datar  dengan  kemiringan  lereng  datar  hingga landai. Nilai kemiringan lereng yang datar tersebut, mempunyai  faktor  kerentanan  yang  tinggi.  Hal tersebut  dipengaruhi  oleh  daerah  yang  datar merupakan  tempat  terakumulasi  aliran  permukaan dan  meresap  kedalam  tanah  (filtrasi).  Daerah dengan  topografi  cekungan  dapat  membentuk genangan  air  permukaan  dan  terinfiltrasi  kedalam tanah (Kalhor et al.,  2019). Derah dengan dengan kemiringan  lereng  yang  semakin  terjal,  akan mempunyai  faktor  kerentanan  yang  kecil,  hal tersebut  disebabkan  oleh  aliran  permukaan  akan bergerak  cepat  mengikuti  alur  lembah,  lereng  dan relatif  sedikit  yang  masuk  kedalam  tanah  (proses infiltrasi).  Daerah  dengan  kemiringan  lereng  yang terjal  akan  menunjukkan  nilai  faktor  yang  relatif kecil.  Hasil  penilaian  kemiringan  lereng  dalam  % disajikan  pada Tabel  9 ,  dan  secara  keuangan (spasial) disajikan pada Gambar 5 .
Tabel 9.
Perhitungan nilai kemiringan lereng (T). Kemiringan lereng (%) Kelas    Nilai  Bobot  Skor 0-2 10 10 1 10 2-6 9 9 1 9 6-12 5 5 1 5 12-18 3 3 1 3 >18 1 1 1 1

## Zona Tak Jenuh (I)
Zona  tak  jenuh  merupakan  bagian  atau material  batuan  yang  berada  dari  permukaan tanah hingga muka airtanah (BSN, 2005; Febriarta &  Purnama,  2020).  Berdasarkan  sifat  batuan dalam mengalirkan aliran airtanah, parameter zona tak jenuh memiliki pengaruh yang paling berpengaruh,  hal  tersebut  dipengaruhi  oleh  laju kecepatan  airtanah  menuju  zona  jenuh  yang dipengaruhi  oleh  sifat  porositas  batuan  (Purnama & Cahyadi, 2019). Bobot untuk parameter zona tak jenuh  adalah  5  poin.  Berdasarkan  kondisi penyusun  litologi  (batuan)  menurut  KESDM, (2015a) terdapat lima (5) litologi yaitu batugamping
Pemetaan Zona Kerentanan Air Tanah di Pulau Bakalan, Kabupaten Banggai Kepulauan, Indonesia .............................. (Febriarta, et al)
terumbu  (perbukitan  karst)  yang  berada  di  bagian utara  dan batugamping terumbu, kategori tersebut memiliki nilai 10 point (tertinggi) ( Tabel 10 ).
Tabel 10. Perhitungan nilai zona tak jenuh (I). Zona Tak Jenuh Kelas Nilai Bobot  Skor Batugamping terumbu (karst) Batugamping, karst 10 5 50 Batugamping bioklastik, fragmen Batugamping 10 5 50 Rombakan batugamping Batugamping rombakan 8 5 40 Pasir Kerikil, lanau Pasir kerikil, lanau 6 5 30 Batugamping terumbu Batugamping 3 5 15
Litologi (batuan) batugamping koral bioklastik, fragmen (bagian selatan) memiliki bobot parameter tinggi,  dikarenakan  media  aliran  airtanah  berupa celah,  rekah,  saluran  hasil  pelarutan,  sehingga memiliki  potensi  kecepatan  sumber  pencemar (polutan)  yang  tinggi  menuju  zona  jenuh  air (akuifer).  Rombakan  batugamping  berada  di morfologi  dataran  bagian  tengah  pulau.  Penilaian nilai  zona tak zenuh  disajikan pada Tabel 10 , dan secara keruangan (spasial) disajikan pada Gambar 5.

## Konduktivitas Hidraulik (C)
Nilai konduktivitas hidraulik atau nilai kelulusan  airtanah,  dipengaruhi  oleh  media  aliran airtanah  dan  jenis  batuan.  Kondisi  jenis  batuan (litologi)  akuifer  didominasi  oleh  batugamping dengan  nilai  rata-rata  konduktivitas  hidraulik sebesar  0,94  m/hari.  Yang  et  al.,  (2019) menyebutkan  bahwa  nilai  konduktivitas  hidraulik akuifer dapat dipengaruhi juga oleh sifat pembawa airtanah,  seperti  pada  daerah  genangan  air (rawa/gambut) di daerah dengan morfologi cekungan  atau  alur  lembah  pada  material  pasir sedang  hingga  kasar  memiliki  kelulusan  airtanah yang  lebih  tinggi  yaitu  5,7  m/hari.  Secara  umum nilai konduktivitas hidraulik akuifer relatif kecil yang dipengaruhi  oleh  sifat  lempung  dengan  nilai pergerakan 0,02 m/hari ( Tabel 11 ). Nilai konduktivitas hidraulik yang kecil tersebut dipengaruhi oleh sifat lempung yaitu akuitard. Sifat akuitard  yaitu  kemampuan  sifat  batuan  dapat menahan dan menyimpan airtanah, tetapi melolosakan  atau  mengalirkan  air  sangat  kecil (White et al., 2017). Distribusi spasial (keruangan) nilai konduktivitas hidraulik disajikan pada Gambar 5 .
Tabel 11. Perhitungan nilai konduktivitas hidraulik (C).
konduktivitas hidraulik (m/hari) Kelas Nilai  Bobot Skor Gambut/rawa (5,7) 4,1-12,3 2 3 6 Lempung pasiran (2,5) 0,04-4,1 1 3 3 Batugamping (0,94) 0,04-4,1 1 3 3 Lempung (0,02) 0,04-4,1 1 3 3
Tabel 12.

## Luas desa terhadap tingkat kerentanan airtanah DRASTIC.
Tingkat Kerentanan Desa Bakalan Desa Bulungkobit Desa Bungin Total Luas (km 2 ) % Luas (km 2 ) % Luas (km 2 ) % Luas (km 2 ) % Sangat Tinggi 8,802 21,58 1,265 3,10 9,988 24,48 20,056 49,16 Tinggi 4,147 10,16 1,007 2,47 0,375 0,92 5,529 13,55 Sedang 4,886 11,98 2,616 6,41 0,041 0,10 7,543 18,49 Rendah 6,712 16,45 0,376 0,92 0,579 1,42 7,667 18,79 Sangat Rendah 0,001 0,003 - - - - 0,001 0,003

## Zona Kerentanan Airtanah DRASCTIC
Penentuan zona kerentanan airtanah, diperoleh  dari  metode  perhitungan  multi-kriteria dari pendekatan sifat hidrogeologi batuan DRASTIC  yang  dianalisis  secara  keruangan (spasial)  dengan sistem informasi geografis (SIG). Penentuan zona kerentanan diperoleh dari perhitungan  penjumlahan  linier  ( Persamaan  2 ) skor setiap parameter. Berdasarkan hasil perhitungan linier tersebut diperoleh nilai total skor (skor  indeks  DRASTIC).  Untuk  mengklasifikasikan zona  kerentanan  airtanah,  mengacu  rentang  nilai dari Tabel  4 .  Berdasarkan  perhitungan,  nilai indekss DRASTIC berkisar 99-193. Zona kerentanan  airtanah  terhadap  pencemaran  sangat rendah,  diperoleh  dari  rentang  nilai  indekss DRASTIC 86-100.
Daerah  dengan  tingkat  kerentanan  sangat rendah,  merupakan  daerah  yang  mempunyai proporsi  <5  %  terjadi  pencemaran  airtanah,  atau dapat  dikatakan  tidak  berpotensi  terkontaminasi dengan  sumber  pencemar  (Putranto  et  al,  2019). Hasil  dari  analisis  zona kerentanan sangat rendah hanya  seluas  0,001  km 2 atau  0,003%  ( Tabel  12 ) dari  luas  Pulau  Bakalan.  Zona  atau  daerah tersebut berada di morfologi dataran bagian utara, yang secara administrasi berada di Desa Bakalan. Karakteristik  dengan  zona  kerentanan  sangat rendah  di  Pulau  Bakalan  memiliki  kondisi kedalaman  muka  airtanah  >10  m,  kondisi  jenis litologi (batuan) akuifer berupa batugamping koral, kondisi tekstur tanah  permukaan berupa lempung, dengan  kondisi  kemiringan  lereng  >18%  atau >15°.  Kondisi  zona  tak  jenuh  dengan  penyusun batuan  berupa  batugamping  terumbu  dengan potensi  sifat  konduktivitas  hidraulis  0,02  m/hari (lempung).
Zona  tingkat  kerentanan  airtanah  rendah, berada  pada  rentang  nilai  indekss  DRASTIC  101- 120. Zona kerentanan airtanah terhadap pencemaran  rendah,  dapat  diartikan  zona  atau daerah  dengan  potensi  kerentanan  airtanah dengan  proporsi  hingga  20%  dari  sumber pencemar  (polutan).  Potensi  pencemaran  pada zona  tersebut  dengan  intensitas  yang  kecil.  Nilai indekss kerentanan airtanah terhadap pencemaran dengan  metode  DRASTIC  rendah  di  Pulau Bakalan, memiliki karakteristik kedalaman airtanah berkisar  5-10  m,  dengan  jenis  litologi  (batuan) akuifer berupa batugamping koral, dengan perkembangan  jenis  tanah  berupa  lempung pasiran  yang  mempunyai  potensi  filtrasi  yang relatif  pelan  dan  permeabilitas/kelulusan  air  (K)
atau kemampuan batuan melululuskan air di dalam rongga-rongga  dalam  batuan  yang  rendah  pada penyusun  batu  lempung  (0,02  m/hari).  Faktor kemiringan lereng dengan dengan rentang 12-18% atau  8-15°.  Faktor  kemiringan  lereng  dan  dan kedalaman  muka  airtanah  merupakan  faktor  yang dominan pada zona kerentanan rendah. Luas zona kerentanan  airtanah  rendah  adalah  seluas  7,667 km 2 atau  seluas  18,794%  ( Tabel  12 )  dari  luas Pulau  Bakalan.  Zona  tingkat  kerentanan  rendah terdistribusi  pada  bagian  dataran  bagian  utara, yang  sebagian  besar  secara  administrasi  masuk wilayah Desa Bakalan ( Gambar 6 ).
Zona  tingkat  kerentanan  airtanah  terhadap pencemaran  sedang, memiliki  rentang  nilai indekss  DRASTIC  121-140.  Zona  atau  daerah dengan  tingkat  kerentanan  sedang,  memiliki potensi  tercemar  (kontaminiasi)  hingga  proporsi 60%  dari  luas  wilayah  (Pacheco  et  al.,  2015). Potensi  pencemaran  dapat  terjadi  pada  intensitas yang  relatif  lama  dan  menerus  dari  sumber pencemaran (polutan). Karakteristik tingkat kerentanan  airtanah  sedang  di  Pulau  Bakan dipengaruhi oleh nilai kedalaman air  berkisar 5 m, dengan  kondisi  perkembangan  tanah  lapisan  atas dengan tekstur lempung pasiran dengan ketebalan hingga 30 cm. Faktor kemiringan lereng berkisar 2- 4% atau 2-3°, dengan morfologi landai. Kemiringan lereng  tersebut  mempunyai  potensi  polutan (sumber  pencemar)  relatif  mudah  menuju  zona jenuh  air.  Kondisi  tersebut  juga  dipengaruhi  oleh kecepatan  aliran  (konduktivitas  hidraulik)  oleh batuan  penyusun  berupa  batugamping,  dengan potensi 0,94 m/hari. Kondisi zona jenuh Karakteristik  zona  vadoze  atau  zona  tak  jenuh tersusun  atas  batuan  batugamping  terumbu, dengan  perkembangan  ruang  rongga  antar  butir atau  media  aliran  air  yang  tidak  terlalu  rapat. Berdasarkan  analisis  keruangan  (spasial),  tingkat zona kerentanan sedang diketahui seluas 7,54 km 2 atau  13,55%  dari  luas  Pulau  Bakalan  ( Tabel  12 ). Secara  keruangan  (spasial),  zona  kerentanan sedang  terdistribusi  di  utara  dan  selatan  Pulau Bakalan ( Gambar 6 ). Secara dominan kerentanan sedang  terdapat  di  morfologi  landai  hingga  datar. Kondisi  imbuhan  airtanah  yang  masuk  kedalam tanah dengan potensi yang tinggi (235 mm/tahun), litologi (batuan) penyusun akuifer berupa batuamping terumbu, yang merupakan bagian dari Karst Banggai Kepulauan.
Zona kerentanan airtanah terhadap pencemaran  tinggi,  diperoleh  dari  nilai  indeks DRASTIC  141-148.  Pada  kondisi  ini,  zona  atau daerah  yang  memiliki  potensi  kerentanan  tinggi
Pemetaan Zona Kerentanan Air Tanah di Pulau Bakalan, Kabupaten Banggai Kepulauan, Indonesia .............................. (Febriarta, et al)
mempunyai  proporsi  atau  peluang  terkontaminiasi hingga  80%  oleh  polutan  (sumber  pencemar), yang  dipengaruhi  oleh  faktor  yang  pendukung berupa  kondisi  airtanah  dangkal  dan  sifat  batuan sangat  cepat  dalam  mengalirkan  air.  Tingkat kerentanan  tinggi  di  Pulau  Bakalan  terdapat  pada kondisi  kedalaman  muka  airtanah  0-3  m  dari permukaan  tanah.  Dangkalnya  kondisi  jenuh  air dapat mempercepat kontaminan (sumber pencemar)  masuk  menuju  zona  jenuh  air (Devianto et al, 2019; Voutchkova et al., 2021).
Nilai  imbuhan  airtanah  yang  tinggi  (235 mm/tahun), termasuk parameter yang mempengaruhi  tingkat  kerentanan  tinggi.  Faktor penyusun  batuan,  yaitu  batugamping  koral  dalam sifat zona tak jenuh, termasuk media dengan aliran yang  memiliki  tingkat  porositas  yang  baik  dalam meresapkan  menuju  akuifer  airtanah.  Anim- Gyampo  et  al.,  (2019)  menyebutkan  semakin besar media rongga antar butir, dapat menurunkan kualitas  airtanah,  hal  tesebut  dikarenakan  media tersebut  dengan  mudah  mengalirkan  polutan  dari atasnya. Hal tersebut memicu terjadianya kerentanan  pencemaan.  Faktor  kemiringan  lereng yang  landai  hingga  dapat  mempercepat  kondisi jenuh  pada  tanah  lapisan  permukaan  yang terfiltrasi menuju zona tak jenuh. Zona atau daerah dengan  kerentanan  tinggi  terdapat  pada  bagian utara  dan  bagian  tengah  disekitar  Danau  Bakalan dan  Danau  Kobit.  Luas  kerentanan  airtanah  tinggi seluas  5,529  km 2 atau  seluas  13,553%  dari  luas Pulau Bakalan ( Gambar 6 ).
Tingkat  kerentanan  tinggi  terdistribusi  di wilayah  Desa  Bulungkopit  dan  Desa  Bakalan. Kerentanan  airtanah  sangat  tinggi,  merupakan kondisi  suatu  zona  yang  mempunyai  potensi proporsi  dapat  tercemar  oleh  semua  polutan (sumber pencemar). Hal tersebut dipengaruhi oleh sifat batuan  yang  sangat  porus  terhadap  airtanah, dengan kedalaman airtanah 0-3 m dari permukaan tanah.  Selisih  kedalaman  tersebut  cukup  dangkal, yang  mengakibatkan  sumber  pencemar  cepat masuk  ke  zona  jenuh  air  (akuifer).  Kondisi kedalaman  airtanah  yang  dangkal  dengan  kondisi topografi  yang  rendah  terdapat  potensi  intrusi  air laut bila  dalam pengunaan airtanah yang berlebihan.  Hal  tersbut  juga  dipengaruhi  oleh  pola perkembangan  permukiman  yang  berada  Faktor penyusun  batuan  berupa  batugamping  terumbu (karst)  mempunyai  potensi  yang  cepat  dalam mengaliran sumber polutan  yang dipengaruhi oleh rongga  antar  butir  media  aliran  berupa  celah, rekah hingga saluran dari perkembangan pelarutan (Febriarta et al, 2020). Nilai konduktivitas hidraulis  pada  batugamping  sebesar  0,94  m/hari. Faktor kemiringan lereng 0-6% atau 0-4° termasuk
landai,  yang  mengakibatkan  terkonsentrasinya aliran  permukan  sehingga  menjadi  jenuh,  dan mudah  teresapkan  (fitrasi)  kedalam  tanah  hingga masuk  dengan  mudah  ke  zona  jenuh  (akuifer) dengan  relatif  cepat.  Zona  tingkat  kerentanan airtanah  sangat  tinggi  terdapat  pada  bagian  utara Pulau  Bakalan,  yang  sebagian  bersar  berada  di perbukitan  karst.  Luas  zona  dengan  tingkat kerentanan  pencemaran  airtanah  sangat  tinggi seluas  20,056  km 2 atau  49,16%.  Secara  umum tingkat  kerentanan  pencemaran  airtanah  di  Pulau Bakalan  didominasi  oleh  kerentanan  sangat  tinggi (49,16%).

## KESIMPULAN
Penilaian kerentanan airtanah terhadap pencemaran dengan metode DRASTIC menghasilkan  nilai  indeks  kerentanan  dengan interval yang panjang, sehingga dapat mengetahui interaksi  parameter  akuifer  yang  dominan  atau memicu  kerentanan.  Penilaian  kerentanan  dapat diterapkan  pada  akuifer  bebas  (tidak  tertekan). Pemodelan  keruangan  (spasial)  dengan  metode DRASTIC  lebih  sesuai  untuk  skala kecil.  Terdapat kekurangan  dengan  metode  ini,  yaitu  diperlukan data yang banyak.
Saran  yang  dapat  dilakukan  pada  zona kerentanan  sangat  rendah  adalah  pemantauan berkala  di  sumber-sumber  airtanah.  Kerentanan rendah mempunyai peluang terkontaminasi sampai  dengan  20%.  Saran  yang  dapat  dilakukan pada zona kerentanan rendah adalah pemantauan berkala  dengan  pengukuran  kualitas  air  sesuai baku  mutu  dan  edukasi  tentang  pengelolaan limbah,  terutama  limbah  domestik.  Kerentanan sedang mempunyai peluang terkontaminasi sampai  dengan  60%.  Saran  yang  dapat  dilakukan pada zona kerentanan sedang adalah pemantauan berkala  dengan  pengukuran  kualitas  baku  mutu, dan  pembuatan  instalasi  pembuangan  limbah (IPAL).  Kerentanan  tinggi  mempunyai  peluang terkontaminasi  sampai  dengan  80%.  Saran  yang dapat  dilakukan  pada  zona  kerentanan  tinggi adalah  dengan  pengelolaan  limbah  dengan wawasan  lingkungan  dengan  pembuatan  instalasi pembuangan limbah (IPAL) dan pembuatan sumur pantau.  Kerentanan  sangat  tinggi  mempunyai peluang  tercemar  oleh  seluruh  polutan  (sumber pencemar).  Saran  yang  dapat  dilakukan  pada zona  kerentanan  sangat  tinggi  dalam  menjaga potensi  pencemaran  airtanah  adalah  dengan pemantaun  berkala  kuantitas  dan  kualitas  air,  dan pembuatan peraturan daerah tentang perlindungan dan pengelolaan limbah.
Gambar 6 . Tingkat kerentanan airtanah terhadap pencemaran dengan metode DRASTIC Pulau Bakalan.

## UCAPAN TERIMA KASIH
Ucapan  terima  kasih  disampaikan  kepada Dinas Lingkungan  Hidup  (DLH) Kabupaten Banggai Kepuluan Provinsi Sulawesi Tengah yang telah mendukung dan memfasilitasi penelitian ini.

## DAFTAR PUSTAKA
Alfiandy,  S.,  Permana,  D.  S.,  Nurjaman,  A. W.,  Kurnia, W.  G.,  Prastika,  L.,  Panggabean,  H.,  Rahman,  M. S. D. P., Wulandari, S., Setiawan, H., Hendrawan, A.,  &  Miranda.  (2020).  Analisis  Iklim  Provinsi Sulawesi  Tengah  berdasarkan  Data  Pemantau Cuaca  Otomatis  BMKG. Buletin  GAW  Bariri , 1 (1),
1 –11. https://doi.org/https://gawpalu.id/bgb/index.php/bg b/article/view/5
Aller, L., Bennett, T., Lehr, J. H., Petty, R. J., & Hackett, G. (1987a). DRASTIC; a Standardized System for Evaluating Ground water Pollution Potential Using Hydrogeologic  Settings .  U.  S.  Environ.  Prot. Agency.
Aller,  L.,  Bennett,  T.,  Lehr,  J.,  Petty,  R.,  &  Hackett,  G. (1987b). DRASTIC;  a  Standardized  System  for Evaluating Ground water Pollution Potential Using
Hydrogeologic Settings . United States
Environmental Protection Agency.
Anim-Gyampo, M., Anornu, G. K., Appiah-Adjei, E. K., &
Agodzo,  S.  K.  (2019).  Quality  and  health  risk assessment  of  shallow  groundwater  aquifers within the Atankwidi basin of Ghana. Groundwater for Sustainable Development , 9 , 100217. https://doi.org/10.1016/j.gsd.2019.100217
BIG.  (2018). DEMNAS  2314-441  dan  DEMNAS2314- 423 . Badan Informasi Geospasial (BIG). BMKG.  (2020a). Analisis  Hujan  Agustus  2020  dan Prakiraan Hujan Oktober, November dan Desember 2020 . Pusat Informasi Perubahan Iklim Kedeputian Bidang Klimatologi Badan Meteorologi Klimatologi dan Geofisika.
BMKG.  (2020b). Analisis  Hujan  April  2020  dan Prakiraan  Hujan  Juni,  Juli  dan  Agustus  2020 . Badan  Meteorologi,  Klimatologi,  dan  Geofisika (BMKG).
BouDagher-Fadel, & Marcelle, K. (2013). Biostratigraphic  and  Geological  Significance  of Planktonic Foraminifera . UCLPress. London. BPLH. (2013). Kajian Potensi Mataair dan Pengelolaannya  Kabupaten  Banggai  Kepulauan (Laporan  Ak).  Badan  Pengelolaan  Lingkungan Hidup (BPLH) Kabupaten Banggai Kepulauan.
BPLH.  (2014). Penyusunan  Ekoregion  dan  Rencana Perlindungan  dan  Pengelolaan  Lingkungan  Hidup (RPPLH) Kabupaten Banggai Kepulauan (Laporan Ak).  Badan  Pengelolaan  Lingkungan  Hidup (BPLH) Kabupaten Banggai Kepulauan.
BPS.  (2020). Kabupaten  Banggai  Kepulauan  Dalam
Pemetaan Zona Kerentanan Air Tanah di Pulau Bakalan, Kabupaten Banggai Kepulauan, Indonesia .............................. (Febriarta, et al)
Angka  2020 .  Badan  Pusat  Statistik  (BPS) Kabupaten Banggai Kepuluan. BSN.  (2002). Standar  Nasional  Indonesia  No.  19- 6728.1-2002  Penyusunan  Neraca  Sumber  Daya  - Bagian  1:  Sumber  Daya  Air  Spasial.  Badan Standardisasi  Nasional .  Badan  Standardisasi Nasional.
BSN.  (2005). Standar  Nasional  Indonesia  (SNI)  no.13- 1712-2005  tentang  penyelidikan  potensi  air  tanah skala  1:100.000  atau  lebih  besar .  Badan Standardisasi Nasional (BSN). BSN.  (2016). Standar  Nasional  Indonesia  (SNI)  Nomor 8283:2016 Metode Pengukuran Kedalaman Menggunakan Alat Perum Gema Untuk Menghasilkan Peta Batimetri. Badan
Standaridisasi Nasional (BSN).
Cahyadi,  A.  (2019).  Analisis  Kerentanan  Airtanah terhadap Pencemaran di Pulau Koral Sangat Kecil dengan  Menggunakan  Metode  GOD. Jurnal
Geografi , 16 (1), 48 –53. https://doi.org/https://doi.org/10.15294/jg.v16i1.184 11.
Devianto,  L.  A.,  Lusiana,  N.,  &  Ramdani,  F.  (2019). Analisis  Kerentanan  Pencemaran  Air  Tanah  di Kota  Batu  Menggunakan  Analisis  Multikriteria Spasial dengan Indekss DRASTIC. Jurnal Wilayah Dan Lingkungan , 7 (2), 90 –104. https://doi.org/10.14710/jwl.7.2.90-104 DPU.  (1996). Standar  Kebutuhan  Air  Rumah  Tangga . Pengembangan  Kawasan  Perkotaan,  Kawasan Perdesaan Dirjen Pekerjaan Umum.
Febriarta, E., & Larasati, A. (2020). Karakteristik Akuifer Air  Tanah  Dangkal  Di  Endapan  Muda  Merapi Yogyakarta. Jurnal Sains Dan Teknologi Lingkungan , 12 (2), 84 –99.
https://doi.org/https://doi.org/10.20885/jstl.vol12.iss 2.art1
Febriarta, E., Marfai, M. A., Hizbaron, D. R., & Larasati, A.  (2020).  Kajian  Spasial  Multi  Kriteria  DRASTIC Kerentanan Air tanah Pesisir Akuifer Batugamping di  Tanjungbumi  Madura. Jurnal  Ilmu  Lingkungan , 18 (3),
476 –487.
https://doi.org/10.14710/jil.18.3.476-487
Febriarta,  E.,  &  Oktama,  R.  (2020).  Pemetaan  Daya Dukung  Lingkungan  Berbasis  Jasa  Ekosistem Penyedia  Pangan  Dan  Air  Bersih  Di  Kota Pekalongan. Jurnal  Ilmu  Lingkungan , 18 (2),  283 – 289. https://doi.org/10.14710/jil.18.2.283-289 Febriarta,  E.,  Oktama,  R.,  &  Purnama,  S.  (2020). Analisis  Daya  Dukung  Lingkungan  Berbasis  Jasa Ekosistem  Penyediaan  Pangan  dan  Air  Bersih  di Kabupaten  Semarang. Geomedia:  Majalah  Ilmiah Dan  Informasi  Kegeografian , 18 (1),  12 –24. https://doi.org/10.21831/gm.v18i1.30612
Febriarta, E., Prabawa, B. A., & Rosaji, F. S. C. (2018). Sumberdaya  Air  di  Pulau  Pelapis  Kepulauan Karimata,  Kabupaten  Kayong  Utara,  Kalimantan Barat. SEMINAR  NASIONAL  IV  Pengelolaan Pesisir  Dan  Daerah  Aliran  Sungai , 4 ,  174 –181. https://doi.org/https://doi.org/10.17605/osf.io/v6nx8 Febriarta,  E.,  &  Purnama,  S.  (2020).  Identifikasi Keterdapatan  Airtanah  Dengan  Electromagnetic Very Low Frequency ( EM-VLF ) di Non Cekungan Airtanah  Kecamatan  Ungaran  Timur. Jurnal
Geosains Dan Teknologi , 3 (2), 52 –62.
https://doi.org/https://doi.org/10.14710/jgt.3.2.2020 .52-62
Febriarta,  E.,  &  Widyastuti,  M.  (2020).  Kajian  Kualitas Air  Tanah  Dampak  Intrusi  Di  Sebagian  Pesisir
Kabupaten  Tuban. Jurnal  Geografi :  Media
Informasi Pengembangan Dan Profesi Kegeografian , 17 (2), 39 –48. https://doi.org/10.15294/jg.v17i2.24143 Ferreira,  J.  P.,  Lobo,  Oliveira,  &  Manuel,  M.  (2004). Groundwater Vulnerability Assessment in Portugal. Geofisica International , 4 , 541 –550.
Fetter,  C.  W.  (2014). Applied  Hydrogeology .  Pearson New Internasional Edition.
Gunawan,  W.  A.  F.,  Sisinggih,  D.,  &  Dermawan,  V.
(2012).  Cekungan  Airtanah  Negara  Kabupaten
Jembrana. Jurnal Pengairan Universitas Brawijaya , 4 (2).
Harter,  T.,  &  Walker,  L.  .  (2001). Booklet:  Assessing Vulnerability of Groundwater . California
Departement of Health Services.
Kalhor, K., Ghasemizadeh, R., Rajic, L., &
Alshawabkeh, A. (2019). Assessment of groundwater  quality  and  remediation  in  karst aquifers:  A  review. Groundwater  for  Sustainable
Development , 8 , 104 –121. https://doi.org/10.1016/j.gsd.2018.10.004 KESDM.  (2012). Geologi .  Kementerian  Energi  dan Sumber Daya Mineral. https://geoportal.esdm.go.id/geologi/ KESDM.  (2015a). Litologi  Akuifer .  Kementerian  Energi dan Sumber Daya Mineral. https://geoportal.esdm.go.id/geologi
KESDM.  (2015b). Produktivitas  Akuifer .  Kementerian Energi dan Sumber Daya Mineral.
KESDM. (2018). Cekungan Air Tanah (CAT) . Kementerian  Energi  dan  Sumber  Daya  Mineral. https://geoportal.esdm.go.id/geologi/
Lathamani,  R.,  Janardhana,  M.  R.,  Mahalingam,  B.,  & Suresha,  S.  (2015).  Evaluation  of  Aquifer Vulnerability Using Drastic Model and GIS: A Case Study  of  Mysore  City,  Karnataka,  India. Aquatic
Procedia , 4 , 1031 –1038. https://doi.org/10.1016/j.aqpro.2015.02.130 Linggasari,  S.,  Cahyadi,  T.  .,  &  Ernawati,  R.  (2019). Overview Metode Perhitungan Kerentanan Airtanah Terhadap Rencana Penambangan.
Prosiding  Nasional  Rekayasa  Teknologi  Industri Dan  Informasi  XIV  Tahun  2019  (ReTII). , 14 ,  123 – 129.
Pacheco, F. A. L., Pires, L. M. G. R., Santos, R. M. B., &  Sanches  Fernandes,  L.  F.  (2015).  Factor weighting  in  DRASTIC  modeling. Science  of  The Total Environment , 505 , 474 –486. https://doi.org/10.1016/j.scitotenv.2014.09.092 Paski, J. A. I. (2018). Analisis Kondisi Atmosfer dengan Memanfaatkan Citra Satelit Cuaca dan
Karakteristik Tanah pada Kejadian Tanah Longsor di  Pesisir  Barat Lampung  sepanjang  Tahun  2014. Jurnal Lingkungan Dan Bencana Geologi , 9 (2), 65. https://doi.org/10.34126/jlbg.v9i2.161
Purnama,  S.,  &  Cahyadi,  A.  (2019).  Groundwater Vulnerability to Pollution in Kasihan District, Bantul Regency, Indonesia. Forum  Geografi , 33 (2), 140 – 152. https://doi.org/10.23917/forgeo.v33i2.7672 Purnama, S., Tivianton, T. A., Cahyadi, A., & Febriarta, E.  (2019).  Kajian  Daerah  Imbuhan  Airtanah  di Kabupaten  Ngawi. Kajian  Daerah  Imbuhan Airtanah  Di  Kabupaten  Ngawi , 16 (1),  54 –59. https://doi.org/10.15294/jg.v16i1.18358
Putranto,  T.  T.  (2019).  Studi  Kerentanan  Airtanah Terhadap  Pencemaran  dengan  Menggunakan Metode  Drastic  Pada  Cekungan  Airtanah  (CAT) Karanganyar-Boyolali,  Provinsi  Jawa  Tengah.
Jurnal Ilmu Lingkungan , 17 (1), 159. https://doi.org/10.14710/jil.17.1.159-171
Putranto,  T.  T.,  Widiarso,  D.  .,  &  Yuslihanu,  F.  (2016). Studi Kerentanan Air Tanah Terhadap Kontaminan Mengginakan Metode Drastic di Kota Pekalongan. Teknik , 37 (1), 26 –31.
Putranto, T. T., Winarno, T., & Susanta, A. P. A. (2020). Risk  Assessment  of  Groundwater  Abstraction Vulnerability Using Spatial Analysis: Case Study at
Salatiga Groundwater Basin, Indonesia. Indonesian Journal on Geoscience , 7 (2). https://doi.org/10.17014/ijog.7.2.215-224
RI. (2010). Peraturan Menteri Energi dan Sumber Daya Mineral  Republik  Indonesia  No.17  Tahun  2012 tentang Penetapan Kawasan Bentang Alam Karst . Menteri  Energi  dan  Sumber  Daya  Mineral Republik Indonesia.
Singhal,  B.  B.  .,  &  Gupta,  R.  .  (2010). Applied Hydogeology of Fracture Rock . Springer Dordrecht Heidelberg London.
Sisultan.  (2017). Peta  Tanah  Skala  1:50.000  Sulawesi Tengah .  Sisultan  Informasi  Sumberdaya  Lahan Pertanian (Sisultan) Kementerian Pertanian.
Sugianti,  K.,  Mulyadi,  D.,  &  Maria,  R.  (2016).  Analisis kerentanan pencemaran air tanah dengan pendekatan DRASTIC di Bandung Selatan. Jurnal Lingkungan  Dan  Bencana  Geologi , 7 (1),  19 –33. https://doi.org/http://dx.doi.org/10.34126/jlbg.v7i1.9 1
Surono, & Sukido. (1993). Peta Geologi Lembar Kolaka, Sulawesi, . Badan Geologi.
Thapa,  R.,  Gupta,  S.,  Guin,  S.,  &  Kaur,  H.  (2018). Sensitivity  analysis  and  mapping  the  potential groundwater  vurnerability  zones  in  Birbhum district,  India:  A  comparative  approach  between vurnerability models. Water Science , 32 , 44 –66. Todd,  D.  K.,  &  Mays,  L.  W.  (2005). Groundwater Hydrology (3rd ed.). John Wiley & Sons, Inc. Trinanda,  T.  C.  (2017).  Pengelolaan  Wilayah  Pesisir Indonesia  dalam  Rangka  Pembangunan  Berbasis Pelestarian  Lingkungan. Matra  Pembaruan ,  75 – 84. https://doi.org/10.21787/mp.1.2.2017.75-84
Vienastra,  S.,  &  Febriarta,  E.  (2018).  Karakteristik Airtanah  di  Pulau  Yeben,  Kabupaten  Raja  Ampat,
Papua Barat. Prosiding Pertemuan Ilmiah Tahunan Ke-3 Perhimpunan Ahli Airtanah Indonesia , 3 (November), 108 –113.
https://doi.org/https://doi.org/10.17605/OSF.IO/EZ HDT
Vienastra,  S.,  &  Febriarta,  E.  (2020).  Penentuan  Zona Kerentanan  Airtanah  Metode  Simple  Vertical Vulnerability Di Pulau Yeben. JURNAL SWARNABHUMI : Jurnal Geografi Dan Pembelajaran Geografi , 5 (2), 58.
https://doi.org/10.31851/swarnabhumi.v5i2.4431
Voutchkova,  D.  D.,  Schullehner,  J.,  Rasmussen,  P.,  &
Hansen,  B.  (2021).  A  high-resolution  nitrate vulnerability assessment of sandy aquifers (DRASTIC-N). Journal of Environmental Management , 277 , 111330.
https://doi.org/10.1016/j.jenvman.2020.111330 Vrba,  J.,  &  Zaporozec,  A.  (1994). Guidebook  on Mapping  Groundwater  Vulnerability .  International Association of Hydrogeologist.
White,  W.  B.,  Herman,  J.  S.,  Herman,  E.  K.,  & Rutigliano. (2017).
Advances in Karst Science:Karst  Groundwater  Contaminiation  and Public . Springer Switzerland. http://www.springer.com/series/15147 Wicaksono,  A.  P.,  Febriarta,  E.,  Nurani,  D.  T.  T.,  & Larasati,  A.  (2020).  Evaluasi  Kebutuhan  Air Persemaian Di Kawasan Karst Nggorang Manggarai  Barat,  Labuan  Bajo,  Nusa  Tenggara Timur. Jurnal  Ilmu  Lingkungan , 18 (3),  572 –581. https://doi.org/10.14710/jil.18.3.572-581
Yang,  S.,  Ge,  W.,  Chen,  H.,  &  Xu,  W.  (2019). Investigation  of  soil  and  groundwater  environment in  urban  area  during  post-industrial  era:  A  case study of brownfield in Zhenjiang, Jiangsu Province, China. China Geology , 3 (4), 504 –514.
https://doi.org/10.31035/cg2018128