Rain infiltration into a soil slope leads to propagation of the wetting front, transport of air in pores and deformation of the soils, in which coupled processes among the solid, liquid and gas phases are typically involved. Most previous studies on the unsaturated flow and its influence on slope stability were based on the singlephase water flow model (i.e., the Richards Equation) or the waterair two-phase flow model. The effects of gas transport and soil deformation on the movement of groundwater and the evolution of slope stability were less examined with a coupled solid-water-air model. In this paper, a numerical model was established based on the principles of the continuum mechanics and the averaging approach of the mixture theory and implemented in an FEM code for analysis of the coupled deformation, water flow and gas transport in porous media. The proposed model and the computer code were validated by the Liakopoulos drainage test over a sand column, and the significant effect of the lateral air boundary condition on the draining process of water was discussed. On this basis, the coupled processes of groundwater flow, gas transport and soil deformation in a homogeneous soil slope under a long heavy rainfall were simulated with the proposed three-phase model, and the numerical results revealed the remarkable delaying effects of gas transport and soil deformation on the propagation of the wetting front and the evolution of the slope stability. The results may provide a helpful reference for hazard assessment and control of rainfall-induced landslides. Rain infiltration into a soil slope leads to propagation of the wetting front, transport of air in pores and deformation of the soils, in which coupled processes among the solid, liquid and gas phases are typically involved. Most previous studies on the unsaturated flow and its influence on slope stability were based on the singlephase water flow model (ie, the Richards Equation) or the waterair two-phase flow model. The effects of gas transport and soil deformation on the movement of groundwater and the evolution of slope stability in this paper, a numerical model was established based on the principles of the continuum mechanics and the averaging approach of the mixture theory and implemented in FEM code for analysis of the coupled deformation, water flow and gas transport in porous media. The proposed model and the computer code were validated by the Liakopoulos drainage test over a sand column, and the significant effect On the basis of the lateral air boundary condition on the draining process of water was discussed. On this basis, the coupled processes of water was discussed. , and the numerical results revealed the remarkable delaying effects of gas transport and soil deformation on the propagation of the wetting front and the evolution of the slope stability.