Shalaho Dina Devy Teknik Pertambangan, Fakultas Teknik, Universitas Mulawarman, Samarinda E-Mail: [email protected]

Shalaho Dina Devy Teknik Pertambangan, Fakultas Teknik, Universitas Mulawarman, Samarinda E-Mail: Shalaho.D2@Gmail.Com

SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

Pemodelan airtanah dan Neraca Airtanah Dampak Penambangan Batubara Open Pit pada Lipatan Sinklin di Daerah Muara Lawa, Kabupaten Kutai Barat, Provinsi Kalimantan Timur

Shalaho Dina Devy Teknik Pertambangan, Fakultas Teknik, Universitas Mulawarman, Samarinda E-mail: [email protected]

Abstract Ch a n g e s in land use impacts on groundwater availability both in quality and quantity. Study of geological, hydrogeological conditions, hydrological, hydrogeological boundary conditions contribute to determine the presence of groundwater in the groundwater basin and determine the natural groundwater flow modeling patterns. Location of the study are in Muara Lawa, including the formation Pulaubalang, Pamaluan, and Balikpapan on Lampanan syncline structure. Exploration drilling results indicate, that the hydrogeological study area into the aquifer system of folded sedimentary rocks consisting of seven aquifer layers alternating between akuitar, aquifers, and the base layer in the form of akuiklud. The area bounded by the limits of the model study of surface water in two major rivers, the Lawa River (east) and the Perak River (west), as well as the groundwater divide with the highest head (north and south. The results of the analysis of the type of regional aquifers aquifer models including the category of semi depressed with aquitards dominance in the surface layer. The pattern of groundwater flow and head height can be predicted by modeling using Visual Modflow. There is an increase zone budget is almost five times higher than the natural condition to the active mining conditions, particularly in the green zone, that of the natural conditions of 1,502 m3 day-1 to 12,930 m3 day-1 in the active mining. This is due, in the mining area there is a change in land use and hidrostratigrafi, namely the formation of pits that reach depths of 70 meters.

Keywords: syncline, open pit, groundwater flow, zone budget.

Abstrak Perubahan tata guna lahan berdampak pada ketersediaan airtanah baik secara kualitas maupun kuantitas. Struktur geologi dan aktivitas penambangan batubara open pit sangat mempengaruhi perubahan pola aliran airtanah. Kondisi seperti ini terjadi pada lokasi penambangan batubara open pit yang berada di kecamatan Muara Lawa yang menjadi daerah model penelitian. Secara metodologis, penelitian ini menggunakan kajian pendekatan dari kondisi hidrogeologi, hidrologi, dan kondisi batas hidrogeologi. Hasil dari penelitian menjelaskan, bahwa Formasi yang berkembang daerah penelitian yaitu, terdiri dari Formasi Pulaubalang, Pamaluan, dan Balikpapan dengan struktur Sinklin Lampanan, dengan sumbu sinklin yang membentang dari timur laut menuju ke barat daya. Hasil kajian hidrogeologi terungkap, bahwa daerah penelitian masuk dalam sistem akuifer batuan sedimen terlipat yang terdiri dari tujuh lapisan akuifer yang berselang-seling antara akuitar, akuifer, dan lapisan dasar berupa akuiklud. Daerah model penelitian dibatasi oleh batas air permukaan, yaitu Sungai Lawa dan Sungai Perak, serta batas pemisah airtanah (utara dan selatan). Pada analisis tipe akuifer, daerah model termasuk kategori akuifer semi tertekan dengan dominasi akuitar. Perubahan morfologi dan geologi akibat aktivitas penambangan batubara open pit berpengaruh terhadap perubahan aliran airtanah, antara lain: daerah pit tambang yang mengalami penurunan elevasi sampai -60 m dpl dan penambahan elevasi hingga 40 m pada daerah disposal. Kondisi ini berdampak perubahan pola arah aliran airtanah ke arah pit dan penurunan muka airtanah piezometrik dan kenaikan neraca airtanah hampir lima kali lipat dari kondisi alami ke kondisi penambangan aktif yang disebabkan perubahan tata guna lahan dan hidrostratigrafi. Perubahan-perubahan dampak aktivitas penambangan dapat diprediksi dengan memodelkan menggunakan Visual Modflow.

Kata Kunci: sinklin, open pit, aliran airtanah, neraca airtanah, void.

1. Pendahuluan Penambangan batubara open pit merupakan kegiatan yang menyebabkan perubahan tata guna lahan, morfologi, geologi, dan hidrogeologi. Perubahan ini berdampak pada perubahan topografi, perlapisan batuan bawah permukaan, dan akuifer. Penambangan open pit ini mendasarkan pada aktivitas penggalian ke arah vertikal hingga mencapai lapisan endapan batubara yang diinginkan. Salah satu perusahaan yang melakukan penambangan metode ini adalah PT. Trubaindo Coal Mining (PT. TCM). PT. TCM merupakan salah satu pemegang Perjanjian Karya Pengusahaan Pertambangan Batubara (PKP2B) generasi kedua. Perusahaan ini berlokasi di Kabupaten Kutai Barat yang berjarak 323 km dari Samarinda, ibukota Provinsi Kalimantan Timur (BPS Kubar, 2011). Daerah model penelitian berada di

69 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

Blok Selatan PT. TCM yang berada di kecamatan Muara Lawa, Kabupaten Kutai Barat, Kalimantan Timur.

2. Metode Penelitian Metode penelitian yang diterapkan adalah metode induktif dengan pendekatan analitik, yang terdiri dari kajian kondisi hidrologi, hidrogeologi, perubahan tata guna lahan, dan pemodelan airtanah, khususnya dampak akibat aktivitas penambangan batubara open pit. Selain itu, penelitian ini menggunakan ekperimen semu (Quasi Experiment Research) yang disebabkan banyaknya data primer di lapangan dan bukan hanya menggabungkan teori-teori untuk menarik suatu kesimpulan tertentu. Adapun urutan tahapan-tahapan dari metode penelitian ini meliputi: (1) persiapan; (2) pelaksanaan dan pekerjaan lapangan; (3) pengolahan dan analisis data, analisis data pada studi kasus ini menggunakan perangkat lunak Visual Madflow; dan (5) penyelesaian. Hasil akhir penelitian ini berupa prediksi model airtanah alami dan ketika operasi penambangan.

3. Hasil dan Pembahasan 3.1. Geologi Daerah Model Penelitian Van Bemmelen (1949), mengungkapkan, secara fisiografis kutai Barat masuk dalam Zona Cekungan Kutai. Sementara itu, Supriatna dkk. (1995) berpendapat, stratigrafi Cekungan Kutai ditandai oleh beberapa pengendapan formasi batuan dengan pemerian yang khas dari lingkungan daratan hingga laut dangkal.

A C 115°37'48"E 115°47'12"E 115°56'36"E PETA GEOLOGI Qa REGIONAL !. !. Samut Batang E Tpkb 4.52.25 0 4.5 U Kilometers D D U Kecamatan : Muara Lawa Kabupaten : Kutai Barat D Propinsi : Kalimantan Timur U n a 20 na K e t e r a n g a n : pa10 15 w Tmpb am a n L L li 0°41'40"S ink . Batas daerah model penelitian S S

k Batas daerah target penelitian Tomp a r e Batas administrasi P 31

. /"!. Kota kabupaten S Kota kecamatan !. H! KotaJaringan kecamatan jalan Suwakalang Jaringan sungaijalan Tmbp Jaringan sungai

!. Dasam

!. Jelimusibak 0°50'0"S

Peta Situasi D B 115°37'48"E 115°47'12"E 115°56'36"E 114°0'0"E 115°30'0"E 117°0'0"E

Formasi dan Struktur Geologi Kapuas Hulu Kutai Timur Bontang Qa (Aluvial) 0°0'0" Murung Raya 0°0'0" Kutai Kartanegara Kutai Barat Tpkb (Formasi Kampungbaru) Lokasi Penelitian Samarinda A B Tmbp (Formasi Balikpapan) 250m 250m Gunung Mas Barito Utara Penajam PaserBalikpapan Utara Tmpb (Formasi Pulaubalang) Kapuas 0m 0m Pasir Tomp (Formasi Pamaluan) Pulang PisauBarito Selatan Selat Makasar - 1°30'0"S Palangka Raya Tabalong 1°30'0"S + Barito Timur -250m -250m Katingan Antiklin 114°0'0"E 115°30'0"E 117°0'0"E -500m -500m Sinklin Sumber Data : Patahan 1. Peta Rupa Bumi Indonesia, skala 1 : 50.000, Bakorsutanal C D 2. Peta Geologi Bersistem, skala 1:50.000 Strike/dip 3. Supriatna dkk., 1995 250m 250m 15 4. Mine Plan Dept. PT. TCM, 2013 0m 0m Strike/dip dari foto udara -250m -250m Sesar Berbalik (Reverse Fault) Di susun Oleh : Shalaho Dina Devy D Sesar Turun (Normal Fault) Jurusan Teknik Geologi -500m -500m U D bagian turun dan U bagian naik Universitas Gadjah Mada Yogyakarta - + + Bidang geser arah ke depan 2 0 1 5 - Bidang geser arah ke belakang Gambar 1: Peta Geologi Regional Daerah Penelitian

Pada Gambar 1, menjelaskan korelasi dari geologi regional yang dihubungkan hasil pemboran eksplorasi. Stratigrafi daerah penelitian dari paling tua hingga paling muda pada daerah model penelitian terdiri dari: (1) Formasi Pamaluan (Miosen Awal-Miosen Bawah); (2) Formasi Pulaubalang (Miosen Tengah-Miosen Akhir); dan (3) Formasi Balikpapan (Miosen Tengah-Miosen Akhir). Kolom stratigrafinya dapat dilihat pada Gambar 2.

70 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

Gambar 2: Kolom Stratigrafi Daerah Penelitian Sumber: Supriatna dkk. (1995)

Secara umum, hasil pemboran ekplorasi lapisan batuan daerah model didominasi oleh batulanau, batupasir, batulempung, dan perselingan batubara dengan ketebalan yang bervariasi. Selain itu, dari ketiga formasi batuan tersebut teridentifikasi satuan-satuan batuan yang mendominasi tiap formasi, yaitu: satuan batulempung (Formasi Pamaluan), satuan perselingan batupasir kuarsa dan batulempung (Formasi Pulaubalang), satuan perselingan batulanau dan batupasir (Formasi Pulaubalang), dan satuan batupasir (Formasi Balikpapan).

Struktur geologi yang berkembang daerah model penelitian berupa lipatan sinklin, yaitu Sinklin Lampanan, yang dipengaruhi oleh fisiografi dari Antiklinorium Samarinda yang terdapat sesar-sesar kecil yang mengakibatkan terjadinya perpotongan perlapisan batuan. Lokasi model penelitian berada ditengah-tengah sumbu lateral sinklin yang membujur diagonal dari arah barat daya ke timur laut. Struktur lipatan sinklin mempunyai arah kemiringan jurus (strike) sebesar N 40° - 50° E dengan kemiringan perlapisan (dip) sebesar 16° - 25° yang berlokasi sebelah utara. Sementara itu, pada sebelah selatan mempunyai arah strike N 180° - 190° E dengan dip 20° - 26°.

3.2. Hidrologi Kondisi hidrologi berperan penting dalam proses penghitungan imbuhan airtanah dan penentu dalam kuantifikasi airtanah. Pemodelan airtanah harus mengenai kondisi daerah imbuhan, konduktivitas akuifer, kondisi batas hidrogeologi. Kapasitas air permukaan sangat dipengaruhi oleh karakteristik kondisi Daerah Aliran Sungai (DAS), tata guna lahan, curah hujan, dan daerah imbuhan. Daerah imbuhan sangat dipengaruhi oleh siklus hidrologi. Faktor penentu besarnya kuantitas imbuhan airtanah

71 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018 adalah intensitas hujan. Menurut Mononobe (Kamiana, 2005), nilai intensitas hujan merupakan jumlah hujan pada suatu daerah tiap satuan waktu.

(1) Besaran nilai intensitas hujan (I) pada Daerah Aliran Sungai (DAS) sangat dipengaruhi oleh waktu I konsentrasi ) dari hujan yang merata di seluruh daerah DAS. Waktu konsentrasi adalah waktu perjalanan yang diperlukan oleh air dari tempat yang paling jauh (hulu DAS) sampai ke titik pengamatan aliran air. SecaraT rinci, rumus waktu konsentrasi ) dapat dilihat di bawah ini (Hammer & Mac Kichan, 1981). T (2) T Keterangan Tc adalah waktu konsentrasi (menit), L atau jarak maksimum aliran DAS (m), V merupakan kecepatan aliran (km jam-1), dan H merupakan beda tinggi daerah hulu dengan lokasi yang ditinjau (km).

