Joint Editorial Invigorating Hydrological Research through Journal Publications

Editors of several journals in the field of hydrology met during the General Assembly of the European Geosciences Union-EGU in Vienna in April 2017. This event was a follow-up of similar meetings held in 2013 and 2015. These meetings enable the group of editors to review the current status of the journals and the publication process, and to share thoughts on future strategies. Journals were represented at the 2017 meeting by their editors, as shown in the list of authors. The main points on invigorating hydrological research through journal publications are communicated in this joint editorial published in the above journals.Water Resource

HESS Opinions: Linking Darcy's equation to the linear reservoir

In groundwater hydrology, two simple linear equations exist describing the relation between groundwater flow and the gradient driving it: Darcy's equation and the linear reservoir. Both equations are empirical and straightforward, but work at different scales: Darcy's equation at the laboratory scale and the linear reservoir at the watershed scale. Although at first sight they appear similar, it is not trivial to upscale Darcy's equation to the watershed scale without detailed knowledge of the structure or shape of the underlying aquifers. This paper shows that these two equations, combined by the water balance, are indeed identical provided there is equal resistance in space for water entering the subsurface network. This implies that groundwater systems make use of an efficient drainage network, a mostly invisible pattern that has evolved over geological timescales. This drainage network provides equally distributed resistance for water to access the system, connecting the active groundwater body to the stream, much like a leaf is organized to provide all stomata access to moisture at equal resistance. As a result, the timescale of the linear reservoir appears to be inversely proportional to Darcy's <q>conductance</q>, the proportionality being the product of the porosity and the resistance to entering the drainage network. The main question remaining is which physical law lies behind pattern formation in groundwater systems, evolving in a way that resistance to drainage is constant in space. But that is a fundamental question that is equally relevant for understanding the hydraulic properties of leaf veins in plants or of blood veins in animals.Water Resource

HESS Opinions "Catchments as meta-organisms – a new blueprint for hydrological modelling"

Catchment-scale hydrological models frequently miss essential characteristics of what determines the functioning of catchments. The most important active agent in catchments is the ecosystem. It manipulates and partitions moisture in a way that supports the essential functions of survival and productivity: infiltration of water, retention of moisture, mobilization and retention of nutrients, and drainage. Ecosystems do this in the most efficient way, establishing a continuous, ever-evolving feedback loop with the landscape and climatic drivers. In brief, hydrological systems are alive and have a strong capacity to adjust themselves to prevailing and changing environmental conditions. Although most models take Newtonian theory at heart, as best they can, what they generally miss is Darwinian theory on how an ecosystem evolves and adjusts its environment to maintain crucial hydrological functions. In addition, catchments, such as many other natural systems, do not only evolve over time, but develop features of spatial organization, including surface or sub-surface drainage patterns, as a by-product of this evolution. Models that fail to account for patterns and the associated feedbacks miss a critical element of how systems at the interface of atmosphere, biosphere and pedosphere function.In contrast to what is widely believed, relatively simple, semi-distributed conceptual models have the potential to accommodate organizational features and their temporal evolution in an efficient way, a reason for that being that because their parameters (and their evolution over time) are effective at the modelling scale, and thus integrate natural heterogeneity within the system, they may be directly inferred from observations at the same scale, reducing the need for calibration and related problems. In particular, the emergence of new and more detailed observation systems from space will lead towards a more robust understanding of spatial organization and its evolution. This will further permit the development of relatively simple time-dynamic functional relationships that can meaningfully represent spatial patterns and their evolution over time, even in poorly gauged environments.Water Resource

Thermodynamics of saline and fresh water mixing in estuaries

The mixing of saline and fresh water is a process of energy dissipation. The freshwater flow that enters an estuary from the river contains potential energy with respect to the saline ocean water. This potential energy is able to perform work. Looking from the ocean to the river, there is a gradual transition from saline to fresh water and an associated rise in the water level in accordance with the increase in potential energy. Alluvial estuaries are systems that are free to adjust dissipation processes to the energy sources that drive them, primarily the kinetic energy of the tide and the potential energy of the river flow and to a minor extent the energy in wind and waves. Mixing is the process that dissipates the potential energy of the fresh water. The maximum power (MP) concept assumes that this dissipation takes place at maximum power, whereby the different mixing mechanisms of the estuary jointly perform the work. In this paper, the power is maximized with respect to the dispersion coefficient that reflects the combined mixing processes. The resulting equation is an additional differential equation that can be solved in combination with the advection-dispersion equation, requiring only two boundary conditions for the salinity and the dispersion. The new equation has been confronted with 52 salinity distributions observed in 23 estuaries in different parts of the world and performs very well.Water Resource

