The variational iteration method (VIM) is applied to find the approximate solutions of the Whitham-Broer-Kaup (WKB) shallow water model. The results obtained are compared with the exact solutions. The results show that the variational iteration method is efficient and practically well suited for solving coupled nonlinear water wave equations.