Topographic delineation is critical to watershed hydrologic modeling, which may significantly influence the accuracy of model simulations. In most traditional delineation methods, however, surface depressions are fully filled and hence, watershed-scale hydrologic modeling is based on depression-less topography. In reality, dynamic filling and spilling of depressions affect hydrologic connectivity and surface runoff processes, especially in depression-dominated areas. Thus, accounting for the internal hydrologic connectivity within a watershed is crucial to such hydrologic simulations. The objective of this study was to improve watershed delineation to further reveal such complex hydrologic connectivity. To achieve this objective, a new algorithm, HUD-DC, was developed for delineation of hydrologic units (HUs) associated with depressions and channels. Unlike the traditional delineation methods, HUD-DC considers both filled and unfilled conditions to identify depressions and their overflow thresholds, as well as all channels. Furthermore, HUs, which include puddle-based units and channel-based units, were identified based on depressions and channels and the detailed connectivity between the HUs was determined. A watershed in North Dakota was selected for testing HUD-DC, and Arc Hydro was also utilized to compare with HUD-DC in depression-oriented delineation. The results highlight the significance of depressions and the complexity of hydrologic connectivity. In addition, HUD-DC was utilized to evaluate the variations in topographic characteristics under different filling conditions, which provided helpful guidance for the identification of filling thresholds to effectively remove artifacts in digital elevation models.Water 2020, 12, 7 2 of 14 new watershed delineation methodologies are needed to capture the detailed connectivity between depressions and reveal the real hydrologic processes.Generally, a typical depression consists of a ponding area and at least one threshold (i.e., pour point). Depressions in an area exhibit a hierarchical characteristic. Following a filling process, the depressions that share the same threshold can merge into a larger, higher-level depression [15,16]. Such a merging process continues until all highest-level depressions are generated. In the existing delineation methods, some algorithms focus on delineation of highest-level depressions only, while others center on characterizing all level depressions. For example, Temme et al. [21] developed an algorithm to search DEM cells for identification of the highest-level depressions and the algorithm was further incorporated into a geomorphological model, LAPSUS (landscape process modelling at multi-dimensions and scales), to investigate the dynamic landscape evolution. Arc Hydro is an ArcGIS-based tool that handles data pre-processing for hydrologic modeling [22]. Differently from the method by Temme et al. [21], in Arc Hydro all highest-level depressions are identified based on the differences between the fully filled DEM and the original ...