the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
21 Sep 2021
21 Sep 2021
Morphological controls on Hortonian surface runoff: An interpretation of steady-state energy patterns, maximum power states and dissipation regimes within a thermodynamic framework
- 1Institute of Water Resources and River Basin Management, Karlsruhe Institute of Technology – KIT, Karlsruhe, Germany
- 2Institute for Hydromechanics, Karlsruhe Institute of Technology – KIT, Karlsruhe, Germany
- 3Max-Planck Institute for Biochemistry, Hans Knöll Str. 10, 07745 Jena, Germany
- 4Engler-Bunte-Institut, Water Chemistry and Water Technology – KIT, Karlsruhe, Germany
- 1Institute of Water Resources and River Basin Management, Karlsruhe Institute of Technology – KIT, Karlsruhe, Germany
- 2Institute for Hydromechanics, Karlsruhe Institute of Technology – KIT, Karlsruhe, Germany
- 3Max-Planck Institute for Biochemistry, Hans Knöll Str. 10, 07745 Jena, Germany
- 4Engler-Bunte-Institut, Water Chemistry and Water Technology – KIT, Karlsruhe, Germany
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 analysed 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 hillslope structures at the macro- and microscale control the free energy balance of Hortonian overland flow. We put special emphasis on the transitions of surface runoff processes at the hillslope scale, as hillslopes energetically behave distinctly different in comparison to fluvial systems. 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 the driver of erosion and sediment transport. We then use this framework at a macroscopic scale to analyse 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. At the microscale, we analyse flow concentration in rills and its influence on the distribution of energy and dissipation in space. Therefore, we develop a new numerical method for the Catflow model, which allows a dynamical separation of Hortonian surface runoff between a rill- and a sheet flow domain. We calibrated the new Catflow-Rill model to rainfall simulation experiments and observed overland flow in the Weiherbach catchment and found evidence that flow accumulation in rills serves as a means to redistribute energy gradients in space, therefore minimizing energy expenditure along the flow path, while also maximizing overall power of the system. Our results indicate that laminar sheet flow and turbulent rill flow on hillslopes develop to a dynamic equilibrium that corresponds to a maximum power state, and that the transition of flow from one domain into the other is marked by an energy maximum in space.
Samuel Schroers et al.
Status: open (until 18 Dec 2021)
-
RC1: 'Comment on hess-2021-479', Keith Beven, 01 Oct 2021
reply
I fully understand that this is a contribution that forms part of an on-going research program to introduce and explore the consequences of thermodynamic principles in hydrological theory. As such the attempt to do this for hillslope runoff is of interest, despite the rather restrictive assumptions under which the analysis is carried out. My main comment is that the discussion should reflect those limitations more, rather than making inferences that are far beyond the results presented.
It would be also useful (I think to this research program as a whole) if some of the concepts could be discussed more clearly. In this case in particular it is the Df component of the energy balance. This is initially the energy dissipation that is the residual of potential and kinetic energy components, and then suddenly becomes equivalent to stream power in a discussion related rivers and sediment transport (even when sediment transport is never actually considered in the analyses). I think a section that discusses how Dfshould be characterised for laminar/transitional/turbulent flows before and after a threshold for sediment mobilisation might be really useful to the reader. As it is, it seems only that a constant Manning n is used to convert to velocity, which is certainly one characterisation but surely might not be considered adequate. Perhaps this suggests a bit more than minor revision.
Further comments are included on the manuscript.
Keith Beven
-
AC1: 'Reply on RC1', Samuel Schroers, 04 Oct 2021
reply
First, we would like to thank Keith Beven for his time, effort and insightful comments. We also very much welcome the upload of a commented manuscript and will gladly incorporate his suggestions to clarify the presentation of our approach in a revised manuscript.
As we understand the principal point of criticism made by KB, we agree that the manuscript will benefit from a clearer definition of our theory and a reflection of the main simplifications. We also agree that the results and inferences made apply to landscape forms that are primarily formed by overland flow and the related forces and energy conversions. We will clarify this in the revised manuscript.
Keith Beven’s suggestion to highlight the different behavior of the Df component of the energy balance for surface runoff with and without sediment load is appreciated. We agree that there is an important difference between viscous dissipation of energy due to friction and energy transfer to the sediment. The former implies that free energy is indeed lost, while the latter means that free energy creates motion of the sediment particles. A water and sediment balance of the rainfall runoff simulation experiments by Gerlinger revealed that suspended sediments may increase the density of the fluid-sediment mix up to 10%. Unfortunately, the available measurements and data are too limited to analyse underlying dynamics and include this into the analysis. In conclusion, we admit that the description of the residual of the energy balance needs some clarification, and it is a good idea to elaborate on the differences of Df between surface runoff with and without sediment transport. This has direct implications for the inference of Maning’s n from such rainfall simulation experiments.
-
AC2: 'Reply on RC1', Samuel Schroers, 12 Oct 2021
reply
We would like to thank KB again for his comment on the Df term (dissipation) and wanted to use this discussion as an opportunity to present a possible extension of the energy scheme which is presented in the current manuscript.
Attached we outline how sediment transport affects the Df term and we distinguish four dissipation regimes, depending on whether sediment is eroded or deposited.If sediment is eroded some energy of the flow will be used for the additionally transported mass, and similarly if sediment is deposited some of the energy of the sediment-water mix becomes available.
This concept could be incorporated as a clarification of the Df term and also be part of the discussion of a revised manuscript.
-
RC2: 'Reply on AC2', Keith Beven, 12 Oct 2021
reply
Well yes ... it seems a straightforward extension in theory, since again it is not too difficult to write down balance equations. But in practice you have now split your significant unknown term Df into a number of unknown terms (even given the simplifying assumptions of steady state flows, well mixed sediment etc) that hydrologists and geomorphologists have struggled with for over a century (and we STILL resort to Manning's n or some assumed constant fraction of kinetic energy). So I think the critical thing in an revision is to really make it clear how this can help frame the problem of closing that energy balance and having some clearer understanding of how to evaluate the Df term - splitting it into component parts might provide some cursory understanding - but what is needed to allow a proper evaluation (remembering that we are actually interested in transient dynamic flows when we often cannot easily close the mass balance equations in practice)?
Just to add - I certainly do not want to discourage this research programme but by staying in the "theoretical" domain it does not yet seem to be really addressing the most significant question.
-
AC3: 'Reply on RC2', Samuel Schroers, 13 Oct 2021
reply
We thank KB for these thoughtful comments. We fully agree that the steady state surface runoff is rarely reached during rainfall runoff events. In fact, we plan to address this point in a follow up study (as explained below). We also share his concerns about the validity of empirical flow laws (in particular Manning's n) in case of very shallow overland flows with flow depths of a few mm. We already stressed that this is particularly delicate upslope of the potential energy maximum, in the current manuscript. Also, we acknowledge that the presented concept does indeed not provide closure of the energy balance. However, we would like to stress that in this study the presentation of the Df term is principally intended to provide a basis for the analysis of the spatial distribution thereof, not it’s magnitude. The question how free energy of flow is dissipated by viscous forces goes beyond the goal of this study, but we agree that this point is important, and that a revision should at least incorporate some perspective.
Yet, solving all these issues is beyond the scope of this study, and more importantly beyond the scope of the field observations at hand. Our work is clearly only a first step, but it underpins that a closer look at the spatial pattern of the steady state energy balance of surface runoff and it’s residual Df is helpful to
- identify distinct differences between hillslopes and rivers,
- single out hot spots, where the flow can perform the largest amount of physical work,
- test thermodynamic concepts (maximum power, least work, etc.) on this dynamical system .
We also like to stress that our work is not purely theoretical. We analyze overland flow measured during sprinkling experiments, which reached steady state conditions according to Scherer et al. (2012). As these experiments provide measured flow velocities in rills and back calculated values of Manning’s n. Hence, we regard these experiments, despite all their short comings (very shallow flow depth and limited validity of all available flow laws) as best suited to test our framework, at least during, albeit rare, steady state conditions. Because the experiments did not provide observed flow depth along the flow path, this approach required a numerical model, as well as expanding the CATFLOW code by a rill domain.
This analysis corroborated the existence of the potential energy maxima and their sensitivity to presence/absence of rills. And we provide first evidence that the transition zone from laminar to mixed flow coincides roughly with the location of the potential energy maximum (when rills are present). We also found that kinetic energy in the sheet and rill domain are equal at the end of the stripe. This is analogue to the maximum power configuration of an electrical circuit and at least an interesting incidence. Of course, two experiments are not sufficient to conclude whether this holds in general.
All In all, we think that our work corroborates that an energetic perspective on overland flow is helpful and has practical implications which should be further explored.
The next logical step is to give up the steady state assumption and work with transient conditions. The mentioned energy residual’s distribution might follow similar thermodynamic concepts in time as in space.
The framework is not restricted to the steady state, but requires numerical simulations of overland flow, ideally with the shallow water equations (as the diffusion wave approach assumes a steady state momentum balance).
Another key challenge is a proper design of experiments to disentangle viscous dissipation and energy transfer to and from the sediments. If this energy transfer is ignored during simulations, this might imply that dissipation is over- or underestimated (see explanation in last reply). And finally, the concept of maximum power i.e., an equal distribution of kinetic energy among rill and sheet flow during steady state, is a testable hypothesis.
-
RC3: 'Reply on AC3', Keith Beven, 13 Oct 2021
reply
An excellent reply (and I stress again that I do not want to dicourage this research programme) .... but: you do have to be careful in avoiding circularity of reasoning when coefficients (like Df or its parts) have to be back-caclulated from experimental data (in the same way that it is always possible to get a value for n given observations of velocities or discharges, which is really the only reason why it is still being used as a convenience not requiring much thought). I can see that testing any hypotheses about Df is going to be highly challenging even at small plot scales, while inferences at larger scales (including inferences about minimisation principles) are going to depend very much on assumptions about the distributed and transient nature of Df in model simulations. So I would suggest that you proceed with caution, remembering that rejecting your favoured hypothesis is often more valuable to improving understanding.
-
AC3: 'Reply on RC2', Samuel Schroers, 13 Oct 2021
reply
-
RC2: 'Reply on AC2', Keith Beven, 12 Oct 2021
reply
-
AC1: 'Reply on RC1', Samuel Schroers, 04 Oct 2021
reply
-
CC1: 'Comment on hess-2021-479', John Ding, 17 Oct 2021
reply
The authors focus the analysis on the rising limbs of overlandflow hydrographs, e.g., Figure 10, Panels (a) and (c) for the sheet and rill flows, and the sheet flow only, respectively. These are essentially the summation or S-curve hydrographs. Not mentioned are their falling limbs after the termination of the effective rainfall, i.e. recession hydrographs.
Because of the absence of the forcing effective rainfall intensity term, Peff, an additional, and simpler, analysis by the authors will likely lead to linearization of the recession hydrograph by an inverse fractional-power transformation of the discharge, 1/Q (1-1/c2). This was derived earlier for a nonlinear storage-discharge function by this writer, e.g., Ding (1974, Eq. 4; second equation after Eq. 11).
For Manning friction law, the flow depth exponent c2 = 5/3 shown in their Table 2.
References
Ding, J.Y., 1974. Variable unit hydrograph. J. Hydrol., 22: 53-69.
-
AC4: 'Reply on CC1', Samuel Schroers, 19 Oct 2021
reply
We thank JD for his comments on our manuscript. The purpose of this study is an evaluation of energy patterns during steady state runoff. While transient floods and runoff are certainly dominant in nature, design of experiments for parameter estimation often departs from a steady state assumption. The results shown in Fig. 10(c) and (d) represent therefore merely the temporal routes of hydraulic variables for reaching a steady state runoff regime and an analysis of the latter, not the routes themselves. However, as mentioned in our previous reply (AC3) to KB’s comment (RC2), we agree that a further study should incorporate temporal patterns.
This study presents therefore just a first step towards an analysis covering distribution of the energy residual in space and time, where we will of course consider JD’s comments on transient dissipation regimes.
-
AC4: 'Reply on CC1', Samuel Schroers, 19 Oct 2021
reply
Samuel Schroers et al.
Samuel Schroers et al.
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 462 | 110 | 24 | 596 | 4 | 4 |
- HTML: 462
- PDF: 110
- XML: 24
- Total: 596
- BibTeX: 4
- EndNote: 4
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1