“…Our approach is based on a Vertex-Centered Finite Volume Method (hereinafter referred to as VCFVM) and sets all unknowns on grid vertices. It should be noted that the vertex-centered finite volume method has been widely studied for solving anisotropic diffusion/parabolic equations on general polygonal/polyhedral meshes in the past decade, see [41,40,45,36,5,29]. The numerical scheme proposed in this paper is partially inherited from [29], but substantial modifications in the spatial and temporal discretization have been made to consider the properties of Richards' equation, including the advection flux involved and the high nonlinearity between the soil hydraulic parameters and the soil water potential.…”