## PERAMBATAN GELOMBANG DAN ARUS RESIDU PASANG SURUT TELUK MAYALIBIT: MODEL HIDRODINAMIKA PASANG SURUT 2D

## WAVE PROPAGATION AND RESIDUAL CURRENT OF TIDE IN MAYALIBIT BAY: 2D HYDRODYNAMICS MODEL OF TIDE
Asep Sandra Budiman 1* , Alan F. Koropitan 2 , dan I Wayan Nurjaya 2 1 Program Studi Ilmu Kelautan, Sekolah Pascasarjana IPB, Bogor
* E-mail: assandra81@gmail.com
2 Departemen Ilmu dan Teknologi Kelautan, FPIK-IPB, Bogor

## ABSTRACT
Study  on  tidal waves and  residual  currents  in Mayalibit Bay was  conducted by constructing the 2D numerical hydrodynamics to find the characteristics of tide and residual currents inside the Bay. The 2D hydrodynamic equations with non-linier terms were solved by finite difference methods explicitly. M2  and K1  wave propagations were  observed  and analyzed at 25  observation points. M2 and K1 residual currents were  obtained and described spatially.  Simulation results of tidal currents were validated  with  field measurements.  The validation result showed  that  the tidal currents  between the model results  and field measurements  were  quite fit.  Generally, simulation  results  showed  the significant differences between the tide inside and outside the Bay. This was strongly explained by the simulation results of M2 and K1 wave propagations. The M2 and K1 amplitudes were 10 times lower inside  than outside the Bay with the phase lags of about 180°(M2) and 160°(K1). This result indicated that the tide inside and outside the Bay were  in the opposite conditions.  The  tide  had an  ebb inside while the flood was outside and vice versa.  M2 residual currents was flowed into the Bay uniformly, while K1 residual currents flowed out to the open sea in the Northern part and it flowed into the Bay in the Southern part of the Bay. Eventually, both of them converged in the middle of the Bay.
Keywords : tide, tide-induced residual current, numerical method

## ABSTRAK
Studi dinamika massa air terkait pasut dan arus residu pasut di Teluk Mayalibit telah dilakukan dengan membangun sebuah model hidrodinamika 2D untuk mengetahui karakteristik pasut dan arus residu di dalam teluk.  Persamaan  Hidrodinamika  2D  dengan suku non-liniernya  diselesaikan  secara eksplisit melalui  metode  beda hingga  ( finite difference ).  Perambatan  gelombang M2  dan K1  diamati dan dianalisis di 25 titik pengamatan. Arus residu M2 dan K1 ditentukan dan digambarkan secara spasial. Validasi model dilakukan terhadap arus pasut. Hasil  validasi  menunjukkan  bahwa arus pasut antara hasil  model dan pengukuran  sudah cukup bersesuaian, baik amplitudo maupun fasanya.    Secara umum, hasil simulasi menunjukkan bahwa terdapat perbedaan yang signifikan antara pasang surut di dalam dan di luar teluk. Hal ini diperkuat dengan hasil simulasi perambatan gelombang pasut M2 dan K1. Amplitudo  M2 dan K1 10 kali lebih rendah  di dalam teluk dengan  keterlambatan fasa  sebesar 180 ° (M2) dan 160 ° (K1). Hasil simulasi model menunjukkan bahwa pasut di dalam dan di luar teluk berada pada kondisi yang berkebalikan dimana kondisi pasut di dalam teluk akan surut ketika di luar pasang dan begitupun sebaliknya. Arus residu M2 bergerak ke dalam teluk dengan pola yang seragam sementara  arus residu K1  mengalir  ke luar  menuju  laut terbuka di bagian Utara namun mengalir masuk  ke dalam teluk di bagian Selatan. Namun, pada akhirnya keduanya bertemu di bagian tengah teluk.
Kata kunci: pasut, arus residu pasut, metode numeric

