## PEMODELAN TINGGI GELOMBANG TSUNAMI MENGGUNAKAN DELFT3D DI PANTAI BANDEALIT JEMBER
Baby Citra Nusantari Wardana a , Retno Utami Agung Wiyono b, * , Gusfan Halik b
a Program Studi S1 Teknik Sipil, Fakultas Teknik, Universitas Jember, Jember
b Jurusan Teknik Sipil, Fakultas Teknik, Universitas Jember, Jember *Corresponding author, email address: retnoutami@unej.ac.id
A R T I C L E     I N F O A B S T R A C T Article History: Bandealit Beach is a tourist attraction that is still very beautiful because it is included in the Meru Betiri National Park area which is located in Tempurejo  District,  Jember  Regency. Apart  from  that,  what  makes  the beach  still  beautiful  is  the  infrequent  number  of  tourists  visiting  these tourist attractions due to very limited access. However, behind the beauty of the beach, quite a lot of people live around the beach and there are also busy activities such as schools. Its location which is located on the south side of the island of Java directly facing the Indian Ocean is of course very possible for a tsunami to occur, considering that this beach was once affected by the Banyuwangi Tsunami in 1994. Therefore, this study aims to simulate a tsunami at Bandealit Beach using Delft3D-Flow to analyze the travel time and height of the tsunami. Tsunami modeling was carried out based on two scenarios, namely scenario 1 as model validation using the characteristics of the 1994 Banyuwangi Tsunami. Model validation calculations  used  the  MAPE  method  to  interpret  the  error  percentage <10%. Scenario 2 models the worst PuSGeN earthquake with a moment magnitude of 8.7 according to BMKG. In modeling, 3 observation points are used. The simulation results show the height of the tsunami waves at each point based on the Tsunami Banyuwangi 1994, namely 8.7 m, 9.1 m, and 10.4 m with  a travel time around 32 -35 minutes. Then for the worst  earthquake  PuSGeN  Mw  8.7,  the  tsunami  heights  were  15.4  m, 16.4 m, and 17.3 m with a travel time of 20 minutes.
Received 20 January 2023 Accepted 16 February 2024 Online 30 Juni 2024

## Keywords:
Tsunami simulation Delft3D Bandealit beach
© 2024 Magister Teknik Sipil Unsyiah. All rights reserved

## 1. PENDAHULUAN
Kabupaten  Jember  merupakan  salah  satu  daerah  di  Jawa  Timur  yang  dikenal  memiliki  keindahan pantainya  yaitu  seperti  Pantai  Bandealit.  Selain  keindahannya,  dalam  sejarah  pantai  tersebut  pernah mengalami tsunami pada Juni 1994 yang merupakan dampak dari tsunami yang terjadi di Dusun Pancer, Desa  Sumberagung,  Kecamatan  Pesanggran,  Banyuwangi.  Keruskaan  yang  ditimbulkan  akibat  bencana tersebut  cukup  parah yaitu  73  rumah rusak  dan  5  orang  meninggal,  dengan  ketinggian  run-up  10  meter untuk daerah Bandealit (Wisyanto, 2003).
Pantai  Bandealit  merupakan  objek  wisata  yang  berlokasi  di  Taman  Nasional  Meru  Betiri,  Desa Andongrejo,  Kecamatan  Tempurejo,  Kabupaten  Jember.  Jumlah  penduduk  Desa  Andongrejo  yaitu sebanyak 5974 jiwa dengan persentase penduduk sebesar 7.38% (Badan Pusat Statistik, 2023).
Potensi terjadinya tsunami tentu saja masih sangat mungkin terjadi, mengingat letak pantai tersebut berada di sisi selatan Kabupaten Jember yang secara langsung berhadapan dengan Samudera Hindia dan dikenal dengan tempat  pertemuan dua lempeng tektonik,  yaitu Lempeng Indo – Australia dan Lempeng Eurasia yang masih sangat aktif secara tektonik dan dapat menjadi sumber bencana salah satunya tsunami.