Pada kondisi alami, daerah resapan air tanah sangat dipengaruhi oleh jumlah limpasan air permukaan. Rumus yang digunakan untuk menghitung volume air limpasan ke Metode Departemen Konservasi Tanah Service (SCS) Amerika Serikat ditunjukkan pada Persamaan 3 (Hammer & Mac Kichan, 1981).

(3) Keterangan, P adalah curah hujan, dan S adalah retensi potensial maksimum, semua dalam satuan mm. Nilai retensi potensial maksimum (S, mm) dapat dikorelasikan dengan Persamaan 4.

(4) Sementara itu, untuk menghitung evapotranspirasiS nyata dengan metode Thornthwaite yang mengacu pada Persamaan 5 (Seiler & Gat, 2007).

(5) ET Keterangan Persamaan 5 sebagai berikut; ETr: evapotranspirasi nyata (mm/tahun); P: curah hujan (mm/tahun); Tm: temperatur rata-rata tahunan (°C). Menurut Lerner (1990), metode perhitungan imbuhan airtanah untuk daerah tropis dengan curah hujan tinggi menggunakan rumus sebagai berikut;

U = P – ETr – Ro (6)

Keterangan Persamaan 6, sebagai berikut; U: imbuhan airtanah (mm/tahun); P: curah hujan tiap tahun (mm/tahun); ETr: evapotranspirasi nyata (mm/tahun); Ro: air limpasan akibat efek dari impermeabialitas tanah (mm/tahun). Hasil perhitungan imbuhan airtanah kondisi alami dapat dilihat pada Tabel 1.

