Studi Perubahan Garis Pantai Tahun 2012-2022 di Pesisir Pantai Bagik Kembar, Kota Mataram, Nusa Tenggara Barat

## Safira Ashilah*, Agus Anugroho Dwi Suryoputro dan Dwi Haryo Ismunarti
Departemen Oseanografi, Fakultas Perikanan dan Ilmu Kelautan, Universitas Diponegoro Jl. Prof. H. Sudarto, SH, Tembalang, Semarang, Jawa Tengah 50275, Indonesia
Email: * safiraashilah@students.undip.ac.id

## Abstrak
Pantai Bagik Kembar merupakan pantai yang berada di Kecamatan Sekarbela, Kota Mataram yang mengalami perubahan garis  pantai  pada  setiap  tahunnya.  Perubahan  ini  menyebabkan  wilayah  pesisir  Pantai  Bagik  Kembar  seringkali mengalami banjir yang dapat merugikan aktivitas dan pemukiman warga sekitar. Tujuan dari penelitian ini adalah untuk menganalisis  luasan  abrasi  dan  akresi  yang  terjadi  pada  tahun  2012  hingga  2022  di  pesisir  Pantai Bagik  Kembar  dan sekitarnya. Metode yang digunakan yaitu analisis DSAS (Digital Shoreline Analysis System) dengan menghitung nilai NSM dan EPR  menggunakan data citra Landsat 7 pada tahun 2012 dan Landsat 8 pada tahun 2017 dan 2022. Analisis ini dikaitkan dengan data kelerengan pantai, sedimen, pasang surut BIG, data angin dan peramalan gelombang serta arus sejajar  pantai  dari  perhitungan  gelombang  pecah.  Hasil  penelitian  menunjukkan  pesisir  Pantai  Bagik  Kembar  dan sekitarnya mengalami perubahan garis pantai yang didominasi oleh terjadinya abrasi dengan luas pada tahun 2012 – 2017 sebesar 2,512 ha dan pada tahun 2017 – 2022 luas abrasi 1,892 ha. Sedangkan akresi yang terjadi hanya 0,296 ha pada tahun  2012  –  2017  dan  0,114  ha  pada  tahun  2017  –  2022.  Faktor  penyebab  perubahan  garis  pantai  didominasi  oleh gelombang yang dibangkitkan oleh angin sehingga menyebabkan terjadinya transpor sedimen oleh arus sejajar pantai.
Kata kunci: Abrasi, DSAS, Gelombang
Abstract
Study of 1 Shoreline Changes in 2012-2022 on the Coast of Bagik Kembar Beach, Mataram City,

## West Nusa Tenggara
Bagik Kembar Beach is  a beach located in  Sekarbela  District, Mataram City, which often experiences changes in  the coastline every year. These changes cause the coastal area of Bagik Kembar Beach to often experience disasters that can harm activities and local residents. The purpose of this study is to analyze the extent of shoreline changes that occurred in 2012 to 2022 on the coast of Bagik Kembar Beach and its surroundings. The method used is DSAS (Digital Shoreline Analysis System) analysis by calculating NSM and EPR values using Landsat 7 image data in  2012 and Landsat  8 in 2017 and 2022. This analysis is associated with coastal slope data, sediment, BIG tides, wind data and wave forecasting and parallel coastal currents from breaking wave calculations. The results showed that the coast of Bagik Kembar Beach and its surroundings experienced changes in the coastline dominated by abrasion with an area in 2012-2017 of 2.512 ha and in 2017-2022 the abrasion area was 1.892 ha. While accretion that occurred was only 0.296 ha in 2012 - 2017 and 0.114 ha in 2017 - 2022. Factors causing shoreline change are dominated by wind-generated waves that cause sediment transport by longshore current.

## Keywords:  Abrasion, DSAS, wave

## PENDAHULUAN
Kota Mataram khususnya di sepanjang pesisir Pantai Bagik Kembar merupakan salah satu pantai yang berada di Kecamatan Sekarbela yang mengalami perubahan garis pantai. Perubahan garis pantai pada lokasi ini didominasi dengan terjadinya abrasi. Menurut Widada et al. (2022), proses abrasi dapat merusak bangunan pantai  dan  mundurnya  garis  pantai  yang menjadi  pintu  masuknya  air  laut  ke  darat  dan  terjadi  banjir  rob di wilayah  pesisir.  Abrasi  pantai  pada  lokasi  penelitian  ini  dominan  terjadi  pada  musim  barat  dimana  sering terjadi angin kencang, yang menimbulkan gelombang besar. Selain terjadi banjir akibat  luapan air laut saat gelombang  besar,  rumah-rumah  nelayan  di  Pantai  Bagik  Kembar  juga  tertimbun  oleh  sedimen  pantai  yang terbawa oleh gelombang (BAPPEDA Kota Mataram, 2019). Oleh karena itu, pada lokasi ini perlu dilakukan
penanganan  dengan  melakukan  mitigasi  bencana  pada  pesisir  pantai  Bagik  Kembar,  Kota  Mataram  dan sekitarnya. Untuk itu diperlukan studi perubahan garis pantai yang bertujuan untuk menganalisis luasan abrasi dan akresi yang terjadi pada tahun 2012-2022. Berdasarkan kajian terdahulu oleh Kusumawardani et al. (2019), pada daerah Pantai Bagik Kembar dan sekitarnya dilakukan pemetaan dan analisis perubahan garis pantai dari tahun 2003 – 2018 dengan citra Landsat 7 dan Landsat 8 dengan memetakan garis pantai dengan transformasi indeks yaitu  NDWI (Normalize  Difference  Water  Index) dan  dilengkapi oleh data  kecepatan  arus  laut pada wilayah  Ampenan.  Tujuan  penelitian  ini  menggunakan  pendekatan  atau  metode  berbeda  yang  belum digunakan  pada  penelitian  sebelumnya  di  daerah  ini  dengan  melakukan  pemetaan  garis  pantai  dengan menggunakan  metode  analisis  DSAS (Digital  Shoreline  Analysis  System) serta  diperbarui  pemetaan  garis pantai  pada  tahun  2022  di  pesisir  Pantai  Bagik  Kembar  dan  sekitarnya.  Dengan  metode  DSAS  ini  dapat menghitung laju dan jarak perubahan garis pantai berdasarkan acuan titik perpotongan dari garis transek yang telah  dibuat  dalam  waktu  tertentu  dan  berbasis  geospasial.  dengan  harapan  dapat  memberikan  informasi mengenai perubahan garis pantai yang signifikan.

