Bottom changes in coastal areas with complex shorelines

A model for the sea-bottom change simulations in coastal areas with complex shorelines is proposed. In deep and intermediate water depths, the hydrodynamic quantities are calculated by numerically integrating the contravariant Boussinesq equations, devoid of Christoffel symbols. In the surf zone, the propagation of the breaking waves is simulated by the nonlinear shallow water equations. The momentum equation is solved inside the turbulent boundary layer in order to calculate intrawave hydrodynamic quantities. An integral formulation for the contravariant suspended sediment advection-diffusion equation is proposed and used for the sea-bottom dynamic simulations. The proposed model is applied to the real case study of Pescara harbor (in Italy).