FUNWAVE Modeling Parameters for Marmara Sea Tsunami

[M.fidan Taşkıran]

Dear FUNWAVE Team,

As part of our undergraduate thesis, we are conducting a regional tsunami simulation for the Sea of Marmara using FUNWAVE-TVD. We processed GEBCO bathymetry and generated the required depth and initial η files. However, we would like to confirm the correctness of several modeling parameters and would greatly appreciate your technical guidance.

Our main model inputs are:
• Bathymetry: GEBCO_2025 (region: 40.0–41.6°N, 26.0–30.0°E)
• Grid size: Nx = 960, Ny = 384
• Resolution: DX = DY = 0.0008333° (≈ 30 m)
• Coordinate system: LAT/LON
• Depth file: GEBCO-based, land masked as 0 m
• Initial sea-surface displacement (η):
• Segment: Adalar Fault
• Center: (40.88°N, 29.12°E)
• Rupture length: ~20 km
• Assumed width: ~10 km
• Slip: 5 m reverse fault
• Method: Simplified elliptical Gaussian uplift distribution

We would like to request your feedback on the following points:
1. Is this resolution (≈30 m) and grid size appropriate for a semi-enclosed basin like the Marmara Sea?
2. Is our use of LAT/LON coordinates with DX/DY specified in degrees correct for FUNWAVE?
3. For fault-generated tsunami initialization, is there a recommended method (e.g., Okada-based deformation) or reference implementation you advise?
4. Do you have example cases, recommended parameter sets, or validation studies for regional tsunami scenarios?

Any guidance or documentation you can provide would be extremely valuable for improving the scientific reliability of our work.

Thank you for your support.
Muhammed Fidan Taşkıran