## MATERI DAN METODE
Materi 1 Penelitian Dalam 1 penelitian  ini  menggunakan  data 1 primer  dimana  data 1 ini  berupa  data  hasil  pengukuran  dan pengambilan sampel di lapangan serta data sekunder 1 yaitu data yang diperoleh dari 1 instansi maupun website resmi  yang  digunakan  sebagai  data  pendukung.  Data  primer  yang  digunakan  dalam  penelitian 1 yaitu  data kelerengan 1 pantai, sedimen dasar, dan 1 data batimetri. Sedangkan data sekunder yang digunakan adalah 1 data pasang  surut  BIG  bulan  Oktober  tahun  2022,  data  angin  ECMWF  per  musim  tahun  2021-2022,  data  citra satelit 1 Landsat 7 tahun 2012 dan Landsat 8 1 tahun 2017 dan 2022 dari https://earthexplorer.usgs.gov serta peta Rupa Bumi Indonesia dari BIG.

## Metode Penelitian
Metode  penelitian  yang 1 digunakan  adalah  metode 1 kuantitatif  karena  data  yang  digunakan  pada penelitian ini berupa angka serta analisis menggunakan 1 statistik (Sugiyono, 2009). Metode kuantitatif 1 yang dimaksud dalam penelitian ini adalah perhitungan kelerengan pantai, ukuran butir sedimen, arah dan kecepatan angin, tinggi dan 1 periode gelombang serta nilai luasan perubahan garis pantai pada lokasi penelitian. Dalam mengambil  data  lapangan  metode  yang  digunakan  adalah  metode purposive 1 sampling yaitu  teknik pengambilan 1 sampel dengan sengaja dan berdasarkan 1 suatu pertimbangan dan tujuan tertentu (Safitri et al., 2019).  Penentuan  titik  sampel  dalam 1 penelitian  ini  dilakukan  dengan  menyesuaikan 1 titik  terjadinya perubahan garis 1 pantai berdasarkan interpretasi citra.

## Metode Pengumpulan dan Analisis Data

## 1. Metode Pengolahan Garis Pantai dengan Citra Satelit 1 Landsat
Garis pantai diolah dengan mengunduh data citra satelit1Landsat 7 untuk tahun 2012 dan citra Landsat 1 8 untuk  tahun  2017  dan  2022  diolah  dengan 1 menggunakan software Envi  5.3  untuk  melakukan  koreksi radiometrik. Setelah itu dilakukan tahap layer stacking untuk menggabungkan band citra. Kemudian dilakukan pemotongan  citra  sesuai  dengan  lokasi  penelitian.  Setelah  pemotongan  citra  dilakukan  penentuan band treshould atau nilai ambang batas. Band yang digunakan untuk Landsat 7 adalah band 5 atau band hijau/SWIR1 sedangkan untuk Landsat 8 digunakan band 6 atau band hijau/SWIR1 kemudian dilakukan perkalian citra biner dengan band treshould serta di ekstrasi ke dalam bentuk vektor. Tahap ini dilakukan 1 agar dapat membedakan daratan dan lautan untuk mengetahui garis pantai dengan mudah. Pengolahan citra selanjutnya dilakukan pada software ArcMap 10.8 untuk melakukan analisis garis pantai menggunakan 1 metode DSAS (Digital Shoreline Analysis System) . Tahapan awal pada software ini yaitu dengan memasukkan citra yang telah diekstrasi pada software Envi  5.3  lalu  dilakukan  pemotongan  garis  pantai  area  penelitian.  Kemudian  citra  pada  tiap  tahun digabungkan  dengan  menggunakan tools  Merge .  Setelah  itu  dilanjutkan  dengan  metode  DSAS  yang menghasilkan perhitungan 1 EPR (End Point Rate) dan NSM 1 (Net Shoreline Movement) agar dapat dilakukan analisis nilai luasan abrasi dan1akresi yang terjadi. Selanjutnya dilakukan layout peta sebagai hasil informasi perubahan garis pantai.

## 2. Metode Pengolahan Pasang Surut
Pengolahan 1 data pasang surut ditujukan untuk mengetahui nilai MSL, HHWL dan LLWL yang didapat dari  website  BIG  kemudian 1 diolah 1 menggunakan  metode Admiralty.  Setelah  itu  hasil  data download dikurangi nilai rata-rata untuk mendapatkan nilai amplitudonya (A)
𝐴 = 𝑁𝑖𝑙𝑎𝑖 𝑂𝑏𝑠𝑒𝑟𝑣𝑎𝑠𝑖 − 𝑀𝑆𝐿
Selain MSL, analisis digunakan untuk mendapatkan konstanta 1 harmonik pasang surut sehingga dapat diperoleh tipe 1 pasang surut menggunakan bilangan formzahl (F) menggunakan persamaan berikut 𝐹 = 𝐴𝑂 1 + 𝐴𝐾 1 𝐴𝑀 2 + 𝐴𝑆 2 Serta dapat diperoleh juga nilai LLWL dan HHWL (Ryanto et al., 2022) 𝐿𝐿𝑊𝐿 = 𝑀𝑆𝐿 − ( 𝐴 ( 𝑀 2 − 𝐾 1 + 𝑂 1 + 𝑃 1 + 𝐾 2)) 𝐻𝐻𝑊𝐿 = 𝑀𝑆𝐿 + ( 𝐴 ( 𝑀 2 + 𝐾 1 + 𝑂 1 + 𝑃 1 + 𝐾 2)) 3. Metode Pengolahan Data Angin dan Peramalan Gelombang
Data 1 angin  dapat  diunduh  dari  website  ECMWF 1 ( European  Centre  for  Medium  Range  Weather Forecasts) https://cds.climate.copernicus.eu/. Data ini mencakup kecepatan 1 dan arah angin komponen u dan v per 1 jam dengan satuan m/s 1 selama 4 musim dalam satu tahun (Desember 2021 – Oktober 2022). Setelah mengunduh data  tersebut,  file  NetCDF  hasil  unduhan  dikonversi  menjadi  file  TXT  menggunakan software ODV  4.0.  Selanjutnya,  data 1 diolah  menggunakan  Microsoft  Excel  untuk  mengelompokkan  data  angin  per musim.  Setelah  itu  data  kecepatan  dan  arah  angin  per  musim  dimasukkan  ke Easywave  Algorithm dimana Easywave ini dapat melakukan peramalan gelombang dengan   metode SMB (Svendrup Munk Bretschneider) berbasis bahasa pemrograman Python (Siregar, et al., 2020). Peramalan gelombang menghasilkan nilai tinggi signifikan 8 (Hs)  dan  periode 8 signifikan  (Ts) 8 gelombang.  Kemudian  dilakukan  analisis  mengenai  kondisi hidro-oseanografi  pada  lokasi  penelitian  yang  dikaitkan  dengan  kecepatan  arah  angin  berhembus  serta pengaruh dari gelombang terhadap perubahan garis pantai yang terjadi.

