Robust　and　reliable　numerical　models　are　vital　to　solve　the　Richards＇　equation,　which　depicts　the　variably　saturated　flows　in　porous　media.　In　this　study,　the　Richards＇　equation　is　discretized　spatially　with　the　numerical　manifold　method　(NMM)　and　temporally　with　the　backward　Euler　scheme,　in　which　the　under-relaxation　and　mass　lumping　techniques　are　introduced　to　keep　the　numerical　stability　and　mass　balance.　Several　examples　are　performed　to　validate　the　correctness　and　accuracy　of　the　proposed　model.　The　numerical　results　demonstrate　the　potential　applicability　of　the　proposed　model　to　solve　saturated-unsaturated　seepage　problems.