We solve one-dimensional Kirchhof transformed Richards equation numerically using finite difference method with various time-stepping schemes, forward in time central in space (FTCS), backward in time central in space (BTCS), Crank–Nicolson (CN), and a stabilized Runge–Kutta–Legendre super time-stepping (RKL), and compare their performances.