## 4. Pengolahan Gelombang Pecah
Untuk mengetahui nilai sudut gelombang pecah (α b ) digunakan nilai tinggi 1 signifikan (Hs) dan periode signifikan 1 (Ts) gelombang yang telah diperoleh dari peramalan gelombang per musim. Menurut Triatmodjo (2012) dalam Kristanto et al., (2022), persamaan yang digunakan adalah sebagai berikut :
1. Menghitung 8 L 0
𝐿 0 =  1,56 𝑇 2
2. Menghitung C 0
𝐶 0 = 𝐿 0 𝑇 ⁄
3. Menghitung 8 d/L 0
𝑑 𝐿 0 ⁄
4. Mencari nilai 8 d/L melalui 8 interpolasi nilai 8 d/L 0 menggunakan 8 tabel C-1 8 pada SPM 1984 5. Mencari 8 nilai panjang 8 gelombang (L)
𝐿 = 𝑑 ∶   𝑑 𝐿 ⁄
6. Mencari 8 nilai cepat 8 rambat gelombang (C)
𝐶 = 𝐿 𝑇 ⁄
7. Menghitung 8 koefisien 8 pendangkalan  (Ks)  melalui 8 interpolasi  nilai 8 d/L 0 menggunakan  tabel  C-1 pada 8 SPM 1984
8. Menghitung 8 H’ 0
𝐻′ 0 = 𝐻 𝑠 𝐾 𝑠 ⁄ 9. Menghitung H’ 0 /gT 2
10. Menghitung Hb
𝐻 𝑏 𝐻′ 0 = 1 3,3 (𝐻′ 0 𝐿 0 ) ⁄ 1 3 ⁄
11. Menghitung db
𝑑 𝑏 𝐻 𝑏 = 1 𝑏 − (𝑎𝐻 𝑏 𝑔𝑇 2 ) ⁄
dimana a dan 8 b merupakan fungsi kemiringan 8 pantai m dengan 1 persamaan:
𝑎 =  43,75 (1 − 𝑒 −19𝑚 ) 𝑏 = 1,56 (1 + 𝑒 −19,5𝑚 ) Triatmodjo (1999)
12. Menghitung panjang 8 gelombang pecah 8 (Lb) 𝐿 𝑏 = 𝑇 𝑠 √𝑔𝑑 𝑏 13. Menghitung sudut 8 gelombang 8 pecah (α b ) sin 𝛼 𝑏 𝐿 𝑏 = sin 𝛼 0 𝐿 0
5. Metode Pengolahan Arus 1 Sejajar Pantai (Longshore 1 Current) Menurut Triatmodjo (1999) dalam Kristanto et al. (2022), untuk menghitung 1 nilai kecepatan longshore current (arus 1 sejajar pantai) berdasarkan nilai tinggi 1 gelombang pecah (Hb) dan sudut 1 gelombang pecah (αb) dapat menggunakan 1 rumus empiris seperti berikut:
𝑉 = 1,17(𝑔 𝑥 𝐻 𝑏 ) 1 2 ⁄ sin 𝛼 𝑏 cos 𝛼 𝑏 Keterangan 1 : V = Kecepatan 1 longshore 8 current g
= Percepatan 8 gravitasi Hb = Tinggi 1 gelombang 8 pecah αb = Sudut 1 gelombang 8 pecah 6. Analisis Ukuran Butir Sedimen
Untuk  melakukan  analisis  ukuran  butir 1 sedimen,  sampel  sedimen 1 sebanyak  200 gram  dapat  diambil pada tiap stasiun yan dapat dilihat pada gambar 1. Pengambilan titik sedimen tersebut didasarkan oleh tujuan penelitian  atau  dengan  metode purpossive  sampling. Setelah  pengambilan  sampel  dikeringkan  kemudian dilakukan  pengayakan.  Selama  proses  pengayakan,  sedimen  dipisah  sesuai 1 ukuran  butirnya  menggunakan sieve  shaker 1 dengan  amplitudo  60 1 dalam  waktu  15  menit.  Nilai  ukuran 1 butir  sedimen  dapat  dihitung pada 1 setiap  lapisan  mesh  yang  berukuran  (2 8 μm,  0.5 8 μm,  0.3 8 μm,  0.150 8 μm  dan  0.063 8 μm).  Pemipetan diperlukan apabila terdapat sedimen yang sangat 1 halus, yaitu 1 <0,063, yang tidak 1 dapat diukur dengan sieve shaker 1 (Holme dan Mclntyre, 1984).
7. Perhitungan Transpor Sedimen
Menurut  Putra et  al. (2022),  transpor  sedimen  yang  diakibatkan  oleh longshore  current diprediksi dengan menggunakan 1 persamaan empiris yang berdasar kepada kondisi 1 gelombang pada daerah penelitian. Persamaan tersebut merupakan korelasi antara 1 transpor sedimen dengan komponen 1 fluks energi gelombang sepanjang pantai. terdapat persamaan empiris yang 1 digunakan untuk memprediksi transpor 1 sedimen akibat longshore current , seperti persamaan yang dikembangkan oleh Komar (1971) sebagai berikut:
1. Menghitung C b
𝐶 𝑏 =   √𝑔𝑑 𝑏 2. Menghitung fluks 8 energi 8 gelombang sepanjang 8 pantai saat pecah
𝑃 1 = 𝜌𝑔 8 𝐻 𝑏 𝐶 𝑏 sin 𝛼 𝑏 cos 𝛼 𝑏
Menurut CERC (1984) dalam Putra et al. (2022), untuk menghitung angkutan sedimen pertahun dengan satuan m3 /tahun adalah sebagai berikut:
𝑄 𝑠 =  1290𝑃 1
untuk menghitung angkutan sedimen perhari adalah sebagai berikut:
𝑄 𝑠 =  3,534𝑃 1
Keterangan  : 𝑃 1  :  Komponen 8 fluks 8 energi 8 gelombang  sepanjang 8 pantai  saat  pecah  (Nm/d/m) 𝜌 :  Massa jenis 1 air  laut  (1025  kg/m 3 ) 𝑔 :  Percepatan 8 gravitasi  (9,81  m/s 2 ) 𝐻𝑏 :  Tinggi 8 gelombang pecah  (m) 𝐶𝑏 : Kecepatan gelombah pecah (m/s) 𝛼𝑏 : Sudut 8 datang 8 gelombang pecah (°) 𝑄𝑠 : Volume angkutan sedimen dalam satuan waktu (m 3 /hari dan m 3 /tahun)

