Abstract. Conventional models of pumping tests in unconfined aquifers often neglect the unsaturated flow process. This study concerns the coupled unsaturatedsaturated flow process induced by vertical, horizontal, and slant wells positioned in an unconfined aquifer. A mathematical model is established with special consideration of the coupled unsaturatedsaturated flow process and the well orientation. Groundwater flow in the saturated zone is described by a three-dimensional governing equation and a linearized three-dimensional Richards' equation in the unsaturated zone. A solution in the Laplace domain is derived by the Laplacefinite-Fourier-transform and the method of separation of variables, and the semi-analytical solutions are obtained using a numerical inverse Laplace method. The solution is verified by a finite-element numerical model. It is found that the effects of the unsaturated zone on the drawdown of a pumping test exist at any angle of inclination of the pumping well, and this impact is more significant in the case of a horizontal well. The effects of the unsaturated zone on the drawdown are independent of the length of the horizontal well screen. The vertical well leads to the largest water volume drained from the unsaturated zone(W) during the early pumping time, and the effects of the well orientation on Wvalues become insignificant at the later time. The screen length of the horizontal well does not affectW for the whole pumping period. The proposed solutions are useful for the parameter identification of pumping tests with a general well orientation (vertical, horizontal, and slant) in unconfined aquifers affected from above by the unsaturated flow process.