One hundred years from its original definition by Fuller (), the probabilistic concept of return period is widely used in hydrology as well as in other disciplines of geosciences to give an indication on critical event rareness. This concept gains its popularity, especially in engineering practice for design and risk assessment, due to its ease of use and understanding; however, return period relies on some basic assumptions that should be satisfied for a correct application of this statistical tool. Indeed, conventional frequency analysis in hydrology is performed by assuming as necessary conditions that extreme events arise from a stationary distribution and are independent of one another. The main objective of this paper is to investigate the properties of return period when the independence condition is omitted; hence, we explore how the different definitions of return period available in literature affect results of frequency analysis for processes correlated in time. We demonstrate that, for stationary processes, the independence condition is not necessary in order to apply the classical equation of return period (i.e., the inverse of exceedance probability). On the other hand, we show that the time‐correlation structure of hydrological processes modifies the shape of the distribution function of which the return period represents the first moment. This implies that, in the context of time‐dependent processes, the return period might not represent an exhaustive measure of the probability of failure, and that its blind application could lead to misleading results. To overcome this problem, we introduce the concept of Equivalent Return Period, which controls the probability of failure still preserving the virtue of effectively communicating the event rareness.

Hydroclimatic variables such as precipitation and temperature are often measured or simulated by climate models at coarser spatiotemporal scales than those needed for operational purposes. This has motivated more than half a century of research in developing disaggregation methods that break down coarse‐scale time series into finer scales, with two primary objectives: (a) reproducing the statistical properties of the fine‐scale process and (b) preserving the original coarse‐scale data. Existing methods either preserve a limited number of statistical moments at the fine scale, which is often insufficient and can lead to an unrepresentative approximation of the actual marginal distribution, or are based on a limited number of a priori distributional assumptions, for example, lognormal. Additionally, they are not able to account for potential nonstationarity in the underlying fine‐scale process. Here we introduce a novel disaggregation method, named Disaggregation Preserving Marginals and Correlations (DiPMaC), that is able to disaggregate a coarse‐scale time series to any finer scale, while reproducing the probability distribution and the linear correlation structure of the fine‐scale process. DiPMaC is also generalized for arbitrary nonstationary scenarios to reproduce time varying marginals. Additionally, we introduce a computationally efficient algorithm, based on Bernoulli trials, to optimize the disaggregation procedure and guarantee preservation of the coarse‐scale values. We focus on temporal disaggregation and demonstrate the method by disaggregating monthly precipitation to hourly, and time series with trends (e.g., climate model projections), while we show its potential to disaggregate based on general nonstationary scenarios. The example applications demonstrate the performance and robustness of DiPMaC.

The need of understanding and modelling the space–time variability of natural processes in hydrological sciences produced a large body of literature over the last thirty years. In this context, a multifractal framework provides parsimonious models which can be applied to a wide-scale range of hydrological processes, and are based on the empirical detection of some patterns in observational data, i.e. a scale invariant mechanism repeating scale after scale. Hence, multifractal analyses heavily rely on available data series and their statistical processing. In such analyses, high order moments are often estimated and used in model identification and fitting as if they were reliable. This paper warns practitioners against the blind use in geophysical time series analyses of classical statistics, which is based upon independent samples typically following distributions of exponential type. Indeed, the study of natural processes reveals scaling behaviours in state (departure from exponential distribution tails) and in time (departure from independence), thus implying dramatic increase of bias and uncertainty in statistical estimation. Surprisingly, all these differences are commonly unaccounted for in most multifractal analyses of hydrological processes, which may result in inappropriate modelling, wrong inferences and false claims about the properties of the processes studied. Using theoretical reasoning and Monte Carlo simulations, we find that the reliability of multifractal methods that use high order moments (>3) is questionable. In particular, we suggest that, because of estimation problems, the use of moments of order higher than two should be avoided, either in justifying or fitting models. Nonetheless, in most problems the first two moments provide enough information for the most important characteristics of the distribution.

An operative methodology for rainfall thresholds definition is illustrated, in order to provide at critical river section optimal flood warnings. Threshold overcoming could produce a critical situation in river sites exposed to alluvial risk and trigger the prevention and emergency system alert. The procedure for the definition of critical rainfall threshold values is based both on the quantitative precipitation observed and the hydrological response of the basin. Thresholds values specify the precipitation amount for a given duration that generates a critical discharge in a given cross section and are estimated by hydrological modelling for several scenarios (e.g.: modifying the soil moisture conditions). Some preliminary results, in terms of reliability analysis (presence of false alarms and missed alarms, evaluated using indicators like hit rate and false alarm rate) for the case study of Mignone River are presented.