## Lokasi Penelitian
Penelitian dilakukan di sepanjang pesisir Pantai Bagik Kembar dan sekitarnya di Kecamatan Sekarbela, Kota  Mataram,  NTB.  Lokasi  penelitian  ini  secara  astronomis  terletak  di  antara  8°36’8,075”LS  - 116°4’23,034”BT  sampai  8  °35’15,043”LS  -  116°4’20,071”BT.  Lokasi  penelitian 1 dapat  dilihat  pada gambar 1 berikut.
Gambar 1. Peta 1 Lokasi Penelitian di sekitar Pantai Bagik Kembar, Kota Mataram, Nusa Tenggara Barat

## HASIL DAN PEMBAHASAN
Kelerengan 1 Pantai Pengukuran 1 kelerengan pantai dilakukan dengan 1 menggunakan alat konvensional berupa busur yang diikat tali yang terdapat pemberat berupa benda kecil untuk dapat melihat sudut  kemiringan dari pantainya. Pengukuran ini dilakukan 1 pada 4 titik pengamatan 1 dengan 3 kali 1 pengukuran yang mencakup daerah dengan tipe  morfologi  pantai  yang  terlihat  berbeda  dari  muara  kecil  hingga  muara  kali  ancar.  Hasil  pengukuran tersebut disajikan pada Tabel 1 yang diklasifikasikan berdasarkan klasifikasi Van Zuidam dan disertai dengan nilai rata-rata sebagai berikut.
Tabel 1 1. Nilai 1 Kelerengan Pantai Stasiun Koordinat Nilai 1 Kelerengan Pantai ( o ) Keterangan 1 Pengukuran 1 Rata-rata (%) Bujur(°) Lintang(°) Ke-1 Ke-2 Ke-3 L1 116.07326 -8.59321 15 16 17 21 – 55 Curam L2 116.07316 -8.5927 10 15 15 21 – 55 Curam L3 116.07294 -8.59143 9 12 12 21 – 55 Curam L4 116.07328 -8.59853 10 21 21 21 – 55 Curam
Berdasarkan 1 hasil  penelitian  yang  telah 1 dilakukan  di  Pantai  Bagik  Kembar  dan  sekitarnya  dapat diketahui  bahwa  pantai  ini  dikategorikan  dengan  kemiringan  pantai  yang  curam  dengan  didominasi  oleh sedimen berjenis pasir yang dimana sampel sedimen diambil berada di sekitar muara sungai Ancar. Hal ini
sesuai dengan penelitian yang dikaji oleh BAPPEDA Kota Mataram (2019) bahwa secara 1 keseluruhan pantai Kota  Mataram  adalah 1 jenis  pantai  berpasir  karena 1 didominasi  oleh  sedimen 1 pasir  hitam.  Batimetri  atau kedalaman laut 1 di depan pantai Kota Mataram juga cenderung landai. Kemiringan pantai di atas garis 1 surut hingga tepi pantai secara umum 1 lebih curam. Hal ini sesuai dengan karakteristik sedimen berpasir, di mana kemiringan 1 dasar pantai tergantung pada bentuk dan ukuran material dasar.  Triatmodjo (1999) menyatakan bahwa  kemiringan  pantai  berpasir  umumnya  berkisar 1 antara  1:20  hingga  1:50.  Bentuk  profil 1 pantai dipengaruhi  oleh 1 serangan  gelombang,  sifat-sifat 1 sedimen,  kondisi  gelombang,  arus,  dan  batimetri  pantai. Penelitian  sebelumnya  juga  menunjukkan  bahwa  terdapat  proses  geomorfologi  yang  terjadi  di  perairan Kabupaten Lombok 1 Barat dan Kota Mataram, termasuk pengendapan 1 oleh sungai, gelombang, 1 dan 1 arus 1 (Kusumawadani et al., 2018). Pasang 8 Surut Data 8 pasang 8 surut  yang  diperoleh  dari  BIG  dianalisis  menggunakan 1 metode  admiralty . Hasil  yang diperoleh berupa komponen harmonik 1 pasang surut, yaitu S 0 , M 2 , S 2 , N 2 , K 1 , O 1 , M 4 , MS 4 , K 2 , dan P 1 . Tabel 2. Komponen 1 Harmonik Pasang Surut Metode Admiralty Constants S 0 8 M 2 8 S 2 8 N 2 8 K 1 8 O 1 8 M 4 8 MS 4 8 K 2 8 P 1 8 A (Cm) 127 23 17 6 31 24 1 2 6 12 g (°) 0 355 172 200 103 98 148 5 172 103
Nilai dari komponen harmonik yang didapatkan dapat digunakan 1 untuk mengetahui nilai elevasi dan tipe  karakteristik  pasang  surut.  Diketahui  nilai  elevasi  pasang  surut  yaitu  MSL (Mean  Sea  Level) ,  HHWL (Highest High-Water Level) , LLWL (Lowest Low Water Level) , MHWL (Mean High Water Level) , MLWL (Mean Low Water Level) sebagai berikut.
Tabel 3. Nilai 8 Elevasi Muka Air Laut Pasang 8 Surut Metode Admiralty Vertical Datum: HHWL =  241 Cm MHWL =  205 Cm MSL =  127 Cm MLWL =  48 Cm LLWL =  12 Cm
Hasil 1 pengolahan  data  pasang  surut  menggunakan  metode 1 Admiralty  ini  diperoleh  nilai 1 Formzahl sebesar 1,37 dimana 0,25 < F≤ 1,5 yang artinya tipe 1 pasang surut pada perairan Kota Mataram ini termasuk pada tipe pasang surut campuran harian ganda. Dari hasil pengolahan data pasang 1 surut tersebut didapatkan bentuk grafik yang sesuai dengan tipe pasang surutnya seperti yang ditunjukkan pada gambar berikut.
Gambar 2. Grafik Pasang Surut Perairan Kota Mataram Bulan Oktober 2022
0 50 100 150 200 250 2 9 -S e p -2 0 2 2 3 0 -S e p -2 0 2 2 1 -O c t- 2 0 2 2 2 -O c t- 2 0 2 2 3 -O c t- 2 0 2 2 4 -O c t- 2 0 2 2 5 -O c t- 2 0 2 2 6 -O c t- 2 0 2 2 7 -O c t- 2 0 2 2 8 -O c t- 2 0 2 2 9 -O c t- 2 0 2 2 1 0 -O c t- 2 0 2 2 1 1 -O c t- 2 0 2 2 1 2 -O c t- 2 0 2 2 1 3 -O c t- 2 0 2 2 1 4 -O c t- 2 0 2 2 1 5 -O c t- 2 0 2 2 1 6 -O c t- 2 0 2 2 1 7 -O c t- 2 0 2 2 1 8 -O c t- 2 0 2 2 1 9 -O c t- 2 0 2 2 2 0 -O c t- 2 0 2 2 2 1 -O c t- 2 0 2 2 2 2 -O c t- 2 0 2 2 2 3 -O c t- 2 0 2 2 2 4 -O c t- 2 0 2 2 2 5 -O c t- 2 0 2 2 2 6 -O c t- 2 0 2 2 2 7 -O c t- 2 0 2 2 2 8 -O c t- 2 0 2 2 2 9 -O c t- 2 0 2 2 3 0 -O c t- 2 0 2 2 3 1 -O c t- 2 0 2 2 1 -N o v -2 0 2 2 T id e (cm ) DD-MM-YYYY

