This paper is the outcome of a community initiative to identify major unsolved scientific problems in hydrology motivated by a need for stronger harmonisation of research efforts. The procedure involved a public consultation through on-line media, followed by two workshops through which a large number of potential science questions were collated, prioritised, and synthesised. In spite of the diversity of the participants (230 scientists in total), the process revealed much about community priorities and the state of our science: a preference for continuity in research questions rather than radical departures or redirections from past and current work. Questions remain focussed on process-based understanding of hydrological variability and causality at all space and time scales. Increased attention to environmental change drives a new emphasis on understanding how change propagates across interfaces within the hydrological system and across disciplinary boundaries. In particular, the expansion of the human footprint raises a new set of questions related to human interactions with nature and water cycle feedbacks in the context of complex water management problems. We hope that this reflection and synthesis of the 23 unsolved problems in hydrology will help guide research efforts for some years to come.

New precipitation (P) datasets are released regularly, following innovations in weather forecasting models, satellite retrieval methods, and multi-source merging techniques. Using the conterminous US as a case study, we evaluated the performance of 26 gridded (sub-)daily P datasets to obtain insight into the merit of these innovations. The evaluation was performed at a daily timescale for the period 2008–2017 using the Kling–Gupta efficiency (KGE), a performance metric combining correlation, bias, and variability. As a reference, we used the high-resolution (4 km) Stage-IV gauge-radar P dataset. Among the three KGE components, the P datasets performed worst overall in terms of correlation (related to event identification). In terms of improving KGE scores for these datasets, improved P totals (affecting the bias score) and improved distribution of P intensity (affecting the variability score) are of secondary importance. Among the 11 gauge-corrected P datasets, the best overall performance was obtained by MSWEP V2.2, underscoring the importance of applying daily gauge corrections and accounting for gauge reporting times. Several uncorrected P datasets outperformed gauge-corrected ones. Among the 15 uncorrected P datasets, the best performance was obtained by the ERA5-HRES fourth-generation reanalysis, reflecting the significant advances in earth system modeling during the last decade. The (re)analyses generally performed better in winter than in summer, while the opposite was the case for the satellite-based datasets. IMERGHH V05 performed substantially better than TMPA-3B42RT V7, attributable to the many improvements implemented in the IMERG satellite P retrieval algorithm. IMERGHH V05 outperformed ERA5-HRES in regions dominated by convective storms, while the opposite was observed in regions of complex terrain. The ERA5-EDA ensemble average exhibited higher correlations than the ERA5-HRES deterministic run, highlighting the value of ensemble modeling. The WRF regional convection-permitting climate model showed considerably more accurate P totals over the mountainous west and performed best among the uncorrected datasets in terms of variability, suggesting there is merit in using high-resolution models to obtain climatological P statistics. Our findings provide some guidance to choose the most suitable P dataset for a particular application.

Model intercomparison studies are carried out to test and compare the simulated outputs of various model setups over the same study domain. The Great Lakes region is such a domain of high public interest as it not only resembles a challenging region to model with its transboundary location, strong lake effects, and regions of strong human impact but is also one of the most densely populated areas in the USA and Canada. This study brought together a wide range of researchers setting up their models of choice in a highly standardized experimental setup using the same geophysical datasets, forcings, common routing product, and locations of performance evaluation across the 1×10⁶ km² study domain. The study comprises 13 models covering a wide range of model types from machine-learning-based, basin-wise, subbasin-based, and gridded models that are either locally or globally calibrated or calibrated for one of each of the six predefined regions of the watershed. Unlike most hydrologically focused model intercomparisons, this study not only compares models regarding their capability to simulate streamflow (Q) but also evaluates the quality of simulated actual evapotranspiration (AET), surface soil moisture (SSM), and snow water equivalent (SWE). The latter three outputs are compared against gridded reference datasets. The comparisons are performed in two ways – either by aggregating model outputs and the reference to basin level or by regridding all model outputs to the reference grid and comparing the model simulations at each grid-cell. The main results of this study are as follows: The comparison of models regarding streamflow reveals the superior quality of the machine-learning-based model in the performance of all experiments; even for the most challenging spatiotemporal validation, the machine learning (ML) model outperforms any other physically based model. While the locally calibrated models lead to good performance in calibration and temporal validation (even outperforming several regionally calibrated models), they lose performance when they are transferred to locations that the model has not been calibrated on. This is likely to be improved with more advanced strategies to transfer these models in space. The regionally calibrated models – while losing less performance in spatial and spatiotemporal validation than locally calibrated models – exhibit low performances in highly regulated and urban areas and agricultural regions in the USA. Comparisons of additional model outputs (AET, SSM, and SWE) against gridded reference datasets show that aggregating model outputs and the reference dataset to the basin scale can lead to different conclusions than a comparison at the native grid scale. The latter is deemed preferable, especially for variables with large spatial variability such as SWE. A multi-objective-based analysis of the model performances across all variables (Q, AET, SSM, and SWE) reveals overall well-performing locally calibrated models (i.e., HYMOD2-lumped) and regionally calibrated models (i.e., MESH-SVS-Raven and GEM-Hydro-Watroute) due to varying reasons. The machine-learning-based model was not included here as it is not set up to simulate AET, SSM, and SWE. All basin-aggregated model outputs and observations for the model variables evaluated in this study are available on an interactive website that enables users to visualize results and download the data and model outputs.

