## ESTIMASI KANDUNGAN DO (DISSOLVED OXYGEN) DI KALI SURABAYA DENGAN METODE KRIGING

## Alan Prahutama
Jurusan Statistika, Fakultas Matematika dan Ilmu Pengetahuan Alam, Universitas Diponegoro, Semarang Alamat e-mail : alan.prahutama@gmail.com

## ABSTRAK
Kota Surabaya merupakan salah satu kota terbesar dengan pemukiman penduduk yang cukup padat. Kali Surabaya merupakan salah satu sungai terbesar di Surabaya. Peningkatan sektor industri, pedatnya pemukiman penduduk menyebabkan pencemaran air di Kali Surabaya. Pengukuran tingkat kebersihan air menggunakan DO ( Dissolved Oxygen ). DO  merupakan oksigen terlarut yang digunakan untuk mengukur kualitas kebersihan air. Semakin besar nilai kandungan DO menunjukan bahwa kualitas air tersebut semakin bagus. Kriging merupakan salah satu metode geostatistika untuk mengestimasi titik yang tidak tersempel dengan menggunakan unsur spasial pada lokasi yang tersempel. Salah satu estimasi titik didalam kriging menggunakan bobot. Penentuan bobot adalah dengan menggunakan model semivariogram.  Model yang digunakan yaitu model Gaussian. Hasil yang diperoleh bahwa kandungan DO Kali Surabaya di titik sesudah outlet PT. Suparma menunjukan kandungan DO sebesar 4.1171.
Kata Kunci :  DO, Kriging, Model Gaussian, Semivariogram.

