Prediction of the movement of water and solutes in the vadose zone requires information on the distribution of spatial trends and heterogeneities in porous media. The present study describes different lithofacies origination mainly from glaciofluvial deposits. Among different lithofacies, hydrological relationships were investigated. By means of a two-dimensional hydrological model it was evaluated how the flow of water and leaching of metribuzin (4-amino-6-tert-butyl-4,5-dihydro-3-methylthio-1,2,4-triazin-5-one) was affected. Two selected large outcrop sections consisting of glacial outwash deposits were used in the modeling study. Eleven different lithofacies were distinguished and described in terms of texture distribution, sorting, bedding style, and external boundaries based on excavated soil profiles from 27 locations representing seven predominantly sandy landforms in Denmark. Undisturbed soil columns were sampled from each of the lithofacies and brought to the laboratory to be analyzed. With respect to their soil hydraulic properties, the different lithofacies formed four different hydrofacies having relatively homogeneous, hydrogeological properties. Two large outcrop sections from one of the locations (a gravel pit) located near the terminal moraine of the former Weichsel glacier were used for the HYDRUS-2D modeling. Modeling results revealed that the spatial distribution of sedimentary bodies affected water flow and the leaching of metribuzin.