This study presents the development of a novel three-dimensional discontinuous Galerkin coastal and estuarine modelling system named DGCEMS. The model has low spurious mixing and second-order convergence of surface water elevation, horizontal velocity, and the tracer field. It has the capability to simulate salt–freshwater interactions in the presence of wetting and drying boundaries.