The notion that landscape features have coevolved over time is well known in the Earth sciences. Hydrologists have recently called for a more rigorous connection between emerging spatial patterns of landscape features and the hydrological response of catchments, and have termed this concept catchment coevolution. In this paper we review recent literature on this subject and attempt to synthesize what we have learned into a general framework that would improve predictions of hydrologic change. We first present empirical evidence of the interaction and feedback of landscape evolution and changes in hydrological response. From this review it is clear that the independent drivers of catchment coevolution are climate, geology, and tectonics. We identify common currency that allows comparing the levels of activity of these independent drivers, such that, at least conceptually, we can quantify the rate of evolution or aging. Knowing the hydrologic age of a catchment by itself is not very meaningful without linking age to hydrologic response. Two avenues of investigation have been used to understand the relationship between (differences in) age and hydrological response: (i) one that is based on relating present landscape features to runoff processes that are hypothesized to be responsible for the current fingerprints in the landscape; and (ii) one that takes advantage of an experimental design known as space‐for‐time substitution. Both methods have yielded significant insights in the hydrologic response of landscapes with different histories. If we want to make accurate predictions of hydrologic change, we will also need to be able to predict how the catchment will further coevolve in association with changes in the activity levels of the drivers (e.g., climate). There is ample evidence in the literature that suggests that whole‐system prediction of catchment coevolution is, at least in principle, plausible. With this imperative we outline a research agenda that implements the concepts of catchment coevolution for building a holistic framework toward improving predictions of hydrologic change.

Droughts are often long-lasting phenomena, without a distinct start or end and with impacts cascading across sectors and systems, creating long-term legacies. Nevertheless, our current perceptions and management of droughts and their impacts are often event-based, which can limit the effective assessment of drought risks and reduction of drought impacts. Here, we advocate for changing this perspective and viewing drought as a hydrological–ecological–social continuum. We take a systems theory perspective and focus on how “memory” causes feedback and interactions between parts of the interconnected systems at different timescales. We first discuss the characteristics of the drought continuum with a focus on the hydrological, ecological, and social systems separately, and then we study the system of systems. Our analysis is based on a review of the literature and a study of five cases: Chile, the Colorado River basin in the USA, northeast Brazil, Kenya, and the Rhine River basin in northwest Europe. We find that the memories of past dry and wet periods, carried by both bio-physical (e.g. groundwater, vegetation) and social systems (e.g. people, governance), influence how future drought risk manifests. We identify four archetypes of drought dynamics: impact and recovery, slow resilience building, gradual collapse, and high resilience–big shock. The interactions between the hydrological, ecological, and social systems result in systems shifting between these types, which plays out differently in the five case studies. We call for more research on drought preconditions and recovery in different systems, on dynamics cascading between systems and triggering system changes, and on dynamic vulnerability and maladaptation. Additionally, we advocate for more continuous monitoring of drought hazards and impacts, modelling tools that better incorporate memories and adaptation responses, and management strategies that increase societal and institutional memory. This will help us to better deal with the complex hydrological–ecological–social drought continuum and identify effective pathways to adaptation and mitigation.