## I. PENDAHULUAN
Teluk Mayalibit merupakan sebuah perairan semi -tertutup (hampir tertutup) di tengah-tengah pulau Waigeo Kabupaten Raja Ampat, Provinsi Papua Barat yang berada di posisi  koordinat  antara 0,01°-0,44°LS dan 130,53°-131,25°BT dengan luas total sekitar 4000 km 2 . Teluk ini dikelilingi oleh pantai- pantai bertebing yang dicirikan dengan relief sedang-tinggi, batu gamping putih, batuan beku, basal, masif, dan keras dengan tinggi tebing 2-100 m serta kemiringan 20% sampai terjal (PEMKAB RA,  2006) dan terhubung langsung dengan perairan lepas melalui suatu jalur yang panjang (38 km) dan sempit (McKenna et al., 2012).
Kepulauan Raja Ampat menjadi se- buah kabupaten maritim otonom (kabupaten bahari) pada tahun 2003. Kondisi ini mem- buka peluang untuk pengelolaan lokal yang lebih besar terhadap sumberdaya laut yang kaya. Di sisi lain, pemanfaatan sumberdaya alam yang tidak terkendali di daratan dapat berdampak negatif pada keberlanjutan eko- sistem laut. Teluk Mayalibit sebagai bagian dari wilayah perairan Raja Ampat mendapat dampak terbesar dengan begitu banyaknya kegiatan pembangunan di pesisir dan sifat teluk yang dangkal dan juga  semi  tertutup (Huffard et al., 2010).
Pola arus yang sangat mempengaruhi distribusi materi-materi di dalam kolom air, khususnya arus residu pasang surut. Arus residu pasang surut memiliki peranan yang sangat penting dalam proses dinamika estuari dan pesisir (Van Manh and  Yanagi,  2000). Besar dan arah dari aliran arus residu pasang surut akan menentukan pertukaran massa air dan proses penyebaran serta pengendapan dari berbagai materi, komposisi sedimen dan polutan di wilayah pesisir dan teluk dalam jangka panjang (Ramming and  Kowalik, 1980).
Penelitian mengenai perambatan ge- lombang pasang surut dan arus residu di wi- layah teluk dan perairan semi tertutup saat ini
masih berupa kajian-kajian yang terpisah atau tidak terintegrasi seperti penelitian me- ngenai dinamika pasang surut Laut Jawa oleh Yusuf and Yanagi (2013) melalui model 3D COHERENS,  kajian pola arus dan corange pasang surut Teluk Benete NTB (Pratama et al., 2012),  dan  model hidrodinamika pasang surut  yang ada di Pulau Baai Bengkulu (Su- piyati,  2005).    Penelitian-penelitian  tersebut masih terpisah dan hanya menitikberatkan pada  satu aspek pembahasan, pasang surut atau  arus residunya  saja.  Penelitian yang le- bih terintegrasi adalah seperti yang dilakukan oleh Fujiie et al. (2002) yang mengkaji pa- sang surut, arus pasang surut, dan transpor sedimen di teluk Manila melalui model hi- drodinamika 2D dan 3D)  sementara kajian serupa di perairan Teluk Mayalibit belum pernah dilakukan.  Kajian mengenai sirkulasi massa air dalam kaitannya dengan arus re- sidu pasang-surut di Teluk Mayalibit diper- lukan untuk memberikan gambaran kemung- kinan pola sebaran atau fate dari materi, polutan, atau zat-zat buangan bilamana ter- jadi pencemaran dimana informasinya sangat bermanfaat untuk pemerintah daerah setem- pat agar dapat mengelola wilayah perairan ini dengan baik. Salah satu metode  yang dapat digunakan adalah dengan aplikasi model numerik melalui pembangunan sebuah model hidrodinamika pasang surut 2D.  Penelitian ini bertujuan mengkaji perambatan gelom- bang  pasang surut dan arus residu pasang surut  di Teluk Mayalibit,  Kabupaten Raja Ampat, Papua Barat.

## II. METODE PENELITIAN

## 2.1. Waktu dan Tempat Penelitian
Wilayah  kajian dan daerah model adalah perairan Teluk Mayalibit, Kabupaten Raja Ampat, Papua yang secara geografis ter- letak pada posisi koordinat  antara 130,52 ° - 131,25 ° BT dan 0,01 ° - 0,44 ° LS dengan luas sebesar  4000  km 2 (lebar) (Gambar  1). Pene- litian yang meliputi pembangunan model dan juga analisis  yang dilakukan pada bulan Maret - Mei 2014.
Gambar 1. Wilayah kajian dan daerah model.

## 2.2. Data
Posisi stasiun ditentukan dengan GPS GARMIN, arus untuk validasi model diukur dengan RCM-7 di koordinat  130,9163 o BT dan 0,3109 o LS dengan penambatan ( moor- ing ) selama 20 jam selama EWIN P2O LIPI tahun 2008.  Penulisan program untuk pem- bangunan model dan simulasi dilakukan di PC OS Windows 7 dengan bahasa pem- rograman Fortran 90. Visualisasi output hasil model menggunakan Surfer 10 dan Matlab R2007b. Data yang digunakan untuk masuk- an model meliputi data batimetri dan pasang surut.  Data batimetri menggunakan peta laut no. 477 tahun 1996 dan no.186 tahun 2001 kepulauan Raja Ampat yang dikeluarkan DISHIDROS TNI AL. Pasang-surut meng- gunakan data elevasi dari NAOTIDE/NAO. 99b ( National Astronomical Observatory , Jepang) (Matsumoto et al., 2000). Suatu model peramalan pasang surut global dengan resolusi ½ o x ½ o dan merupakan data asimila- si dari TOPEX/Poseidon selama 5 tahun de- ngan 16 konstanta harmonik (M2, S2, K1, O1, N2, P1, K2, Q1, M1, J1, OO1, 2N2, Mu2, Nu2, L2, T2). Data lain yang dimasuk- kan dalam model selengkapnya disajikan pada Tabel 1.

## 2.3. Model

