This paper focuses on the application of an efficient technique, namely, the fractional natural decomposition method (FNDM). The numerical solutions of the model containing the water transport in unsaturated porous media, called Richards equation, are extracted. This model is used to describe the non-locality behaviors which cannot be modeled under the framework of classical calculus. To demonstrate the effectiveness and efficiency of the scheme used, two cases with time-fractional problems are considered in detail. The numerical stimulation is presented with results accessible in the literature, and corresponding consequences are captured with different values of parameters of fractional order. The attained consequences confirm that the projected algorithm is easy to implement and very effective to examine the behavior of nonlinear models. The reliable algorithm applied in this paper can be used to generate easily computable solutions for the considered problems in the form of rapidly convergent series.