Mathematical and Numerical Analysis of Tsunami Generation and Propagation: The Initial Tsunami Condition
Abstract
This work investigates the appropriate initial condition for tsunami propagation that allows replacing a kinematic system with its corresponding static formulation, where the bottom is fixed at its final earthquake-induced deformation. Using the two-dimensional (2D) nonlinear shallow-water equations (NSW) and a 2D dispersive abcd-Boussinesq system, we prove that for the NSW system, the appropriate initial free-surface condition is exactly the final bottom deformation. For the Boussinesq model, this condition is obtained by first applying a \textit{ Boussinesq filtering operator} to the bottom profile. Furthermore, for both models, we establish an error estimate between the kinematic and static solutions that depends linearly on the earthquake rise time. This result provides a rigorous mathematical explanation for the errors observed in previous numerical tsunami simulations. As a direct consequence, the error vanishes as the bottom motion becomes instantaneous. We validate the theoretical results using an Okada solution for the seabed deformation. In particular, both the NSW model and the passive approach, which initializes the free surface with the unfiltered bottom deformation, overestimate the amplitudes of the tsunami waves relative to the linear potential theory for the configurations considered.
Related
- A numerical method for two-dimensional Boussinesq-Full Dispersion systems modelling internal wave propagation
- Direct Vertex-Derivative Coefficient Propagation for Homogeneous Numerical Integration (HNI)
- Numerical and Machine Learning Prediction of Blast and Fragment Propagation in RC Frames Under Internal Explosion
- A Mathematical Study of an SIS Epidemic Model: Global Asymptotic Stability Analysis and Construction of Nonstandard Numerical Schemes
- Convergence analysis of a numerical scheme for the Cahn-Hilliard-Navier-Stokes system with dynamical boundary condition and its application to moving contact line problem
BibTeX
@misc{lamsahel2026mathematical,
title = {{Mathematical and Numerical Analysis of Tsunami Generation and Propagation: The Initial Tsunami Condition}},
author = {Noureddine Lamsahel and Min Chen},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.02708},
eprint = {2611.02708},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.02708},
}
RIS
TY - UNPB
TI - Mathematical and Numerical Analysis of Tsunami Generation and Propagation: The Initial Tsunami Condition
AU - Noureddine Lamsahel
AU - Min Chen
PY - 2026
DA - 2026/11/02
PB - arXiv
DO - 10.48550/arxiv.2611.02708
UR - https://arxiv.org/abs/2611.02708
AB - This work investigates the appropriate initial condition for tsunami propagation that allows replacing a kinematic system with its corresponding static formulation, where the bottom is fixed at its final earthquake-induced deformation. Using the two-dimensional (2D) nonlinear shallow-water equations (NSW) and a 2D dispersive abcd-Boussinesq system, we prove that for the NSW system, the appropriate initial free-surface condition is exactly the final bottom deformation. For the Boussinesq model, this condition is obtained by first applying a \textit{ Boussinesq filtering operator} to the bottom profile. Furthermore, for both models, we establish an error estimate between the kinematic and static solutions that depends linearly on the earthquake rise time. This result provides a rigorous mathematical explanation for the errors observed in previous numerical tsunami simulations. As a direct consequence, the error vanishes as the bottom motion becomes instantaneous. We validate the theoretical results using an Okada solution for the seabed deformation. In particular, both the NSW model and the passive approach, which initializes the free surface with the unfiltered bottom deformation, overestimate the amplitudes of the tsunami waves relative to the linear potential theory for the configurations considered.
ER -