The hierarchy of channel networks in landscapes displays features that are characteristic of nonequilibrium complex systems. Here we show that a sequence of increasingly complex ridge and valley networks is produced by a system of partial differential equations coupling landscape evolution dynamics with a specific catchment area equation. By means of a linear stability analysis we identify the critical conditions triggering channel formation and the emergence of characteristic valley spacing. The ensuing channelization cascade, described by a dimensionless number accounting for diffusive soil creep, runoff erosion, and tectonic uplift, is reminiscent of the subsequent instabilities in fluid turbulence, while the structure of the simulated patterns is indicative of a tendency to evolve toward optimal configurations, with anomalies similar to dislocation defects observed in pattern-forming systems. The choice of specific geomorphic transport laws and boundary conditions strongly influences the channelization cascade, underlying the nonlocal and nonlinear character of its dynamics.
Keywords: detachment limited; drainage area; landscape evolution model; ridge and valley patterns; river networks.