This study attempts to develop a numerical scheme for 2-D transient analysis under unsaturated conditions. First, the unsaturated groundwater flow was described using the mass conservation law. Then, the Finite Difference Method and Backward Euler approximation were applied for space and time discretization, respectively. |