“…In numerical modeling, a single-phase flow model based on the Richards equation was utilized widely to simulate the unsaturated flow in porous media [5,15]. The change in pore water pressure might be affected by the interstitial air, i.e., the Lisse effect, which is why the two-phase flow model was developed to simulate the dynamic behavior of pore air and its impacts on the pore water pressure [6,11,12,18,24,29]. To investigate the impacts of changes in pore water/air pressure on the slope stability, the limit equilibrium method evolved in the single-/two-phase flow model [18,30,31].…”