Abstract: A hybrid finite-difference/finite-volume scheme is developed to solve the 2D fully nonlinear Boussinesq equations. Finite volume method, in conjunction with the MUSTA scheme, is used to evaluate the flux terms while finite difference method is used to approximate the rest terms in the conservative governing equations. The third order Runge-Kutta method with TVD property is adopted for time marching. The proposed model has the advantages of shock capturing, easy coding and strong stability preserving, and contains less adjustable parameters. Two typical numerical tests are conducted for model validation, and the calculated results are in good agreement with experimental data.