Tabel 1: Prediksi Imbuhan Airtanah Daerah Model Penelitian Pada Kondisi Alami Parameter Total Satuan Curah hujan (P) 2.665,1 mm th-1 -1 Evapotranspirasi aktual (ETr) 1.180,5 mm th Limpasan air permukaan (Ro) 549,33 mm th-1 Imbuhan airtanah (R) 935,27 mm th-1

72 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

3.3. Hidrogeologi Daerah Model Penelitian Hidrogeologi daerah penelitian ditentukan oleh litologi, keterdapatan sumber air, dan karakteristik akuifer. Berdasarkan klasifikasi Mandel dan Shiftan (1981) dan oleh Puradimadja (1993) disesuaikan dengan kondisi geomorfologi dan geologi indonesia, maka daerah penelitian termasuk dalam Sistem Akuifer Batuan Sedimen Terlipat (Irawan & Puradimaja, 2013). Sedimen terlipat merupakan lapisan sedimen yang dipengaruhi struktur perlipatan.

D A' Sungai Lawa B'

C' Sungai Perak 200 m E 100 m 0 m -100 m Sungai Lawa -200 m Sungai Perak -300 m D' A

Sungai Perak

C E'

D A' B' K e t e r a n g a n: PETA HIDROGEOLOGI Batas daerah model penelitian C' Batas daeragdaerah target penelitian E Jaringan sungai E 00.5 1 2 D' Simbol akuifer: A Kilometer Akuifer : Muara Lawa B Akuitar Kecamatan Kabupaten : Kutai Barat Akuiklud Propinsi : Kalimantan Timur C E' Gambar 3: Hidrogeologi Daerah Penelitian

Daerah model penelitian memiliki dua jenis media akuifer, yaitu akuifer media rekahan dan pori. Akuifer media rekahan dijumpai pada lapisan batubara yang tersisipi pengotor yang membentuk retakan-retakan batubara (cleat) yang berfungsi sebagai tempat aliran airtanah. Retakan-retakan ini membentuk rangkaian retakan yang sejajar yang biasanya berorientasi tegak lurus perlapisan. Pengotor yang menyisipi rekahan terdiri dari unsur mineral atau karbonat, lempung, jenis sulfida atau sulfat yang secara umum nampak pada permukaan batubara yang mengelupas. Sementara itu, akuifer media pori banyak dijumpai pada batupasir dengan rongga-rongga antar butir yang rapat yang berfungsi sebagi media aliran airtanah.

Pengamatan karakteristik airtanah dapat dilakukan berdasarkan pengamatan pada lokasi kemunculannya di permukaan yang berupa munculnya mata air (artesis) dari bekas sumur bor. Tipe mata air yang berkembang di daerah penelitian adalah mata air rekahan (fracture artesian springs) pada perpotongan lapisan batupasir dengan batulanau/batulempung/batubara, serta lapisan batulanau pasiran/batupasir lanauan.

Deformasi yang berupa lipatan sinklin, membentuk akuifer batuan sedimen terlipat yang berakibat kecilnya potensi airtanah di daerah penelitian. Di sisi lain, batuan penyusun yang didominasi oleh batulanau, batulempung pasiran, dan batubara yang mempunyai sifat semi kedap air, menambah makin kecilnya potensi airtanah yang ada. Hasil interpretasi dan analisis log pemboran, uji pemompaan disimpulkan, bahwa lapisan akuifer, yang berfungsi sebagai lapisan pembawa air, dibatasi dan ditutupi

73 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018 oleh lapisan penyekat pada bagian atas oleh batuan yang relatif bersifat semi kedap/semi permeabel, yaitu lapisan batulanau. Keterdapatan batupasir yang mayoritas berupa sisipan yang sangat kompak dan berumur tua sebagi akibat dari proses tektonik menyebabkan lapisan batupasir ini kemungkinan kecil berfungsi sebagai akuifer yang potensial (tinggi).

Kedalaman akuifer semi tertekan mempunyai kedalaman yang bervariasi, namun akuifer yang berpotensi sebagai akuifer produktif dikelompokkan menjadi dua, yaitu (1) akuifer kedalaman kurang dari 40 m, dan (2) akuifer dengan kedalaman lebih dari 40 m. Kedalaman akuifer kurang dari 40 m berada di daerah akuifer atas yang tersebar dari utara pada dataran dan punggung/sayap sinklin sampai ke selatan daerah penelitian. Sementara itu, akuifer kedalaman di atas 40 m menempati wilayah tengah daerah penelitian pada bagian bawah poros/sumbu sinklin. Ketebalan akuifer yang relatif besar pada bagian tengah lipatan dan mempunyai penyebaran yang luas memberikan cadangan airtanah yang baik. Walaupun demikian, hal ini sangat dipengaruhi juga oleh jumlah imbuhan airtanah yang masuk ke dalam akuifer.

