Including hillslope processes in models of fluvial landform evolution is still challenging. Since applying the respective models for fluvial and hillslope processes to the entire domain causes scaling problems and makes the results dependent on the spatial resolution, the domain is explicitly subdivided into channels and hillslopes in some models. The transition from hillslopes to channels is typically attributed to a given threshold catchment size as a proxy for a minimum required discharge. Here we propose a complementary approach for delineating channels based on the discrete representation of the topography. We assume that sites with only one lower neighbor are channelized. In combination with a suitable model for hillslope processes, this concept initiates the self-organization of channels and hillslopes. A numerical analysis with a simple model for hillslope dynamics reveals no scaling issues, so the results appear to be independent of the spatial resolution. The approach predicts a break in slope in the sense that all channels are distinctly less steep than hillslopes. On a regular lattice, the simple D8 flow-routing scheme (steepest descent among the eight nearest and diagonal neighbors) harmonizes well with the concept proposed here. The D8 scheme works well even when applied to the hillslopes. This property simplifies the numerical implementation and increases its efficiency.</p
