Quantifying parametric uncertainty is important since the reliability of simulation results strongly depends on the parametrization of a model (Beven, 1995;Reinecke et al., 2019). Diagnostic model analyses are helpful to investigate parameter and process dynamics as well as to better understand process representations in hydrological models (Guse et al., 2016;Pfannerstill et al., 2015). The problem of equifinal parameter sets is well known in hydrology (Beven, 2006;Her & Chaubey, 2015;Kelleher et al., 2017) and time series of different variables (e.g., discharge, groundwater levels, and snow depth) can be used to constrain model parameters according to the physical processes that they describe (