Hydrologic models such as the USEPA Stormwater Management Model (SWMM) are commonly used to assess the design and performance of green infrastructure (GI). To accurately represent GI performance models used in design need to be able to address both the hydrology/hydraulics of the catchment and the GI unsaturated (vadose) zone hydrology. While hydrologic models, such as SWMM, address the need for catchment hydrology/hydraulics, they often simplify the unsaturated zone hydrology. This paper presents a methodology utilizing existing components of SWMM to represent unsaturated zone hydrology in an accessible format that does not require adjustments to the SWMM source code. The methodology simulated the unsaturated soil water movement by considering flow caused by differences of soil matric head and flow caused by gravity between soil layers with finite depth/length. The flow flux related to the soil matric head is a function of soil water diffusivity (D) and the soil moisture gradient, where D can be represented by a pump curve in SWMM. The flow flux related to gravity was controlled by unsaturated hydraulic conductivity (K) only and was also simulated by a pump. The methodology was compared to another variably saturated model, HYDRUS, with theoretical soils (with single layers of sand, loam, silt, and clay, as well as dual-layer scenarios). Field data was used to compare the methodology to HYDRUS and the SWMM LID (Low Impact Development) module. In all comparisons the presented methodology and HYDRUS delivered similar results for the vadose zone response to a storm event, while the LID module of SWMM exhibited slower water movement. The results showed that under natural conditions, the approximation of the presented methodology yielded satisfactory results to simulate flow through the unsaturated vadose zone.