In many rainfall-runoff models, at least some calibration of model parameters has to take place. Especially for ungauged or poorly gauged basins this can be problematic, because there is little or no data available for calibration. A possible solution to overcome the problems caused by data scarcity is to set up a measurement campaign for a limited time period. In this study, we determine the minimum amount of data required to determine robust parameter values for a simple model with two parameters. The model is constructed such that the parameters can be determined not only with automatic calibration, but also by recession analysis and a priori from Boussinesq theory. The model has been applied to a research catchment in Switzerland. For automatic calibration and recession analysis, one season (5 months) is found to be sufficient to give robust parameters for simulation of high flows over the full observation period. For automatic calibration, this should be the season with the highest precipitation, for recession analysis the season with least evapotranspiration. The Boussinesq equation is able to give good parameter estimates for modeling high flows, but detailed in situ knowledge of the catchment is required. Automatic calibration outperforms recession analysis and Boussinesq theory by far when it comes to parameter estimation with a focus on prediction of low flows. It was shown that a single set of parameters cannot simultaneously describe high and low flows with a reasonable accuracy, suggesting that more than two parameters are needed to characterize subsurface properties. Key Points Automatic calibration, recession analysis, and Boussinesq theory are compared Automatic calibration leads to the highest model efficiencies One season of data is needed for robust parameter estimation@en

The Wageningen Lowland Runoff Simulator (WALRUS) is a new parametric (conceptual) rainfall-runoff model which accounts explicitly for processes that are important in lowland areas, such as groundwater-unsaturated zone coupling, wetness-dependent flowroutes, groundwater-surface water feedbacks, and seepage and surface water supply (see companion paper by Brauer et al., 2014). Lowland catchments can be divided into slightly sloping, freely draining catchments and flat polders with controlled water levels. Here, we apply WALRUS to two contrasting Dutch catchments: the Hupsel Brook catchment and the Cabauw polder. In both catchments, WALRUS performs well: Nash-Sutcliffe efficiencies obtained after calibration on 1 year of discharge observations are 0.87 for the Hupsel Brook catchment and 0.83 for the Cabauw polder, with values of 0.74 and 0.76 for validation. The model also performs well during floods and droughts and can forecast the effect of control operations. Through the dynamic division between quick and slow flowroutes controlled by a wetness index, temporal and spatial variability in groundwater depths can be accounted for, which results in adequate simulation of discharge peaks as well as low flows. The performance of WALRUS is most sensitive to the parameter controlling the wetness index and the groundwater reservoir constant, and to a lesser extent to the quickflow reservoir constant. The effects of these three parameters can be identified in the discharge time series, which indicates that the model is not overparameterised (parsimonious). Forcing uncertainty was found to have a larger effect on modelled discharge than parameter uncertainty and uncertainty in initial conditions.@en

The transfer of parameter sets over different temporal and spatial resolutions is common practice in many large-domain hydrological modelling studies. The degree to which parameters are transferable across temporal and spatial resolutions is an indicator of how well spatial and temporal variability is represented in the models. A large degree of transferability may well indicate a poor representation of such variability in the employed models. To investigate parameter transferability over resolution in time and space we have set up a study in which the Variable Infiltration Capacity (VIC) model for the Thur basin in Switzerland was run with four different spatial resolutions (1 km × 1 km, 5 km × 5 km, 10 km × 10 km, lumped) and evaluated for three relevant temporal resolutions (hour, day, month), both applied with uniform and distributed forcing. The model was run 3150 times using the Hierarchical Latin Hypercube Sample and the best 1 % of the runs was selected as behavioural. The overlap in behavioural sets for different spatial and temporal resolutions was used as an indicator of parameter transferability. A key result from this study is that the overlap in parameter sets for different spatial resolutions was much larger than for different temporal resolutions, also when the forcing was applied in a distributed fashion. This result suggests that it is easier to transfer parameters across different spatial resolutions than across different temporal resolutions. However, the result also indicates a substantial underestimation in the spatial variability represented in the hydrological simulations, suggesting that the high spatial transferability may occur because the current generation of large-domain models has an inadequate representation of spatial variability and hydrologic connectivity. The results presented in this paper provide a strong motivation to further investigate and substantially improve the representation of spatial and temporal variability in large-domain hydrological models.@en

