Streamflow is often the only variable used to evaluate hydrological models. In a previous international comparison study, eight research groups followed an identical protocol to calibrate 12 hydrological models using observed streamflow of catchments within the Meuse basin. In the current study, we quantify the differences in five states and fluxes of these 12 process-based models with similar streamflow performance, in a systematic and comprehensive way. Next, we assess model behavior plausibility by ranking the models for a set of criteria using streamflow and remote-sensing data of evaporation, snow cover, soil moisture and total storage anomalies. We found substantial dissimilarities between models for annual interception and seasonal evaporation rates, the annual number of days with water stored as snow, the mean annual maximum snow storage and the size of the root-zone storage capacity. These differences in internal process representation imply that these models cannot all simultaneously be close to reality. Modeled annual evaporation rates are consistent with Global Land Evaporation Amsterdam Model (GLEAM) estimates. However, there is a large uncertainty in modeled and remote-sensing annual interception. Substantial differences are also found between Moderate Resolution Imaging Spectroradiometer (MODIS) and modeled number of days with snow storage. Models with relatively small root-zone storage capacities and without root water uptake reduction under dry conditions tend to have an empty root-zone storage for several days each summer, while this is not suggested by remote-sensing data of evaporation, soil moisture and vegetation indices. On the other hand, models with relatively large root-zone storage capacities tend to overestimate very dry total storage anomalies of the Gravity Recovery and Climate Experiment (GRACE). None of the models is systematically consistent with the information available from all different (remote-sensing) data sources. Yet we did not reject models given the uncertainties in these data sources and their changing relevance for the system under investigation.

Hydrological models are being applied for impact assessment across a wide range of resolutions. In this study, we quantify the effect of model resolution on the simulated hydrological response in five mesoscale basins in the Swiss Alps using the distributed hydrological model Spatial Processes in Hydrology (SPHY). We introduce a new metric to compare a range of values resulting from a distributed model with a single value: the density-weighted distance (DWD). Model simulations are performed at two different spatial resolutions, matching common practices in hydrology: 500m×500m matching regional-scale models, and 40km×40km matching global-scale modeling. We investigate both the intra-basin response in seasonal streamflow and evapotranspiration from the high-resolution model and the difference induced by the two different spatial resolutions, with a focus on four seasonal extremes, selected based on temperature and precipitation. Results from the high-resolution model show that the intra-basin response covers a surprisingly large range of anomalies and show that it is not uncommon to have both extreme positive and negative flux anomalies occurring simultaneously within a catchment. The intra-basin response was grouped by land cover, where different dominant runoff-generating processes are driving the differences between these groups. The low-resolution model failed to capture the diverse and contrasting response from the high-resolution model, since neither the complex topography nor land cover classes were properly represented. DWD values show that, locally, the hydrological response simulated with a high-resolution model can be a lot more extreme than a low-resolution model might indicate, which has important implications for global or continental scale assessments carried out at coarse grids of 0:5°×0:5° or 0:25°× 0:25° resolution.