## Pola Hidrodinamika di Perairan Nunukan sebagai Usulan pada Permasalahan Abrasi Pulau-pulau Kecil
WULAN FITRY HARDHIYANTI 1 , YESSI NIRWANA KURNIADI 1 , EVA MUSTIKASARI 2 , YOGI NOVIADI 3
1 Jurusan Teknik Sipil Institut Teknologi Nasional
2 Pusat Penelitian dan Pengembangan Sumber Daya Laut dan Pesisir, Kementerian
Kelautan dan Perikanan
3 Pusat Penelitian dan Pengembangan Geologi Kelautan, Kementerian Energi dan Sumber Daya Mineral Email: wulanfitryh@gmail.com

## ABSTRAK
Perairan Nunukan yang berada di Provinsi Kalimantan Utara, merupakan salah satu wilayah terluar  Indonesia  yang  berbatasan  langsung  dengan  Negara  Malaysia.  Berada  dalam wilayah  terluar, Kabupaten  Nunukan  menyimpan  banyak  potensi dalam sektor kelautan dan perikanan. Sangat penting untuk menjaga zona ekonomi eksklusif pulau-pulau kecil terluar Indonesia dari abrasi. Besar abrasi yang mungkin terjadi dapat diketahui dengan melakukan  analisis  pemodelan  mengenai  hidrodinamika  wilayah  tersebut.  Pemodelan numerik dilakukan menggunakan perangkat lunak MIKE 21 dengan modul Hydrodynamic (HD).  Pemodelan  yang  dilakukan  terdiri  dari  grid  persegi  dan  mesh  terstruktur.  Hasil penelitian  menunjukkan  terdapat  beberapa  daerah  di  Kabupaten  Nunukan  terjadi peningkatan kecepatan arus antara 0,1625 m/s sampai 0,3 m/s yang berpotensi terjadi abrasi. Abrasi terjadi di barat daya Pulau Sebatik, Perairan Nunukan sepanjang 1.014,26 meter. Penelitian sebelumnya menunjukkan, ukuran butir di wilayah tersebut sebesar pasir kelanauan sebesar 0,01 mm. Terdapat hubungan antara kecepatan dan ukuran butir yang menyatakan wilayah tersebut memiliki potensi abrasi.
Kata kunci: hidrodinamika, pemodelan, arus, pasang surut, abrasi

## ABSTRACT
Sebatik Island is located in Nunukan sea at North Kalimantan Province. This island is one of the outermost area in Indonesia which is in border line with Malaysia. These outlying islands keep a lot of potential in marine, and fishery. It is important to keep Indonesia exclusive economic zone by protecting this outlying island from abrasion. Abrasion can be found by analyzing hydrodynamic modelling in the area. Hydrodynamic modelling is using software  MIKE  21  with  Hydrodynamic  modules  (HD).  The  numerical  hydrodynamic modelling was consists of rectangular grid and unstructured mesh modelling. The modelling showed there are some area in Kabupaten Nunukan are has increase of current velocity between 0.1625 m/s and 0.3 m/s that potentially to abrasion. Abrasion occurred in the southwest of Sebatik Island, Nunukan sea along 1,014.26 meters. Based on the observation before, the grain size in the area was 0.01 mm sand silt. There is correlation between the current velocity and grain size in the area has potentially to abrasion.
Keywords: hydrodynamic, modelling, current, tides, abrasion

## 1. PENDAHULUAN
Langkah  mendasar  dalam  mewujudkan  pengamanan  garis  pantai  di  Kabupaten  Nunukan adalah  perlu  adanya  informasi  mengenai  hidro-oseanografi  di  wilayah  tersebut.  Menurut Pariwono (2006) salah satu aspek penting untuk mengetahui dan memahami agar pengelolaan sumber  daya  dapat  dilaksanakan  dengan  tepat  yaitu  mengetahui  dinamika  perairannya. Hidrodinamika  di  suatu  perairan  dapat  diketahui  setelah  dilakukan  beberapa  pengamatan seperti pasang surut, batimetri dan pola arus di wilayah tersebut. Pengambilan data dilakukan pada wilayah di Kabupaten Nunukan, Kalimantan Utara tepatnya pada koordinat Utara 4,09615 Timur  117,71347.  Pemodelan  dilakukan  dengan  menggunakan rectangular  grid  dan unstructured mesh. Hasil pemodelan yang sudah di validasi dapat digunakan untuk analisis selanjutnya sehingga diketahui wilayah yang terdampak abrasi dilihat dari adanya peningkatan kecepatan arus.

