A 2D numerical method with shock capturing properties was developed for Boussinesq-type equations in an attempt to describe wave-breaking phenomena. The model is based on hybridization of the 4th-order finite difference scheme and TVD method. The model tests were performed through comparison of model outputs with experimental data on two different problems; 2D solitary wave deformation over uneven bottom and wave deformations over submerged breakwaters. It was found that the model can simulate a wide range of 2D wave deformations without instability although it had a tendency to underestimate wave dissipations associated with wave breaking.