3.4. Model Konseptual Daerah Model Penelitian Model konseptual adalah representasi sederhana dari model alami dari sistem hidrogeologi fisik dan perilaku hidrolik. Model ini merupakan dasar dari model analitik dan numerik yang diformulasikan untuk meniru kondisi asli di lapangan. Model mendasarkan pada kondisi alami daerah penelitian, dan faktor-faktor yang mempengaruhi aliran airtanah, baik disebabkan oleh keadaan alami, seperti: daerah imbuhan, luahan, presipitasi, evapotranspirasi, pola aliran air (permukaan dan airtanah), dan akibat adanya kegiatan yang dapat merubah kondisi alami tersebut, misalnya dampak aktivitas penambangan open pit.

Daerah imbuhan airtanah daerah penelitian meliputi lapisan akuifer yang tersingkap di permukaan. Litologi lapisan akuifer didominasi oleh batupasir yang disisipi oleh batulanau pasiran, pasir, dan batubara. Fluktuasi imbuhan airtanah sangat dipengaruhi oleh besaran presipitasi di daerah penelitian. Sumber airtanah berasal dari air permukaan pada daerah rawa, jaringan sungai-sungai yang berada di daerah lapisan akuifer.

Pada Gambar 4, merupakan deskripsi daerah imbuhan dan daerah luahan airtanah daerah model penelitian. Daerah imbuhan digambarkan dengan poligon warna kuning dan daerah luahan poligon selain warna kuning. Berdasarkan sayatan A–A’ yang memotong tegak lurus dari struktur sinklin memperlihatkan bentuk dari lapisan imbuhan mengikuti lapisan akuifer. Begitu pula daerah imbuhan pada lapisan akuifer yang terlihat di sayatan B–B’.

74 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

PETA KENAMPAKAN ATAS DAN SAYATAN DAERAH IMBUHAN E 01 2 4 6 S. Jeluk B' Kilometer Kecamatan : Muara Lawa S. Perak Kabupaten : Kutai Barat S. Tunau Propinsi : Kalimantan Timur

S. Lawa K e t e r a n g a n: Batas daerah model penelitian Daerah imbuhan airtanah Sungai

Peta Situasi

114°0'0"E 115°30'0"E 117°0'0"E

Kapuas Hulu Kutai Timur A' Bontang 0°0'0" Murung Raya 0°0'0" Kutai Kartanegara Kutai Barat Lokasi Penelitian Samarinda

Gunung Mas Barito Utara Penajam PaserBalikpapan Utara Kapuas A Daerah imbuhan Daerah imbuhan Daerah imbuhan A' Pasir Pulang PisauBarito Selatan Selat Makasar 1°30'0"S Palangka Raya Barito TimurTabalong 1°30'0"S 200 Katingan 114°0'0"E 115°30'0"E 117°0'0"E 100 Sumber Data : 1. Peta Rupa Bumi Indonesia Lembar Penarong (1815-13), 0 3 x eksagerasi vertikal skala 1 : 50.000, Bakosurtanal, 1991 2. Peta Topografi PT. TCM, 2013 3. Peta Geomorfologi Lembar Longiram, skala 1 : 100.000 Pusat Penelitian dan Pengembangan Geologi, Bandung, 1996 B B' 4. Peta Geologi Lembar Longiram (1815), skala 1 : 250. 000, Pusat Penelitian dan Pengembangan Geologi, Bandung, 1994 S. Perak Daerah imbuhan Daerah imbuhan 200 S. Lawa 5. Mine Plan Dept. PT. TCM, 2013 50 6. Peta SRTM 100 Disusun oleh : 0 Shalaho Dina Devy 3 x eksagerasi vertikal Jurusan Teknik Geologi Universitas Gadjah Mada Yogyakarta 2 0 1 5 Gambar 4: Daerah Imbuhan dan Luahan Daerah Penelitian

KONDISI ALAMI KONDISI ALAMI DAN MODEL KONSEPTUAL A A' S. Perak S. Tunau S. Lawa SEBELUM PENAMBANGAN 200! DAERAH PENELITIAN 150! 100! Kecamatan : Muara Lawa 50! Kabupaten : Kutai Barat 00.5 1 2 Propinsi 0! : Kalimantan Timur Km -50! Keterangan: Geologi: 3 X Eksagerasi Vertikal Konduktivitas hidrolika (m dtk-1) Satuan batupasir Balikpapan Satuaan perselinagnbatupasir kuarsa Simbol batuan: K1 (Akuifer atas) 1,1 × 10-4 Satuan perselingan batulanau dan dan batulempung Pulaubalang Batupasir Batulempung K2 (Akuitar atas) 6,56 × 10-7 batupasir Pulaubalang Satuan batulempung Pamaluan Batulanau K3 (Akuifer tengah) 1,09 × 10-4 K4 (Akuitar tengah) 2,54 × 10-7 MODEL KONSEPTUAL KONDISI ALAMI K5 (Akuifer bawah) 1,2 × 10-4 A A' K6 (Akuitar bawah) 5,8 × 10-8 u S. Tunau S. Lawa -11 S. Perak a K7 (Akuiklud) 2 × 10 200! n u ! T 150 . -1 K1 S Curah hujan rata-rata (mm th ) 2.665,1 ! K6 100 K4 Evapotraspirasi (mm th-1) 1.180,5 K2 K3 K5 K7 50! Limpasan air permukaan (mm th-1) 549,33 00.5 1 2 -1 0! Imbuhan airtanah (mm th ) 935,27 Km -50! 3 X Eksagerasi Vertikal Keterangan: Batas pemisah airtanah (no flow boundary) Evapotranspirasi Batas tanpa aliran (no flow boundary) Limpasan air permukaan Batas muka air permukaan (river boundary) Muka piezometrik Arah aliran airtanah K Konduktivitas hidrolika Presipitasi Batas daerah target penelitian

115°36'50"E 115°40'40"E 115°44'30"E a aw Peta Situasi K7 S. L 114°0'0"E 115°30'0"E 117°0'0"E K5 Kapuas Hulu Kutai Timur Bontang K6 0°0'0" Murung Raya 0°0'0" 02 4E 8 12 Kutai Kartanegara K4 Kutai Barat Lokasi Penelitian Samarinda Km K3 Gunung Mas Barito Utara Penajam PaserBalikpapan Utara K e t e r a n g a n: Kapuas 0°42'40"S 0°42'40"S K2 Pasir Pulang PisauBarito Selatan Selat Makasar 1°30'0"S