## Jurnal Arsip Rekayasa Sipil dan Perencanaan (JARSP) Journal of Archive in Civil Engineering and Planning
E-ISSN: 2615-1340; P-ISSN: 2620-7567
Journal homepage: http://www.jurnal.unsyiah.ac.id/JARSP/index
Simulasi  tsunami  telah  dilakukan  pada  beberapa  penelitian  di  pantai  Jember  dengan  menggunakan parameter  gempa tektonik yang berbeda. Armono dkk. (2021),  memodelkan tsunami  dengan parameter gempa berkekuatan M 9.1 yang merupakan skenario terburuk  dipantai selatan Jawa Timur  dengan hasil ketinggian  tsunami  hingga  20  meter.  Pratama  (2017)  melakukan  simulasi  tsunami  Pantai  Watu  Ulo menggunakan  parameter  gempa  dari  pemodelan  Monte  Carlo  dengan  kekuatan  8,9  SR.  Ketinggian penyetelan  tertinggi  dari  titik  pengamatan  adalah  5,52  meter,  dan  waktu  tempuh  34  menit.  Sedangkan Rikarda  dkk.  (2020)  memodelkan  tsunami  Pantai  Puger  menggunakan  Delft3D  dengan  kombinasi karakteristik sumber gempa di Jawa Selatan dan menghasilkan tinggi tsunami 22,34 meter dalam waktu 40 menit.  Penelitian  -  penelitian  lain  menggunakan  Delft3D  juga  pernah  dilakukan  oleh  MA  ddk.  (2020), Kurniawan dkk. (2020), dan Kartini dkk. (2020).
Berdasarkan  latar  belakang  dan  penelitian  sebelumnya,  belum  ada  yang  melakukan  pemodelan tsunami di Pantai Bandealit dengan skrenario Tsunami Banyuwangi 1994 dan PuSGeN Mw 8.7. di sekitar lokasi  tersebut  masih  banyak  warga  dan  komoditas  yang  berlangsung  sampai  saat  ini.  Oleh  karena  itu, penelitian  ini  bertujuan  untuk  mensimulasikan  gelombang  tsunami  menggunakan  karakteristik  gempa terburuk  PuSGeN  di  Pantai  Bandealit  Kabupaten  Jember  dengan  memanfaatkan model DelftDashboard dan Delft3D . Validasi model dilakukan dengan menggunakan data observasi Tsunami Banyuwangi 1994 dari  penelitian  sebelumnya.  Model  tersebut  akan  menghasilkan  beberapa  komponen  seperti  tinggi gelombang dan waktu tiba. Diharapkan dengan mengetahui komponen tersebut dapat berkontribusi dalam mitigasi bencana tsunami.

## 2. KAJIAN PUSTAKA

## 2.1 Definisi Tsunami
Tsunami berasal dari bahasa Jepang, Tsu berarti pelabuhan dan Nami berarti gelombang. Jadi, tsunami dapat diartikan sebagai peristiwa gelombang laut yang tinggi dan besar yang tiba di dasar laut di wilayah pesisir beberapa saat setelah terjadi gempa bumi, letusan gunung berapi, dan longsor di dasar laut. Tsunami adalah gelombang air yang sangat besar yang disebabkan oleh berbagai gangguan di dasar laut. Tsunami tidak terlihat ketika berada jauh di tengah lautan. Namun ketika mencapai daerah yang dangkal, gelombang bergerak lebih cepat dan semakin besar (Hartuti, 2009).

## 2.2 Proses Pembangkitan Gelombang Tsunami
Tsunami  dapat  terjadi  ketika  gempa  bumi,  letusan  gunung  berapi,  dan  longsor  di  dasar  laut,  atau tumbukan meteorit. Pergerakan air laut dan deformasi vertikal dasar laut menyebabkan perbedaan tinggi muka  air  laut.  Proses  penyeimbangan  ini  menimbulkan  gelombang  laut  yang  diperkirakan  setinggi  50 meter. Pada saat yang sama, panjang gelombang mencapai ribuan kilometer, dengan kecepatan gelombang yang  bergerak  500  hingga  1.000  Km  per  jam.  Gelombang  laut  biasanya  akan  menghantam  pantai  atau pelabuhan terdekat dalam waktu 10 sampai 30 menit setelah tsunami terjadi (Hartuti, 2009).