## 2. TINJAUAN PUSTAKA

## 2.1. Abrasi
Abrasi merupakan erosi pada dasarnya. Erosi sendiri adalah proses pengikisan batuan, tanah, maupun  padatan  lainnya  yang  disebabkan  oleh  gerakan  air, es, atau  angin. Selain terkait masalah kerusakan lingkungan, kondisi di suatu wilayah juga turut mempengaruhi tingginya gelombang  air  laut  serta  tanah  sehingga  mengalami  erosi.  Abrasi  sendiri  memiliki  5  kelas dalam kerusakannya yaitu ringan, sedang, berat, amat berat, dan amat sangat berat.

## 2.2. Arus
Arus laut adalah suatu pergerakan massa air secara vertikal serta horizontal sehingga menuju suatu keseimbangannya dengan wilayah yang sangat luas. Arus dipengaruhi oleh beberapa faktor yaitu perbedaan permukaan air laut, angin, perbedaan temperatur, perbedaan salinitas dan kepadatan air, pasang surut, dan bentuk pantai. Arus salah satunya dipengaruhi oleh angin dan dipengaruhi oleh musim. Di Indonesia sendiri dikenal 4 musim, yaitu musim barat yang terjadi selama periode bulan Desember-Februari, musim peralihan barat ke timur yang terjadi selama periode bulan Maret-Mei, musim timur yaitu terjadi selama periode bulan Juni- Agustus, dan musim peralihan timur ke barat yaitu terjadi selama periode bulan September- November.

## 2.3. Pasang Surut
Pasang surut adalah fluktuasi muka air laut sebagai fungsi waktu karena adanya gaya tarik benda-benda  di  langit,  terutama  matahari  dan  bulan  terhadap  massa  air  laut  di  bumi. Perbedaan vertikal antara pasang tinggi dan pasang rendah disebut rentang pasang surut (tidal range). Jenis pasang surut berbeda-beda, tergantung letak wilayah tersebut. Terdapat empat jenis pasang surut yang ditunjukkan pada dimana bentuk kurva berbeda-beda dari besarnya pasang surut terhadap fungsi waktu. J enis-jenis pasang surut anatar lain pasang surut harian ganda (semi diurnal tide), pasang surut harian tunggal (diurnal tide), pasang surut campuran condong ke harian ganda (mixed tide prevailing semidiurnal), pasang surut campuran condong ke harian tunggal (mixed tide prevailing diurnal). Pada umumnya sifat pasang  surut  di  perairan  dapat  ditentukan  dengan  menggunakan  rumus  Formzahl  seperti pada Persamaan 1 berikut:
𝐹 = 𝐾 1 + 𝑂 1 𝑀 2 + 𝑆 2 …(1)
dimana:
𝑂 1 = unsur pasang surut tunggal utama yang disebabkan oleh gaya tarik bulan,
𝐾 1 = unsur pasang surut tunggal yang disebabkan oleh gaya tarik matahari,
𝑀 2 = unsur pasang surut ganda utama yang disebabkan oleh gaya tarik bulan,
𝑆 2 = unsur pasang surut ganda utama yang disebabkan oleh gaya tarik matahari.
Nilai Formzahl ditunjukkan pada Tabel 1 berikut:

## Tabel 1. Nilai Formzahl
(Sumber: Pratikto et al, 1996)

