The HYDRUS unsaturated flow and transport model was modified to simulate the effects of non-linear air-water interfacial (AWI) adsorption, solution surface tension-induced flow, and variable solution viscosity on the unsaturated transport of per- and polyfluoroalkyl substances (PFAS) within the vadose zone. These modifications were made and completed between March 2019 and May 2019, and were implemented into both the one-dimensional (1D) and two-dimensional (2D) versions of HYDRUS. Herein, the model modifications are described and validated against the available literature-derived PFAS transport data (i.e., 1D experimental column transport data). The results of both 1D and 2D example simulations are presented to highlight the function and utility of the model to capture the dynamic and transient nature of the temporally and spatially variable interfacial area of the AWI (Aaw) as it changes with soil moisture content (Θw) and how it affects PFAS unsaturated transport. Specifically, the simulated examples show that while AWI adsorption of PFAS can be a significant source of retention within the vadose zone, it is not always the dominant source of retention. The contribution of solid-phase sorption can be considerable in many PFAS-contaminated vadose zones. How the selection of an appropriate AawΘw function can impact PFAS transport and how both mechanisms contribute to PFAS mass flux to an underlying groundwater source is also demonstrated. Finally, the effects of soil textural heterogeneities on PFAS unsaturated transport are demonstrated in the results of both 1D and 2D example simulations.