## 2.3.1. Persamaan Pembangun
Terdapat 4 jenis gaya yang bekerja pada massa air laut, yaitu gaya pada gradien tekanan, gaya coriolis, gaya gravitasi, dan
gaya friksi per unit massa (Ramming and Kowalik, 1980; Pond and Pickard, 1983; Stewart, 2002). Persamaan hidrodinamika diturunkan dari Hukum II Newton yang disebut hukum kekekalan momentum yang menyatakan bahwa perubahan momentum terhadap waktu sama dengan total gaya yang bekerja. Hukum ini dijabarkan dalam bentuk persamaan Hidrodinamika 2D sebagai berikut (Ramming and Kowalik, 1980): komponen-x:
u A z u k fv x p y u v x u u t u ∆ + ∂ ∂ + + ∂ ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 1 ρ ...........(1)
komponen-y:
v A z v k fu y p y v v x v u t v ∆ + ∂ ∂ + − ∂ ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 1 ρ ...........(2)
dimana: ∆ =  Laplace operator 2 dimensi
      ∂ ∂ + ∂ ∂ ∧ ∧ 2 2 2 2 y j x i = Kecepatan arah-x (m/det) = Kecepatan arah–y (m/det)
T =  variabel waktu (det)
f =  parameter Coriolis ( φ ϖ sin 2 = f )
^ ^ , j i =  unit vektor pada sumbu x dan y = koefisien viskositas Eddy lateral atau koefisien pertukaran momentum pada arah  horizontal
Tabel 1.  Komponen dan nilai hidrodinamika yang digunakan dalam pemodelan.
Komponen Hidrodinamika Nilai Jumlah sel/grid x 1450 Jumlah sel/grid y 890 Lebar sel/grid x ( ) 50 m Lebar sel/grid y ( ) 50 m Langkah Waktu 6 detik Waktu Simulasi 2,592,618 detik (+/- 30hari) Percepatan Gravitasi ( g ) 10 m/s Koefisien gesekan dasar ( r ) 0,003 Phi 3,141592654 Batas terbuka Selatan dan Timur (elevasi NAOTIDE/NAO99b) -16 Konstanta harmonik pasang surut untuk kajian perambatan gelombang pasang surut - 4 Konstanta harmonik pasang surut (M2, S2, K1, O1) untuk kajian arus residu M2, K1, dan M4 Batimetri Peta laut DISHIDROS 1996 dan 2001 = tekanan permukaan air laut (kg/m.det 2 ) =  percepatan gravitasi bumi (m/ det 2 ) ρ =  densitas air laut (kg/m 3 ) = koefisien viskositas eddy vertikal atau koefisien pertukaran momentum arah vertikal
Air laut diasumsikan sebagai fluida incompressible sehingga ditambahkan persa- maan kontinuitas ke dalam sistem persamaan yang ada di atas dalam bentuk (Ramming and Kowalik, 1980):
0 = ∂ ∂ + ∂ ∂ y v x u
.....................................(3)

## 2.3.2. Asumsi Model Hidrodinamika Pasang Surut 2D
Asumsi-asumsi yang diterapkan dalam persamaan model hidrodinamika pasang surut 2D ini antara lain (Pond and Pickard, 1983):
• Tekanan atmosfer permukaan ( a P ) adalah konstan, sehingga pada  turunan parsialnya
terhadap x dan y yaitu sama dengan nol.
      = ∂ ∂ = ∂ ∂ 0 y P x P a a ρ ρ …..…………………..(4)
