Time discretization along with space discretization is important in the numerical simulation of subsurface flow applications for long run. In this paper, we derive theoretical convergence error estimates in discrete-time ...