1°30'0"S Batas daerah target penelitian Palangka Raya Barito TimurTabalong A' Katingan Jaringan sungai A K1 114°0'0"E 115°30'0"E 117°0'0"E k Sumber Data : Daerah imbuhan airtanah ra e 1. Peta Rupa Bumi Indonesia, skala 1 : 50.000, Bakosurtanal P 2. Mine Plan Dept. PT. TCM, 2013 Arah aliran air sungai . S 3. Peta Topografi PT. TCM, 2013 Akuifer

Akuitar 0°46'5"S Disusun oleh : 0°46'5"S Shalaho Dina Devy Akuiklud Jurusan Teknik Geologi Universitas Gadjah Mada Yogyakarta 115°36'50"E 115°40'40"E 115°44'30"E 2 0 1 5 Gambar 5: Model Konseptual Sebelum Penambangan (Kondisi Alami) Daerah Penelitian

Gambar 5 merupakan model konseptual kondisi alami daerah model penelitian. Kondisi batas hidrolika yang menjadi batas hidrogeologis daerah model penelitian bagian permukaan meliputi: (1) batas pemisah airtanah (groundwater divide) berada di daerah pegunungan bagian utara, dan selatan, (2) batas muka air permukaan internal (internal head controlled boundary) yang ditandai dengan keberadaan

75 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018 sungai Lawa di timur dan sungai Perak di sebelah barat daerah model penelitian, dan (3) batas aliran airtanah ke luar (outflow boundary) di daerah barat dari daerah model penelitian yang dibatasi oleh sungai Perak. Selain batas permukaan, terdapat kondisi batas hidrolika secara vertikal antara lain: (1) batas tanpa aliran (internal zero-flow boundary) yang berada di bawah lapisan akuiklud dari satuan batulempung formasi Pamaluan, (2) batas sisi timur berupa batas air permukaan internal (internal head controlled boundary) yang ditandai dengan keberadaan sungai seperti sungai Lawa, Perak, Jeluk dan Tunai, dan (3) batas sisi utara dan selatan berupa batas pemisah airtanah (groundwater divide) yang berada pada puncak pegunungan.

3.5. Penambangan Batubara Open pit Daerah Model Penelitian Metode penambangan open pit merupakan kegiatan pengambilan bijih/batubara melalui pengupasan tanah/batuan penutup yang berada di atasnya dan berhubungan langsung dengan udara terbuka. Perubahan litologi akibat penambangan open pit berdampak terganggunya akuifer, misalnya terjadi perubahan muka airtanah, pola aliran airtanah, kualitas-kuantitas airtanah dan sumber airtanah. Di samping itu, terganggunya akuifer mengakibatkan terbentuknya void atau pit lakes akibat lubang bekas pit yang terisi oleh air (permukaan dan airtanah).

Pada perencanaan penambangan, pit tambang mempunyai elevasi terendah mencapai -60 m. Selain itu, terdapat kenaikan elevasi pada daerah disposal dengan kenaikan hingga 30 m. Penurunan elevasi dari pit menyebabkan terpotongnya lapisan pembawa air potensial (akuifer). Kondisi ini makin terlihat apabila elevasi bottom layer dari void menembus lapisan akuifer yang menyebabkan masuknya aliran airtanah ke dalam pit. Semakin dalam pit, maka beresiko terhadap masuknya airtanah yang berkonsekuensi tidak ditambangnya batubara potensial karena tenggelam atau membentuk void.

3.6. Pemodelan Airtanah dan Zona Neraca Airtanah pada Kondisi Alami Daerah Model Penelitian Pemodelan airtanah digunakan untuk mensimulasi dan memprediksi kondisi akuifer dan perubahannya. Manfaat dari pemodelan ini adalah untuk mempresentasikan keadaan airtanah, seperti pola aliran airtanah, kuantitas, dan kualitas airtanah, pada waktu sekarang dan akan datang, khususnya dampak penambangan batubara open pit. Pemodelan airtanah daerah penelitian sangat dipengaruhi oleh kondisi: litologi/akuifer dan geometrinya, head hidrolika, batas hidrogeologis, serta diskretisasi model. Pemodelan pada penelitian ini menggunakan perangkat lunak Visual Modflow.

Pemodelan airtanah diawali dengan pengumpulan data-data mengenai kondisi geologi dan hidrologi pada cekungan airtanah. Pada penelitian ini digunakan metode numerik dengan pendekatan metode beda hingga (finite difference). Metode ini didasarkan dalam membagi daerah model dengan wilayah-wilayah domain dengan grid seragam yang disebut sebagai diskretisasi model. Aspek mendasar model numerik adalah representasi sistem hidrogeologi dengan cara mendiskretisasi model dengan dimensi tertentu. Diskretisasi daerah model penelitian menggunakan bentuk ortogonal grid dengan metode beda hingga (finite difference). Keakuratan dalam analisis dipengaruhi oleh ketelitian dalam penentuan dimensi grid yang mewakili tiap kondisi batas hidrolika/hidrogeologi dan geologi, seperti batas muka air dan struktur batuan.

Pola aliran airtanah sangat dipengaruhi oleh lapisan litologi dan struktur dari batuan. Data-data mengenai batas-batas hidrologi, seperti head air permukaan, lapisan kedap (impermeable), daerah tanpa aliran (no flow); tata guna lahan atau pemanfaatan lahan untuk kegiatan tambang, data meteorologi dan kontur permukaan merupakan data-data pelengkap dalam memodelkan pola aliran permukaan secara konseptual.

Zone groundwater budget atau zona neraca airtanah merupakan kapasitas/kuantitas air yang terdapat dalam akuifer yang berasal dari neraca airtanah pada suatu zona/daerah yang mencakup komponen aliran dari daerah berdekatan dan lapisan-lapisan akuifernya. Pada hasil pemodelan komponen aliran ini merupakan aliran masuk dan keluar dalam suatu zona, seperti constant head, drains, recharge, evapotranspirasi, river leakage, stream leakage, dan general head. Zona neraca airtanah kondisi alami

76 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018 digunakan sebagai referensi perubahan zona neraca airtanah pada perubahan tata guna lahan dan hidrogeologi daerah penelitian dampak penambangan open pit.

