This paper proposes two efficient numerical schemes to deal with time-fractional integro-differential reaction-convection-diffusion equations of the Volterra type in a rectangular domain. The fractional derivative is taken in the Caputo sense of order \( \alpha \in (0,1).\) In order to construct the computational semi-discrete schemes, we discretize the fractional operator using finite difference techniques, namely the L1 and the L1-2, on a uniform mesh in the time direction. A composite trapezoidal rule approximates the integral part. The semi-discrete problems become a set of two-point boundary value problems (BVPs) in the spatial domain. The cubic B-spline collocation method is employed to solve the corresponding BVPs. The optimal error estimation and convergence analysis are established under suitable regularity assumptions on the initial data and the true solution of the considered model problem. It is shown in the theoretical establishment that the L1 based scheme achieves \( (2-\alpha ) \) order of accuracy, while the rate of convergence of the L1-2 based scheme is quadratic. Finally, a numerical experiment is performed in support of the theoretical findings.