## Angin
Data angin diperoleh dari ECMWF diolah menjadi grafik windrose berupa data per musim dari musim barat tahun 2021 – 2022 hingga musim timur pada tahun 2022. Pada musim barat didominasi angin menuju arah barat daya dengan didominasi kecepatan 11,08-17,11 knot. Musim peralihan I didominasi angin menuju arah tenggara dengan kecepatan 7,0-11,08 knot. Musim timur didominasi angin menuju arah tenggara dengan kecepatan 5,7-8,8 knot. Pada musim peralihan II dominasi angin menuju arah tenggara dengan kecepatan 3,6- 5,7 knot.
(a). Windrose Musim 8 Barat  arah 8 utama angin  240 o (Desember 8 2021  – 8 Februari 2022) (b). Windrose Musim  Peralihan  I  arah utama angin 135 o (Maret – Mei 2022)
(c). Windrose Musim 8 Timur 8 arah 8 utama angin 135 o (Juni 8 – Agustus 2022)
(d). Windrose Musim  Peralihan  II  arah utama  angin  135 o (September  –  Oktober 2022)
Gambar 3. Windrose Musim Barat – Musim Peralihan II Tahun 2021-2022

## Peramalan Gelombang
Data  tinggi  gelombang 1 signifikan  (Hs)  dan  periode  gelombang 1 signifikan  (Ts)  diperoleh  dari  data angin  per  musim  selama  1  tahun  (Desember  2021  –  Oktober  2022)  yang  diunduh  melalui  ECMWF.  Data peramalan  gelombang  ini  diolah  dengan  metode  SMB (Sverdurp  Munk  Bretsceneider) menggunakan Easywave Algorithm.
Tabel 4. Hasil Tinggi dan Periode Gelombang Signifikan dengan metode SMB Variabel Musim Barat Musim Peralihan I Musim Timur Musim Peralihan II Hs (m) 0.65 0.14 0.08 0.1 Ts (s) 4.15 1.79 0.91 1.3 Arah datang Gelombang ( o ) 240 (Barat Daya) 135 (Tenggara) 135 (Tenggara) 135 (Tenggara)
Dari analisis data angin dan peramalan gelombang yang diperoleh dari Easywave Algorithm dihasilkan windrose atau mawar angin yang menghasilkan arah angin bervariasi pada setiap musim. Tinggi 8 gelombang signifikan 8 (Hs) dan periode gelombang 8 signifikan (Ts) juga memiliki nilai yang bervariasi pada tiap musim. Menurut Triatmodjo (1999), tinggi dan periode gelombang yang dibangkitkan dipengaruhi oleh angin. Angin yang 1 berhembus  di  atas  permukaan  air  akan  memindahkan  energinya  ke  air.  Kecepatan 1 angin  akan menimbulkan 1 tegangan pada permukaan 1 laut, sehingga air yang semula tenang 1 akan terganggu dan timbul riak 1 gelombang  kecil  di  atas  permukaan 1 air.  Apabila  kecepatan  angin  bertambah,  riak  tersebut  menjadi semakin besar, dan apabila angin berhembus terus 1 akhirnya akan terbentuk gelombang. Tinggi 8 gelombang 8 signifikan  maksimum  terdapat  pada  musim 8 barat  dengan  nilai 8 sebesar  0,65  m  dengan periode gelombang 8 signifikan sebesar 4,15 s. Sedangkan tinggi 8 gelombang 8 signifikan minimum terdapat pada musim 8 timur dengan nilai sebesar 0.08 m dan periode 8 gelombang signifikannya sebesar 0.91 s.

## Gelombang Pecah
Nilai  gelombang  pecah  dapat  diketahui  dari  hasil  perhitungan  gelombang  signifikan  (Hs)  yang  telah diperoleh pada Tabel 8 . Berikut ini merupakan tabel nilai gelombang pecah pada tiap musim.
Tabel 5. Nilai 8 Gelombang 8 Pecah No.  Musim α 0 ( o ) 8 H b (m) 8 d b (m) 8 α b ( o ) 8 1. Musim 8 Barat 240 0,65 0,7 22,48 2. Musim 8 Peralihan I 135 0,15 0,16 18,88 3. Musim Timur 135 0,06 0,07 24,52 4. Musim Peralihan II 135 0,09 0,10 14,31
Berdasarkan  hasil  perhitungan  pada  Tabel 1 5  diketahui  pada  musim  barat  memiliki  nilai  arah  datang gelombang 240° dimana gelombang menuju arah barat 1 daya dan pecah pada kedalaman 0,7 m dengan tinggi gelombang sebesar 0,65 m kemudian menghasilkan sudut gelombang pecah pada arah 22,48°. Pada musim peralihan I gelombang 1 datang menuju arah tenggara dengan 1 sudut 135° dan pecah pada kedalaman 0,16 m dengan tinggi gelombang 1 pecah sebesar 0,15 m kemudian menghasilkan sudut 1 gelombang pecah pada arah 18,88°.  Untuk  musim  timur  gelombang  datang  menuju  arah  tenggara  dengan  sudut  135°  dan  pecah  pada kedalaman 0,107 m dengan tinggi gelombang pecah sebesar 0,06 m kemudian menghasilkan sudut gelombang pecah  pada  arah  24,52°.  Selanjutnya  pada  musim  peralihan  II 1 gelombang 1 datang  ke  arah  tenggara  dengan sudut  135°  dan  pecah  pada  kedalaman  0,1  m  dengan  tinggi  gelombang 1 pecah  sebesar  0,09  m  kemudian menghasilkan sudut gelombang 1 pecah pada arah 1 14,31°.
Arus 8 Sejajar 8 Pantai 8 (Longshore 8 Current) Nilai 8 kecepatan 8 arus 8 sejajar 8 pantai (longshore 8 current) diperoleh  dari  hasil  perhitungan 8 gelombang pecah. Berdasarkan hasil perhitungan tersebut diketahui nilai kecepatan arus sejajar pantai berkisar 8 antara 0,26 – 1,04 m/s. Nilai kecepatan arus 1 terbesar terjadi pada musim barat dengan arah arus menuju utara menyusuri pantai.
Tabel 6. Nilai 8 Kecepatan 8 Arus 8 Sejajar 8 Pantai No.  Musim V(m/s) 1. Musim 1 Barat 1,04 2. Musim 1 Peralihan I 0,43 3. Musim Timur 0,34 4. Musim Peralihan II 0,36