Setelah ditentukan kondisi batas-batas hidrolika, yang telah diuraikan sebelumnya, kemudian ditentukan nilai dari data-data masukan untuk pemrograman yang didasarkan pada batas konseptual model. Konduktivitas hidrolika batuan (K) dihasilkan dari uji akuifer dengan metode uji pemompaan dan uji slug yang dilakukan pada 9 lokasi bekas lubang pemboran. Selain itu, dilakukan pula pengamatan muka airtanah piezometrik beserta head hidrolika dari 9 lubang bor bekas pemboran tersebut.

Data-data masukan yang diperlukan untuk pemodelan, yaitu batas daerah imbuhan, data sifat lapisan batuan/akuifer, data storativitas batuan/lapisan akuifer, data konduktivitas hidrolika (K), dan ketinggian head hidrolika. Secara singkat, nilai K batuan untuk data masukan dapat dilihat pada Tabel 2, serta kedalaman muka airtanah piezometrik dan head hidrolika observasi/kalibrasi pada Tabel 3.

Tabel 2: Data Masukan Pemodelan Aliran Airtanah Daerah Model Penelitian Kondisi Alami (Kalibrasi) Data Unit Nilai Lokasi K (konduktivitas hidrolika) m det-1 1,1 × 10-4 Akuifer atas (K1) 6,56 × 10-7 Akuitar atas (K2) 1,09 × 10-4 Akuifer tengah (K3) 2,54 × 10-7 Akuitar tengah (K4) 1,2 × 10-4 Akuifer bawah (K5) 5,8 × 10-8 Akuitar bawah (K6) 2 × 10-11 Akuiklud bawah (K7) - Constant head m 170 – 210 Batas pemisah airtanah - River, stream m 60 – 70 Sungai Lawa 40 – 70 Sungai Perak 60 – 70 Rawa, Sungai Jelau dan Tunau - Imbuhan airtanah mm th-1 935,27

Tabel 3: Nilai Head Hidrolika Observasi dan Kalibrasi Pemodelan Airtanah Kondisi Alami Sumur bor Elevasi Kedalaman muka Head- Head- No bekas permukaan airtanah observasi pemodelan pemboran (m dpl) piezometrik (m) (m dpl) (m dpl) 1 OW-1 168 58,8 107 109,2 2 OW-2 136,5 37,4 97,5 99,1 3 DD_13G 143,2 47,5 93 95,7 4 OW-4 139,2 38,6 99 100,6 5 OW-6 139,9 40,8 97,2 99,1 6 OW-8 140,8 45,8 92 95 7 OW-9 141,6 45,5 98 96,1 8 OW-10 127,1 33,5 95,1 93,6 9 OW-13 109,2 19,2 94,2 90

Hasil pemodelan mengenai perubahan ekuipotensial head dan arah aliran airtanah sebelum penambangan (kondisi alami) hasil kalibrasi dapat dilihat pada Gambar 6.

77 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

115°36'36"E 115°40'2"E 115°43'28"E u a a n w PETA ARAH ALIRAN AIRTANAH u a 180 T L . . KONTUR HEAD DAN ZONE BUDGET 170 S S 0 16 160 KONDISI ALAMI (KALIBRASI) 0 150 3 140 1 0 0 4 0 11 3 0 1 2 0 1 5 1 1 0 0 120 9 9

0 E

0 01.25 2.5 5 7

11 0

!. OW-1 0 1 !. OW-2 1 Kilometer 00 1 2

!. OW-4 0°42'9"S 0 Kecamatan: Muara Lawa 100 0 3 k 1 Kabupaten : Kutai Barat ra !. DD-13G 90 e !. OW-8 170 Propinsi : Kalimantan Timur P . 0 S 6 1 100 !. OW-13 K e t e r a n g a n : A A' !. Sumur bor bekas pemboran !. OW-9 Batas daerah target penelitian 8 !. OW-10 0 1 00 Arah aliran airtanah 6 120

0 90 !. OW_6 Ekuipotensial head hidrolika (kalibrasi) 120 Ekuipotensial head hidrolika lap. ke-3 1 40 10 Jaringan sungai 40 Inactive cells 60 1 20 Batas muka airtanah piezometrik 0°45'36"S 130

80 40 1 1 0 0 Neraca airtanah Zona 3 -1 160 (m hari ) 0 Peta Situasi 150 17 139.240 114°0'0"E 115°30'0"E 117°0'0"E

80.082 Kapuas Hulu Kutai Timur 5.870 Bontang 0°0'0" Murung Raya 0°0'0" Kutai Kartanegara Kutai Barat Lokasi Penelitian Samarinda 115°36'36"E 115°40'2"E 115°43'28"E Gunung Mas Barito Utara Penajam PaserBalikpapan Utara Kapuas

A A' Pasir Pulang PisauBarito Selatan Selat Makasar 1°30'0"S Palangka Raya 1°30'0"S ! S. Perak Barito TimurTabalong 200 Katingan 114°0'0"E 115°30'0"E 117°0'0"E 100150!

0 50 0

5 0 Sumber Data : 0 ! 6

80 1

0 3 1

0 100 1

1 1

0 0 1. Peta Rupa Bumi Indonesia, skala 1 : 50.000, Bakosurtanal 7

70 4

60 0 ! 1

2 2. Mine Plan Dept. PT. TCM, 2013 -10050 90 1 3. Visual Modflow (lisensi Teknik Geologi UGM) -2000!

-50! Disusun oleh : 3 X Eksagerasi Vertikal Shalaho Dina Devy Jurusan Teknik Geologi Universitas Gadjah Mada Yogyakarta 2 0 1 5 Gambar 6: Pola Aliran Airtanah dan Zona Neraca Airtanah Kondisi Alami (Kalibrasi)

Hasil pemodelan dengan Visual Modflow, seperti terlihat pada Gambar 6, dapat disimpulkan, nilai head hidrolika kondisi alami daerah penelitian mempunyai hasil bervariasi dengan tinggi maksimum 210 m dpl dan minimum 40 m dpl. Head hidrolika tinggi berada di daerah batas pemisah airtanah yang berupa perbukitan utara dan selatan serta head hidrolika rendah berada di rawa, Sungai Tunau, Jelau, Perak, dan Lawa. Penyebaran ekuipotensial head hidrolika secara rata-rata dari elevasi tinggi di daerah utara, selatan, dan timur, kemudian menuju ke elevasi rendah daerah barat, di Sungai Perak. Pola sebaran ekuipotensial head hidrolika pada lapisan akuifer teratas menggambarkan, ekuipotensial head hidrolika makin turun dan rapat ketika mencapai jaringan sungai, seperti Sungai Tunau dan Jelau. Kerapatan menurun ini menyebabkan aliran airtanah pada lapisan akuifer teratas mengalir menuju ke sungai yang mempunyai head hidrolika rendah dengan kuantitas terbatas. Jadi, akibat kondisi tersebut dapat disimpulkan, airtanah menjadi pengisi/pemberi air pada aliran air sungai atau sebagai gaining stream/effluent.