• Pengaruh gaya coriolis terhadap gerak  pa- da    massa air diabaikan mengingat daerah model yang relatif kecil dan berada dekat dengan khatulistiwa.
• Diasumsikan tidak ada stratifikasi densitas air laut ( ρ konstan).
• Tidak ada sumber ( source ) dan kebocoran ( sink ) air laut yang terjadi di dalam daerah model, artinya evaporasi  dan presipitasi diabaikan serta dasar laut bersifat imper- meable.
• Tidak ada sumber momentum (gaya-gaya luar) yang terjadi pada area, seperti gerakan kapal, tsunami dan gempa.
• Batas tertutup tidak bergeser dengan naik- turunnya permukaan air laut.
• Pasang surut merupakan gaya pembangkit utama di Teluk Mayalibit. Faktor angin di- abaikan karena pengaruhnya relatif kecil.
Persamaan hidrodinamika laut yang digunakan berdasarkan asumsi-asumsi di atas adalah (Ramming and Kowalik, 1980):
komponen x:
u A H H x g fv y u v x u u t u x b x s ∆ + + − + + ∂ ∂ − = − ∂ ∂ + ∂ ∂ + ∂ ∂ ζ τ ζ τ ζ ) ( ) ( ............. (5) komponen y: v A H H y g fu y v v x v u t v y b y s ∆ + + − + + ∂ ∂ − = + ∂ ∂ + ∂ ∂ + ∂ ∂ ζ τ ζ τ ζ ) ( ) ( ..............(6) (5)
dimana ζ = elevasi (m), s τ = gaya gesekan permukaan, dan b τ = gaya gesekan dasar.
Persamaan (7) dan (8)  selanjutnya  di- integrasikan terhadap kedalaman  secara ver- tikal dari dasar (z = -H) sampai ke permuka- an (z = ζ ) untuk mendapatkan persamaan transpor massa sehingga diperoleh persaman kecepatan rata-rata dalam bentuk transpor pada  massa air 2D adalah sebagai berikut (Ramming and Kowalik, 1980).
komponen-x:
U A H V U rU W W W dx gH y U V x U U t U y x x ∆ + + − + + ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 / 1 2 2 2 2 ) ( λ ζ ................(7) komponen-y: V A H V U rV W W W dx gH y V V x V U t V y x y ∆ + + − + + ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 / 1 2 2 2 2 ) ( λ ζ ..............(8)
persamaan kontinuitas memiliki bentuk:
0 = ∂ ∂ + ∂ ∂ + ∂ ∂ t y V x
U ζ ................................... (9)
Suku gaya tekanan angin dalam kajian mengenai pasang surut dapat diabaikan (Bowden,  1983). Persamaan Hidrodinamika 2D pasang surut pada akhirnya memiliki bentuk:
komponen x:
U A H V U rU dx gH y U V x U U t U ∆ + + − ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 / 1 2 2 ) ( ζ ................(10) komponen y: V A H V U rV dy gH y V V x V U t V ∆ + + − ∂ − = ∂ ∂ + ∂ ∂ + ∂ ∂ 2 2 / 1 2 2 ) ( ζ ..........................(11)
persamaan kontinuitas:
0 = ∂ ∂ + ∂ ∂ + ∂ ∂ t y V x U ζ ...............................(12)
dimana: ∫ − ∂ = ζ 0 H z u U ; ∫ ∂ = ζ 0 H z v V =
kecepatan tranpor arah sumbu-x dan y ζ : elevasi muka air laut dari muka air laut rata-rata (MSL) H : kedalaman total perairan
r : Koefisien gesekan dasar Diskritisasi Persamaan Kontinuitas adalah  :
        ∆ − − + ∆ − − ∆ − = + y n j i V n j i V x n j i U n j i U t n j i n j i , 1 , 1 , , , 1 , ζ ζ .....................(13)
dimana U = H ; V = H adalah kecepatan transport arah-x,y (m 2 /dt). Suku-suku Non- Linier didiskritisasi dengan metode finite dif- ference menggunakan skema 2 langkah wa- ktu ( two-time level scheme ) (model GF2 oleh Crean et al., 2008 dalam Kowalik and Murty, 1993). Seluruh proses diskritisasi numerik persamaan hidrodinamika 2 dimensi secara eksplisit tersebut harus memenuhi kriteria stabilitas Courant-Freiderichs-Lewy (CFL) sebagai berikut:
max 2 gH L t ∆ ≤ ∆
…………………………(14)
Dimana: ∆ L =  min  [ ∆ x, ∆ y] dan H max =  max [d+ ς ].  Kriteria kestabilan dipenuhi dengan membuat model hidrodinamika teluk Maya- libit ini memiliki ukuran sel 55  m x 55  m dengan langkah waktu sebesar 6 detik.

## 2.4. Analisis Data
Penelitian ini secara garis besar terba- gi atas beberapa tahap yaitu tahap uji kesta- bilan dan validasi model, tahap menggambar- kan perambatan gelombang pasut secara spasial, dan tahap penentuan arus residu pa-
sut serta menggambarkannya secara spasial (Gambar 2).
Validasi dilakukan di titik validasi dengan koordinat 130,92 o BT dan 0,31 o LS (bulatan hitam pada Gambar 3) dengan mem- bandingkan komponen-komponen harmonik arus pasang surut antara hasil model dan pe- ngukuran  yang  melalui metode least square dengan bantuan yaitu  perangkat lunak Ttide (Pawlowicz et al ., 2002) .
Gambar 2.  Diagram alir pembangunan model.
Gambar 3. Titik validasi model.
Arah arus dominan  saat pasang dan surut dibandingkan dengan bantuan grafik pencar arus untuk melihat inklinasi atau sudut kemiringan antara hasil model dan pengukuran.  Perambatan gelombang pasang surut pada Teluk Mayalibit dilakukan dengan menggambarkan amplitudo dan juga fasa dari komponen harmonik semidiurnal M2 sebagai komponen yang dominan di wilayah ini be- serta komponen diurnal K1 yaitu secara spa- sial pada 25 titik grid di dalam domain model (Gambar  4).  Komponen harmonik pasang
surut M2, K1, M4, M6, M8, dan M10 di- peroleh melalui analisis Ttide terhadap hasil simulasi elevasi selama 20 hari.  Amplitudo dan fasa masing-masing komponen harmonik pada 25 titik grid kemudian diinterpolasi dan digambarkan secara spasial  sehingga dapat diketahui pola perambatannya.
Pola pergerakan arus residu M2 dan K1 diperoleh dengan mensimulasikan model se- lama 30 hari dimana yang menjadi tenaga penggerak arus ( forcing) adalah 4 komponen pasut yaitu M2, S2, K1, dan O1. Arus residu ditentukan dengan merata-ratakan hasil ke- luaran model yang berupa arus pada selang 10 periode yang terakhir yaitu untuk masing- masing komponen M2 dan K1. Arus residu M2 dan K1 yang telah diperoleh kemudian digambarkan secara spasial  dengan bantuan Surfer 10 agar dapat diketahui pola-pola per- gerakannya. Pola arus residu yang terjadi akan menentukan fate materi-materi di dalam teluk.
Arus residu dari 4 komponen pasang surut (M2, S2, K1, dan O1) ditentukan mela- lui persamaan-persamaan  (Ramming and Kowalik, 1980):
Gambar 4. Titik-titik grid untuk pengamatan perambatan gelombang pasang surut.
…………………..(15)
atau
…………….……(16)
dimana: adalah arus residu (m/s),  u dan v adalah komponen kecepatan arah- x dan arah- y (m/s), T adalah periode suatu kompon- en pasut (detik), H adalah kedalaman  per- airan (m), ζ adalah elevasi (m), dan t adalah waktu (detik).

## III. HASIL DAN PEMBAHASAN

## 3.1. Hasil Validasi Arus Pasang Surut
Hasil validasi model pada umumnya masih belum mendekati hasil pengukuran di lapangan. Hal ini dapat dilihat terutama dari selisih kecepatan arus pasang surut antara ke- duanya yang masih relatif besar yaitu antara 5,95-54,45 cm/s dengan beda fasa  berkisar antara 1,53-88,51 ° atau memiliki beda waktu sebesar 0,6 menit - 3 jam (Tabel 2). Hasil ini menunjukkan bahwa model belum cukup mendekati kondisi di lapangan.
Perbedaan kecepatan yang relatif besar ini dapat diakibatkan karena posisi titik validasi yang berada di kanal  atau jalur sempit de- ngan topografi yang rumit dimana faktor- faktor non-linier seperti turbulen atau gesek- an sangat intens sehingga mengakibatkan perubahan kecepatan arus yang cepat dalam ruang dan waktu (Pond and Pickard, 1983). Perbedaan kecepatan arus pasang surut yang relatif besar juga dapat diakibatkan karena perbedaan interval waktu yang digunakan dalam analisis Ttide dimana hasil simulasi model yang dianalisis adalah arus pasang surut selama 20 hari sedangkan hasil pengu- kuran hanya 20 jam saja karena keterbatasan data yang dimiliki.  Validasi hasil model dan pengukuran melalui analisis harmonik sudah sangat mencukupi karena proses  ini merupa- kan bagian terpenting dalam validasi sebuah model dan karena hasilnya dapat menunjuk- kan presisi suatu model (Bernardes, 2007). Titik validasi yang berada di dalam jalur teluk (Gambar  2)  yang  menyebabkan sulit memperoleh hasil model yang benar-benar mendekati kondisi di lapangan karena efek- efek non-linier yang relatif intens di titik ini. Validasi dilanjutkan dengan membuat grafik pencar ( Scatter plot ) arus pasang surut hasil model dan pengukuran (Gambar  5).
Tabel 2. Perbandingan komponen harmonik arus pasut hasil model dan pengukuran
Konstanta Arus (cm/s) Fasa (derajat) Selisih  Arus (cm/s) Beda Fasa(°) Beda Waktu Pengukuran *M2 81,29 270,21 54,45 88,51 3 jam 3 menit Model 26,84 181,7 Pengukuran *M4 16,6 178,84 15,99 49,02 51 menit Model 0,61 129,82 Pengukuran *M6 19,15 332,1 18,99 12,87 9 menit Model 0,16 319,23 Pengukuran *M8 8,16 241,37 8,04 1,53 0.6 menit Model 0,12 242,9 Pengukuran *M10 6,19 63,81 5,95 21,1 9 menit Model 0,24 42,71 Pengukuran *MS4 16,76 341,97 16,31 9,44 9 menit 36 detik Model 0,45 332,53
Gambar 5. Scatter plot arus pasang surut hasil model dan pengukuran.
Arus bergerak ke Barat Laut sampai Utara saat pasang atau masuk ke dalam teluk dan bergerak ke Tenggara sampai Selatan saat surut. Perbedaan kemiringan grafik yang dihasilkan antara model dan pengukuran ter- jadi karena perbedaan komponen arus yang digunakan. Scatter plot arus hasil model di- plot  yang  berdasarkan komponen-v  ( meri- dional: Utara-Selatan ) dan komponen-v ( zo- nal:  Timur-Barat )  sedangkan  s catter plot arus  pada  pengukuran diplot berdasarkan komponen–v saja karena hanya data  dalam
komponen tersebut saja yang tersedia. Grafik dibuat dalam orientasi Utara-selatan karena posisi titik validasi tepat berada di dalam kanal dalam arah tersebut (Gambar 6).
Perbedaan kemiringan pada scatter plot arus pasang surut antara hasil model dan pengukuran di titik validasi ini masih dapat ditolerir karena nilainya masih relatif kecil yaitu sebesar 12,85°. Selisih kecepatan arus saat pasang dan surut antara hasil model dan pengukuran masih relatif besar (Tabel  3).
Gambar 6. Titik validasi yang berada di kanal dengan orientasi utara-selatan.
Tabel 3. Perbandingan arus hasil model dan pengukuran di titik validasi.
Pasang Surut (dalam cm/s) Model Pengukuran Model Pengukuran Min 10 8,07 10 2,26 Max 57 151,92 59 54,28 Arah Utara Barat Laut-Utara Tenggara-Selatan Selatan
Hal ini dapat diakibatkan oleh posisi titik validasi yang berada di kanal  dimana dapat terjadi perubahan yang cepat terhadap nilai kecepatan arus disamping data hasil pengukuran yang diplot pada kondisi pasang purnama (10-11 Nopember 2011), sedangkan hasil model diplot selama 20 hari dimana kondisi pasang purnama dan pasang perbani masuk dalam perhitungan. Validasi kemudi- an dilanjutkan dengan menggambarkan arus pasang surut hasil model dan pengukuran dalam bentuk stick plot dimana waktu) dan interval dibuat sama (10-11 Nopember 2011, interval data 20 jam per-30 menit)  dengan orientasi Utara-Selatan (Gambar 7) . Arus pasang surut pada  hasil model memiliki kecepatan antara 5 –  89,4  cm/s (saat pasang) dan antara 4,5  –  75 cm/s (saat
surut) dengan keterlambatan sekitar 3 jam dari hasil pengukuran.  Arus pengukuran me- miliki kecepatan antara 8,07  –  151,92 cm/s (saat pasang) dan antara 2,26  –  54,28 cm/s. Hal ini menunjukkan bahwa hasil model le- bih mendekati hasil pengukuran ketika diplot pada waktu  yang sama meskipun masih memiliki keterlambatan fasa. Namun demi- kian, hasil ini dinilai yang terbaik yang dapat diperoleh setelah melalui berbagai proses penyesuaian komponen-komponen hidrodi- namika (koefisien gesekan dasar, batimetri smoothing) .  Vektor kecepatan hasil model tampak lebih miring dalam arah Barat Laut dan arah Selatan karena adanya komponen-u sedangkan pada pengukuran tidak digunakan karena  hanya pada komponen-v saja yang ada.
Gambar 7. Stick plot arus hasil model dan pengukuran di titik validasi selama 20 jam.
3.2. Pola Perambatan Gelombang M2 dan K1
Hasil  analisis Ttide terhadap data pa- sang surut perairan Teluk Mayalibit menun- jukkan bahwa pada  komponen semidiurnal M2 merupakan komponen yang paling  do- minan (amplitudo paling tinggi) di wilayah
ini. Komponen M2 dan K1 diolah dan Diana- lisis  lebih lanjut untuk dapat diketahui  peri- laku atau pola perambatannya di Teluk Mayalibit. Pola perambatan gelombang pasut M2 dan K1 dapat dilihat pada Gambar 8. dan 9 di bawah ini.
Gambar 8. Hasil simulasi co-amplitude M2.
Gambar 9. Hasil simulasi co-phase M2.
Pada umumnya, hasil simulasi co- amplitude M2    yang  menunjukkan bahwa amplitudo  M2 di dalam teluk mengalami reduksi sekitar  10 kali  (Gambar  8). Gelom- bang M2 mengalami reduksi atau pengurang- an ketika masuk ke dalam teluk akibat kanal atau jalur teluk yang panjang, sempit dan berkelok. Gelombang pasut yang datang dari luar teluk tersebut mula-mula tertahan dan terakumulasi di titik  4 yang memiliki topo- grafi berupa belokan atau tikungan yang tajam sehingga mengakibatkan tingginya ele- vasi dititik tersebut. Gelombang terus me- rambat  hingga masuk ke dalam teluk namun amplitudonya kini lebih rendah dari se- belumnya.  Gelombang M2 di dalam teluk mengalami keterlambatan fasa  sebesar 180 ° atau sekitar 6 jam 12 menit (Gambar 9). Hal ini menunjukkan bahwa terjadi kondisi yang berkebalikan antara pasut di dalam dan di luar teluk, ketika elevasi di luar teluk sedang pasang maka pada saat yang sama di dalam teluk sedang surut dan begitupun sebaliknya. Hal ini akan mengakibatkan perbedaan juga pada pola arus di dalam dan juga  di luar teluk, terutama arahnya. Di luar teluk, arus akan bergerak masuk ketika pasang, namun pada saat yang sama di dalam teluk arus akan bergerak ke luar karena saat itu sedang surut
dan begitupun sebaliknya.
Gelombang K1 yang dapat merambat masuk ke dalam teluk memiliki pola yang re- latif sama dengan gelombang M2 dimana amplitudonya menjadi 10 kali lebih rendah di dalam teluk dan mengalami keterlambatan fasa  sebesar 160 ° atau  juga  sekitar 12 jam (Gambar 10 dan 11). Hasil simulasi fasa pe- rambatan gelombang K1 (Gambar  11)  se- makin memperkuat terjadinya kondisi pasut beserta pola arus yang berkebalikan antara di dalam dan di luar teluk. Keberadaan jalur teluk yang panjang dan sempit diduga kuat memiliki peranan yang besar dalam mem- pengaruhi sirkulasi massa air di dalam teluk terkait  kondisi  pasang surut dan arus.  Hasil simulasi ini belum dapat dibuktikan kebena- rannya karena tidak adanya data pengukuran lapangan untuk elevasi dan arus di dalam dan di luar teluk disamping kajian mengenai ele- vasi dan arus Teluk Mayalibit ini terbilang baru karena belum adanya penelitian yang serupa sehingga tidak ada referensi yang dapat dijadikan acuan untuk menjelaskan ke- benaran  hasil simulasi model. Data elevasi dan arus di dalam dan di luar teluk dengan interval waktu yang lebih panjang (minimal 2-3 minggu) diperlukan untuk validasi model yang lebih baik.
Gambar 10. Hasil simulasi co-amplitudo K1.
Gambar 11. Hasil simulasi co-amplitudo K1.