## 2.3 Patahan
Patahan atau sesar adalah zona patahan antara dua blok batuan. Patahan  memungkinkan blok untuk bergerak relatif satu sama lain. Pergerakan ini dapat terjadi dengan cepat, dalam bentuk gempa bumi atau dapat  terjadi  secara  perlahan,  dalam  bentuk  rambatan.  Panjang  patahan  dapat  berkisar  dari  beberapa milimeter hingga ribuan kilometer. Sebagian besar patahan menghasilkan perpindahan berulang sepanjang waktu. Permukaan patahan bisa horizontal atau vertikal atau beberapa sudut yang berubah – ubah. (USGS, 2022).

## 3. METODE PENELITIAN
Simulasi  tsunami  dilakukan  dengan  melakukan  pemodelan  hidrodinamik  menggunakan DelftDashboard  dan  Delft3D-Flow.  Penelitian  ini  dimulai  dengan  mengumpulkan  data  yang  diperlukan
untuk  simulasi.  Pertama,  grid  dari  domain  simulasi  numerik  dibuat  menggunakan  DelftDashboard. Penelitian  ini  menggunakan  ukuran  grid  yang  berbeda  di  wilayah  pesisir  dan  lepas  pantai  dengan menyesuaikan penggunaan dua data batimetri dengan resolusi yang berbeda. Resolusi grid pada dasarnya menentukan  keakuratan  dan  efisiensi  hasil.  Sejumlah  besar  sel  grid  dengan  resolusi  grid  tinggi  akan membutuhkan  waktu  dan  daya  komputasi  yang  tinggi.  Dengan  demikian,  resolusi  grid  yang  tepat diperlukan untuk mendapatkan hasil yang stabil dan waktu komputasi yang efisien. Setelah domain simulasi disiapkan, parameter patahan dari gempa yang ditentukan dimasukkan ke dalam toolbox DelftDashboard untuk menghasilkan tsunami. Proses interpolasi, pengaturan parameter, model running, dan tampilan hasil model menggunakan Delft3D.
Pemodelan tsunami dalam penelitian ini menggunakan dua skenario simulasi. Skenario 1 merupakan simulasi  gelombang  tsunami  berdasarkan  karakteristik  Tsunami  Banyuwangi  1994.  Skenario  2  adalah simulasi tsunami menggunakan gempa terburuk PuSGeN mw 8.7 dari BMKG. Simulasi kedua skenario ini menghasilkan  elevasi  muka  air  yang  dianalisis  untuk  mendapatkan  tinggi  dan  waktu  tiba  gelombang tsunami di setiap titik pengamatan yang telah ditentukan.

## 3.1 Skenario 1
Skenario  ini  adalah  model  validasi.  Simulasi  tsunami  menggunakan  karakteristik  Tsunami Banyuwangi tahun 1994. Data karakteristik gempa yang digunakan dalam skenario ini diperoleh dari USGS dan disajikan pada Tabel 1.
Tabel 1 . Karakteristik gempa Skenario Long. (  ) Lat. (  ) Mag. (Mw) Depth (Km) Strike (  ) Dip (  ) Slip (  ) Tsunami Banyuwangi 1994 112.835 -10.477 7.8 11.5 97 83 91 PuSGeN Mw 8,7 112.689 -10.193 8.7 12 275 7 90