The physics behind Van der Burgh's empirical equation, providing a new predictive equation for salinity intrusion in estuaries

The practical value of the surprisingly simple Van der Burgh equation in predicting saline water intrusion in alluvial estuaries is well documented, but the physical foundation of the equation is still weak. In this paper we provide a connection between the empirical equation and the theoretical literature, leading to a theoretical range of Van der Burgh's coefficient of 1ĝ•2 < K < 2ĝ•3 for density-driven mixing which falls within the feasible range of 0 < K < 1. In addition, we developed a one-dimensional predictive equation for the dispersion of salinity as a function of local hydraulic parameters that can vary along the estuary axis, including mixing due to tide-driven residual circulation. This type of mixing is relevant in the wider part of alluvial estuaries where preferential ebb and flood channels appear. Subsequently, this dispersion equation is combined with the salt balance equation to obtain a new predictive analytical equation for the longitudinal salinity distribution. Finally, the new equation was tested and applied to a large database of observations in alluvial estuaries, whereby the calibrated K values appeared to correspond well to the theoretical range.Water Resource

Maximum power of saline and fresh water mixing in estuaries

According to Kleidon (2016), natural systems evolve towards a state of maximum power, leading to higher levels of entropy production by different mechanisms, including gravitational circulation in alluvial estuaries. Gravitational circulation is driven by the potential energy of fresh water. Due to the density difference between seawater and river water, the water level on the riverside is higher. The hydrostatic forces on both sides are equal but have different lines of action. This triggers an angular moment, providing rotational kinetic energy to the system, part of which drives mixing by gravitational circulation, lifting up heavier saline water from the bottom and pushing down relatively fresh water from the surface against gravity; the remainder is dissipated by friction while mixing.With a constant freshwater discharge over a tidal cycle, it is assumed that the gravitational circulation in the estuarine system performs work at maximum power. This rotational flow causes the spread of salinity inland, which is mathematically represented by the dispersion coefficient. In this paper, a new equation is derived for the dispersion coefficient related to density-driven mixing, also called gravitational circulation. Together with the steady-state advection-dispersion equation, this results in a new analytical model for densitydriven salinity intrusion. The simulated longitudinal salinity profiles have been confronted with observations in a myriad of estuaries worldwide. It shows that the performance is promising in 18 out of 23 estuaries that have relatively large convergence length. Finally, a predictive equation is presented to estimate the dispersion coefficient at the downstream boundary. Overall, the maximum power concept has provided a new physically based alternative for existing empirical descriptions of the dispersion coefficient for gravitational circulation in alluvial estuaries.Water Resource

An analytical approach to determining resonance in semi-closed convergent tidal channels

An analytical model is used to investigate the resonant behavior in a semi-closed channel. The main integral quantities of the tidal wave are obtained by means of a linearized one-dimensional model as a function of three dimensionless parameters, representing cross-section convergence, friction and distance to the closed boundary. Arbitrary along-channel variations of width and depth are accounted for by using a multi-reach approach, whereby the main tidal dynamics are reconstructed by solving a set of linear equations satisfying the continuity conditions of water level and discharge at the junctions of the sub-reaches. We highlight the importance of depth variation in the momentum equation, which is not considered in the classical tidal theory. The model allows for a direct characterization of the resonant response and for the understanding of the relative importance of the controlling parameters, highlighting the role of convergence and friction. Subsequently, the analytical model is applied to the Bristol Channel and the Guadalquivir estuary. The proposed analytical relations provide direct insights into the tidal resonance in terms of tidal forcing, geometry and friction, which will be useful for the study of semi-closed tidal channels that experience relatively large tidal ranges at the closed end.Water Resource

