Modeling tsunamis from earthquake sources near Gorringe Bank southwest of Portugal

• Published in: Journal of Geophysical Research, Washington, DC Washington, DC, 1997-12, Vol.102 (C13), p.27931-27949
• Main Authors: Gjevik, B., Pedersen, G., Dybesland, E., Harbitz, C. B., Miranda, P. M. A., Baptista, M. A., Mendes‐Victor, L., Heinrich, P., Roche, R., Guesmia, M.
• Format: Article
• Language: English
• Subjects: Earth, ocean, space; Exact sciences and technology; External geophysics; Marine; Other topics; Physics of the oceans
• ISSN: 0148-0227; 2169-9275; 2156-2202; 2169-9291
• DOI: 10.1029/97JC02179

The Azores‐Gibraltar fracture zone with the huge bathymetric reliefs in the area southwest of Portugal is believed to have been the source of large historic tsunami events. This report describes simulations of tsunami generation and propagation from sources near the Gorringe Bank. The well‐documented 1969 tsunami event is examined both with a ray‐tracing technique and with finite difference models based on various shallow water equations. Both methods show that the most likely source location is southeast of the Gorringe Bank near the epicenter location determined from seismic data. The tsunami source is calculated by formulas given by Okada [1985] for surface deformation of an elastic half‐space caused by faulting. Observed wave amplitude and travel time and values computed from an initial wave field according to Okada [1985] formulas show acceptable agreement for most stations along the coast of Portugal and Spain. However, in order to explain a large primary wave with downward displacement observed on the coast of Morocco, an alternative source model with a larger area of downward displacement has been introduced. This also leads to a better overall fit with observed travel time. Implications for disastrous events, as the one in 1755, are also discussed. Linear hydrostatic shallow water models are used for most of the simulations, but the importance of nonlinearity and dispersion is examined with the Boussinesq equations. The sensitivity of the solution to changes in the location and the strength of the source is discussed, and a series of grid refinement studies are performed in order to assess the accuracy of the simulations.