## Sumber: USGS dan BMKG
Tabel  1  menunjukkan  karakteristik  gempa  yang  digunakan  sebagai  parameter  input  pada  menu Tsunami DelftDashboard untuk  membangkitkan  tsunami  dalam  simulasi.  Bujur  dan  lintang  adalah koordinat pusat gempa. Magnitudo mewakili kekuatan gempa, kedalaman mewakili kedalaman hiposenter, sedangkan strike, dip, dan slip adalah arah gerak patahan.
Magnitudo gempa adalah salah satu aspek terpenting dari parameter sumber gempa yang dapat diukur dan  dikuantifikasi  secara  langsung.  Umumnya,  pada  permulaan  gempa  bumi,  magnitudo  merupakan parameter  sumber  pertama  yang  dikuantifikasi  untuk  menentukan  kekuatan  gempa.  Magnitudo  tersebut dapat dimanfaatkan untuk mempelajari pola kejadian gempa. Hubungan empiris antara luas sesar, panjang sesar, dan lebar sesar dapat diperoleh dengan menggunakan pendekatan regresi, yaitu persamaan regresi empiris  yang  menghubungkan  dimensi  sesar  dengan  magnitudo  gempa.  Oleh  karena  itu,  perhitungan panjang patahan (L), lebar patahan (W), dan dislokasi (D) pada model ini menggunakan persamaan regresi dari Wells dan Coppersmith (1994). Pada simulasi ini digunakan nilai L, W, dan D masing – masing sebesar 144.6 Km, 22.2 Km, dan 20 m untuk skenario 1. Kemudian  untuk scenario 2 nilai L, W, dan D diperoleh dari BMKG yaitu  masing – masing 141 Km, 91 Km, dan 27 m.
Untuk  Skenario  1,  tinggi  gelombang  Tsunami  Banyuwangi  1994 yang  disimulasikan  dibandingkan dengan tinggi gelombang tsunami dari pengamatan lapangan Maramai dan Tinti (1997). Dengan demikian, penentuan titik pengamatan Skenario 1 dipilih berdasarkan pengamatan lapangan Maramai dan Tinti yang meliputi Pantai Watu Ulo, Pantai Puger, dan Tanjung Pelindu, Pantai Sukamade, Pantai Rajegwesi, Pantai Pancer,  dan  Pantai  Pulau  Merah.  Validasi  ini  dihitung  dengan  menggunakan  MAPE  (Mean  Absolute Percentage Error), yang ditunjukkan pada Persamaan. 1
= 1 − 100% (1)
Dengan n adalah  jumlah  observasi, x oi dan x pi adalah  tinggi  gelombang  observasi  dan  simulasi. Interpretasi nilai MAPE diklasifikasikan menjadi 4 kategori: <10% sangat akurat, 10-20% baik, 20-50% wajar,  dan  >50%  tidak  akurat.  Dengan  demikian,  persentase  validasi  yang  ditentukan  dalam  model  ini kurang dari 10%.

## 3.2  Skenario 2
Skenario 2 merupakan simulasi tsunami dengan kemungkinan gempa terburuk. Pada simulasi ini titik pengamatan difokuskan di Pantai bandealit. Simulasi pada Skenario 2 dilakukan dengan mengubah input data  parameter  gempa tektonik.  Parameter  gempa  yang  digunakan  pada  skenario  ini  adalah  data  gempa PuSGeN yang diperoleh dari  BMKG yang dapat dilihat pada Tabel 1. Berdasarkan Tabel  1, Magnitudo pada skenario ini adalah 8.7 Mw yang juga dapat disebut sebagai gempa yang mungkin terjadi. Skenario gempa terburuk tentunya berpotensi menimbulkan dampak yang cukup besar, terutama di wilayah pesisir.

