In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.