A model for tsunamis propagation is derived by coupling a weakly compressible liquid layer, standing for the ocean, to a viscoelastic solid layer of constant depth, representing the upper layer of the Earth. Below this layer, the Earth is assumed to be perfectly rigid. The resulting equations are averaged over the depth of the water layer for the liquid part and are integrated over the thickness of the solid layer for the solid part. The obtained system of equations is hyperbolic and admits an exact equation of energy conservation. The model is dispersive and includes an elastic branch, an acoustic branch and a gravity branch. The elasticity of the solid layer entails: 1) a reverse dispersion effect i.e. the phase velocity of the gravity branch decreases for small wave numbers if the wave number decreases; 2) an arrival time delay of the tsunami and 3) a leading negative phase i.e. a negative water elevation before the arrival of the main wave. All these effects are absent if the sea floor is assumed to be perfectly rigid. The arrival time delay due to elasticity comes in addition to the delay due to seawater compressibility. These effects are in agreement with previous works on tsunamis propagation. This preliminary model is a promising alternative approach for tsunamis propagation due to its ability to capture the main physical effects with a relatively simple mathematical structure of equations, which is easy to solve numerically.