## 3.3   Pemodelan Hidrodinamik
Model Delft3D-Flow ,  program  simulasi  hidrodinamika  multidimensi  (2D  atau  3D)  yang dikembangkan oleh Deltares (2014) digunakan untuk mensimulasikan tsunami berdasarkan penetuan dua skenario.  Penelitian  ini  melakukan  komputasi  dua  dimensi  dengan  komponen  horizontal  yang  sama (kalkulasi  rata-rata  kedalaman  dua  dimensi)  karena  untuk  simulasi  tsunami,  skala  panjang  dan  waktu horizontal  jauh  lebih  signifikan  daripada  skala  vertikal.  Persamaan  yang  mengatur Delft3D-Flow ditunjukkan  pada  Persamaan.  (2),  (3),  dan  (4).  Persamaan  momentum  untuk  x  dan  y  ditunjukkan  pada Persamaan.  (2)  dan  (3)  sedangkan  persamaan  kontinuitas  rata-rata  kedalaman  2D  ditunjukkan  pada Persamaan. (4), seperti yang dijelaskan dalam Deflt3D Flow Manual. Setelah file input disiapkan, kondisi batas  dan  awal  ditentukan.  Persamaan.  (2),  (3),  dan  (4)  digunakan  pada  langkah  awal  perhitungan hidrodinamika  hingga  simulasi  selesai.  Setelah  menghitung  persamaan,  tinggi  gelombang  tsunami  pada setiap node untuk semua langkah waktu dapat diperoleh.
+ + + + − + −
= − 1
+ + 1 (   +   ) +
(2) + + + + − − +
= − 1 + + 1 (   +   ) + (3) + 1 (   +   ) + 1 (   +   ) = (   +   ) (4)
Dalam  Persamaan.  (2),  (3),  dan  (4),  ada  beberapa  istilah  yang  dijelaskan  di  bawah  ini.  U  dan  V didefinisikan sebagai kecepatan rata-rata kedalaman.
= 1 + = (5)
= 1 + = (6)
Q merupakan kontribusi tiap satuan luas karena pemborosan air, presipitasi, dan evaporasi.
= ( − ) +  – (7)
dan η merepresentasikan  gradien  tekanan, adalah  koefisien  viskositas  eddy  vertikal  (m 2 /s). Sementara itu, dan η mewakili ketidakseimbangan horizontal dari tekanan Reynold (m/s 2 ), dan η mewakili  kontribusi  sumber  daya  eksternal  atau  inundasi  momentum.    adalah  parameter  Coriolis (frekuensi inersia) (1/s),   adalah rembesan aliran menuju x (m/s),   kecepatan fluida menuju y (m/s), dan kecepatan menuju  z (m/s). Selanjutnya ζ adalah elevasi  permukaan air dengan t  sebagai  waktu dan d sebagai kedalaman (m), Gηη =  R dan Gξξ = sebagai ϕ adalah garis lintang dan R adalah jari- jari  Bumi  (6378.137  km,  WGS  84),  juga    dan    adalah  koordinat  kurva  horizontal,  q in dan  q out merepresentasikan sumber lokal dan resapan air per satuan volume [1/s], P adalah air hujan, dan E adalah genangan non-lokal karena penguapan.
Pemodelan numerik dimulai dengan input data dan pengaturan parameter pada DelftDashboard . Data yang digunakan dalam pemodelan ini adalah data Batimetri Nasional dengan resolusi 6 arc-second (kurang lebih 180 meter) yang diperoleh dari Badan Informasi Geospasial untuk wilayah yang dekat dengan garis pantai. Sedangkan di laut lepas pantai digunakan data GEBCO_2021 dengan resolusi 15 arc-second (kurang lebih  450 meter) yang diperoleh dari BODC (British Oceanographic Data Center). Berdasarkan resolusi batimetri, domain pemodelan yang digunakan dalam penelitian ini adalah grid persegi berukuran 0,005 x 0,005 derajat setara dengan 500x500 meter untuk area lepas pantai, dan 0,0025 derajat setara dengan 250 meter  untuk  area  pantai.  Selanjutnya,  parameter  sesar  gempa  ditentukan  pada  toolbox  Tsunami DelftDashboard  berdasarkan  data  gempa  pada  Tabel  1  untuk  mendapatkan  elevasi  permukaan  awal. Simulasi dijalankan berdasarkan dua skenario yang telah ditentukan. Keluaran model dari hasil simulasi dapat divisualisasikan menggunakan Quickplot di Delft3D, yang akan menampilkan elevasi permukaan air untuk analisis lebih lanjut.

## 4. HASIL DAN PEMBAHASAN

