Infiltration and drainage processes in multi-layered soils are complicated by contrasting hydraulic properties. The objective of this study was to evaluate the performances of the hysteretic and non- hysteretic models to simulate the infiltration and drainage processes from three different natural soil profiles containing as manyas20texturallydifferentlayers.Hydraulicpropertieswereestimatedfromsoiltexturesusingpedotransferfunctions and were calibrated and validated using measured water contents during infiltration and drainage phases, respectively. The results supported the use of the Arya-Paris pedotransfer function to estimate the wetting curve when contact angles are incorporated. The unique Kozeny-Carmen equation parameter was evaluated by optimizing the estimated saturated hydraulic conductivity. The calibrated numerical model (Hydrus-1D) accurately simulated soil water content profiles and water volumes during the infiltration and drainage phases. The mean error of prediction (MEP) between the measured and estimated soil water contents varied from 0.030 to 0.010 cm3 cm 3, and the standard deviation of prediction (SDP) from 0.003 to 0.057 cm3 cm 3. The simulation was improved for more heterogeneous soil profiles when hysteresis was taken into account. The measured and simulated results indicated that the soil profile with vertical heterogeneity in soil texture can store more water than the similar textured vertically homogeneous soils under drained conditions.