Globally, many rivers are experiencing declining water quality, for example, with altered levels of sediments, salts, and nutrients. Effective water quality management requires a sound understanding of how and why water quality differs across space, both within and between river catchments. Land cover, land use, land management, atmospheric deposition, geology and soil type, climate, topography, and catchment hydrology are the key features of a catchment that affect: (1) the amount of suspended sediment, nutrient, and salt concentrations in catchments (i.e., the source), (2) the mobilization ,and (3) the delivery of these constituents to receiving waters. There are, however, complexities in the relationship between landscape characteristics and stream water quality. The strength of this relationship can be influenced by the distance and spatial arrangement of constituent sources within the catchment, cross correlations between landscape characteristics, and seasonality. A knowledge gap that should be addressed in future studies is that of interactions and cross correlations between landscape characteristics. There is currently limited understanding of how the relationships between landscape characteristics and water quality responses can shift based on the other characteristics of the catchment. Understanding the many forces driving stream water quality and the complexities and interactions in these forces is necessary for the development of successful water quality management strategies. This knowledge could be used to develop predictive models, which would aid in forecasting of riverine water quality.WIREs Water2018, 5:e1260. doi: 10.1002/wat2.1260This article is categorized under: Science of Water > Hydrological Processes Science of Water > Water Quality

Groundwater recharge is one of the important factors determining the groundwater development potential of an area. Even though recharge plays a key role in controlling groundwater system dynamics, much uncertainty remains regarding the relationships between groundwater recharge and its governing factors at a large scale. Therefore, this study aims to identify the most influential factors of groundwater recharge, and to develop an empirical model to estimate diffuse rainfall recharge at a global scale. Recharge estimates reported in the literature from various parts of the world (715 sites) were compiled and used in model building and testing exercises. Unlike conventional recharge estimates from water balance, this study used a multimodel inference approach and information theory to explain the relationship between groundwater recharge and influential factors, and to predict groundwater recharge at 0.5° resolution. The results show that meteorological factors (precipitation and potential evapotranspiration) and vegetation factors (land use and land cover) had the most predictive power for recharge. According to the model, long-term global average annual recharge (1981–2014) was 134 mm yr⁻¹ with a prediction error ranging from −8 to 10 mm yr⁻¹ for 97.2 % of cases. The recharge estimates presented in this study are unique and more reliable than the existing global groundwater recharge estimates because of the extensive validation carried out using both independent local estimates collated from the literature and national statistics from the Food and Agriculture Organization (FAO). In a water-scarce future driven by increased anthropogenic development, the results from this study will aid in making informed decisions about groundwater potential at a large scale.

Uncertainty analysis is an integral part of any scientific modeling, particularly within the domain of hydrological sciences given the various types and sources of uncertainty. At the center of uncertainty rests the concept of equifinality, that is, reaching a given endpoint (finality) through different pathways. The operational definition of equifinality in hydrological modeling is that various model structures and/or parameter sets (i.e., equal pathways) are equally capable of reproducing a similar (not necessarily identical) hydrological outcome (i.e., finality). Here we argue that there is more to model equifinality than model structures/parameters, that is, other model components can give rise to model equifinality and/or could be used to explore equifinality within model space. We identified six facets of model equifinality, namely, model structure, parameters, performance metrics, initial and boundary conditions, inputs, and internal fluxes. Focusing on model internal fluxes, we developed a methodology called flux mapping that has fundamental implications in understanding and evaluating model process representation within the paradigm of multiple working hypotheses. To illustrate this, we examine the equifinality of runoff fluxes of a conceptual rainfall‐runoff model for a number of different Australian catchments. We demonstrate how flux maps can give new insights into the model behavior that cannot be captured by conventional model evaluation methods. We discuss the advantages of flux space, as a subspace of the model space not usually examined, over parameter space. We further discuss the utility of flux mapping in hypothesis generation and testing, extendable to any field of scientific modeling of open complex systems under uncertainty.