## 2.4. Perangkat Lunak MIKE 21
Mike  21  adalah  pemodelan  terhadap  pesisir  dan  perairan  laut.  Memasukkan  data-data pendukung seperti data batimetri, data arus, data sungai, dan besar sedimen yang terdapat di wilayah tersebut, maka akan didapatkan hasil tinggi muka air yang terjadi, kecepatan arus, debit, dan sedimentasi yang terjadi dalam jangka waktu tertentu. Modul yang digunakan pada MIKE 21 adalah hydrodinamic module (HD). Pada hydrodynamic module terdapat persamaan yang menggambarkan aliran dan perbedaan muka air dengan Persamaan 2, Persamaan 3, Persamaan 4 berikut:
𝜕𝜁 𝜕𝑡 + 𝜕𝑝 𝜕𝑥 + 𝜕𝑞 𝜕𝑦 = 𝜕𝑑 𝜕𝑡 …(2)
𝜕𝑝 𝜕𝑡 + 𝜕 𝜕𝑥 ( 𝑝 2 ℎ ) + 𝜕 𝜕𝑦 ( 𝑝𝑞 ℎ ) + 𝑔ℎ 𝜕𝜁 𝜕𝑥 + 𝑔𝑝√𝑝 2 +𝑞 2 𝐶 2 .ℎ 2 − 1 𝜌 𝑤 [ 𝜕 𝜕𝑥 (ℎ𝜏 𝑥𝑥 ) + 𝜕 𝜕𝑦 (ℎ𝜏 𝑥𝑦 )] − Ω 𝑞 − 𝑓𝑉𝑉 𝑥 + ℎ 𝜌 𝑤 𝜕 𝜕𝑥 (𝑝 𝑎 ) = 0 …(3)
𝜕𝑞 𝜕𝑡 + 𝜕 𝜕𝑦 ( 𝑞 2 ℎ ) + 𝜕 𝜕𝑥 ( 𝑝𝑞 ℎ ) + 𝑔ℎ 𝜕𝜁 𝜕𝑦 + 𝑔𝑝√𝑝 2 +𝑞 2 𝐶 2 .ℎ 2 − 1 𝜌 𝑤 [ 𝜕 𝜕𝑦 (ℎ𝜏 𝑥𝑥 ) + 𝜕 𝜕𝑥 (ℎ𝜏 𝑥𝑦 )] − Ω 𝑝 − 𝑓𝑉𝑉 𝑦 + ℎ 𝜌 𝑤 𝜕 𝜕𝑥𝑦 (𝑝 𝑎 ) = 0 …(4)
dimana:
ℎ(𝑥, 𝑦, 𝑡) = kedalaman air [m], 𝑑(𝑥, 𝑦, 𝑡) = kedalaman air dalam berbagai waktu [m], 𝜁(𝑥, 𝑦, 𝑡) = elevasi permukaan [m], 𝑝, 𝑞(𝑥, 𝑦, 𝑡) = flux density dalam arah 𝑥 dan 𝑦 [m 3 /s/m], = (uh,vh); (u,v), = depth average velocity dalam arah 𝑥 dan 𝑦 ,
𝐶(𝑥, 𝑦)
= tahanan Chezy [m 1/2 /s], 𝑔 = kecepatan gravitasi [m/s 2 ], 𝑓(𝑉) = faktor gesekan angin,
𝑉, 𝑉 𝑥 , 𝑉 𝑦 (𝑥, 𝑦, 𝑡) = kecepatan angin dalam arah 𝑥 dan 𝑦 [m/s],
Ω(𝑥, 𝑦) = parameter Coriolis [s -1 ], 𝑝 𝑎 (𝑥, 𝑦, 𝑡) = tekanan atmosfer [kg/m/s 2 ],
No.  Bilangan Formzahl Jenis Pasang Surut 1 0,00 – 0,25 Ganda (semi diurnal) 2 0,26 – 1,50 Campuran dengan tipe ganda yang menonjol ( mixed, mainly semi diurnal) 3 1,51 – 3,00 Campuran dengan tipe tunggal yang menonjol ( mixed, mainly diurnal) 4 >3,00 Tunggal ( diurnal)
𝜌 𝑤 = berat jenis air [kg/m 3 ], (𝑥, 𝑦)
= kordinat ruang [m],
𝑡
= waktu [s],
𝜏 𝑥𝑥 , 𝜏 𝑥𝑦 , 𝜏 𝑦𝑦 = komponen effective shear stress [N/m 2 ].

## 3. METODOLOGI PENELITIAN

## 3.1. Lokasi Pengamatan
Lokasi penelitian berada di Kabupaten Nunukan, Kalimantan Utara dengan koordinat utara 4,09615 timur 117,71347 yang dilihat pada Gambar 1 .

## Gambar 1. Lokasi pengamatan
(Sumber: Google Earth, 2016)

