A new numerical tool is presented which models the two-dimensional contaminant transport through saturated porous media using a meshfree method called the radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes is needed. An advection-dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The Galerkin weak form of the governing equation is formulated using two-dimensional meshfree shape functions constructed using thin plate spline radial basis functions. A computer program is developed for the implementation of the RPIM procedure. Three numerical examples are presented and the results are compared with those obtained from the analytical solution and finite element method. The experimental results are also used to validate the approach. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. © Springer-Verlag 2008.