## 3.3. Arus Residu
Arus residu komponen pasut M2 dan K1  ditentukan dengan mensimulasikan mo- del selama 30  hari dengan menggunakan 4 komponen pasut di batas terbuka sebagai forcing atau penggerak arusnya, yaitu M2,
S2, K1, dan O1. Arus residu M2 dan K1 ditentukan dengan merata-ratakan arus hasil simulasi selama 10 periode terakhir M2 dan K1. Hasil simulasi arus residu untuk masing- masing komponen M2 dan K1  dapat dilihat pada Gambar 12 dan 13 di bawah ini.
Gambar 12. Hasil simulasi arus residu M2 (cm/s).
Gambar 13. Hasil simulasi arus residu K1 (cm/s).
Arus residu M2 di dalam teluk ber- kisar antara 0,01 – 0,2 cm/s dan pada umum- nya lebih mengarah ke dalam  teluk  dengan pola yang seragam (Gambar  12). Pola arus residu M2 seperti  ini  dalam jangka waktu lama akan membawa masuk materi-materi di perairan semakin jauh ke dalam teluk.  Arus residu K1 di dalam teluk berkisar antara 0,01 – 0,15 cm/s dengan pola yang berbeda antara bagian Utara dan Selatan (Gambar 13). Arus residu K1 mengalir ke luar teluk menuju laut terbuka di sebelah Utara namun bergerak masuk  ke dalam teluk di sebelah Selatan hingga keduanya bertemu di bagian tengah teluk.  Pola arus K1  ini akan semakin men- jaga materi-materi  tetap berada di dalam teluk dan terakumulasi terutama di bagian te- ngah  teluk.  Keberadaan arus residu M2 dan K1  ini  harus diwaspadai bilamana terjadi pembuangan materi atau polutan yang ber- lebihan akibat aktifitas manusia di sekitar teluk karena akan membuat pencemaran yang semakin parah di dalam teluk. Hal ini di- karenakan besar dan arah dari aliran arus re- sidu pasang surut akan menentukan pertu- karan massa air dan proses penyebaran serta pengendapan dari berbagai materi, komposisi
sedimen dan polutan di wilayah pesisir dan teluk dalam jangka panjang (Ramming and Kowalik, 1980).