Zona neraca airtanah alami daerah model penelitian (Gambar 6) dibagi menjadi tiga zona, yaitu pertama, zona putih yang meliputi batas Sub DAS Tunau dan Jelau dengan batas hidrogeologis (daerah perbukitan utara dan selatan serta daerah Sungai Perak dan Lawa), dengan besar neraca airtanah 139.240 m3 hari-1. Kedua, zona kuning yang dibatasi oleh daerah Sub DAS Tunau dan Jelau dengan neraca airtanah sebesar 80.082 m3 hari-1. Ketiga, zona hijau yang dibatasi oleh kawasan tambang dengan neraca airtanah sebesar 5.870 m3 hari-1. Ketiga zona tersebut pada kondisi alami masih berupa kawasan hutan yang mempunyai nilai limpasan air permukaan kecil dan imbuhan airtanah yang besar.

Hasil kalibrasi head hidrolika observasi dan pemodelan, kemudian diplot pada diagram scatter, memperoleh nilai RMS sebesar 8,875 persen. Nilai ini menjelaskan, terdapat penyimpangan atau deviasi dari hasil model dengan kondisi alami yang masih wajar, yaitu kurang dari 10 persen. Penyimpangan ini disebabkan heterogenitas litologi/akuifer dan variasi konduktivitas hidrolika tiap lapisan akuifer yang mempengaruhi selisih/perbedaan head hidrolika observasi (observed head) dengan pemodelan (calculated head). Diagram scatter head observasi dan pemodelan dapat dilihat pada Gambar 7.

78 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

Gambar 7: Diagram Scatter Head Observasi (Observed Head) dan Pemodelan (Calculated Head)

3.7. Pemodelan Airtanah dan Zona Neraca Airtanah pada Penambangan Aktif Daerah Model Penelitian Tabel 4 menjelaskan, terdapat perubahan kedalaman muka airtanah piezometrik yang berkorelasi dengan perubahan head hidrolika, seperti pada sumur observasi bekas lubang bor OW-8 dari 95 m dpl pada kondisi alami menjadi 58 m dpl pada penambangan aktif. Tingginya perubahan kedalaman airtanah piezometrik OW-8 pada aktivitas penambangan dibandingkan pada kondisi sebelum penambangan disebabkan berada kawasan tambang aktif.

Tabel 4: Selisih Nilai Head Hidrolika Pemodelan Airtanah Alami dan Penambangan Aktif Head- Sumur bor Elevasi Kedalaman muka Head- tambang Selisih No bekas permukaan airtanah awal/alami aktif head (m) pemboran (m dpl) piezometrik (m) (m dpl) (m dpl) 1 OW-1 168 60,4 109,2 107,6 1,6 2 OW-2 136,5 39,5 99,1 97 2,1 3 DD_13G 143,2 53,2 95,7 90 5,7 4 OW-4 139,2 42 100,6 97,2 3,4 5 OW-6** 146,9 79,9 99,1 67 32,1 6 OW-8* 133,8 75,8 95 58 37 7 OW-9* 147,1 82,1 96,1 65 31,1 8 OW-10* 115,1 48,1 93,6 67 26,6 9 OW-13** 139,2 66,2 90 73 17 Keterangan: *) lokasi di daerah pit; **) lokasi di daerah disposal

Hasil pemodelan airtanah dari Gambar 8, memperlihatkan, terdapat tiga kawasan zona neraca airtanah pada penambangan aktif, yang pembagiannya seperti pada zona kondisi alami, yaitu zona putih mempunyai luas 145.500 m3 hari-1, zona kuning sebesar 88.241 m3 hari-1, dan zona hijau sebesar 31.241 3 -1 m hari .

79 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018

Berdasarkan hubungan antara ekuipotensial head hidrolika dan letak perlapisan akuifer, yang mengacu hasil pemodelan airtanah pada Gambar 8, dapat diambil kesimpulan, terdapat kecenderungan pola arah aliran airtanah yang mengarah ke pit, yaitu dari ekuipotensial head hidrolika teratas pada lapisan ke-1 hingga ke lapisan ke-4. Kerapatan ekuipotensial head hidrolika yang menurun pada daerah pit mengakibatkan arah aliran airtanah terfokus menuju ke pit. Namun, arah aliran airtanah ke pit kembali pulih pada pada lapisan ke-5 hingga ke-7 akibat adanya lapisan akuitar (lapisan ke-4) yang membatasi/menyekat lapisan ke-3 dan lapisan ke-5. Secara umum, arah aliran airtanah di luar daerah pit yang berasal dari lapisan akuifer pertama hingga terakhir mengalir dari timur ke barat.

115°36'27"E 115°39'54"E 115°43'21"E u a a n w 180 u a PETA ARAH ALIRAN AIRTANAH T L . . S S KONTUR HEAD DAN ZONE BUDGET PENAMBANGAN AKTIF 0 12

1 20 E 1 0 01.25 2.5 5 0 12 100 0 !. OW-1 !. OW-2 Kilometer

!. OW-4 0°42'9"S 0 Kecamatan: Muara Lawa

k 0 Kabupaten : Kutai Barat ra 6 e !. DD-13G Propinsi : Kalimantan Timur P !. OW-8 . S 60 0 40 K e t e r a n g a n : A !. OW-13 A' !. Sumur bor bekas pemboran 70 Batas daerah target penelitian !. OW-9 !. OW-10 6 Arah aliran airtanah 0 200 Ekuipotensial head hidrolika lap. ke-3 !. OW_6 40 Jaringan sungai 0 14 40 Inactive cells 60 80 Batas muka airtanah piezometrik Saluran terbuka 8 0 0 8 7 100 0 0°45'36"S

0 1 12 0 Neraca airtanah 0 Zona (m3 hari-1) 145.500 Peta Situasi 114°0'0"E 115°30'0"E 117°0'0"E