## Gambar 4. Peta Arah 1 Arus Sejajar Pantai pada Musim Barat

## Analisis Ukuran Butir Sedimen
Data  ukuran  butir  sedimen 1 diperoleh  dari  hasil  pengolahan  sampel  sedimen  dasar  yang  diambil  di lapangan dengan menggunakan alat berupa shieve shaker dan neraca digital. Hasil yang diperoleh disajikan pada 8 tabel 8 berikut.
Tabel 8 7. Hasil Pengolahan Ukuran 8 Butir 8 Sedimen Stasiun Koordinat Kedalaman air (m) Jenis Sedimen Bujur Lintang 1 116,074683° -8,587944° 0,9 Pasir 2 116,073389° -8,587631° 0,8 Pasir 3 116,071828° -8,587208° 0,6 Pasir 4 116,070839° -8,585422° 3,9 Pasir 5 116,070581° -8,586656° 4,5 Pasir 6 116,070611° -8,587372° 5,2 Pasir 7 116,070781° -8,588269° 5,1 Pasir 8 116,069297° -8,589522° 6,8 Pasir 9 116,068733° -8,587431° 6,5 Pasir 10 116,068547° -8,585650° 6,6 Pasir 11 116,069225° -8,599478° 6,5 Pasir 12 116,066892° -8,584822° 7,8 Pasir Lanauan 13 116,066886° -8,587514° 8,4 Pasir Lanauan 14 116,068172° -8,590406° 7,9 Pasir Lanauan
Berdasarkan hasil yang diperoleh dari hasil laboratorium dan pengolahan data diperoleh jenis sedimen yang tersebar di Pantai Bagik Kembar dan sekitarnya yaitu sedimen berupa pasir dan pasir lanauan. Perairan ini didominasi oleh sedimen berjenis pasir. Titik pengambilan sampel dapat dilihat pada gambar 1.

## Nilai Transpor Sedimen Sejajar Pantai
Dari  hasil  pengolahan  jenis  sedimen  dapat  dihitung  nilai  transport  sedimen  perhari  dan  pertahunnya. Transpor sedimen tersebar dari garis 1 pantai sampai pada daerah 1 gelombang pecah. Sehingga didapatkan nilai transpor 1 sedimen permusim yang disajikan pada tabel berikut.
Tabel 8. Nilai Transpor 1 Sedimen Sejajar Pantai No.  Musim 8 Q s 8 (m 3 /tahun) Q s 8 (m 3 /hari)
1. Musim 8 Barat 973.267,59 2.666,32 2. Musim Peralihan I 94.097,50 257,783 3. Musim Timur 31.288,95 85,71 4. Musim Peralihan II 35.898,84 98,35
Berdasarkan  hasil  yang  diperoleh  dapat  diketahui  bahwa  potensi  transpor  sedimen  pertahun  berkisar antara 31.288,95 – 973.267,59 m 3 /tahun dimana dalam sehari potensi transpor sedimen berkisar antara 85,71 – 2.666,32 m 3 /hari. Potensi transpor sedimen terbesar terjadi pada musim barat karena dipengaruhi oleh nilai kecepatan arus sejajar pantai yang paling tinggi dari musim lainnya.
Triatmodjo (1999) menyatakan bahwa arus sejajar pantai dapat mengangkut sedimen yang telah  yang telah tererosi oleh gelombang, dan terus membawanya sepanjang pantai. Proses ini dikenal sebagai transpor sedimen  sepanjang  pantai.  Untuk  memprediksi  transpor  sedimen  yang  disebabkan  oleh  longshore  current, digunakan  rumus  empiris  yang  berdasarkan  kondisi  gelombang  di  lokasi  penelitian.  Rumus  tersebut merupakan  hubungan  sederhana  antara  transpor  sedimen  dan  fluks  energi  gelombang  sepanjang  pantai. Menurut BAPPEDA Kota Mataram (2019), kondisi muara 1 sungai di Kota Mataram, terutama dengan adanya pengendapan 1 lidah pasir, menunjukkan adanya 1 pola angkutan sedimen 1 sejajar pantai (longshore transport) . Secara umum, orientasi garis pantai Kota Mataram menghadap 1 ke arah barat, sehingga arah gelombang 1 yang membentuk  sudut  adalah  dari  barat  laut 1 atau  barat  daya.  Dengan demikian, pola  angkutan 1 sedimen  sejajar pantai di Kota 1 Mataram adalah ke 1 arah utara.

