Multi-year drought frequency analysis at multiple sites by operational hydrology - A comparison of methods

This paper compares two generators of yearly water availabilities from sources located at multiple sites with regard to their ability to reproduce the characteristics of historical critical periods and to provide reliable results in terms of the return period of critical sequences of different length. The two models are a novel multi-site Markov mixture model explicitly accounting for drought occurrences and a multivariate ARMA. In the case of the multisite Markov mixture model parameter estimation is limited to a search in the parameter space guided by the value of parameter k to show the sensitivity of the model to this parameter. Application to two of the longest time series of streamflows available in Sicily (Italy) shows that the models can provide quite different results in terms of estimated return periods of historic droughts, although they seem to perform more uniformly when it comes to simulate drought-related statistics such as drought length, severity and intensity. The role of parameter selection for the multisite Markov mixture model and of the marginal probability of generated flows in providing results consistent with the characteristics of the observed series is discussed. Both models are applied to the system of sources supplying the city of Palermo (Sicily) and its environs showing the applicability of the newly developed multisite Markov mixture model to medium-to-large scale water resources systems.