The propagation and associated run-up process of nearshore tsunamis in the vicinity of shorelines have been analyzed by using a two-dimensional numerical model. The governing equations of the model are the nonlinear shallow-water equations. They are discretized explicitly by using a finite volume method and the numerical fluxes are reconstructed with a HLLC approximate Riemann solver and weighted averaged flux method. The model is applied to two problems; The first problem deals with water surface oscillations, while the second one simulates the propagation and subsequent run-up process of nearshore tsunamis. Predicted results have been compared to available analytical solutions and laboratory measurements. A very good agreement has been observed.