A meta-analysis on 192 peer-reviewed articles reporting on applications of the variable infiltration capacity (VIC) model in a distributed way reveals that the spatial resolution at which the model is applied has increased over the years, while the calibration and validation time interval has remained unchanged. We argue that the calibration and validation time interval should keep pace with the increase in spatial resolution in order to resolve the processes that are relevant at the applied spatial resolution. We identified six time concepts in hydrological models, which all impact the model results and conclusions. Process-based model evaluation is particularly relevant when models are applied at hyper-resolution, where stakeholders expect credible results both at a high spatial and temporal resolution.@en

In order to investigate to what extent rainfall fluctuations observed with different types of instruments reflect the properties of the rainfall process itself and to what extent they are merely instrumental artefacts we are in the process of developing a stochastic model of rainfall. The starting point for the development of the model has been the notion that at the spatial and temporal scales associated with many types of surface rainfall measurements, rainfall is a discrete process describing the arrival of raindrops of different sizes at the ground. A fundamental question is whether this raindrop arrival process can be considered a homogeneous (Poisson) process or whether it behaves as a clustering (or possibly even scaling) process, as has recently been proposed in the literature. We have tested the classical Poisson homogeneity hypothesis in rainfall on a 35 min time series of 10 s raindrop size spectra collected with a 50 cm2 optical disdrometer. The rain rates calculated from the spectra indicated roughly uncorrelated fluctuations around a constant mean rain rate of about 3.5 mm h-1. Two types of analysis of the drop counts were carried out, a global analysis taking into account all drops regardless of their size and an analysis considering the drop counts in the 16 0.21 mm diameter intervals separately. The first type of analysis revealed that even for the more or less stationary time series under consideration the total raindrop arrival rate was overdispersed with respect to the homogeneous Poisson process. The second type of analysis demonstrated that this rejection of the homogeneity hypothesis could be attributed entirely to raindrops with diameters smaller than 1.14 mm. Although these drops account for 66% of the raindrop concentration in the air and 55% of the raindrop arrival rate at the ground, they only account for 14% of the rain rate and 2% of the radar reflectivity factor (on the basis of the mean drop size distribution during the experiment). In other words, although clustering may be a significant phenomenon for the smallest raindrops, the analyzed data seem to indicate that for moderate rain rates the arrival rate fluctuations of the raindrops which contribute most to rain rate and radar reflectivity factor behave according to Poisson statistics.@en

A meta-analysis on 192 peer-reviewed articles reporting applications of the Variable Infiltration Capacity (VIC) model in a distributed way reveals that the spatial resolution at which the model is applied has increased over the years, while the calibration and validation time interval has remained unchanged. We argue that the calibration and validation time interval should keep pace with the increase in spatial resolution in order to resolve the processes that are relevant at the applied spatial resolution. We identified six time concepts in hydrological models, which all impact the model results and conclusions. Process-based model evaluation is particularly relevant when models are applied at hyper-resolution, where stakeholders expect credible results both at a high spatial and temporal resolution.@en

Existing tidal input reduction approaches applied in accelerated morphodynamic simulations aim to capture the dominant tidal forces in a single or double representative tidal cycle, often referred to as a “morphological tide.” These strongly simplified tidal signals fail to represent the tidal extremes and hence poorly allow to represent hydrodynamics in the intertidal areas. Here, a generic method is developed to construct a synthetic spring-neap tidal cycle that (a) represents the original signal; (b) is exactly periodic; and (c) is derived directly from tidal time series or harmonic constituents. The starting point is a fortnightly modulation of the semidiurnal tide to represent spring-neap variations, while conserving periodicity. Diurnal tides and higher harmonics of the semidiurnal tide are included to represent the asymmetry of the tide. The amplitudes and phases of the synthetic signal are then fitted to histograms of water levels and water level gradients derived from the original sea surface elevation time series. A depth-averaged model of the Ems estuary (The Netherlands) demonstrates the effects of alternative tidal input reduction techniques. Adopting the new approach, the along-estuary variation in tidal wave shape is well-represented, leading to an improved representation of extreme tidal conditions. Especially the more realistic representation of intertidal dynamics improves the overall hydrodynamics and residual sand transport patterns, approaching nonschematized tidal dynamics.@en