## 4.1 Pemodelan Numerik Skenario 1
Simulasi  dalam  skenario  ini  bertujuan  untuk memvalidasi  model  yang digunakan  dalam  penelitian. Hasil simulasi akan menampilkan elevasi permukaan air yang digunakan sebagai parameter analisis untuk mendapatkan tinggi gelombang di setiap titik pengamatan. Gambaran lokasi titik pengamatan pada model ditunjukkan pada Gambar 1.
Tanjung Pelindu Puger Watu Ulo Bandealit Sukamade Rajegwesi Pancer Pulau Merah
Gambar  1  menunjukkan  lokasi  pengamatan  pada  model  yang  diilustrasikan  pada  DelftDashboard. Pembuatan  domain  model  untuk  simulasi  tsunami  meliputi  lokasi  –  lokasi  tersebut.  Parameter  tsunami untuk  Skenario  1  menggunakan  peristiwa  Tsunami  Banyuwangi  1994,  menghasilkan  elevasi  muka  air spasial yang ditunjukkan pada Gambar 2 (a-e).
Gambar 2. Elevasi muka air spasial skenario 1: a) elevasi awal, b) 10 menit, c) 20 menit, d) 30 menit, e) 40 menit.
Gambar  2  (a-e)  menunjukkan  elevasi  muka  air  maksimum  secara  keseluruhan  dalam  domain  yang dimodelkan.  Gambaran  ketinggian  air  dalam  penelitian  ini  diatur  untuk  ditampilkan  setiap  10  menit. Dimana  x  pada  gambar  menunjukkan  daerah  koordinat  bujur  dan  y  adalah  koordinat  lintang.  Gempa tersebut memicu tinggi muka air awal tsunami, ditunjukkan pada Gambar 2a. Berdasarkan legenda pada gambar, data gempa yang diperoleh dari USGS menghasilkan gelombang awal setinggi 9 meter. Setelah 10 menit (Gambar 2b), tinggi gelombang menjadi 15 meter. Tinggi gelombang terus merambat ke arah daratan setelah 20 menit. Akibatnya, tinggi gelombang pada Gambar 2c berkurang menjadi 5 meter. Semakin dekat ke  daratan  dalam  30 menit,  gelombang  tsunami  pada Gambar  2d meningkat  menjadi  10,5 meter  karena semakin  dekat  ke  pantai,  permukaan  air  di  laut  semakin  dangkal,  sehingga  gelombang  semakin  tinggi.
b)
c) d) e) a)
Akibatnya, dalam waktu 40 menit gelombang mencapai pantai, menghasilkan tinggi gelombang maksimum sekitar 9,5 meter, seperti terlihat pada Gambar 2e.
Berdasarkan elevasi muka air secara keseluruhan, dilakukan analisis untuk mendapatkan waktu tiba dan  tinggi  gelombang  maksimum  pada  setiap  titik  pengamatan  yang  telah  ditentukan.  Kedua  parameter tersebut  diperoleh  dari  hasil  perhitungan  hidrodinamika Delft3D-Flow melalui  Persamaan.  2  dan  3. Terakhir,  tinggi  gelombang  maksimum  tersebut  digunakan  untuk  memvalidasi  model  dengan membandingkannya dengan tinggi gelombang terukur di observasi lapangan. Persentase kesalahan dihitung menggunakan MAPE (Persamaan 1). Hasil perhitungan validasi ditunjukkan pada Tabel 2.
Tabel 2. Hasil validasi pemodelan skenario 1
Lokasi Ketingian Tsunami (m) (Maramai dan Tinti, 1997) Ketinggian Tsunami (m) (Hasil Simulasi) Error (%) Tanjung Pelindu Puger Watu Ulo Sukamade Rajegwesi Pancer Pulau Merah 3.2 4.88-5.85 6.5-7.5 6.1-6.25 6.3-9.5 4.55-4.7 4.35-6.15 3.12 4.99 7.28 7.11 7.60 5.03 6.12 2.5 14.7 2.93 13.76 20 7.02 0.48 MAPE 8.77
Sumber: Maramai dan Tinti (1997)
Berdasarkan Tabel 2, persentase kesalahan rata-rata dari seluruh titik pengamatan adalah 8,77%. Nilai ini  memenuhi  ketentuan  validasi  nilai  MAPE,  yaitu  kurang  dari  10%.  Oleh  karena  itu,  pemodelan gelombang tsunami ini dapat diterima.