Evidence suggests that catchment state variables such as groundwater can exhibit multiyear trends. This means that their state may reflect not only recent climatic conditions but also climatic conditions in past years or even decades. Here we demonstrate that five commonly used conceptual “bucket” rainfall‐runoff models are unable to replicate multiyear trends exhibited by natural systems during the “Millennium Drought” in south‐east Australia. This causes an inability to extrapolate to different climatic conditions, leading to poor performance in split sample tests. Simulations are examined from five models applied in 38 catchments, then compared with groundwater data from 19 bores and Gravity Recovery and Climate Experiment data for two geographic regions. Whereas the groundwater and Gravity Recovery and Climate Experiment data decrease from high to low values gradually over the duration of the 13‐year drought, the model storages go from high to low values in a typical seasonal cycle. This is particularly the case in the drier, flatter catchments. Once the drought begins, there is little room for decline in the simulated storage, because the model “buckets” are already “emptying” on a seasonal basis. Since the effects of sustained dry conditions cannot accumulate within these models, we argue that they should not be used for runoff projections in a drying climate. Further research is required to (a) improve conceptual rainfall‐runoff models, (b) better understand circumstances in which multiyear trends in state variables occur, and (c) investigate links between these multiyear trends and changes in rainfall‐runoff relationships in the context of a changing climate.

It has been widely shown that rainfall‐runoff models often provide poor and biased simulations after a change in climate, but evidence suggests existing models may be capable of better simulations if calibration strategies are improved. Common practice is to use “least squares”‐type objective functions, which focus on hydrological behavior during high flows. However, simulation of a drying climate may require a more balanced consideration of other parts of the flow regime, including mid‐low flows and drier years in the calibration period, as a closer analogue of future conditions. Here we systematically test eight objective functions over 86 catchments and five conceptual model structures in southern and eastern Australia. We focus on performance when evaluated over multiyear droughts. The results show significant improvements are possible compared to least squares calibration. In particular, the Refined Index of Agreement (based on sum of absolute error, not sum of squared error) and a new objective function called the Split KGE (which gives equal weight to each year in the calibration series) give significantly better split‐sample results than least squares approaches. This improvement held for all five model structures, regardless of basin characteristics such as slope, vegetation, and across a range of climatic conditions (e.g., mean precipitation between 500 and 1,500 mm/yr). We recommend future studies to avoid least squares approaches (e.g., optimizing NSE or KGE with no prior transformation on streamflow) and adopt these alternative methods, wherever simulations in a drying climate are required.

While the majority of hydrological prediction methods assume that observed interannual variability explores the full range of catchment response dynamics, recent cases of prolonged climate drying suggest otherwise. During the ∼decade‐long Millennium drought in south‐eastern Australia significant shifts in hydrologic behavior were reported. Catchment rainfall‐runoff partitioning changed from what was previously encountered during shorter droughts, with significantly less runoff than expected occurring in many catchments. In this article, we investigate the variability in the magnitude of shift in rainfall‐runoff partitioning observed during the Millennium drought. We re‐evaluate a large range of factors suggested to be responsible for the additional runoff reductions. Our results suggest that the shifts were mostly influenced by catchment characteristics related to predrought climate (aridity index and rainfall seasonality) and soil and groundwater storage dynamics (predrought interannual variability of groundwater storage and mean solum thickness). The shifts were amplified by seasonal rainfall changes during the drought (spring rainfall deficits). We discuss the physical mechanisms that are likely to be associated with these factors. Our results confirm that shifts in the annual rainfall‐runoff relationship represent changes in internal catchment functioning, and emphasize the importance of cumulative multiyear changes in the catchment storage for runoff generation. Prolonged drying in some regions can be expected in the future, and our results provide an indication of which catchments characteristics are associated with catchments more susceptible to a shift in their runoff response behavior.