## 3.2. Data Penelitian
Data  yang  digunakan  penelitian  didapat  dari  beberapa  sumber  yaitu  batimetri  hasil pengukuran oleh Tim Sebatik dari Pusat Penelitian dan Pengembangan Geologi Kelautan pada tahun 2005 dan data arus, pasang surut yang dilakukan Pusat Riset Kementerian Kelautan dan  Perikanan  pada  tahun  2016.  Untuk  data  perbandingan  batimetri,  digunakan  data batimetri yang telah di digitasi oleh GEBCO ( General Bathymetric Chart of the Oceans) yang dapat diunduh dengan data pada interval waktu tertentu.

## 3.3. Batas Penelitian
Batasan-batasan  dalam  penelitian  ini  antara  lain  tidak  memperhitungkan tidal  bore,  tidak memperhitungkan perubahan garis pantai, dan tidak memperhitungkan angkutan sedimen yang terjadi.

## 3.4. Alur Penelitian
Analisis  pendahuluan  dilakukan  dengan  mengolah  data-data  yang  diperlukan  seperti batimetri. Batimetri yang digunakan pada penelitian ini terdiri dari grid series (rectangular) dan mesh generator. Wilayah yang ditinjau merupakan selat yang dikelilingi oleh pulau dan berdekatan dengan wilayah Malaysia. Data batimetri berbentuk grid series atau rectangular harus di-digitasi terlebih dahulu dengan membuat grid dengan jarak antar grid sebesar 50 meter. Setelah data batimetri sudah ter-digitasi, data tersebut di plot pada grid series MIKE
21. Pemodelan menggunakan mesh generator, data batimetri dalam bentuk .xyz dimasukkan pada modul mesh Generator, lalu data tersebut dihubungkan antar noda menggunakan draw arc. Data batimetri telah memiliki bentuk baru dalam MIKE, selanjutnya membuat pemodelan menggunakan flow model sehingga dapat diketahui hidrodinamika di wilayah penelitian. Data kuat arus dan pasang surut hasil pengukuran digunakan untuk koreksi terhadap data yang dihasilkan  dari  pemodelan,  dengan  cara  membandingkannya.  Validasi  dilakukan  dengan metode  RMSE  ( Root  Mean  Square  Error)  dan  Correlation.  Apabila  data  pemodelan  telah tervalidasi  oleh  data  pengukuran,  maka  pemodelan  bisa  digunakan  untuk  dasar  analisis selanjutnya  yaitu  mengetahui  wilayah  yang  berpotensi  abrasi  dengan  melihat  perubahan kecepatannya.

## 4. PENGOLAHAN DATA DAN PEMBAHASAN

## 4.1. Data Validasi
Pasang  surut  dan  kecepatan  arus  didapat  dari  hasil  pengukuran  lapangan  menggunakan tidemaster dan ADCP. Data hasil pengukuran pasang surut oleh tidemaster dilakukan dengan pengamatan  15  menit  selama  22  hari.  Kecepatan  rata-rata  perairan  Nunukan  tergambar dalam grafik sebaran kecepatan  arus  hasil  pengukuran  ADCP ( Acoustic Doppler Current Profiler)   yang  ditunjukan  pada Gambar 2 . Terlihat bahwa kecepatan arus di perairan nunukan mengikuti pola geometri perairan yang mengarah ke timur laut-tenggara. Rata-rata kecepatan arus berkisar 40 cm/detik. Arah pergerakan kecepatan arus secara visual mengikuti geometri pantai.

## Gambar 2. Kecepatan arus rata-rata
(Sumber: Badan Penelitian dan Pengembangan Kelautan dan Perikanan, 2016)
Pada  pengukuran  yang  dilakukan  di  lokasi  Pengamatan  oleh  KKP  selama  22  hari  dengan interval waktu 15 menit didapat hasil dengan tenggang pasang surut diantara 2 hingga -2 yang diukur berdasarkan mean sea level yang dapat dilihat pada Gambar 3 . Jenis pasang surut yang terdapat pada wilayah tinjauan adalah pasang surut condong harian ganda (semi diurnal tide) dengan periode pasang surut rata-rata 12 jam.
-60 -40 -20 0 20 40 60 -60 -40 -20 0 20 40 60 K e ce p a ta n A ru s A ra h U ta ra , N [ cm /s ] Kecepatan Arus Arah Timur, E [cm/s] Diagram Sebaran Kecepatan Arus Rata - rata Nunukan
Gambar 3. Pasang surut pengamatan tidemaster (Sumber: Badan Penelitian dan Pengembangan Kelautan dan Perikanan, 2016)

