in this work, we introduce an efficient meshless technique for solving the two-dimensional variable-order time-fractional mobile/immobile advection-diffusion model with Dirichlet boundary conditions. The main advantage of this scheme is to obtain a global approximation for this problem which reduces such problems to a system of algebraic equations. To approximate the first and fractional variable-order against the time, we use the finite difference relations. The proposed method is based on the Moving Kriging (MK) interpolation shape functions. To discretization this model in space variables, we use the MK interpolation. Duo to the fact that the shape functions of MK have Kronecker’s delta property, boundary conditions are imposed directly and easily. To illustrate the capability of the proposed technique on regular and irregular domains, several examples are presented in different kinds of domains and with uniform and nonuniform nodes. Also, we use this scheme to simulating anomalous contaminant diffusion in underground reservoirs.