## Perubahan Garis Pantai
Data  hasil  analisis  perubahan 8 garis 8 pantai  diperoleh  dari  hasil 8 analisis 8 DSAS (Digital 8 Shoreline Analysis 8 System) yang dibagi menjadi 2 yaitu analisis perubahan 8 garis 8 pantai pada tahun 2012 – 2017 dan analisis  perubahan 8 garis 8 pantai  pada  tahun  2017  –  2022.  Citra  yang 8 digunakan  pada  tahun  2012  adalah citra 8 Landsat 7 dan pada tahun 2017 sampai 2022 menggunakan citra Landsat 8. Berdasarkan analisis dengan metode DSAS diperoleh nilai EPR (End Point Rate) untuk 1 mengetahui laju 8 perubahan 8 garis pantai dan nilai NSM (Net 8 Shoreline 8 Movement) untuk mengetahui panjang jarak 1 perubahan garis pantai. Pada tabel 9. dapat diketahui bahwa nilai EPR tertinggi pada tahun 2012 – 2017 untuk akresi mencapai 4,57 m/tahun dan terendah mencapai 0,5 m/tahun dengan rata-rata 2,15 m/tahun. Sedangkan untuk abrasi tertinggi mencapai 9,19 m dan terendah 0,13 dengan rata-rata 3,23 m. Nilai  EPR tertinggi pada tahun 2017  – 2022 untuk akresi mencapai 1,72  m/tahun  dan  terendah  mencapai  0,31  m/tahun  dengan  rata-rata  0,86  m/tahun. Sedangkan untuk  abrasi tertinggi mencapai 10,94 m/tahun dan terendah 0,10 m/tahun dengan rata-rata abrasi 2,94 m/tahun.
Tabel 8 9. Hasil EPR 8 Pengolahan 8 DSAS Tahun EPR (End Point Rate) (m/th) Akresi 8 Akresi 8 Abrasi 8 Abrasi 8 Rata - Rata 8 Min 8 Max 8 Min 8 Max 8 Akresi 8 Abrasi 8 2012 - 2017 0,5 4,57 0,13 9,19 2,15 3,23 2017 - 2022 0,31 1,72 0,10 10,94 0,86 2,94
Pada tabel  10. dapat  diketahui  bahwa  nilai  NSM  tertinggi  pada  tahun  2012  –  2017  untuk  akresi mencapai 25,32 m dan terendah mencapai 2,78 m dengan rata – rata 11,91 m. Sedangkan untuk abrasi tertinggi mencapai 50,93 m dan terendah sebesar 0,70 dengan rata-rata 17,91 m. Kemudian nilai NSM tertinggi pada tahun  2017 –  2022  untuk  akresi  mencapai  8,81  m  dan  terendah  mencapai  1,57 m  dengan  rata-rata  4,40  m. Sedangkan untuk abrasi tertinggi mencapai 56,08 m/tahun dan terendah 0,51 m dengan rata-rata 15,09 m.
Tabel 10. Hasil 8 NSM 8 Pengolahan 8 DSAS Tahun NSM 8 (Net 8 Shoreline 8 Movement) (m) Akresi 8 Akresi 8 Abrasi 8 Abrasi 8 Rata - Rata 8 Min 8 Max 8 Min 8 Max 8 Akresi 8 Abrasi 2012 – 2017 2,78 25,32 0,70 50,93 11,91 17,91 2017 – 2022 1,57 8,81 0,51 56,08 4,40 15,09
Hasil  yang  diperoleh  dari  perhitungan  metode  DSAS  diketahui  perubahan  garis  pantai  yang  terjadi berdasarkan nilai total luasan 8 area yang terjadi abrasi 8 dan 8 akresi. Pada 1 tabel 11. dapat dilihat 1 bahwa pada
tahun 2012 – 2017 terjadi abrasi dengan nilai luas sebesar 2,512 ha 8 dan 8 akresi 8 sebesar 0,296 ha. Sedangkan pada 8 tahun 2017 – 2022 terjadi abrasi dengan nilai luas sebesar 1,892 ha dan akresi sebesar 0,114 ha.
Tabel 11. Total 8 Luas 8 Area Abrasi 8 dan 8 Akresi Tahun Total 8 Luas Area 8 (ha) Abrasi Akresi 2012 – 2017 2,512 0,296 2017 – 2022 1,892 0,114
Gambar 5(a). Peta Analisis DSAS di Pantai Bagik Kembar Tahun 2012 – 2017
Gambar 5(b). Peta Analisis DSAS di Pantai Bagik Kembar Tahun 2012 – 2017 Berdasarkan  hasil  analisis  DSAS  perubahan  garis  pantai  di  Pantai  Bagik  Kembar  dan  sekitarnya didominasi oleh terjadinya erosi atau abrasi pantai. Pada tahun 2012 – 2017 nilai luas area abrasi sebesar 2,512 ha dan akresi sebesar 0,296 ha. Pada tahun 2017 – 2022 mengalami hal yang sama dimana pantai didominasi
terjadinya abrasi. Nilai luas area terjadinya abrasi pada tahun ini sebesar 1,892 ha dan akresi sebesar 0,114 ha. Dengan  nilai  tersebut  dapat  diketahui  bahwa  daerah  ini  telah  terjadi  perubahan 1 bentuk  lahan  garis  pantai mundur ke 1 arah darat. Faktor penyebab utama terjadinya perubahan garis pantai di lokasi penelitian ini adalah gelombang laut yang dibangkitkan oleh angin dan pasang surut. Angin yang bertiup di atas permukaan laut dapat  mendorong  air  laut,  menciptakan  gelombang  yang  kemudian  mencapai  pantai.  Menurut  Putra et  al. (2022),  gerak  air  yang  disebabkan  gelombang  akan  membentuk  suatu  lintasan  yang  juga  ikut  membawa partikel-partikel  sedimen  dan  mengakibatkan  terjadinya  transpor  sedimen  dalam  arah  sejajar  garis  pantai sehingga  terjadi  arus  sejajar 1 pantai  atau longshore  current. Arus  ini  mengakibatkan 1 terangkutnya  sedimen dasar dari pesisir pantai 1 dan akan berpotensi menjadi 1 penyebab terjadinya erosi 1 maupun sedimentasi. Selain itu, pasang surut dapat menjadi faktor yang berperan dalam perubahan garis pantai. Pasang surut merupakan naik-turunnya permukaan air laut secara periodik akibat gaya gravitasi antara bumi, bulan, dan matahari. Saat pasang, air laut naik dan mencapai pesisir, sedangkan saat surut, air laut turun dan menarik kembali sedimen pantai.  Dari  hasil  analisis  pasang  surut  diperoleh  muka  air  laut  tertinggi  sebesar  2,4  m  dan  muka  air  laut terendah sebesar 0,12 m. Kondisi pasang surut perairan Kota Mataram memiliki tipe pasang surut campuran condong  harian  ganda,  dimana  dalam  satu  hari  dapat  terjadi  dua  kali  pasang  dan  surut,  dengan  tinggi  dan periode yang berbeda. Menurut Setyawan et al. (2021), pasang naik dapat menyebabkan sedimen mendekati pantai, sementara saat surut, sedimen cenderung maju ke arah laut lepas. Perbedaan dalam ketinggian pasang surut dapat mempengaruhi perubahan garis pantai, baik dalam bentuk akresi (penambahan material sedimen) maupun abrasi (erosi material sedimen). Tetapi perbedaan rentang ketinggian muka air laut dapat dipengaruhi oleh  faktor  geomorfologi  juga,  termasuk  kemiringan  pantai.  Kemiringan  pantai  yang  curam  cenderung menghasilkan perbedaan rentang ketinggian yang lebih besar. Namun, faktor lain seperti topografi dasar laut dan pengaruh alam sekitar juga dapat berperan dalam menentukan besar kecilnya perbedaan rentang ketinggian muka air laut di daerah penelitian.
Proses ini dapat terjadi secara alami dan merupakan karakteristik umum dari dinamika garis pantai di berbagai lokasi pesisir. Berdasarkan 1 penelitian yang dilakukan oleh BAPPEDA Kota Mataram (2019), garis pantai Kota Mataram, termasuk di daerah Pantai Bagek Kembar dan Ampenan, mengalami kemunduran secara bertahap antara tahun 2006 hingga 2017. Abrasi pantai dapat terjadi dalam waktu singkat saat gelombang besar dengan arah tegak lurus terhadap pantai. Sedimen yang pergi ke laut dalam 1 diperkirakan tidak mungkin terjadi karena 1 kondisi bathimetri 1 perairan Kota Mataram yang landai. Namun sedimen 1 yang pergi ke daerah pantai lebih  tinggi  karena 1 overwash terjadi  pada  beberapa  tempat  seperti  di 1 Pantai  Bagek  Kembar.  Hal  ini  dapat menyebabkan banjir akibat luapan air laut dan juga dapat menimbun rumah-rumah nelayan di Pantai Bagek Kembar.