## IV. KESIMPULAN
Pada umumnya validasi model terha- dap pengukuran melalui perbandingan kons- tanta harmonik dan scatter plot belum me- nunjukkan hasil yang cukup baik dan men- dekati kondisi di lapangan karena perbedaan interval data yang digunakan, namun validasi melalui stick  plot lebih menunjukkan bahwa model cukup baik dan  mendekati pengukur- an karena selisih rentang kecepatan antara keduanya yang sangat relatif kecil (model: 5- 89,4 cm/s; pengukuran:  8,07-151,92 cm/s) meskipun masih terdapat keterlambatan fasa sekitar 3 jam.
Hasil simulasi menunjukkan perbeda- an signifikan antara kondisi pasang surut beserta arus di dalam dan di luar teluk. Jalur teluk yang panjang, sempit, dan berkelok di- duga kuat menjadi penyebab perbedaan ini. Kesimpulan ini diperkuat dengan hasil simu- lasi perambatan gelombang pasang surut M2 dan K1 yang datang dari perairan terbuka
dimana amplitudonya menjadi 10 kali lebih rendah di dalam teluk dan mengalami keter- lambatan fasa  sebesar  180 ° (M2) dan  juga 160 ° (K1). Hasil simulasi fasa gelombang M2 dan K1 tersebut menunjukkan bahwa terjadi kondisi  yang berkebalikan antara pasut di dalam dan di luar teluk, yaitu  ketika elevasi di luar teluk pasang maka pada saat keadaan yang sama di dalam teluk surut dan begitu- pun sebaliknya. Hal ini mengakibatkan kon- disi  yang berkebalikan juga  antara  pola arus di dalam dan di luar teluk.
Hasil simulasi arus residu M2 pada umumnya memperlihatkan pola arus yang bergerak  seragam  ke  dalam  teluk  dan dalam jangka waktu yang lama akan membawa ma- suk materi-materi semakin jauh ke dalam te- luk. Arus residu K1 di dalam teluk memiliki pola yang berbeda antara bagian Utara dan Selatan. Arus residu K1 mengalir ke luar te- luk menuju laut terbuka di sebelah Utara namun bergerak masuk ke dalam teluk di sebelah Selatan hingga keduanya bertemu di bagian tengah teluk. Pola arus K1 ini akan semakin menjaga materi-materi tetap berada di dalam teluk dan terakumulasi terutama di bagian tengah teluk. Keberadaan arus residu M2 dan K1 harus diwaspadai sebagai faktor yang dapat memperparah kondisi pencemar- an di dalam teluk bilamana terjadi pem- buangan materi atau polutan yang berlebihan akibat aktifitas manusia.