It is generally acknowledged in the environmental sciences that the choice of a computational model impacts the research results. In this study of a flood and drought event in the Swiss Thur basin, we show that modeling decisions during the model configuration, beyond the model choice, also impact the model results. In our carefully designed experiment we investigated four modeling decisions in ten nested basins: the spatial resolution of the model, the spatial representation of the forcing data, the calibration period, and the performance metric. The flood characteristics were mainly affected by the performance metric, whereas the drought characteristics were mainly affected by the calibration period. The results could be related to the processes that triggered the particular events studied. The impact of the modeling decisions on the simulations did, however, vary among the investigated sub-basins. In spite of the limitations of this study, our findings have important implications for the understanding and quantification of uncertainty in any hydrological or even environmental model. Modeling decisions during model configuration introduce subjectivity from the modeler. Multiple working hypotheses during model configuration can provide insights on the impact of such subjective modeling decisions.@en

Using 30. years (1979-2009) of data from 33 automatic rain gauges in the Netherlands, a study of the space-time variability of rainfall is performed. This study uses 90-day averaged semi-variograms to find seasonal signals in the fitted spherical semi-variogram parameters of rainrates in the Netherlands, for accumulation intervals between 1 and 24. h. These signals can be well-described by simple cosine functions. The dependence of these cosine functions on the accumulation interval is modeled in two different ways: (1) power-law relations between the variogram parameters and the accumulation interval, and (2) power-law relations between the parameters of the cosine functions and the accumulation interval. For the first method the cosine function at the 24-h accumulation interval is also needed. The second of these methods has more parameters, but is shown to model the temporal scaling best. The space-time scaling relations found in this paper are compared to those found by others for similar and contrasting climates. Seasonality is shown to play an important role in determining rainfall spatial variability.@en

The relation between storage and discharge is an essential characteristic of many rainfall-runoff models. The simple dynamical systems approach, in which a rainfall-runoff model is constructed from a single storage-discharge relation, has been successfully applied to humid catchments. Here we investigate (1) if and when the less humid lowland Hupsel Brook catchment also behaves like a simple dynamical system by hydrograph fitting, and (2) if system parameters can be inferred from streamflow recession rates or more directly from soil moisture storage observations. Only 39% of the fitted monthly hydrographs yielded Nash-Sutcliffe efficiencies above 0.5, from which we can conclude that the Hupsel Brook catchment does not always behave like a simple dynamical system. Model results were especially poor in summer, when evapotranspiration is high and the thick unsaturated zone attenuates the rainfall input. Using soil moisture data to obtain system parameters is not trivial, mainly because there is a discrepancy between local and catchment storage. Parameters obtained with direct storage-discharge fitting led to a strong underestimation of the response of runoff to rainfall, while recession analysis leads to an overestimation.@en

This paper presents a novel approach to estimate the vertical profile of reflectivity (VPR) from volumetric weather radar data using both a traditional Eulerian as well as a newly proposed Lagrangian implementation. For this latter implementation, the recently developed Rotational Carpenter Square Cluster Algorithm (RoCaSCA) is used to delineate precipitation regions at different reflectivity levels. A piecewise linear VPR is estimated for either stratiform or neither stratiform/convective precipitation. As a second aspect of this paper, a novel approach is presented which is able to account for the impact of VPR uncertainty on the estimated radar rainfall variability. Results show that implementation of the VPR identification and correction procedure has a positive impact on quantitative precipitation estimates from radar. Unfortunately, visibility problems severely limit the impact of the Lagrangian implementation beyond distances of 100 km. However, by combining this procedure with the global Eulerian VPR estimation procedure for a given rainfall type (stratiform and neither stratiform/convective), the quality of the quantitative precipitation estimates increases up to a distance of 150 km. Analyses of the impact of VPR uncertainty shows that this aspect accounts for a large fraction of the differences between weather radar rainfall estimates and rain gauge measurements. Key Points A new Lagrangian based VPR identification methodRadar rainfall uncertainty estimation from spatial precipitation variabilityImproved weather radar based surface rainfall estimation@en

We discuss two definitions of reproducibility, and question if both definitions are required to be met in computational hydrological studies.@en

Borneo's diverse ecosystems, which are typical humid tropical conditions, are deteriorating rapidly, as the area is experiencing recurrent large-scale wildfires, affecting atmospheric composition and influencing regional climate processes. Studies suggest that climate-driven drought regulates wildfires, but these overlook subsurface processes leading to hydrological drought, an important driver. Here, we show that models which include hydrological processes better predict area burnt than those solely based on climate data. We report that the Borneo landscape has experienced a substantial hydrological drying trend since the early twentieth century, leading to progressive tree mortality, more severe than in other tropical regions. This has caused massive wildfires in lowland Borneo during the past two decades, which we show are clustered in years with large areas of hydrological drought coinciding with strong El Niño events. Statistical modelling evidence shows amplifying wildfires and greater area burnt in response to El Niño/Southern Oscillation (ENSO) strength, when hydrology is considered. These results highlight the importance of considering hydrological drought for wildfire prediction, and we recommend that hydrology should be considered in future studies of the impact of projected ENSO strength, including effects on tropical ecosystems, and biodiversity conservation.@en