## 4.2 Pemodelan Numerik Skenario 2
Skenario  ini  merepresentasikan  hasil  simulasi  tsunami  dengan  menggunakan  parameter  gempa terburuk  PuSGeN  dari  BMKG.  Keluaran  yang  dihasilkan  berupa  elevasi  muka  air  dianalisis  untuk mendapatkan  waktu  tiba  tinggi  gelombang  tsunami  pada  setiap  titik  pengamatan  yang  telah  ditentukan terutama  difokuskan  di  Pantai  Bandealit  .  Gambar  3  dan  Gambar  4(a-d)  menunjukkan  kondisi  spasial elevasi permukaan air dari simulasi gelombang tsunami Skenario 2.
Gambar 3. Elevasi awal muka air spasial skenario 2
Gambar 4. Elevasi muka air spasial skenario 2: a) 10 menit, b) 20 menit, c) 30 menit, d) 40 menit
Berdasarkan  legenda  masing-masing  Gambar  3  dan  Gambar  4,  terlihat  bahwa  elevasi  muka  air maksimum  keseluruhan  pada  model  mencapai  10  hingga  25  meter,  karena  visualisasi  ini  diatur  untuk ditampilkan setiap 10 menit  di mana x adalah area koordinat bujur  dan y adalah area koordinat lintang. Gambar 3 menunjukkan elevasi muka air awal 10 meter. Setelah 10 menit, gelombang maksimum yang ditunjukkan  pada  Gambar  4a  menghasilkan  tinggi  gelombang  10.5  m.  Gempa  bumi  pada  skenario  ini bermagnitudo  besar  yaitu  8.7  Mw,  sehingga  dalam  waktu  20  menit  pada  Gambar  4b  menimbulkan gelombang yang tingginya mencapai 25 meter. Seiring berjalannya waktu, gelombang tsunami merambat ke arah daratan, menghasilkan gelombang dengan tinggi sebesar 20 meter pada Gambar  4c, begitu pula dalam waktu 40 menit gelombang mencapai daratan dengan ketinggian sekitar 20 meter pada Gambar 4d. Waktu dan tinggi gelombang di setiap titik pengamatan dianalisis berdasarkan elevasi muka air maksimum keseluruhan  yang  diperoleh  dalam  simulasi.  Hasil  analisis  tinggi  gelombang  dan  waktu  tiba  di  Pantai Bandealit dapat dilihat pada Tabel 3.
Tabel 3 . Ketinggian dan waktu tiba tsunami Skenario 1 Skenario 2 Lokasi Tinggi Tsunami (m) Waktu (menit) Tinggi Tsunami (m) Waktu (menit) Bandealit 1 8.68 35 15.39 20 Bandealit 2 9.14 32 16.35 20 Bandealit 3 10.38 32 17.29 20
Berdasarkan Tabel 3 pada skenario 1 didapat ketinggtian maksimum pada titik pengamatan Bandealit 1, Bandealit 2, dan Bandealit 3 masing – masing adalah adalah 8.68 meter, 9.14 meter, dan 10.38 meter dengan waktu tiba 32 – 35 menit. Sedangkan untuk skenario 2, waktu tiba gelombang sama yaitu 20 menit dengan  ketinggian  15.39  meter,  16.35  meter,  dan  17.29  meter.  Mengacu  pada  jurnal  Wisyanto  (2003),
a)
b) c) d)
ketinggian tsunami di Pantai Bandealit adalah 10 meter untuk kasus Tsunami Banyuwangi 1994, maka dari ketiga  titik  pengamatan  tersebut  ketinggian  tsunami  yang  sesuai  yaitu  pada  lokasi  Bandealit  3. Perbandingan ketinggian dan waktu tiba tsunami untuk Skenario 1 dan 2 ditunjukkan pada Gambar 5.
Gambar 5. Hasil Simulasi Skenario 1 dan Skenario 2
Berdasarkan grafik yang ditunjukan pada Gambar 5 pada skenario 1 dan skenario 2, waktu tempuh gelombang tsunami untuk mencapai  tepi pantai lebih cepat  skenario 2 dibandingkan dengan skenario 1. Hasil pemodelan tsunami pada 3 titik pengamatan di Pantai Bandealit pada skenario 2 menunjukkan bahwa gelombang  tsunami  memiliki  waktu  tempuh  yang  sama  saat  mencapai  ketinggian  maksimum  yaitu  20 menit. Sedangkan pada skenario 1 waktu tempuh pada keempat titik pengamatan yaitu 32 – 35 menit.
Faktor yang mempengaruhi perbedaan tinggi gelombang dan waktu tempuh gelombang tsunami untuk mencapai pesisir adalah karakteristik dari masing – masing gempa dan juga jarak lokasi sumber patahan ke titik observasi. Pada skenario 2, magnitudo yang digunakan sangat besar dibandingkan dengan skenario 1, lebar patahan dan dislokasi yang digunakan pada skenario 2 juga jauh lebih besar. Tentu saja hal tersebut dapat  menimbulkan  tsunami  yang  sangat  dahsyat,  sehingga  waktu  yang  dibutuhkan  gelombang  untuk mencapai pesisir cukup singkat.

