Hortonian Overland Flow, Hillslope Morphology and Stream PowerI: Spatial Energy Distributions and Steady-state Power Maxima

Abstract. Recent developments in hydrology have led to a new perspective on runoff processes, extending beyond the classical mass dynamics of water in a catchment. For instance, stream flow has been analyzed in a thermodynamic framework, which allows the incorporation of two additional physical laws and enhances our understanding of catchments as open environmental systems. Related investigations suggested that energetic extremal principles might constrain hydrological processes, because the latter are associated with conversions and dissipation of free energy. Here we expand this thermodynamic perspective by exploring how macro and micro hillslope structures control the free energy balance of Hortonian overland flow. This may ultimately help understanding why these structures have evolved to their present shape. To this end, we develop a general theory of surface runoff and of the related conversion of geopotential energy gradients into other forms of energy, particularly kinetic energy as driver of erosion and sediment transport. We then use this framework to analyze how combinations of typical hillslopes profiles and width distributions control the spatial patterns of steady state stream power and energy dissipation along the flow path. Additionally, we provide a first order estimate whether and when rills reduce the overall energy dissipation compared to sheet flow. Finally, we relate accumulated stream power of linear hillslopes to slope angles, closing the loop to Horton's original formulation of erosion force. The analytical analysis of stream power reveals that the common formulation, a function of the depth-discharge product is a reduced version of the more general equations if we neglect changes in velocity and discharge in space. The full equations of stream power result in maximum energy fluxes in space for sinusoidal and exponential hillslope profiles, while linear and negative exponential forms unlimitedly increase these fluxes in the downstream direction. Depending on geometry, rill flow increases or decreases kinetic energy fluxes downslope, effectively counteracting or increasing the dissipation of potential energy. For accumulated power in space for steady state runoff, we find that on linear hillslopes a slope angle of 45° maximizes the conversion of potential energy into dissipation and an angle of 35° maximizes the conversion of potential energy into kinetic energy.