This paper proposes a local meshless radial basis function (RBF) method to obtain the solution of the two-dimensional time-fractional Sobolev equation. The model is formulated with the Caputo fractional derivative. The method uses the RBF to approximate the spatial operator, and a finite-difference algorithm as the time-stepping approach for the solution in time. The stability of the technique is examined by using the matrix method. Finally, two numerical examples are given to verify the numerical performance and efficiency of the method.