## 4.2. Perbandingan Pemodelan dengan Data Validasi
Hasil output pemodelan harus tervalidasi terlebih dahulu sebelum mendapatkan pasang surut terprediksi agar data pemodelan sesuai data yang ada di lapangan. Validasi dilakukan dengan melakukan perbandingan hasil yang didapat dari pengukuran terhadap pemodelan dengan dua  metode,  yakni  RMSE  dan Correlation.  Pemodelan  1  merupakan  penamaan  dari pemodelan  yang  menggunakan  batimetri  dengan  modul grid  series  atau  rectangular, sedangkan untuk pemodelan 2 adalah penamaan untuk pemodelan dengan menggunakan batimetri dengan modul mesh generator. Berikut hasil perbandingan antara hasil pengukuran dan pemodelan pada Tabel 2 .
Tabel 2. Nilai RMSE dan Correlation
Dari hasil pada Tabel 2 , digunakan hasil perbandingan pemodelan 2 ( mesh generator) dengan nilai Manning 32 m 1/2 /s dan Eddy Viscosity 0,22 m 2 /s, karena memiliki nilai RMSE yang rendah sebesar  0,3888  dan  nilai correlation  sebesar  0,9873.  Pada Gambar  4 terlihat  data  hasil pengukuran dan data hasil pemodelan dengan pasang surut, Gambar 5 data hasil pengukuran dan data hasil pemodelan untuk nilai kecepatan arus terhadap arah U (barat-timur) dan arus terhadap arah V (utara-selatan).
Pemodelan Manning [m 1/2 /s] Eddy Viscosity [m 2 /s] RMSE Correlation Pasang Surut [m] Kecepatan Arah U [m/s] Kecepatan Arah V [m/s] Pasang Surut [m] Kecepatan Arah U [m/s] Kecepatan Arah V [m/s] 1 23 0 1,1119 0,1641 0,1758 0,0021 0,0085 0,0247 1 28 0 1,1073 0,1638 0,1758 0,6018 0,0037 0,0223 1 32 0 1,0569 0,1579 0,1749 0,0542 0,1198 0,1087 1 23 0,22 1,0514 0,1641 0,1758 0,0652 0,0078 0,0206 1 28 0,22 1,0549 0,1638 0,1758 0,0584 0,0044 0,0234 1 32 0,22 1,0568 0,1638 0,1749 0,0548 0,0022 0,0266 2 23 0 0,4232 0,1883 0,1995 0,9852 0,4597 0,4691 2 28 0 0,3398 0,1913 0,1986 0,9869 0,5514 0,5451 2 32 0 0,3857 0,1941 0,1985 0,9871 0,6071 0,6025 2 23 0,22 0,4230 0,1882 0,1994 0,9852 0,4697 0,4603 2 28 0,22 0,3996 0,1911 0,1984 0,9869 0,5522 0,5460 2 32 0,22 0,3888 0,1939 0,1982 0,9873 0,6082 0,6036

## Wulan Fitry Hardhiyanti, Yessi Nirwana Kurniadi, Eva Mustikasari, Yogi Noviadi

## Gambar 4. Perbandingan pasang surut hasil pengamatan dan pemodelan

## Gambar 5. Perbandingan kecepatan arus pemodelan dan pengamatan

