In this study, a parameter-uniform numerical scheme is built and analyzed to treat a singularly perturbed parabolic differential equation involving large spatial delay. The solution of the considered problem has two strong boundary layers due to the effect of the perturbation parameter, and the large delay causes a strong interior layer. The behavior of the layers makes it difficult to solve such problem analytically. To treat the problem, we developed a numerical scheme using the weighted average ($$\theta$$
θ
-method) difference approximation on a uniform time mesh and the central difference method on a piece-wise uniform spatial mesh. We established the Stability and convergence analysis for the proposed scheme and obtained that the method is uniformly convergent of order two in the temporal direction and almost second order in the spatial direction. To validate the applicability of the proposed numerical scheme, two model examples are treated and confirmed with the theoretical findings.