1. Introduction
Floods are major weather-related events that continue to cause high economic and human losses all over the Globe, reaching, on average, tens of billions of $US and thousands of fatalities per year. Hence, an understanding of mechanisms of temporal change and variability of floods and flood-related variables (such as precipitation, river discharge, flood loss) is of vast theoretical and practical interest and importance.
The search for ubiquitous, spatially coherent trends in long time-series of hydro-climatic observation records has been very successful for temperature, less successful for precipitation (including intense precipitation), and even less so for river flow. Actually, no persuading, spatially-organized long-term trends could be detected in records of maximum river discharge.
Analyses of the time series of annual maximum (or peak-over-threshold) values of river discharge at a given station (cross-section) and of flood indices do not show a convincing ubiquitous upward trend. In the absence of a gradual trend, the structure of the time series is dominated by strong, and seemingly rather irregular, inter-annual and inter-decadal variability. An intriguing question emerges—what is the mechanism of variability of flood-related variables? Perhaps there exist some links with recognized climate variability and its indices. Therefore, it is quite natural to examine whether the climate variability expressed via atmosphere–ocean oscillation indices can drive the variability of flood-related indices.
Better understanding of spatially- and temporally-organized links between climate variability patterns and flood hazard and flood risk at a range of spatial and temporal scales would be a considerable asset. Anticipating regional predisposition to river floods in a particular phase of oscillation in the atmosphere–ocean system, measured by a range of climate variability indices, could help in improving preparedness to emergency.
The aim of this review paper is to assemble a selective summary of the existing information on links of various climate-variability drivers with variability of characteristics of destructive water abundance at a range of scales. We discuss the modes of climatic variability that may have links to flood-related variables. We review publications reporting on links between variability of intense precipitation, high river discharge and flood indices and natural oscillations in the ocean–atmosphere system, starting from global scale to continental, river-basin, national, regional, to local scale.
2. Modes of Climate Variability
Many indices of inter-annual and inter-decadal climate variability have been proposed in the scientific literature in order to quantify oscillation in the atmosphere–ocean system. Some of them have been studied in the context of links with hydrological variables. Time series of characteristics of climate variability indices of concern may reflect volcanic eruptions and indices of solar irradiance, illustrated by sunspot numbers, as well as various classes of ocean–atmosphere oscillations. The first two indices from the above list are likely the drivers of climate variability. In contrast, ocean–atmosphere oscillations, being manifestations of climate variability, are likely the drivers of climate impact variability, that are the focus of this paper.
Accounting for the heat and mass exchange between the Earth’s atmosphere and the oceans is absolutely necessary to understand large-scale meteorological processes and to interpret the climate variability. The latter consists of quasi-periodic oscillations in the climate system, i.e., the alternating occurrence of positive and negative anomalies of meteorological variables such as temperature or air pressure in a specific reference location, or of differences between values of these variables in two specific reference locations. These spatially defined oscillations in the ocean–atmosphere system have been found to affect the climate and various climate impact variables over large areas, including remote regions, via teleconnections.
There are four principal modes of climate variability that, according to a plethora of references, are of particular relevance for links with abundance of water—i.e., with characteristics of intense precipitation, high river discharge, and floods. They are oscillations of quasi-periodicity of several years: ENSO (El Niño-Southern Oscillation) and NAO (North Atlantic Oscillation), as well as of several decades: AMO (Atlantic Multi-decadal Oscillation) and PDO (Pacific Decadal Oscillation). Brief information about these four modes of climate variability is given below.
The ENSO is the sea surface temperature oscillation over the tropical Pacific Ocean. During the warm phase, El Niño, water off the South American coast is warmer than average (in the absence of cold upwelling), while during the cold phase, La Niña, with increased cold upwelling, warm water is further west off the South American coast than usual. Between the warm phase, El Niño, and the cold phase, La Niña, there is a neutral phase. The Southern Oscillation is the accompanying atmospheric (air pressure) component, coupled with the sea surface temperature (SST) change. Many relevant indices have been proposed and some of them were used in references reviewed in this paper. Among the indices are: ONI (Oceanic Niño Index, 3-month running mean of SST anomalies in the Niño 3.4 region, stretching from the 120° W to 170° W astride the equator, five degrees of latitude on either side; for declaration of El Niño or La Niña, the anomalies must exceed ±0.5 °C for at least five consecutive months), Niño 3.4 (5-month running mean, and SSTs must exceed ±0.4 °C for a period of six months at least), Niño 1 + 2 (0–10° S, 90° W–80° W), Niño 3 (5° N–5° S, 150° W–90° W), Niño 4 (5° N–5° S, 160° E–150° W), and the Trans-Niño Index (TNI), which is defined as the difference in normalized SST anomalies between the Niño 1 + 2 and Niño 4 regions and the event is classified as a “central Pacific El Niño” or “El Niño Modoki”. Also, the Southern Oscillation Index (SOI,
https://www.ncdc.noaa.gov/teleconnections/enso/indicators/soi/) has been used, expressing fluctuations in air pressure between the western and the eastern tropical Pacific, calculated based on the differences in air pressure anomaly between Tahiti and Darwin, Australia, [
1,
2].
The term NAO refers to cyclic air pressure changes in the Icelandic Low and Azores High as well as cyclic displacements of the center of the Azores High. If the (seasonal or monthly) value of the NAO index is positive, there is a stronger air pressure contrast between the Icelandic Low and the Azores High, i.e., the pressure value in the Azores Region is higher than the long-term average, while the corresponding pressure in the Icelandic Low is lower than the long-term average. Several indices were proposed to measure the strength of NAO, related to air pressure in various stations, such as: the Hurrell index—between Lisbon (Portugal) and Stykkishólmur (Iceland), the Rogers index—between Ponta Delgada (Azores) and Akureyri (Iceland), and the Jones index—between Gibraltar and Stykkishólmur. More detailed information on the concept and studies of NAO can be found in Wanner et al. [
3].
The PDO is the decadal-scale temperature oscillation of the Pacific Ocean. Sometimes it is called a long-lived ENSO-like Pacific climate oscillation pattern. Two quasi-periodicities have been determined: of 15–25 and 50–70 years. Temperature anomalies (departures from average) in the East Pacific (off the west coast of North America) have an opposite sign to those in the West Pacific, therefore there is some compensation effect and the impact on the global temperature is small [
4,
5].
The term AMO refers to a coherent pattern of oscillatory changes in the North Atlantic sea-surface temperature (SST) with a period of about 60–90 years. Traditionally, indices of the AMO have been based on average annual SST anomalies in the North Atlantic. Temperature anomalies in the North Atlantic have an opposite sign to those in the South Atlantic, therefore there is a kind of compensation effect and the impact on the global temperature is small. The nature and origin of the AMO is uncertain, and it remains unknown whether it represents a persistent, quasi-periodic driver in the climate system, or merely a transient feature. A deterministic mechanism relying on atmosphere–ocean–sea ice interactions for the AMO was proposed by Dima and Lohmann [
6].
There are several data sets of climate variability indices that have been collected at various sources and made publicly available (e.g.,
Table 1).
3. In Search for Links between Variability of Flood-Related Variables and Climate Oscillations
The structure of the time series of annual maximum river discharge and flood damage worldwide, does not show a persuading ubiquitous upward trend [
7], even if the Clausius–Clapeyron law demonstrates that more water vapor can be stored in the warmer atmosphere, suggesting the increase of the potential for more intense precipitation (that may lead to floods) in the warmer climate. What the time series of annual maximum river discharge shows is a dominating inter-annual and inter-decadal variability that may overshadow gradual (if any) trends. High river discharges tend to be clustered in time [
8].
There is a strong multi-annual and multi-decadal variability in flood-related variables. Perhaps part of this variability is random or chaotic, but there are serious arguments supporting the thesis that climate variability plays, generally, an important role in driving flood hazard over many areas of the world. Prudhomme et al. [
9] examined whether the same few circulation types are systematically associated with floods at the European scale. Their results showed that, indeed, at the river basin scale and at the continental, European scale, some circulation types have significant positive frequency anomalies with flood occurrence, i.e., they are associated with a predisposition to floods. However, the results were not robust—they were shown to depend on various technicalities. Burn and Arnell [
10] investigated the spatial and temporal pattern of high river flow, globally, using annual maximum discharge data for 200 gauging stations. The results of the investigation revealed that diverse regions of the world exhibit synchronous flood responses which can be related to various global climatic conditions, such as atmosphere–ocean oscillations. They also examined global temporal patterns in flooding and identified distinct periods corresponding to increased global flooding activity.
Many research works have been undertaken and reported in literature, aimed at identification of links between the climate variability and various variables related to floods. Climate variability indices can quantify some climate phenomena which can be considered as drivers of variability of such flood-related characteristics, as intense precipitation, high river discharge, flood severity and magnitude (in the sense of Brakenridge, see explanation in [
11,
12]), as well as flood loss (material damage, population, number of fatalities, and a combined index). Clearly, there are several climatic and non-climatic factors influencing flood risk (understood as probability times consequences, that can be disaggregated into three driving factors: hazard, exposure, and vulnerability). Non-climatic drivers are of terrestrial/hydrological (e.g., land-use and land-cover change, river regulation, and structural flood defenses), as well as social/economic nature. For instance, changes in rivers and their drainage basins, and changes in flood damage potential depend on the socio-economy (demography, wealth, development, GDP). However, the climate variability track is likely to play an important role in the interpretation of the variability in the time series of flood-related characteristics globally and, even more so, in many regions. Hence, it seems opportune to examine variability in recorded time series of flood-related characteristics and climate variability indices and seek spatially and temporally organized links between them. Improved understanding of climate-variability impacts on variability of flood-related characteristics would be of considerable value.
The expectation of discovery of absolute deterministic rules would be naïve, yet existence of probabilistic, frequency-type rules of spatially-organized validity is likely. The global spatial scale can be used, but particular focus can be put on areas where links have been found to be strong and where robust regional laws can be formulated. In many papers, authors claim to have unveiled new teleconnections (i.e., co-variability of variables at distant locations), possibly in asynchronous (e.g., lagged) mode, and/or to have found an important original interpretation.
This paper aims at building an inventory of selected information on links between climate variability indices and variability of flood-related characteristics. It is virtually impossible to embrace all the many thousands of related references. Search for the logical sum of “ENSO” and “flood” gives 32,500 hits in Google Scholar, and 733 hits in Web of Science. Search for “NAO” and “flood” gives 24,200 hits in Google Scholar, and 147 hits in Web of Science, while search for “PDO” and “flood” gives 8,020 hits in Google Scholar, and 94 hits in Web of Science. Search for “AMO” and “flood” gives 19,900 hits in Google Scholar, and 36 hits in Web of Science (updated on 13 May 2019). The gross orientation numbers given above are clearly overestimating the count, for several reasons, e.g., in some entries, Flood is an author’s name, while Nao and Amo refer to geographic names and NAO—to National Audit Office. Nevertheless, the number of relevant source items is very high. In brief, the audience is inundated by information on the links between climate variability and variables related to floods, of varying strength and trustworthiness. Hence, some transparent rules of selection of entries to be included in this paper are needed. We give priority to reviewing papers that have gathered many citations, according to bibliometric indices (related to citations count) reported in Web of Science. However, we also include several papers that have not gathered numerous citations yet (some of them being quite recent), but deserve being considered. In result, bibliography contains 120 source items. In this review paper, there is usually reference to one or two papers authored by a single scientist or co-authored by a group of scientists, even if the same authors may have possibly published many other papers in the field. An attempt was made to select the appropriate paper(s) out of a larger set.
Meta-analysis (analysis of analyses) of results from selected studies reported in scientific literature is offered. The bibliography is collected and classified into categories (after regions, climate variability modes, and flood-related variables: precipitation, river flow, and flood losses). Reported results are assessed and a synthesis of findings is presented. It could improve understanding of the links of four variability modes with abundance of water and inform further studies in the area. Information is given on who studied what and whether approaches turned out to be successes or failures (though, definitely, the latter are less likely to be well represented in the scientific literature).
From the viewpoint of flood variability, the dominating mode of inter-annual climate oscillation is undoubtedly the El Niño-Southern Oscillation (ENSO), with impacts felt most strongly in the Pacific region but also in much of the world via teleconnections. Colloquially speaking, once the occurrence of a strong El Niño is declared, a range of extreme weather events worldwide, including floods, are blamed on this phenomenon in the media. Other, possibly relevant, climate variability indices examined in many references are: AMO (Atlantic Multi-decadal Oscillation), PDO (Pacific Decadal Oscillation), and NAO (North Atlantic Oscillation). Some papers tackle other variability indices, such as Northern Hemispheric Teleconnection Indices (Arctic Oscillation—AO, Polar/Eurasia Pattern, Siberian High-Intensity, Scandinavia Pattern, East Atlantic/Western Russia Pattern, East Atlantic Pattern, Pacific/North American Pattern—PNA, and Northern Hemispheric Circum-global Teleconnection—CGT), as well as Monsoon Indices (Webster–Yang Monsoon Index, South Asian Monsoon Index, East Asian–Western North Pacific Monsoon Index), see [
13], Indian Ocean Dipole, see [
14], and AAO (Antarctic Oscillation), see [
15].
Some studies show that the overall impact of climate variability modes (e.g., ENSO and PDO) depends on the combination of ENSO and PDO. For instance, Wang et al. [
16] investigated the inter-decadal modulation of the PDO upon the impact of ENSO on the EAWM (East Asian Winter Monsoon). No robust relationship was found between ENSO and EAWM on the inter-annual timescale, when the PDO is in its high phase. However, when the PDO is in its low phase, ENSO was found to exert a strong impact on the EAWM. Chan and Zhou [
17] found that inter-decadal variations in the early (May–June) summer monsoon rainfall over South China are linked to the ENSO and the PDO (Pacific Decadal Oscillation), with wetter monsoon years during the periods of low PDO index.
Statistical relationships have been established between some quasi-periodic ocean–atmosphere oscillations (e.g., ENSO, NAO) and general wetness (e.g., river runoff) indices in various regions of the world (cf. [
18,
19]). Some research studies [
19,
20] reviewed connections between high river flow (in winter and spring) and the NAO index in winter over many areas. This could be particularly important e.g., in regions with snowmelt as the dominating flood generation mechanism [
21].
4. Review of Links between Climate Oscillations and Flood-Related Variables—Global Studies
The El Niño-Southern Oscillation (ENSO) is perhaps the most commonly studied mode of inter-annual climate variability that was found to have an impact on river flow and flooding. Indeed, many dramatic floods have occurred in some areas during strong warm phases of ENSO (El Niño) and in other areas during strong cold phases of ENSO (La Niña). This climate variability mode largely influences the Pacific region, but its impacts can also be felt in much of the world via teleconnections, possibly of asynchronous (lagged) nature. Chiew and McMahon [
18] examined global teleconnections between ENSO and river flow. They found strong and regionally consistent ENSO-streamflow teleconnections in Australia and New Zealand, South and Central America, and weaker signals in some parts of Africa and North America. They also proposed a hypothesis that the streamflow–ENSO relationship can be stronger than the rainfall–ENSO relationship (so that an amplification effect takes place), due to integrative properties of the streamflow process.
It can be stated that the occurrence of a specific stage of climate variability (e.g., El Niño) can open the door to a visit of abundant water, yet floods may or may not come. In fact, floods can happen any time, regardless of the climate variability phase. However, in the area of interest, they can be much more likely during one phase (e.g., El Niño) and perhaps more unlikely during another phase (e.g., La Niña). Sun et al. [
22] postulated that the effect of ENSO on extreme precipitation is asymmetric, in the sense that a strong, significant effect is only for a single ENSO stage, and much weaker (or none) for another.
There is no doubt that the random component of climate variability as well as of variability of flood-related variables is strong. Nevertheless, Ward et al. [
23,
24] examined statistical links between ENSO and flooding and arrived at interesting and promising frequency-type results globally. They demonstrated that the ENSO quasi-cycle of natural variability likely has serious impacts on flood hazard in many regions of the world and called for further research to reduce uncertainty. In contrast, more recently, Emerton et al. [
25] showed that the interpretation of likelihood of ENSO-driven flood hazard is more complex than perceived and reported in literature.
Ward et al. [
23,
24] arrived at an updated global map of prevalence of wetness (intense precipitation, flooding) for particular phases of the ENSO oscillation. It is clear that various flood-related variables may react differently to climate variability indices in different areas. The behavior of the time series of river discharge may largely differ from the time series of intense precipitation, as per suggestion of integrative properties due to Chiew and McMahon [
18]. Large rivers have basins covered by multiple grid cells and the transformation of rainfall into river runoff modulates the signal and introduces a transport lag, not to forget human influence (anthropogenic effects in drainage basins and river beds) that can be substantial. However, it is not only the scale problem. In many river basins, the flow regime of the main stream differs from the regime of tributaries (e.g., pluvial versus glacial, snowmelt, or rain on snow).
There is a rich body of research aimed at connecting floods with climatic variability at the continental scale. Many relevant studies have been undertaken in Asia, North America, South America, Australia, and Europe, using ENSO, NAO, AMO, and PDO indices. Hodgkins et al. [
26] examined climate-driven variability in the occurrence of major floods across two continents (North America and Europe), only in minimally altered catchments, to focus on climate-driven changes rather than changes due to catchment alterations. Temporal changes in the occurrence of major floods were shown to be dominated by multi-decadal variability (AMO) rather than by long-term trends.
The question whether one individual flood event (or a train of floods in a flood season) can be meaningfully associated with the occurrence of a particular phase of climate variability is not properly posed. A more appropriate question is whether the occurrence of a particular phase of climate variability can cause the risk of experiencing floods (predisposition to abundance of wetness) to be significantly higher or lower than average. The reported strength of correlation may vary with time (strengthening or weakening over various periods), with possible asymmetry [
22]. Even the sign of correlations may change with time. This is the meaning of the so called “epochal behavior” recognized for the Mekong by Räsänen and Kummu [
27]. Ward et al. [
23,
24] speculated that climate change may possibly be modifying the strengths of teleconnections.
Table 2 presents the essential information extracted from selected references, about the links between climate variability and floods, resulting from global studies.
5. Review of Links between Climate Oscillations and Flood-Related Variables—Continental, Regional, and Local Studies
Table 3,
Table 4,
Table 5,
Table 6,
Table 7 and
Table 8 present the most essential information extracted about the links between climate variability and floods, resulting from continental (and bi-continental) as well as regional or local studies, as reported in scientific literature. The tables are organized according to continental coverage of studies.
Table 3 refers to Asia,
Table 4—to North America,
Table 5—to South America,
Table 6—to Australia,
Table 7—to Europe, and
Table 8—to Africa. Large regions within continents are placed first, then large river basins.
A few studies examine more than one mode of climate variability. This is explicitly noted in the tables. If a study refers to more than one continent, but does not cover the Globe, it is mentioned in several tables.
Figure 1 illustrates regional predisposition to abundance of water for various phases of four modes of climatic variability. The map in
Figure 1 summarizes the findings of continental, regional, to local studies reported in
Table 3,
Table 4,
Table 5,
Table 6,
Table 7 and
Table 8 of this section.
6. Concluding Remarks
No persuading, and spatially-coherent, trend could be detected in long time series of flood-related variables, such as high river discharge, flood severity and magnitude, as well as flood loss. The structure of time series of river discharge is dominated by strong, and seemingly rather irregular, inter-annual and inter-decadal variability. Perhaps part of this variability can be random or chaotic, but there are serious arguments supporting the thesis that climate variability (quasi-periodic oscillations in the atmosphere–ocean system) plays, generally, an important role in driving variables related to floods over many areas of the world. Therefore, it is quite natural to examine whether there are links between the climate variability expressed via atmosphere–ocean oscillation indices and the variability of flood-related indices.
Spatially-defined oscillations in the ocean–atmosphere system have been found to affect the climate and various climate impacts over large areas, covering the oceans and adjacent continents, but also more remote (sometimes—very remote) regions, via so-called teleconnections. The strength of the link can vary from region to region and one can identify zones of influence where a climate variability mode (e.g., ENSO, PDO, NAO, or AMO) affects variables related to floods.
Many research works have been undertaken, with results reported in the literature, aimed at identification of links between the climate variability and various variables related to floods. Studying variability of properties of high river flows (e.g., timing, intensity, volume, duration, frequency) and their links with climate variability is a cutting-edge research activity, attracting broad attention due to scientific interests and practical importance. Anticipating regional predisposition to river floods in a particular phase of oscillation in the atmosphere–ocean system, measured by a range of climate variability indices, could help in improving the system of early preparedness to emergency.
The present paper reviewed published results of studies indicating links between climate variability indices and characteristics of abundant wetness (precipitation, river runoff, and flood indices) in different areas. Many studies have been found that report on the existence of such links, at a range of scales—from global (for ENSO and PDO) to local. Some references indicate the necessity to consider a combination of oscillation modes, e.g., one oscillation mode can modulate another one. Some studies combine analysis of change and variability, or seek for change points in system behavior. The list of references considered in this paper consists of 120 source items. The continents ordered after the number of source items, from highest to lowest, are: North America, Asia, Europe, Africa, South America, and Australia. Some reports on the existence of links are of a rather speculative nature, as a signal is sought among strong noise. Some findings are too complex and do not fit into the simple story of this paper.
As shown in
Figure 1, the reviewed references found predisposition to abundance of water during a particular phase of one or more modes of climate variability. More consistent results are reported for South America and El Niño signals wetness there. In Africa, also ENSO is found important, while predisposition to abundance of water is related to La Niña for most areas and El Niño for some. In Europe, it is mostly NAO+ in the Northern and Western parts and NAO- in the Southern part. In Asia, links with ENSO (El Niño for some areas and La Niña for other areas) and PDO were reported. In Australia, it is La Niña and PDO-.
The situation in North and Central America looks most complex (
Figure 1), because all modes: ENSO, NAO, AMO, and PDO may play a role. However, Jiang et al. [
118] showed that no single index of those four, derived from the Pacific and Atlantic oceans, can explain the multi-scale temporal variability and spatial distribution of heavy precipitation in the western USA. Yet, characterization of their integrated effect, including modulation and reflecting the combined physical oceanic-atmospheric processes, allows useful links to be deciphered.
Reported links have spatial, temporal, and seasonal validity. As expected, different source items may give differing representation of links between climate variability and variables related to floods, for instance due to different definitions of indices, or due to different data windows (start-year and end-year) being considered. In brief, there are still considerable uncertainties in the understanding of the climate variability modes and their impact mechanisms on flood hazard and flood risk. This adds to the recognized general uncertainty of climate impacts on water resources [
119].
Reviewed references follow various methodologies to investigate the relationship between climate oscillations and flood-related variables. Correlation and cross-correlation techniques (including time lags) are commonly used. Many papers use principal component or empirical orthogonal function (EOF) analysis (e.g., [
107,
118]), as well as wavelet power spectrum analysis (e.g., [
107]).
Some publications report on convincing physical explanations of mechanisms behind the relationships. One of examples is in the study for North America, by Gershunov and Barnett [
120], who found that the North Pacific oscillation (NPO) modulates the consistency and strength of El Niño anomalies. They hypothesized that a deeper Aleutian low (manifestation of the NPO) shifts the storm track southward while El Niño provides an enhanced eastern tropical Pacific moisture source for the storms to tap. They also demonstrated that seasonal climate anomalies over North America exhibit rather large variability between years characterized by the same ENSO phase. Interdecadal climatic anomalies in the North Pacific (the North Pacific oscillation, NPO) were shown to exert a modulating effect on ENSO teleconnections. Typical El Niño patterns (e.g., low pressure over the northeastern Pacific, dry northwest, and wet southwest, etc.) are strong and consistent only during the high phase of the NPO, which is associated with an anomalously cold northwestern Pacific.
It is hoped that the present review paper could play an important and useful role, if only because of the systematization of scattered knowledge. It could improve understanding of the links of four climate variability modes with abundance of water and inform further studies in the area.