## 4.3. Analisa Kecepatan Arus dan Pasang Surut Pemodelan
Setelah dilakukan validasi dan memilih pemodelan dengan koefisien Manning 32 m 1/2 /s dan Eddy Viscosity 0,22 m 2 /s dilakukan dalam interval waktu 60 hari dimulai pada tanggal 15 Mei hingga 15 Juli (dua bulan). Hasil pemodelan 2 (menggunakan batimetri hasil mesh generator) menunjukan terjadi peningkatan kecepatan arus di beberapa titik hingga time step ke 720. Pada Gambar 6 terlihat hasil pemodelan pasang surut dengan satuan m dengan pergerakan yang ditunjukkan oleh vektor. Pasang surut pada wilayah tinjauan  memiliki interval 0,8 m sampai 1,5 m. Gambar 7 terlihat kecepatan arus dalam arah U (barat-timur) memiliki besar antara -0,5 m/s sampai 0,8 m/s. Pada Gambar 8 ditunjukkan kecepatan arus dalam arah V (utara-selatan) dengan -0,25 m/s sampai 0,8 m/s.
-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 4/30/2016 5/2/2016 5/4/2016 5/6/2016 5/8/2016 5/10/2016 5/12/2016 5/14/2016 5/16/2016 Ele vasi M u ka Air (met er ) Waktu (Jam) Pemodelan Pengukuran -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4
Kecep at an Aru s Ara h V (m/s) Kecepatan Arus Arah U (m/s) Pemodelan Pengukuran
Gambar 6. Hasil output pasang surut
Gambar 7. Hasil kecepatan arus arah U
Gambar 8. Hasil kecepatan arus arah V

## 4.4. Pasang Surut Purnama dan Pasang Surut Perbani
Kondisi  pasang  surut  purnama  terjadi  pada  tanggal  7  Mei  2016,  sedangkan  pasang  surut perbani  terjadi  pada  tanggal  14  Mei  2016.  Pada Gambar  9 saat  kondisi  pasang  surut purnama, pasang surut bergerak dari arah tenggara dan timur menuju barat laut (selat). Saat kondisi surut purnama terlihat pada Gambar 10 terjadi perbedaan arah pasang surut yaitu menuju tenggara dan timur. Tinggi pasang surut pada kondisi purnama berkisar 1,5 meter.

## Wulan Fitry Hardhiyanti, Yessi Nirwana Kurniadi, Eva Mustikasari, Yogi Noviadi
Kondisi  pasang  surut  perbani  terjadi  pada  tanggal  14  Mei  dimana  terlihat  pasang  perbani pada Gambar 11 dan surut perbani pada Gambar 12 . Pada saat kondisi menuju pasang perbani, arah pasang surut bergerak dari arah tenggara dan timur menuju barat laut, namun saat kondisi surut perbani, arah pasang surut bergerak ke menuju timur dan tenggara. Tinggi pasang surut sekitar 0,5 meter. Jenis pasang surut pada perairan ini adalah (semi diurnal tide) dengan periode pasang surut rata-rata 12 jam.
Gambar 9. Kondisi pasang purnama
Gambar 10. Kondisi surut purnama
Gambar 11. Kondisi pasang perbani
Pola Hidrodinamika di Perairan Nunukan sebagai Usulan pada Permasalahan Abrasi Pulau-Pulau Kecil

## Gambar 12. Kondisi surut perbani

## 4.5. Wilayah Berpotensi Abrasi
Dilihat  dari output  hasil  pemodelan  yang  dilakukan  terhadap  wilayah  tinjauan,  beberapa wilayah di Kabupaten Nunukan mengalami peningkatan kecepatan arus. Kecapatan arus pada pemodelan dengan time step ke-0 memperlihatkan kecepatan arus berkisar antara -0,2500 m/s sampai -0,1125 m/s yang terlihat pada Gambar 13 . Peningkatan terlihat pada Gambar 14 dengan time step ke-360 kecepatan arus meningkat berkisar antara 0,0250 m/s sampai 0,4375 m/s, sedangkan pada Gambar 15 dengan time step ke-720, kecepatan relatif sama terhadap time step sebelumnya yakni antara 0,0250 m/s sampai 0,4375 m/s. Dalam hal ini, perlu adanya peningkatan pengamanan terhadap garis pantai pada wilayah yang mengalami peningkatan kecepatan arus yang dapat dilihat pada Gambar 16 . Abrasi terlihat terjadi di Barat Daya Pulau Sebatik, Perairan Nunukan pada titik X 117°42’36’’ Y 4°9’14,4’’ hingga X 117°42’28,8’’  Y  4°8’38,4’’  sepanjang  1014,26  meter.  Pada  wilayah  tinjauan,  berdasarkan penelitian dari Tim Sebatik Pusat Penelitian dan Pengembangan Geologi Kelautan, wilayah tersebut memiliki jenis tanah lanau kepasiran. Ukuran butir lanau kepasiran dapat diketahui ukuran  butirnya  berdasarkan  jenis  tanahnya  menurut  AASHTO (American  Association  of Highway and Trasportation Officials) dalam buku M. Das (1995). Diketahui ukuran butir tanah lanau kepasiran memiliki besar butir 0,01 mm atau 10 -5 m. Menurut Dronkers (2005) dalam bukunya  berjudul Dynamics  of  Coastal  Systems,  terdapat  hubungan  antara  besar  butir dengan  kecepatan  arus  sehingga  dapat  diketahui  kemungkinan  terjadinya  erosi  ataupun abrasi. Pada wilayah tinjauan diketahui kecepatan arus sebesar antara 0,0250 m/s sampai 0,4375 m/s m/s dengan ukuran butir 10 -5 m, pada grafik terlihat wilayah tersebut berpotensi abrasi.
Gambar 14. Kecepatan arus pada time step 360
Gambar 15. Kecepatan arus pada time step 720
Gambar 16. Wilayah berpotensi abrasi

