Self-similarity and vanishing diffusion in fluvial landscapes.

Complex topographies exhibit universal properties when fluvial erosion dominates landscape evolution over other geomorphological processes. Similarly, we show that the solutions of a minimalist landscape evolution model display invariant behavior as the impact of soil diffusion diminishes compared to fluvial erosion at the landscape scale, yielding complete self-similarity with respect to a dimensionless channelization index. Approaching its zero limit, soil diffusion becomes confined to a region of vanishing area and large concavity or convexity, corresponding to the locus of the ridge and valley network. We demonstrate these results using one dimensional analytical solutions and two dimensional numerical simulations, supported by real-world topographic observations. Our findings on the landscape self-similarity and the localized diffusion resemble the self-similarity of turbulent flows and the role of viscous dissipation. Topographic singularities in the vanishing diffusion limit are suggestive of shock waves and singularities observed in nonlinear complex systems.

A minimalist model for coevolving supply and drainage networks.

Numerous complex systems, both natural and artificial, are characterized by the presence of intertwined supply and/or drainage networks. Here, we present a minimalist model of such coevolving networks in a spatially continuous domain, where the obtained networks can be interpreted as a part of either the counter-flowing drainage or co-flowing supply and drainage mechanisms. The model consists of three coupled, nonlinear partial differential equations that describe spatial density patterns of input and output materials by modifying a mediating scalar field, on which supply and drainage networks are carved. In the two-dimensional case, the scalar field can be viewed as the elevation of a hypothetical landscape, of which supply and drainage networks are ridges and valleys, respectively. In the three-dimensional case, the scalar field serves the role of a chemical signal, according to which vascularization of the supply and drainage networks occurs above a critical 'erosion' strength. The steady-state solutions are presented as a function of non-dimensional channelization indices for both materials. The spatial patterns of the emerging networks are classified within the branched and congested extreme regimes, within which the resulting networks are characterized based on the absolute as well as the relative values of two non-dimensional indices.