The HBV.IANIGLA Hydrological Model

Over the past 40 years, the HBV (Hydrologiska Byrans Vattenbalansavdelning) hydrological model has been one of the most used worldwide due to its robustness, simplicity, and reliable results. Despite these advantages, the available versions impose some limitations for research studies in mountain watersheds dominated by ice-snow melt runoff (i.e., no glacier module, a limited number of elevation bands, among other constraints). Here we present HBV.IANIGLA, a tool for hydroclimatic studies in regions with steep topography and/or cryospheric processes which provides a modular and extended implementation of the HBV model as an R package. To our knowledge, this is the first modular version of the original HBV model. This feature can be very useful for teaching hydrological modeling, as it offers the possibility to build a customized, open-source model that can be adjusted to different requirements of students and users.Fil: Toum, Jorge Ezequiel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; ArgentinaFil: Masiokas, Mariano Hugo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; ArgentinaFil: Villalba, Ricardo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; ArgentinaFil: Pitte, Pierre. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; ArgentinaFil: Ruiz, Lucas Ernesto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; Argentin

Hydrology modelling R packages: a unified analysis of models and practicalities from a user perspective

Following the rise of R as a scientific programming language, the increasing requirement for more transferable research, and the growth of data availability in hydrology, R packages containing hydrological models are becoming more and more available to hydrologists. Corresponding to the core of the hydrological studies workflow, their value is increasingly meaningful regarding the reliability of methods and results. Despite package and model distinctiveness, no study has ever 5 provided a comparison of R packages for conceptual rainfall-runoff modelling from a user perspective, contrasting their philosophy, model characteristics and ease of use. We have selected eight packages based on our ability to consistently run their models on simple hydrology modelling examples. We have uniformly analysed the exact structure of seven of the hydrological models integrated in these R packages in terms of conceptual storages and fluxes, spatial discretisation, data requirements and output provided. The analysis showed that very different modelling choices are associated with these packages, which emphasises various hydrological concepts. These specificities are not always sufficiently well explained by the package documentation. Therefore a synthesis of the package functionalities was performed from a user perspective. This synthesis helps inform selection of what packages could/should be used depending on the problem at hand. In this regard, technical features, documentation, R implementations and computational times were investigated. Moreover, by providing a framework for package comparison, this study is a step forward towards supporting more transferable and reusable methods and results for hydrological modelling in R

Benchmarking ensemble streamflow prediction skill in the UK