Lowlands are vulnerable to flooding due to their mild topography in often densely populated areas with high social and economic value. Moreover, multiple physical processes coincide in lowland areas, such as those involved in river-sea interactions and in merging rivers at confluences. Simultaneous occurrence of such processes can result in amplifying or attenuating effects on water levels. Our aim is to understand the mechanisms behind simultaneous occurrence of discharge waves in a river and its lowland tributaries. Here, we introduce a new way of analyzing lowland discharge and water level dynamics, by tracing individual flood waves based on dynamic time warping. We take the confluence of the Meuse River (∼ 33000km2) with the joining tributaries of the Dommel and Aa rivers as an example, especially because the January 1995 flood at this confluence was the result of the simultaneous occurrence of discharge peaks in the main stream and the tributaries and because independent observations of water levels and discharge are available for a longer period. The analysis shows that the exact timing of the arrival of discharge peaks is of little relevance because of the long duration of the average discharge wave compared to typical time lags between peaks. The discharge waves last on average 9 days, whereas the lag time between discharge peaks in the main river and the tributaries is typically 3 days. This results in backwaters that can rise up to 1.5 m over a distance of 4 km from the confluence. Thus, local measures to reduce the impact of flooding around the confluence should account for the long duration of flood peaks in the main system.@en

Understanding the long-term (interannual-decadal) variability of water availability in river basins is paramount for water resources management. Here, the authors analyze time series of simulated terrestrial water storage components, observed precipitation, and discharge spanning 74 yr in the Colorado River basin and relate them to climate indices that describe variability of sea surface temperature and sea level pressure in the tropical and extratropical Pacific. El Niño-Southern Oscillation (ENSO) indices in winter [January-March (JFM)] are related to winter precipitation as well as to soil moisture and discharge in the lower Colorado River basin. The low-frequency mode of the Pacific decadal oscillation (PDO) appears to be strongly correlated with deep soil moisture. During the negative PDO phase, saturated storage anomalies tend to be negative and the "amplitudes" (mean absolute anomalies) of shallow soil moisture, snow, and discharge are slightly lower compared to periods of positive PDO phases. Predicting interannual variability, therefore, strongly depends on the capability of predicting PDO regime shifts. If indeed a shift to a cool PDO phase occurred in the mid-1990s, as data suggest, the current dry conditions in the Colorado River basin may persist.@en

The Meuse is an important river in Western Europe, which is almost exclusively rain-fed. Projected changes in precipitation characteristics due to climate change, therefore, are expected to have a considerable effect on the hydrological regime of the river Meuse. We focus on an important tributary of the Meuse, the Ourthe, measuring about 1600 km2. The well-known hydrological model HBV is forced with three high-resolution (0.088°) regional climate scenarios, each based on one of three different IPCC CO2 emission scenarios: A1B, A2 and B1. To represent the current climate, a reference model run at the same resolution is used. Prior to running the hydrological model, the biases in the climate model output are investigated and corrected for. Different approaches to correct the distributed climate model output using single-site observations are compared. Correcting the spatially averaged temperature and precipitation is found to give the best results, but still large differences exist between observations and simulations. The bias corrected data are then used to force HBV. Results indicate a small increase in overall discharge, especially for the B1 scenario during the beginning of the 21st century. Towards the end of the century, all scenarios show a decrease in summer discharge, partially because of the diminished buffering effect by the snow pack, and an increased discharge in winter. It should be stressed, however, that we used results from only one GCM (the only one available at such a high resolution). It would be interesting to repeat the analysis with multiple models.@en