## KESIMPULAN
Berdasarkan hasil analisis data pada penelitian ini, dapat disimpulkan bahwa Pantai Bagik Kembar dan sekitarnya mengalami perubahan garis pantai yang dinamis selama rentang waktu 2012 hingga 2022. Abrasi merupakan fenomena yang dominan, dengan luas abrasi mencapai 2,512 ha pada tahun 2012-2017 dan 1,892 ha pada tahun 2017-2022. Sementara itu, akresi yang terjadi relatif lebih kecil, yaitu sebesar 0,296 ha pada tahun 2012-2017 dan 0,114 ha pada tahun 2017-2022. Faktor yang mendominasi terjadinya perubahan garis pantai  pada  lokasi  penelitian  ini  adalah  gelombang  yang  dibangkitkan  oleh  angin  sehingga  menyebabkan terjadinya pengangkutan sedimen oleh arus sejajar pantai.
Informasi  ini  memberikan  pemahaman  tentang  dinamika  garis  pantai  di  Pantai  Bagik  Kembar  dan sekitarnya serta faktor-faktor yang berkontribusi pada perubahan garis pantai tersebut. Penelitian ini penting dalam mengembangkan strategi pengelolaan wilayah pesisir yang rentan terhadap abrasi dan meminimalkan risiko bencana bagi pemukiman dan aktivitas manusia di sekitar pantai tersebut.

## DAFTAR PUSTAKA
BAPPEDA  Kota  Mataram.  2019.  Pengelolaan  Garis  Pantai  Kota  Mataram  sebagai  Upaya  Adaptasi  dan Mitigasi  Bencana  Pantai  dan  Kenaikan  Muka  Air  Laut.  PEER  USAID  Project  #5-395.  Kerjasama TDMRC Univ. Syiah Kuala -Universitas Mataram-BAPPEDA Mataram.
Holme, & McIntyre. 1984. An Introduction to Coastal. New York: Harper and Row Publisher. Kristanto, O., Ismanto, A., Satriadi, A., Setiyono, H., dan Atmodjo, W. 2022. Analisis Longshore Current dan Pengaruhnya  Terhadap  Transpor  Sedimen  Dasar  di  Perairan  Teluk  Awur,  Jawa  Tengah. Indonesian Journal of Oceanography , 4(1): 59-68.
Kusumawardani,  K.  P.,  Cahya,  Z.  I.,  Ananto,  W.  H.  G.,  dan  Asri,  G.  H.  M.  2019.  Pemetaan  Dan  Analisis Perubahan Garis Pantai Di Sebagian Pesisir Barat Lombok Barat Menggunakan Normalized Difference Water Index Pada Citra Landsat. Seminar Nasional Geomatika, 3: 911-918.
Pratiwi, A. A., Setiyono, H., Suryoputro, A. A. D., Marwoto, J., & Satriadi, A. 2020. Perubahan Garis Pantai Pada  Morfologi  Gisik  Kantung  di  Pantai  Baron,  Kabupaten  Gunungkidul  Daerah  Istimewa Yogyakarta. Indonesian Journal of Oceanography , 2 (4): 363-369.
Putra, A. N., Handoyo, G., Ismanto, A., Satriadi, A., dan Setiyono, H. 2022. Studi Pengaruh Longshore Current Terhadap  Transpor  Sedimen  Dasar  di  Pantai  Slamaran,  Kota  Pekalongan,  Jawa  Tengah. Indonesian Journal of Oceanography , 4 (1): 36-46.
Ryanto, N. A., Wiyono, R. U. A., dan Hidayah, E. 2022. Studi Peramalan Pasang Surut di Pesisir Pantai Pancer Kecamatan Puger, Kabupaten Jember. Maspari Journal Marine Science Research , 14(1): 1-13.
Setyawan, F.O., Sari, W.K., dan Aliviyanti, D. 2021. Analisis Perubahan Garis Pantai Menggunakan Digital Shoreline  Analysis  System  Di  Kecamatan  Kuala  Pesisir,  Kabupaten  Nagan  Raya,  Aceh. Journal  of Fisheries and Marine Research, 5(2): 368-377. Sugiyono. 2009. Metode Penelitian Kuantitatif Kualitatif dan R&D. ALFABETA, Bandung. Sulaiman, D.M. 2021. Teknologi Pegar untuk Penanggulangan Erosi dan Abrasi Pantai. Deepublish Publisher, Yogyakarta.
Triatmodjo, B. 1999. Teknik Pantai. Beta Offset, Yogyakarta. Widada, S., Ismanto, A., Priambodo, I. B., dan Siagian, H. S. R. 2022. Perubahan Garis Pantai dan Dampaknya Terhadap Banjir Rob di Kecamatan Pekalongan Utara, Kota Pekalongan, Provinsi Jawa Tengah.