SUMMARYA new numerical method for Nwogu's (ASCE Journal of Waterway, Port, Coastal and Ocean Engineering 1993; 119:618) two-dimensional extended Boussinesq equations is presented using a linear triangular ÿnite element spatial discretization coupled with a sophisticated adaptive time integration package. The authors have previously presented a ÿnite element method for the one-dimensional form of these equations (M. Walkley and M. Berzins (International Journal for Numerical Methods in Fluids 1999; 29 (2):143)) and this paper describes the extension of these ideas to the two-dimensional equations and the application of the method to complex geometries using unstructured triangular grids. Computational results are presented for two standard test problems and a realistic harbour model. Copyright ?