Travel time distributions are often used to characterize catchment discharge behavior, catchment vulnerability to pollution and pollutant loads from catchments to downstream waters. However, these distributions vary with time because they are a function of rainfall and evapotranspiration. It is important to account for these variations when the time scale of interest is smaller than the typical time-scale over which average travel time distributions can be derived. Recent studies have suggested that subsurface mixing controls how rainfall and evapotranspiration affect the variability in travel time distributions of discharge. To quantify this relation between subsurface mixing and dynamics of travel time distributions, we propose a new transformation of travel time that yields transformed travel time distributions, which we call Storage Outflow Probability (STOP) functions. STOP functions quantify the probability for water parcels in storage to leave a catchment via discharge or evapotranspiration. We show that this is equal to quantifying mixing within a catchment. Compared to the similar Age function introduced by Botter et al. (2011), we show that STOP functions are more constant in time, have a clearer physical meaning and are easier to parameterize. Catchment-scale STOP functions can be approximated by a two-parameter beta distribution. One parameter quantifies the catchment preference for discharging young water; the other parameter quantifies the preference for discharging old water from storage. Because of this simple parameterization, the STOP function is an innovative tool to explore the effects of catchment mixing behavior, seasonality and climate change on travel time distributions and the related catchment vulnerability to pollution spreading.@en

Because of global warming, the hydrologic behavior of the Rhine basin is expected to shift from a combined snowmelt-and rainfall-driven regime to a more rainfall-dominated regime. Previous impact assessments have indicated that this leads, on average, to increasing streamflow by ~30% in winter and spring and decreasing streamflow by a similar value in summer. In this study, high-resolution (0.088°) regional climate scenarios conducted with the regional climate model REMO (REgional MOdel) for the Rhine basin are used to force a macroscale hydrological model. These climate scenarios are based on model output from the ECHAM5-Max Planck Institute Ocean Model (MPI-OM) global climate model, which is in turn forced by three Special Report on Emissions Scenarios (SRES) emission scenarios: A2, A1B, and B1. The Variable Infiltration Capacity model (VIC; version 4.0.5) is used to examine changes in streamflow at various locations throughout the Rhine basin. Average streamflow, peak flows, low flows, and several water balance terms are evaluated for both the first and second half of the twenty-first century. The results reveal a distinct contrast between those periods. The first half is dominated by increased precipitation, causing increased streamflow throughout the year. During the second half of the century, a streamflow increase in winter/spring and a decrease in summer is found, similar to previous studies. This is caused by 1) temperature and evapotranspiration, which are considerably higher during the second half of the century; 2) decreased precipitation in summer; and 3) an earlier start of the snowmelt season. Magnitudes of peak flows increase during both periods, and the magnitudes of streamflow droughts increase only during the second half of the century.@en

In many climate impact studies hydrological models are forced with meteorological data without an attempt to assess the quality of these data. The objective of this study was to compare downscaled ERA15 (ECMWF-reanalysis data) precipitation and temperature with observed precipitation and temperature and apply a bias correction to these forcing variables. Precipitation is corrected by fitting it to the mean and coefficient of variation (CV) of the observations. Temperature is corrected by fitting it to the mean and standard deviation of the observations. It appears that the uncorrected ERA15 is too warm and too wet for most of the Rhine basin. The bias correction leads to satisfactory results, precipitation and temperature differences decreased significantly, although there are a few years for which the correction of precipitation is less satisfying. Corrections were largest during summer for both precipitation and temperature. For precipitation alone large corrections were applied during September and October as well. Besides the statistics the correction method was intended to correct for, it is also found to improve the correlations for the fraction of wet days and lag-1 autocorrelations between ERA15 and the observations. For the validation period temperature is corrected very well, but for precipitation the RMSE of the daily difference between modeled and observed precipitation has increased for the corrected situation. When taking random years for calibration, and the remaining years for validation, the spread in the mean bias error (MBE) becomes larger for the corrected precipitation during validation, but the overal average MBE has decreased.@en

Hydrologic projections are of vital socio-economic importance. However, they are also prone to uncertainty. In order to establish a meaningful range of storylines to support water managers in decision making, we need to reveal the relevant sources of uncertainty. Here, we systematically and extensively investigate uncertainty in hydrologic projections for 605 basins throughout the contiguous US. We show that in the majority of the basins, the sign of change in average annual runoff and discharge timing for the period 2070-2100 compared to 1985-2008 differs among combinations of climate models, hydrologic models, and parameters. Mapping the results revealed that different sources of uncertainty dominate in different regions. Hydrologic model induced uncertainty in the sign of change in mean runoff was related to snow processes and aridity, whereas uncertainty in both mean runoff and discharge timing induced by the climate models was related to disagreement among the models regarding the change in precipitation. Overall, disagreement on the sign of change was more widespread for the mean runoff than for the discharge timing. The results demonstrate the need to define a wide range of quantitative hydrologic storylines, including parameter, hydrologic model, and climate model forcing uncertainty, to support water resource planning.@en