An evaluation of approaches for modelling hydrological processes in high-elevation, glacierized Andean watersheds

We use two hydrological models of varying complexity to study the Juncal River Basin in the Central Andes of Chile with the aim to understand the degree of conceptualization and the spatial structure that are needed to model present and future streamflows. We use a conceptual semi-distributed model based on elevation bands [Water Evaluation and Planning (WEAP)], frequently used for water management, and a physically oriented, fully distributed model [Topographic Kinematic Wave Approximation and Integration ETH Zurich (TOPKAPI-ETH)] developed for research purposes mainly. We evaluate the ability of the two models to reproduce the key hydrological processes in the basin with emphasis on snow accumulation and melt, streamflow and the relationships between internal processes. Both models are capable of reproducing observed runoff and the evolution of Moderate-resolution Imaging Spectroradiometer snow cover adequately. In spite of WEAP's simple and conceptual approach for modelling snowmelt and its lack of glacier representation and snow gravitational redistribution as well as a proper routing algorithm, this model can reproduce historical data with a similar goodness of fit as the more complex TOPKAPI-ETH. We show that the performance of both models can be improved by using measured precipitation gradients of higher temporal resolution. In contrast to the good performance of the conceptual model for the present climate, however, we demonstrate that the simplifications in WEAP lead to error compensation, which results in different predictions in simulated melt and runoff for a potentially warmer future climate. TOPKAPI-ETH, using a more physical representation of processes, depends less on calibration and thus is less subject to a compensation of errors through different model components. Our results show that data obtained locally in ad hoc short-term field campaigns are needed to complement data extrapolated from long-term records for simulating changes in the water cycle of high-elevation catchments but that these data can only be efficiently used by a model applying a spatially distributed physical representation of hydrological processes.