An isogeometric analysis (IGA) is introduced to obtain a head-based solution to Richards equation for unsaturated flow in porous media. IGA uses Non-Uniform Rational B-Spline (NURBS) as shape functions, which provide a higher level of inter-element continuity in comparison with Lagrange shape functions. The semi-discrete nonlinear algebraic equations are solved using a combination of implicit backward-Euler time-integration and Newton-Raphson scheme. The time-step size is adaptively controlled based on the rate of changes in the pore pressure. The results from the proposed formulation are compared and verified against an analytical solution for one-dimensional transient unsaturated flow in a homogenous soil column. The proposed method is then applied to four more complex problems including two-dimensional unsaturated flow in a two-layered soil and a semi-circular furrow. The test cases in two-layered soil system involve sharp variations in the pressure gradient at the intersection of the two media, where the pore water pressure abruptly changes. It is shown that the proposed head-based IGA is able to properly simulate changes in pore pressure at the soils interface using fewer degrees of freedom and higher orders of approximation in comparison with the conventional finite element method.