## PENDAHULUAN
Air merupakan sumber daya alam yang menjadi kebutuhan penting bagi manusia. Oleh karena itu kebersihan air perlu dijaga. Pencemaran sungai menjadi faktor penting dalam permasalahan lingkungan. Sungai yang menjadi sumber air disuatu pemukiman, terkadang tercemar oleh limbah pabrik, limbah rumah tangga dan zat-zat berbahaya lainnya. Kota Surabaya merupakan ibukota provinsi Jawa Timur dengan luas wilayah 326.81 km 2 dengan kepadatan 7.568 jiwa per km 2 [2]. Tingginya kepadatan penduduk mengakibatkan pencemaran sungai menjadi perhatian. Salah satu sungai di Surabaya adalah Kali Surabaya, merupakan bagian dari Daerah Aliran Sungai Brantas yang mengalir sepanjang  41 km dari DAM Mlirip di Mojokerto melewati Gresik,
Sidoarjo dan berakhir di DAM Jagir Surabaya.
Pengukuran tingkat kualitas air dilihat dari oksigen terlarut (Dissolved Oxygen). Semakin tinggi kandungan Dissolved Oxygen (DO) semakin bagus kualitas air tersebut [8]. Pemantuan kualitas sungai di Kali Surabaya dilakukan di delapan titik sungai. Pengukuran kadar DO dilakukan hanya dibeberapa titik, sehingga mengakibatkan terbatasnya data mengenai kadar DO di sepanjang Kali Surabaya. Kadar DO antara titik lokasi yang satu dengan titik lokasi lainnya bisa berbeda. Oleh karena itu diperlukan metode statistika untuk mengestimasi titik lokasi  yang ditentukan. Salah satu metode yang digunakan untuk estimasi titik lokasi adalah Kriging [4].
Penenlitian-penelitian mengenai kriging sudah banyak dilakukan antara mengenai kandungan air dibawah
permukaan bumi [5], model kriging sebagai pendekatan global untuk simulasi [9], penaksiran kandungan cadangan bauksit di daerah Mempawah [10], pendekatan ordinary kriging pada kasus curah hujan di Kabupaten Karangasem, Bali [1].
Oleh karena itu peneliti tertarik melakukan penelitian yaitu mengestimasi kandungan DO di Kali Surabaya dengan metode Kriging. Titik lokasi yang akan diestimasi yaitu Kali Surabaya di sesudah outlet PT. Suparma. Hal ini dikarenakan lokasi tersebut dekat dengan industri.
Prinsip dasar jarak Euclid adalah mengukur jarak antar dua titik. Diberikan
2 titik , ( i i ) A x y dan , ( j j ) B x y dengan
merupakan jarak dari h
A ke B . Jarak euclid didapat dari h
2 2 ) j y − ( ) i j x = − ( + i h x y [3].
Kriging merupakan metode geostatistika yang digunakan untuk menduga atribut dari sebuah titik lokasi sebagai kombinasi linier dari nilai sampel yang terdapat di sekitar titik lokasi yang akan diduga [4]. Metode kriging menggunakan pembobotan untuk mengestimasi lokasi yang diduga. Bobot kriging diperoleh dari hasil variansi estimasi minimum dengan memperluas penggunaan semivariogram. Kriging memberikan lebih banyak bobot pada sampel dengan jarak terdekat dibandingkan dengan sampel dengan jarak lebih jauh. Nilai bobot Ordinary kriging dapat diperoleh melalui Persamaan sebagai berikut:
( ) ( ) ( ) ( ) ( ) ( ) ⎟⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜⎜ ⎜ ⎜ ⎜ ⎝ ⎛ ⎟⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜⎜ ⎜ ⎜ ⎜ ⎝ ⎛ = ⎟⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜⎜ ⎜ ⎜ ⎜ ⎝ ⎛ − 1 , , 0 1 1 1 , , 1 , , 0 0 1 1 1 1 1 1 1 s s s s s s s s s s s s n n n n n n γ γ γ γ γ γ
μ λ λ M
L L M M O M L M (1)  (1)
( , ) i j s s γ merupakan semivariogram dari lokasi ke-i ke lokasi ke-j, dengan
. 1, i j = = 2, , n K 0 ( , ) i s s γ merupakan semivariogram dari lokasi ke-i ke lokasi
yang diduga [7]. Pembobotan disetiap lokasi adalah i λ dengan μ merupakan pembatas atau Lagrange Multiplier.
Untuk mencari lokasi yang diduga
didapat dari
0 ˆ( ) z s ( ) ( ) ( ) ⎟ ⎟ ⎟ ⎠ ⎞ n s s M 1 ) ⎜ ⎜ ⎜ ⎝ ⎛ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎝ ⎛ = T n z z s z M 1 0 ˆ λ λ
(2)
dengan merupakan nilai atribut pada lokasi ke-i.
( i z s
Pada persamaan 2, untuk mendapatkan nilai semivariogramnya maka ditentukan model semivariogramnya. Model-model semivariogram adalah sebagai berikut:
1. Model Gaussian
( ) ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ − 1 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − + = 2 2 0 exp a h C C h γ 2. Model Spherical
( ) ⎪ ⎩ ⎪ ⎨ ⎧ ≥ ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ − + = h C h a h a h C C h ; ; 2 2 3 3 3 0 γ < a a 3. Model Exponential ( ) ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ − − + = a h C C h exp 1 0 γ 4. Model Linier ( ) Ch h = γ dengan merupakan nugget effect yaitu pendekatan semivariogram pada titik disekitar nol.  Sill ( 0 C ) 0 C C + merupakan Sebuah nilai tertentu yang konstan yang dimiliki oleh semivariogram untuk jarak tertentu sampai dengan jarak yang tidak terhingga. Range adalah jarak maksimum dimana masih terdapat korelasi antar data. Jarak euclid merupakan jarak antara lokasi yang satu dengan lokasi yang lain yang didasarkan pada garis lintang dan garis bujur.
( ) a
h
Sifat-sifat semivariogram adalah sebagai berikut [7]:
1. Semivariogram dari dua data yang
berjarak nol sama dengan nol
2. Nilai semivariogram selalu positif
3. Semivariogram adalah fungsi genap
Proses didalam kringing ada dua, yaitu proses  anisotropi dan isotropi. Proses anisotropi adalah proses spasial dimana dalam penyusunan struktur variogram bergantung pada jarak lurus antara suatu pasangan titik dan arah dari vektor yang menghubungkan titik-titik tersebut. Sedangkan  proses isotropi hanya bergantung pada jarak antara pasangan vektor dan struktur korelasi spasialnya sama untuk semua arah [4].
Oksigen terlarut atau Dissolved Oxygen (DO) merupakan salah satu parameter mengenai kualitas air. Tersedianya oksigen terlarut didalam air sangat menentukan kehidupan di perairan tersebut [6]. Menurut PP No. 82 Tahun 2001, baku mutu kandungan DO disungai adalah 6 Mg/L.

## METODE PENELITIAN

## Sumber Data dan Variabel Penelitian
Penulis melakukan analisis data terhadap data sekunder yaitu data kandungan DO di tujuh titik lokasi sungai Kali Surabaya. Data diambil dari Balai Lingkungan Hidup (BLH) Kota Surabaya tahun 2010. Variabel penelitian yang digunakan pada penelitian ini
ialah data garis lintang, data garis bujur, dan data nilai DO. Ketiga variabel tersebut diukur sebanyak 7 titik lokasi. Garis lintang dan garis bujur yang didapat dari lokasi-lokasi tersebut disajikan dalam Tabel 1.

## Metode Analisis
Langkah-langkah yang dilakukan di dalam pengolahan data sekunder untuk memprediksi nilai DO suatu titik lokasi adalah sebagai berikut :
1. Menentuan nilai nugget effect, sill , dan range menggunakan software easykrig . Proses kriging yang digunakan adalah proses isotropi. Model semivariogram yang digunakan adalah Model Gaussian.
2. Mengukur jarak antara dua titik menggunakan prisip dasar jarak Euclid
3. Berdasarkan nilai nugget effect, sill , dan range yang telah didapat maka menentukan nilai semivariogram antar lokasi.
4. Menentukan nilai bobot Ordinary Kriging .
5. Menentukan nilai DO pada lokasi yang diprediksi atau diduga menggunakan persamaan 2.

## Tabel 1 Garis lintang dan Garis Bujur di setiap Lokasi
No LOKASI GARIS LINTANG GARIS BUJUR 0 Kali Surabaya di sesudah outlet PT.Suparma 07 o 21’072” 112 o 40’461” 1 Kali Surabaya di Intake Jagir 07 o 18’015” 112 o 44’0” 2 Kali Surabaya di Dam Gunungsari 07 o 18’504” 112 o 43’097” 3 Kali Surabaya di intake Karangpilang 07 o 20’864” 112 o 40’931”
4 Kali Surabaya di sesudah pertemuan dengan Kali Tengah 07 o 21’073” 112 o 39’751”
5 Kali Surabaya di sesudah outlet PT. Miwon 07 o 22’921” 112 o 36’693”
6 Kali Surabaya di sebelum outlet PT. SAK 07 o 22’028” 112 o 36’232”

## HAS
D gamb koor kuali Salah air y kand peng tabel di Ka Tab Ta kand adala devia kand menu kand 82  T DO d
Sel
effec softw dari mode Gaus

## SIL PENEL
Data yang baran dat dinat lokasi itas kebersi h satu param yaitu kandun dungan D gestimasian.
l statistika d ali Surabaya
bel 2 Statistika
Banyak Mea Varia Standar D Med Minim Maksim abel 2 dungan DO ah 3,345 asi 1,247. B dungan D unjukan tid dungan DO,
Tahun 2001 di sungai ad anjutnya  m ct, sill , da ware easykr program  a el semivar ssian isotrop

## LITIAN
diperoleh  a ta yang i dari sebua han air di K meter  kuali ngan DO. O
DO dijadi
. Di bawah deskriptif  k a di tujuh ti a Deskriptif D DO k data an ansi Deviasi ian
mum mum menunjuk O di tujuh Mg/L den Berdasarkan DO pada dak memenu
, karena me
1  baku mu
dalah 6 Mg/
menentuan an range rig . Berdasa aplikasi dip riogram pi berikut in
mivariogram M
adalah seb menjelas ah pemanta Kali Suraba itas kebersi Oleh karena ikan uku h ini disaji kandungan itik lokasi.
Data Kandung 7 3,345 1,554 1,247 3,495 1,51 4,67 kan rata-r h titik lok ngan stand n nilai rata-r Tabel uhi  baku  m enurut PP N
utu kandun
/L.. nilai nug mengguna arkan kelua peroleh gra untuk mo ni : Model Gaussia
uah kan auan aya. han a itu uran ikan DO
an rata kasi dard rata 4.2 mutu No. ngan gget akan aran afik odel
Ber nug 1.4 sem terl lok das ada dan den
h Ha lok mo effe ran mo lok Sill yan me lok sem dal di dah Kri dar nila lok
1 2 3 4 5 6 λ λ λ λ λ λ μ ⎡ ⎤ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ = ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎣ ⎦
rdasarkan gget effect 4757, dan ra Sebelum mivariogram lebih dahu kasi denga sar jarak alah misalka n , ( j j ) B x y m
ngan B adal
( ) i j x x = −
asil perhitun kasi disajika
Berdasarka odel Gaussi fect bernilai nge adalah odel Gaussia
( ) = 4757 , 1 h γ
dengan h a kasi. Berdas l , Range da ng sudah enghitung  n kasinya. H mivariogram lam Tabel 4 Sebelum m suatu titik l hulu menen iging berda ri penyelesa ai bobot kasi ke- i ( i λ
875.478 848.457 848.457 1644.20 205.132 322.38 815.498 1078.32 − − − − − − = 161.773 473.27 951.797 1263.43 6.626 27.85 ⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ − − ⎣
Gambar 1 t adalah 0 ange sebesa menent m antara lok ulu jarak  a an menggu Euclid. Pr an  diberika
maka jarak lah
2 ) ( ) i j y y + −
ngan jarak an dalam Ta an grafik ian diperol 0, Sill sebe h 0.95 atau an adalah se
⎜⎜ ⎝ ⎛ ⎜⎜ ⎝ ⎛ − − 95 , 0 exp 1 7 h adalah jarak sarkan nila an jarak Eu didapat, nilai semiv Hasil per m antar l 4. mengestimas lokasi yang ntukan nilai asarkan pers aian persam Ordinary λ ) sebagai b
205.132 815.498 16 322.38 1078.32 47 202.28 338.01 19 338.0 − − − 1 949.98 3
195.42 342.42 2 196.34 939.43 2 1.96 11.76 7 − − − − − − −
diperoleh 0, Sill be ar 0.95. tukan kasi, mengh antara dua unakan pr rinsip das an 2 titik ( A
Euclid anta
. 2
Euclidian abel 3.
semivario leh nilai n esar 1.4757 u ditulis d ebagai berik
⎟⎟ ⎠ ⎞ ⎟⎟ ⎠ ⎞ 2 2 5 h k Euclidian ai Nugget  e uclid antar l maka variogram rhitungan lokasi disa si kandunga g diduga, ter bobot Ord samaan 1. maan 1 dipe Kriging u erikut :
61.773 951.797 6.62 73.27 1263.43 27.8 95.42 196.34 1.9 − − − 342.42 939.43 11.76 245.94 242.11 7.62 242.11 1230.01 18.0 7.622 18.05 0.7 − − − −
nilai ernilai
nilai hitung titik rinsip arnya , ( ) i i x y ara A antar ogram ugget
7, dan dalam kut :
antar effect , lokasi dapat antar nilai ajikan an DO rlebih dinary
Hasil eroleh untuk
26 85 96 0.01319 0.012233 0.09245 6 2
5 73 ⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦ 1 0.006659 0.10066 0.027524 1 ⎡ ⎤ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎣ ⎦
Setelah bobot didapatkan, selanjutnya mengestmasi nilai kandungan DO pada lokasi yang diduga. Lokasi yang diduga mempunyai garis lintang 07 o 21’072” dan garis bujur 112 o 40’461”. Nilai kandungan DO pada lokasi yang diduga merupakan perkalian antara bobot dengan nilai kandungan DO di lokasi yang lainnya. Untuk mencari Nilai kandungan 0 ˆ( ) z s ˆ( ) z s DO pada lokasi yang diduga sesuai dengan persamaan 2. Sehingga diperoleh nilai lokasi yang diduga adalah sebesar 1.58 Mg/L.
0
ˆ( ) z s 0 [2

## KESIMPULAN
Estimasi kandungan DO di titik Kali Surabaya sesudah outlet PT. Suparma menggunakan kriging adalah 1.58 Mg/L. Hasil estimasi kandungan DO, nilainya disekitar rata-rata data. Kandungan DO di titik tersebut dibawah baku mutu kandungan DO di sungai. Hal ini dikarenakan banyaknya industri disekitar Kali Surabaya yang menyebabkan Kali Surabaya tercemar limbah domestik.

## DAFTAR PUSTAKA
[1] Ayuni, N. W. D., 2010, Perbandingan Metode Ordinary Kriging dan Inverse Distance Weighted pada Interpolasi Spasial (Studi kasus: Curah Hujan di Kabupaten Karangasem) , Tugas Akhir Jurusan Matematika FMIPA Udayana, Bali.
] Badan Pusat Statistika (BPS), 2013, “Sosial dan Kependudukan”, www.surabayakota.bps.go.id diakses pada 19 November 2013.
[3] Cressie, N. A. C. (1993), Statistics for Spatial Data: Revised Edition , Inc. John Wiley and Sons, Canada.
[4] Issaks, E. H dan Srivastava, R.M.
(1989), Applied Geostatistics , Oxford University Press, New York.
[5] Kumar, V. dan Remadevi. (2006), “Kriging of Groundwater Levels- A case Study ”, Jornal of Spatial Hydrology Vol. 6, No. 1 Spring 2006.
Tabel 4 Nilai Semivariogram antar lokasi
Lokasi 0 1 2 3 4 5 6 0 0  0.01319 0.012233 0.092451 0.006659  0.10066 0.027524 1  0.01319 0 0.030032 0.168311 0.032447 0.161839 0.015751 2 0.012233 0.030032 0  0.08074 0.029279 0.057229 0.017174 3 0.092451 0.168311  0.08074 0 0.072827 0.032418 0.165843 4 0.006659 0.032447 0.029279 0.072827 0 0.107083 0.060403 5  0.10066 0.161839 0.057229 0.032418 0.107083 0 0.122078 6 0.027524 0.015751 0.017174 0.165843 0.060403 0.122078 0 Tabel 3 Jarak Euclidian antar lokasi Lokasi 0 1 2 3 4 5 6 0 0 0.090015 0.086674 0.241639 0.063889  0.25251 0.130354 1 0.090015 0 0.136222 0.330595 0.141652 0.323785 0.098411 2 0.086674 0.136222 0 0.225345 0.134485 0.188935 0.102786 3 0.241639 0.330595 0.225345 0 0.213719 0.141588  0.32801 4 0.063889 0.141652 0.134485 0.213719 0 0.260743 0.194211 5  0.25251 0.323785 0.188935 0.141588 0.260743 0 0.279159 6 0.130354 0.098411 0.102786  0.32801 0.194211 0.279159 0
[6] Mubarak, Satyari dan Kusdarwati, 2010, Korelasi antara Kosentrasi Oksigen Terlarut pada Kepadatan yang Berbeda dengan Skoring Warna
Daphnia Spp, Jurnal Ilmiah Perikanandan Kelautan, Vol.2, No.1
[7] Oliver, M. A and Webster, R. (2007), Geostatistic for Environmental
Statistic , Second Edition . Wiley,
United Kingdom.
[8] Simanjutak, M. (2007), Oksige Terlarut dan Apparent Oxygen Utilization di Perairan Teluk Klabat
Pulau Bangka. Jurnal ILMU KELAUTAN, Juni 2007. Vol 12 (2): 59-66 .
[9] Simpson et. al. (2001), “Kriging Models for Global Approximation in
Simulation-Based Multidisciplinary
Design Optimization”, AIAA Journal
Vol. 39, No. 12 December 2001 .
[10] Widhita, P. J. A. (2008),
Penaksiran Kandungan Cadangan Bauksit di Daerah Mempawah
Menggunakan Ordinary Kriging dengan Semivariogram Anisotropik ,
Tugas Akhir Departemen Matematika FMIPA Universitas Indonesia,
Depok.