Drainage networks delineated from Digital Elevation Models (DEMs), are the basis for the modelling of geomorphological and hydrological processes, biogeochemical cycling, and water resources management. Besides providing effective models of water flows, automatically extracted drainage networks based on topography can diverge from reality to varying degrees. The variability of such disagreement within catchments has rarely been examined as a function of the heterogeneity of land cover, soil type, and slope in the catchment of interest. This research gap might not only substantially limit our knowledge of the uncertainty of hydrological prediction, but can also cause problems for users attempting to use the data at a local scale. Using 1:100000 scale land cover maps, Quaternary deposits maps, and 2 m resolution DEMs, it is found that the accuracy of delineated drainage networks tends to be lower in areas with denser vegetation, lower hydraulic conductivity, and higher erodibility. The findings of this study could serve as a guide for the more thoughtful usage of delineated drainage networks in environmental planning, and in the uncertainty analysis of hydrological and biochemical predictions. Therefore, this study makes a first attempt at filling the knowledge gap described above.