Skilful hydrological forecasts at sub-seasonal to seasonal lead times would be extremely beneficial for decision-making in water resources management, hydropower operations, and agriculture, especially during drought conditions. Ensemble streamflow prediction (ESP) is a well-established method for generating an ensemble of streamflow forecasts in the absence of skilful future meteorological predictions, instead using initial hydrologic conditions (IHCs), such as soil moisture, groundwater, and snow, as the source of skill. We benchmark when and where the ESP method is skilful across a diverse sample of 314 catchments in the UK and explore the relationship between catchment storage and ESP skill. The GR4J hydrological model was forced with historic climate sequences to produce a 51-member ensemble of streamflow hindcasts. We evaluated forecast skill seamlessly from lead times of 1 day to 12 months initialized at the first of each month over a 50-year hindcast period from 1965 to 2015. Results showed ESP was skilful against a climatology benchmark forecast in the majority of catchments across all lead times up to a year ahead, but the degree of skill was strongly conditional on lead time, forecast initialization month, and individual catchment location and storage properties. UK-wide mean ESP skill decayed exponentially as a function of lead time with continuous ranked probability skill scores across the year of 0.75, 0.20, and 0.11 for 1-day, 1-month, and 3-month lead times, respectively. However, skill was not uniform across all initialization months. For lead times up to 1 month, ESP skill was higher than average when initialized in summer and lower in winter months, whereas for longer seasonal and annual lead times skill was higher when initialized in autumn and winter months and lowest in spring. ESP was most skilful in the south and east of the UK, where slower responding catchments with higher soil moisture and groundwater storage are mainly located; correlation between catchment base flow index (BFI) and ESP skill was very strong (Spearman's rank correlation coefficient = 0.90 at 1-month lead time). This was in contrast to the more highly responsive catchments in the north and west which were generally not skilful at seasonal lead times. Overall, this work provides scientific justification for when and where use of such a relatively simple forecasting approach is appropriate in the UK. This study, furthermore, creates a low cost benchmark against which potential skill improvements from more sophisticated hydro-meteorological ensemble prediction systems can be judged

Characterising urban catchments for explaining storm runoff and application in UK flood estimation

The impacts of urbanisation on catchment hydrology have been the focus of investigation over the last few decades, but quantifying and predicting the impacts remains an ongoing area of active research. One such area has been improving characterisation of urban land cover to predict urbanisation impacts whereby lumped catchment characterisation of urban land cover limits the ability of attribution and modelling methods to consider the spatial role of land cover in runoff response. This thesis evaluates the potential for spatially explicit characterisations of urban land cover based on landscape metrics, commonly employed in landscape ecology, to explain storm runoff in urban catchments and their application in UK flood estimation methods. Rainfall and channel flow monitoring across two towns containing 18 variably urbanised sub-catchments were used to provide high-resolution time-series of rainfall and runoff and to identify storm events which were quantified using a range of hydrological metrics. Analysing storm runoff along a rural-urban gradient showed a lumped measure of urban extent can generally explain differences in the hydrological response between rural and urban catchments but not between more urbanised catchments in which soil moisture does not play a contributing role. Using high resolution geospatial data can improve the representation of the urban environment and landscape metrics can better represent the form and function of urban land cover, improving estimates of the index flood QMED over lumped catchment descriptors. Regression analysis of hydrological metrics showed the potential of landscape metrics for explaining inter-catchment differences in rainfall-runoff and point to the importance of considering the location and connectivity of urban surfaces. Landscape metrics provide a workable means of overcoming the limitations inherent in using lumped characterisation of complex urban land cover and their ability to express connectivity, size and location of urban land cover promises potential applications in hydrological applications such as UK design flood estimation methods

How uncertain are precipitation and peak flow estimates for the July 2021 flooding event?

The disastrous July 2021 flooding event made us question the ability of current hydrometeorological tools in providing timely and reliable flood forecasts for unprecedented events. This is an urgent concern since extreme events are increasing due to global warming, and existing methods are usually limited to more frequently observed events with the usual flood generation processes. For the July 2021 event, we simulated the hourly streamflows of seven catchments located in western Germany by combining seven partly polarimetric, radar-based quantitative precipitation estimates (QPEs) with two hydrological models: a conceptual lumped model (GR4H) and a physically based, 3D distributed model (ParFlowCLM). GR4H parameters were calibrated with an emphasis on high flows using historical discharge observations, whereas ParFlowCLM parameters were estimated based on landscape and soil properties. The key results are as follows. (1) With no correction of the vertical profiles of radar variables, radar-based QPE products underestimated the total precipitation depth relative to rain gauges due to intense collision–coalescence processes near the surface, i.e., below the height levels monitored by the radars. (2) Correcting the vertical profiles of radar variables led to substantial improvements. (3) The probability of exceeding the highest measured peak flow before July 2021 was highly impacted by the QPE product, and this impact depended on the catchment for both models. (4) The estimation of model parameters had a larger impact than the choice of QPE product, but simulated peak flows of ParFlowCLM agreed with those of GR4H for five of the seven catchments. This study highlights the need for the correction of vertical profiles of reflectivity and other polarimetric variables near the surface to improve radar-based QPEs for extreme flooding events. It also underlines the large uncertainty in peak flow estimates due to model parameter estimation.</p

Conceptual hydrological model calibration using multi-objective optimization techniques over the transboundary Komadugu-Yobe basin, Lake Chad Area, West Africa

Study Area: The discharge of the transboundary Komadugu-Yobe Basin, Lake Chad Area, West Africa is calibrated using multi-objective optimization techniques. Study focus: The GR5J hydrological model parameters are calibrated using six optimization methods i.e. Local Optimization-Multi Start (LOMS), the Differential Evolution (DE), the Multiobjective Particle the Swarm Optimization (MPSO), the Memetic Algorithm with Local Search Chains (MALS), the Shuffled Complex Evolution-Rosenbrock’s function (SCE-R), and the Bayesian Markov Chain Monte Carlo (MCMC) approach. Three combined objective functions i.e. Root Mean Square Error, Nash- Sutcliffe efficiency, Kling-Gupta efficiency are applied. The calibration process is divided into two separate episodes (1974–2000 and 1980–1995) so as to ascertain the robustness of the calibration approaches. Runoff simulation results are analysed with a timefrequency wavelet transform. New hydrological insights for the region: For calibration and validation stages, all optimization methods simulate the base flow and high flow spells with a satisfactory level of accuracy. For calibration period, MCMC underestimate it by -0.07 mm/day. The performance evaluation shows that MCMC has the highest values of mean absolute error (0.28) and mean square error (0.40) while LOMS and MCMC record a low volumetric efficiency of 0.56. In all cases, the DE and the SCE-R methods perform better than others. The combination of multi-objective functions and multi-optimization techniques improve the model’s parameters stability and the algorithms’ optimization to represent the runoff in the basin