## 5. KESIMPULAN DAN SARAN

## 5.1. Kesimpulan
Berdasarkan hasil penelitian yang telah dilakukan pada Perairan Nunukan, Kalimantan Utara, dihasilkan  kesimpulan  bahwa  hidrodinamika  Perairan  Nunukan  dibangkitkan  oleh  pasang surut karena wilayah tinjauan berada di selat, pasang Surut pemodelan yang sebelumnya telah dilakukan validasi menggunakan RMSE dan Correlation memiliki ketinggian antara -1,5 meter hingga 1,5 meter, besar kecepatan arus pada titik tinjauan yang bergerak timur laut-
tenggara, dan beberapa wilayah di Perairan Nunukan mengalami peningkatan kecepatan arus diantara 0,0250 m/s sampai 0,4375 m/s. Wilayah yang terdampak abrasi memiliki jenis tanah lanau  pasiran  dengan  ukuran  butiran  sebesar  0,1  mm  sehingga  didapat  bahwa  wilayah tinjauan berpotensi abrasi.

## 5.2. Saran
Data  batimetri  yang  digunakan  sebaiknya  memiliki  skala  yang  lebih  kecil,  sehingga  dapat diperoleh grid  yang  lebih  rapat  agar  hasil  pemodelan  grid  menjadi  lebih  teliti, mempertimbangkan pengaruh sedimentasi pada penelitian selanjutnya, mempertimbangkan pengaruh  sungai-sungai  yang  berada  di  sekitar  perairan  dan  menganalis  perubahan  garis pantai  dengan  penggunaan  bangunan  pelindung  pantai  untuk  mengatasi  abrasi  yang mungkin terjadi pada wilayah yang berpotensi terdampak abrasi.

## UCAPAN TERIMA KASIH
Pusat Penelitian dan Pengembangan Sumberdaya Laut dan Pesisir (P3SDLP), Badan Penelitian dan Pengembangan Kelautan Perikanan, Kementerian Kelautan dan Perikanan,
Pusat Penelitian dan Pengembangan Geologi Kelautan, Badan Penelitian dan Pengembangan Energi Sumber Daya Mineral, Kementrian Energi Sumber Daya dan Mineral.

## DAFTAR RUJUKAN
Badan Penelitian dan Pengembangan Kelautan dan Perikanan. (2016). Laporan Akhir Nunukan. Jakarta: Kementerian Kelautan dan Perikanan.
Dronkers, Job. (2005). Dynamics of Coastal System. Singapore: World Scientific Publishing Co. Pte. Ltd.
Google, US Dept of State Geographer, AfriGis (Pty). Ltd., Digital Globe. (2017, August 25). Google         Earth         Pro.          Dipetik         dari          Google          Earth: https://www.google.com/maps/@4.1530244,117.67254,12793m/data=!3m1!1e3 Das, B.M. (1995). Mekanika Tanah Jilid I. Jakarta: Erlangga. Pusat Penelitian dan Pengembangan Geologi Kelautan. (2005). Laporan Akhir Perairan Sebatik Kabupaten  Nunukan  Kalimantan  Timur.  Bandung:  Badan  Geologi  Kementerian  Energi dan Sumber Daya Mineral.
Pratikto,  W.A;  Armono,  H.D;  Suntoyo.  (1996). Perencanaan  Fasilitas  Pantai  dan  Laut. Yogyakarta: BPFE.