88.241 Kapuas Hulu Kutai Timur 31.241 Bontang 0°0'0" Murung Raya 0°0'0" Kutai Kartanegara Kutai Barat 115°36'27"E 115°39'54"E 115°43'21"E Lokasi Penelitian Samarinda

Gunung Mas Barito Utara Penajam PaserBalikpapan Utara A' Kapuas A Pit DA Pasir 200! S. Perak Pit Pit Pulang PisauBarito Selatan Selat Makasar 1°30'0"S Palangka Raya 1°30'0"S Barito TimurTabalong Katingan 100 ! 150 114°0'0"E 115°30'0"E 117°0'0"E

! 1000 0

0 0 Sumber Data :

1 5 50

0 0

50 50 ! 30 0 1 1. Peta Rupa Bumi Indonesia, skala 1 : 50.000, Bakosurtanal -100 7 50 7 7 10 1 2. Mine Plan Dept. PT. TCM, 2013 -2000! 3. Visual Modflow (lisensi Teknik Geologi UGM) -50 ! Disusun oleh : 3 X Eksagerasi Vertikal Shalaho Dina Devy Jurusan Teknik Geologi Universitas Gadjah Mada Yogyakarta 2 0 1 5 Gambar 8. Pola Aliran Airtanah dan Zona Neraca Airtanah Penambangan Aktif

4. Kesimpulan Hasil pengamatan di lapangan yang kemudian diplotkan dalam peta geologi, memperlihatkan, daerah model penelitian berada di struktur lipatan, yaitu Sinklin Lampanan. Posisi sumbu sinklin berada di antara lapisan batuan utara dan selatan yang mempunyai kemiringan saling berlawanan.

Kecenderungan pola arah aliran airtanah mengarah ke pit, yaitu dari ekuipotensial head hidrolika teratas pada lapisan ke-1 hingga ke lapisan ke-4. Terdapat perubahan kedalaman muka airtanah piezometrik yang berkorelasi dengan perubahan head hidrolika, seperti pada sumur observasi bekas lubang bor OW- 8 dari 95 m dpl pada kondisi alami menjadi 58 m dpl pada penambangan aktif. Tingginya perubahan kedalaman airtanah piezometrik OW-8 pada aktivitas penambangan dibandingkan pada kondisi sebelum penambangan disebabkan berada kawasan tambang aktif.

Selain itu, terjadi kenaikan neraca airtanah hampir lima kali lipat dari kondisi alami ke kondisi penambangan aktif, terutama pada zona hijau, yaitu dari kondisi alami 5.870 m3 hari-1 menjadi 31.241 m3 hari-1 pada penambangan aktif. Hal ini disebabkan, zona hijau adalah kawasan tambang terjadi perubahan tata guna lahan dan hidrostratigrafi, yaitu terbentuknya pit yang mencapai kedalaman 70 meter, yang mengakibatkan terpotongnya lapisan akuifer tengah, yang mempunyai penyebaran luas

80 SPECTA Journal of Technology, Vol. 2, No. 2, July - August 2018 dengan ketebalan hingga 21 meter. Lapisan akuifer tengah mayoritas berupa batupasir yang berasal dari satuan perselingan batulanau dan batupasir Pulaubalang. Akibat terpotongnya lapisan akuifer ini mengakibatkan turunnya tekanan hidrolika dan pengaruh dari gaya gravitasi yang menyebabkan debit aliran airtanah menjadi besar. Sementara itu, pada zona kuning dan putih terjadi kenaikan neraca airtanah secara wajar, yaitu rata-rata kenaikan dua kali lipat, yang disebabkan penurunan imbuhan airtanah dari kondisi alami sebesar 935,27 mm th-1 menjadi 601,71 mm th-1 pada penambangan aktif.

Ucapan Terima Kasih Penulis menyampaikan terimakasih kepada PT. Trubaindo Coal Mining tbk daerah operasi penambangan Muara Lawa di Kutai Barat, Provinsi Kalimantan Timur, yang telah memberikan bantuan teknis, supervisi, dan nonteknis untuk operasional penelitian lapangan, sehingga penelitian ini selesai sesuai dengan target rencana yang diharapkan.

Daftar Pustaka Arsyad, S., 1989, Konservasi Tanah dan Air, IPB Press, Bogor, hlm 85-86 Asdak, C., 1995, Hidrologi dan Pengelolaan Daerah Aliran Sungai, Gadjah Mada University Press, Yogyakarta, hlm. 154 – 210. Hammer, M.J. and Mac-Kichen, K.A., 1981, Hydrology and Quality of Water Resources, John Wiley & Sons, Ltd., New York, p. 41, 140, 214 Healy, R., W. and Cook Peter, G., 2002, Using Groundwater Levels to Estimate Recharge, Hydrogeology Journal, Vol. 10, No. 1, p. 91 – 107. Irawan, D.E., Puradimaja, D.J., 2013, Lembar Kerja Hidrogeologi Umum, Kelompok Keahlian Geologi Terapan Fakultas Ilmu dan Teknologi Kebumian Institut Teknologi Bandung Kabupaten Kutai Barat dalam Angka, 2011, BPS Kabupaten Kutai Barat Kamiana, I., 2010, Teknik Perhitungan Debit Rencana Bangunan Air, Graha Ilmu, Yogyakarta, hlm. 28-30, 203 Lerner, D.N., 1990, Groundwater recharge in Urban Areas, Hydrological Processes and Water Management in Urban Area, Proceedings of the Duisberg Symposium, April 1990, IAHS Publ., No. 198 Mandel, S., Shiftan, Z.L., 1981, Groundwater Resources: Investigation and Development, Academic Press. Inc, USA Notodarmojo, S., 2005, Pencemaran Tanah dan Airtanah, Penerbit ITB, Bandung. Puradimaja, D.J., 1993, Penyusunan Tipologi Paket Penelitian Sumber Daya Air, LAPI ITB-Departmen Transmigrasi, Bandung Seiler, K.P., Gat, J.R., 2007, Groundwater Recharge from Run-Off, Infiltration and Percolation, Springer, The Netherland, p. 75 Supriatna, S., Sukardi, Rustandi, 1995, Peta Geologi Bersistem, Lembar Samarinda, Kalimantan sekala 1:250.000, Pusat Penelitian dan Pengembangan Geologi Bandung. Van Bemmelen, R.W., 1949, The Geology of Indonesia, Government Printing Office, the Hague, Netherland, p. 328 – 360.