## 5. KESIMPULAN DAN SARAN
Dari hasil simulasi menggunakan 3 titik pengamatan didapat tinggi gelombang tsunami untuk Pantai Bandealit  masing  –  masing  sebesar  8.7  m  dan  9.1  m,  dan  10.4  m  berdasarkan  skenario  Tsunami Banyuwangi 1994, lalu untuk skrenario PuSGeN Mw 8,7 didapat ketinggian gelombang tsunami masing – masing sebesar 15.4 m, 16.4 m, dan 17.3 m. Sedangkan waktu yang dibutuhkan gelombang tsunami untuk mencapai  Pantai  Bandealit  berdasarkan  skenario  Tsunami  Banyuwangi  1994  yaitu  32  -  35  menit. Sedangkan berdasarkan skrenario PuSGeN Mw 8,7 waktu tempuh gelombang tsunami yaitu 20 menit.
Pada penelitian ini digunakan grid dengan ukuran 500 x 250 m untuk wilayah dekat pantai dan 500 x 500  m  untuk  wilayah  lepas  Pantai.  Pada  penelitian  selanjutnya  diharapkan  menggunakan  grid  dengan
ukuran yang lebih kecil serta menggunakan data batimetri hanya dari BATNAS karena memiliki resolusi yang  lebih  rapat,  sehingga  didapakan  hasil  pemodelan  yang  lebih  akurat  dan  melakukan  pengecekan koordinat serta elevasi data batimetri terlebih dahulu pada lokasi penelitian.

## DAFTAR PUSTAKA
Armono, H. D., Putra, A. R., & Wahyudi. 2021. Analysis of Tsunami Wave Height, Run-up, and Inundation in  the  Coastal  of  Blitar  and  Malang  Regency. IOP  Conference  Series:  Earth  and  Environmental Science . 936(1). https://doi.org/10.1088/1755-1315/936/1/012013
Badan Pusat Statistik. 2023. Kecamatan Tempurejo Dalam Angka Tahun 2023. Statistik . 1102001. 85. Deltares. 2014. 3D/2D modelling suite for integral water solutions: Hydro-Morphodynamics . User Manual . 710.
Hartuti,
E. R. 2009.
Buku Pintar Gempa. Yogyakarta:
Diva Press . 89. https://doi.org/10.1073/pnas.0703993104
Kartini,  L.,  Utami  Agung  Wiyono,  R.,  Denisio  Edwin  Rikarda,  R.,  Bagus  Kurniawan,  F.,  Halik,  G., Hidayah, E., & Bintang Pratama, M. 2020. Two-dimensional numerical simulations of tsunamis on Puger  beach,  Jember. IOP  Conference  Series:  Earth  and  Environmental  Science .  437(1). https://doi.org/10.1088/1755-1315/437/1/012022
Kurniawan, F. B., Wiyono, R. U. A., Irawan, J. F., Halik, G., Hidayah, E., Widiarti, W. Y., & Pratama, M. B.  2020.  Numerical  simulation  of  groyne  in  Puger  Beach  Jember  District. AIP  Conference Proceedings . 2278. https://doi.org/10.1063/5.0014686
Maramai, A., & Tinti, S. 1997. The 3 June 1994 Java Tsunami: A post-event survey of the coastal effects. Natural Hazards . 15(1), 31–49. https://doi.org/10.1023/A:1007957224367
Pratama, W. A. 2017. Simulasi penjalaran gelombang tsunami akibat gempa tektonik di pantai jember . 91. http://repository.its.ac.id/46390/1/4313100075-Undergraduate_Theses.pdf
Rikarda, R. D. E., Wiyono, R. U. A., Halik, G., Hidayah, E., & Pratama, M. B. 2020. Tsunami simulation in  Puger  Beach  considering  the  combination  of  earthquake  source  in  South  Java. AIP  Conference Proceedings . 2278. https://doi.org/10.1063/5.0014684
USGS. 2022. Moment magnitude, Richter scale . U.S. Department of the Interior. https://www.usgs.gov/faqs/moment-magnitude-richter-scale-what-are-different-magnitude-scales- and-why-are-there-so-many
Wells, D. L., & Coppersmith, Kevin, J. 1994. New empical relationship between magnitude, rupture length, rupture  width,  rupture  area,  and  surface  displacement. Bulletin  of  the  Seismological  Society  of America . 84(4), 974–1002.
Wisyanto.  2003.  Analisis  Potensi  Bahaya  Tsunami  di  Pantai  Selatan  Jawa  Timur  Melalui  Identifikasi Kemiripan Bentuk Tapchan. Alami: Jurnal Teknologi Reduksi Risiko Bencana . 8(2), 33–38).