This paper presents a new numerical tool to model the two-dimensional contaminant transport through saturated porous media using a meshfree method, called 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. The 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 2D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Three numerical examples are presented and the results are compared with those obtained from the finite element method and analytical solutions. In order to test the practical applicability and performance of the RPIM, two case studies of contaminant transport through landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data. © 2009 Springer-Verlag.