Abstract. Historical records are an important source of information on extreme and rare floods and fundamental to establish a reliable flood return frequency. The use of long historical records for flood frequency analysis brings in the question of flood stationarity, since climatic and land-use conditions can affect the relevance of past flooding as a predictor of future flooding. In this paper, a detailed 400 yr flood record from the Tagus River in Aranjuez (central Spain) was analysed under stationary and non-stationary flood frequency approaches, to assess their contribution within hazard studies. Historical flood records in Aranjuez were obtained from documents (Proceedings of the City Council, diaries, chronicles, memoirs, etc.), epigraphic marks, and indirect historical sources and reports. The water levels associated with different floods (derived from descriptions or epigraphic marks) were computed into discharge values using a one-dimensional hydraulic model. Secular variations in flood magnitude and frequency, found to respond to climate and environmental drivers, showed a good correlation between high values of historical flood discharges and a negative mode of the North Atlantic Oscillation (NAO) index. Over the systematic gauge record , an abrupt change on flood magnitude was produced in 1957 due to constructions of three major reservoirs in the Tagus headwaters (Bolarque, Entrepeñas and Buendia) controlling 80 % of the watershed surface draining to Aranjuez. Two different models were used for the flood frequency analysis: (a) a stationary model estimating statistical distributions incorporating imprecise and categorical data based on maximum likelihood estimators, and (b) a time-varying model based on "generalized additive models for location, scale and shape" (GAMLSS) modelling, which incorporates external covariates related to climate variability (NAO index) and catchment hydrology factors (in this paper a reservoir index; RI). Flood frequency analysis using documentary data (plus gauged records) improved the estimates of the probabilities of rare floods (return intervals of 100 yr and higher). Under nonstationary modelling flood occurrence associated with an exceedance probability of 0.01 (i.e. return period of 100 yr) has changed over the last 500 yr due to decadal and multidecadal variability of the NAO. Yet, frequency analysis under stationary models was successful in providing an average discharge around which value flood quantiles estimated by non-stationary models fluctuate through time.

Colombia using regionalization and scaling concepts with the long-term flow field as the scaling variable. All data sets, methods, and results are included in HidroSIG, an available interactive hydrologic atlas of Colombia.

A long-term flood record from the Buffels River, the largest ephemeral river of NW South Africa (9250 km2), was reconstructed based on interpretation of palaeoflood, documentary and instrumental rainfall data. Palaeoflood data were obtained at three study reaches, with preserved sedimentary evidence indicating at least 25 large floods during the last 700 yr. Geochronological control for the palaeoflood record was provided by radiocarbon and optically stimulated luminescence (OSL) dating. Annual resolution was obtained since the 19th century using the overlapping documentary and instrumental records. Large floods coincided in the past within three main hydroclimatic settings: (1) periods of regular large flood occurrence (1 large flood/~30 yr) under wetter and cooler prevailing climatic conditions (AD 1600–1800), (2) decreasing occurrence of large floods (1 large flood/~100 yr) during warmer conditions (e.g., AD 1425–1600 and after 1925), and (3) periods of high frequency of large floods (~ 4–5 large floods in 20–30 yr) coinciding with wetter conditions of decadal duration, namely at AD 1390–1425, 1800–1825 and 1915–1925. These decadal-scale periods of the highest flood frequency seem to correspond in time with changes in atmospheric circulation patterns, as inferred when comparing their onset and distribution with temperature proxies in southern Africa.