## UCAPAN TERIMA KASIH
Penulis mengucapkan terima kasih ke- pada Dr. Gentio Harsono, M.Si (DISHID- ROS TNI AL) atas penyediaan data batimetri Teluk Mayalibit, Kab. Raja Ampat Provinsi Papua Barat.

## DAFTAR PUSTAKA
Bowden,  K.F. 1983. Physical oceanography of coastal waters.  Ellis Horwood Li- mited. New York. 302p.
Dinas Hidro Oseanografi TNI AL. 1996. Peta
laut Indonesia No.  477. Jakarta. 186hlm.
Fuji-Ie, W., T. Yanagi,  and  F.P. Siringan. 2002. Tide, tidal current, and sedi- ment transport in Manila bay. J. La Mer., 40(1):137-145.
Huffard, C.L., J. Wilson, C. Hitipeuw, C. Ro- tinsulu, S. Mangubhai, M.V. Erdman, W. Adnyana,  P.  Barber,  J.  Manu- putty,  M.  Mondong,  G.  Purba,  K. Rhodes, dan H. Toha. 2010. Pengelo- laan berbasis ekosistem di bentang laut kepala burung Indonesia: mengu- bah ilmu pengetahuan menjadi tin- dakan. Ecosystem based management program:  conservation  international, the  nature  conservancy, and WWF Indonesia. 10hlm.
Kowalik,  Z. and T.S.  Murty. 1993. Numeri- cal  modeling of ocean  dynamics. World Scientific Publishing Co.  Pte. Ltd.  London. 481p.
Matsumoto, K., T. Takanezawa, and M. Ooe. 2000. Ocean tide models developed by assimilating TOPEX/POSEIDON altimeter data into hydrodynamical model: a global model and a regional model around Japan. Oceanography , 56:567-581.
McKenna,  S.A.,  GR.  Allen,  and  S.  Suryadi . 2002. A marine  rapid  assesment of the Raja Ampat island, Papua Provin- ce, Indonesia. RAP Bulletin of Biolo- gical Assesment 22 . Conservation In- ternational, Washington. 18p.
Nugraha,  R.B.H. dan  H. Surbakti.  2004.
Simulasi  pola  arus  dua dimensi di perairan teluk pelabuhan ratu pada bulan  September 2004. J.  Kelautan Nasional , 40(1):48-55.
Pawlowicz,  R.,  B.  Beardsley,  and  S.  Lentz.
2002. Classical tidal  harmonic  analy- sis including error estimates in MAT- LAB using T_TIDE. J.  of  Computa- tional Geosciences, 28:929-937. Pemerintah Kabupaten Raja Ampat. 2006. Atlas  sumberdaya  pesisir  kabupaten Raja Ampat. Papua Barat. 6hlm.
Pokavanich, T., I. Polikarpov, A. Lennox, F. Al-Hulail, T. Al-Said, E. Al-Enezi, F. Al-Yamani, N. Stokozov, dan B. Shu- haibar. 2013. Comprehensive  inves- tigation of summer hydrodynamics and water quality characteristic of desertic shallow water body: Kuwait bay. J.  of  Coastal Dynamics ,  12(2): 1253-1264.
Pond, S. and G.L. Pickard. 1983. Introdutory dynamical oceanography. 2 nd ed. Per- gamon Press Ltd. Oxford. 329p. Pratama, T.R., E. Indrayanti, dan I.B. Praset- yawan . 2012.  Kajian pola arus dan c o-range pasang surut di teluk Benete Nusa Tenggara Barat. Oceanography , 1(1):111-120. Ramming, H.G. and  Z.  Kowalik.  1980. Nu- merical modelling of marine  hydro- dynamics: applications to dynamics
physical processes. Oceanography Series 26. Amsterdam. 381p.
Stewart, R. 2003. Introduction to physical oceanography.  Department of Oce- anography. A and M University. Texas. 352p.
Supiyati. 2005. Model hidrodinamika pasang surut di perairan Pulau Baai  Bengku- lu. J. Gradien , 1(2):51-55.
Van Manh, D. and T. Yanagi. 2000. A study on residual flow in the gulf of tong- king. Oceanography , 56:59-68.
Yusuf,  M.  dan  T. Yanagi. 2013. Numerical modelling of tidal dynamics in the Java Sea. J. of Coastal Marine Scien-
ce , 36(1):1-12.
Diterima : 29 Maret 2015 Direview : 25 Mei 2015 Disetujui : 13 Juni 2015