HESS Opinions: Are soils overrated in hydrology?

Traditional hydrological theories are based on the assumption that soil is key in determining water's fate in the hydrological cycle. According to these theories, soil hydraulic properties determine water movement in both saturated and unsaturated zones, described by matrix flow formulas such as the Darcy–Richards equations. They also determine plant-available moisture and thereby control transpiration. Here we argue that these theories are founded on a wrong assumption. Instead, we advocate the reverse: the terrestrial ecosystem manipulates the soil to satisfy specific water management strategies, which are primarily controlled by the ecosystem's reaction to climatic drivers and by prescribed boundary conditions such as topography and lithology. According to this assumption, soil hydraulic properties are an effect rather than a cause of water movement. We further argue that the integrated hydrological behaviour of an ecosystem can be inferred from considerations about ecosystem survival and growth without relying on internal-process descriptions. An important and favourable consequence of this climate- and ecosystem-driven approach is that it provides a physical justification for catchment models that do not rely on soil information and on the complexity associated with the description of soil water dynamics. Another consequence is that modelling water movement in the soil, if required, can benefit from the constraints that are imposed by the embedding ecosystem. Here we illustrate our ecosystem perspective of hydrological processes and the arguments that support it. We suggest that advancing our understanding of ecosystem water management strategies is key to building more realistic hydrological theories and catchment models that are predictive in the context of environmental change.Water Resource

Rainfall-runoff modelling using river-stage time series in the absence of reliable discharge information: A case study in the semi-arid Mara River basin

Hydrological models play an important role in water resources management. These models generally rely on discharge data for calibration. Discharge time series are normally derived from observed water levels by using a rating curve. However, this method suffers from many uncertainties due to insufficient observations, inadequate rating curve fitting procedures, rating curve extrapolation, and temporal changes in the river geometry. Unfortunately, this problem is prominent in many African river basins. In this study, an alternative calibration method is presented using water-level time series instead of discharge, applied to a semi-distributed rainfall-runoff model for the semi-arid and poorly gauged Mara River basin in Kenya. The modelled discharges were converted into water levels using the Strickler-Manning formula. This method produces an additional model output; this is a geometric rating curve equation that relates the modelled discharge to the observed water level using the Strickler-Manning formula and a calibrated slope-roughness parameter. This procedure resulted in good and consistent model results during calibration and validation. The hydrological model was able to reproduce the water levels for the entire basin as well as for the Nyangores sub-catchment in the north. The newly derived geometric rating curves were subsequently compared to the existing rating curves. At the catchment outlet of the Mara, these differed significantly, most likely due to uncertainties in the recorded discharge time series. However, at the Nyangores sub-catchment, the geometric and recorded discharge were almost identical. In conclusion, the results obtained for the Mara River basin illustrate that with the proposed calibration method, the water-level time series can be simulated well, and that the discharge-water-level relation can also be derived, even in catchments with uncertain or lacking rating curve information.Water Resource

Satellite-based drought analysis in the Zambezi River Basin: Was the 2019 drought the most extreme in several decades as locally perceived?

Study region: The study area is the river basin upstream of the Kariba dam located in the Zambezi River at the border of Zambia and Zimbabwe. Study focus: During the dry season of 2019 in Sub-Saharan Africa, extremely low water levels occurred in the Zambezi. According to news media, locals perceived this drought as the worst in several decades. We analyzed the 2019 drought in the Zambezi River Basin upstream of the Kariba dam to determine whether it indeed was the longest, most intense, and severe drought, in terms of precipitation, total water storage and reservoir water level observations over recent decades. New hydrological insights for the region: Data analysis indicates that the 2019 drought indeed had the lowest basin-averaged annual rainfall, most severe local rainfall deficit in the north of the basin, and lowest reservoir level since 1995. However, the rainfall deficit was more severe in 2002, both basin-wide and locally in the south of the basin. The total storage deficit was more severe in 2004, both basin-wide and locally in the central part of the basin. However, as the available storage data did not cover the entire deficit for 2019, its final duration and severity remain unknown. Therefore, it depends on the drought characteristic, hydrological variable, and location within the basin, whether the 2019 drought was indeed the most extreme over recent decades.Water Resource