The interaction between physical drivers from oceanographic, hydrological, and meteorological processes in coastal areas can result in compound flooding. Compound flood events, like Cyclone Idai and Hurricane Harvey, have revealed the devastating consequences of the co-occurrence of coastal and river floods. A number of studies have recently investigated the likelihood of compound flooding at the continental scale based on simulated variables of flood drivers, such as storm surge, precipitation, and river discharges. At the global scale, this has only been performed based on observations, thereby excluding a large extent of the global coastline. The purpose of this study is to fill this gap and identify regions with a high compound flooding potential from river discharge and storm surge extremes in river mouths globally. To do so, we use daily time series of river discharge and storm surge from state-of-the-art global models driven with consistent meteorological forcing from reanalysis datasets. We measure the compound flood potential by analysing both variables with respect to their timing, joint statistical dependence, and joint return period. Our analysis indicates many regions that deviate from statistical independence and could not be identified in previous global studies based on observations alone, such as Madagascar, northern Morocco, Vietnam, and Taiwan. We report possible causal mechanisms for the observed spatial patterns based on existing literature. Finally, we provide preliminary insights on the implications of the bivariate dependence behaviour on the flood hazard characterisation using Madagascar as a case study. Our global and local analyses show that the dependence structure between flood drivers can be complex and can significantly impact the joint probability of discharge and storm surge extremes. These emphasise the need to refine global flood risk assessments and emergency planning to account for these potential interactions.
Flooding in deltas and estuaries is driven by the interactions of oceanographic, hydrological, and meteorological phenomena such as extreme
rainfall, river discharge, storm surge, and wave action. When these co-occur in space and time, they can exacerbate the flood extent, depth, and duration
locally, resulting in a so-called compound flood event (Zscheischler et al., 2018). These events have the potential to cause large social and economic impacts and can directly or indirectly impact flood emergency response and infrastructure (Leonard et al., 2014; Zscheischler et al., 2018). The potentially extreme impacts of compound flood events were
recently highlighted by Cyclone Idai in March 2019. The long-lived tropical
cyclone travelled within the Mozambique Channel, causing catastrophic damage
along its path in Madagascar, Malawi, Zimbabwe, and, most importantly, in
Mozambique, where it made landfall twice. Its second landfall near Beira on
14 March devastated the city and surrounding villages. The combination of
extreme winds of more than 160 km h
Classical methodologies for large-scale flood hazard and flood risk studies mainly consider univariate flooding mechanisms and do not include compound flood events (Ward et al., 2015). These assessments therefore focus on either river (e.g. Alfieri et al., 2014; Dottori et al., 2016; Hirabayashi et al., 2013; Ward et al., 2013, 2017; Winsemius et al., 2013, 2016) or coastal floods (e.g. Brown et al., 2016; Hinkel et al., 2014; Muis et al., 2016; Vousdoukas et al., 2018), neglecting riverine and sea level interactions. Yet, these interactions can significantly influence the magnitude of simulated water levels (Santiago-Collazo et al., 2019). At the global scale, Ikeuchi et al. (2017) found the annual maximum river water level for 2005 in low-lying flat areas to be underestimated by more than 0.5 m when ignoring sea level interactions. Local studies have shown that ignoring the dependence between river discharge and storm surge can underestimate the return period of a given water level within a river mouth (Bevacqua et al., 2017; Couasnon et al., 2018; Moftakhari et al., 2019; Serafin et al., 2019).
Compound flood events can occur due to synoptic weather systems (Seneviratne et al., 2012). Clearly, tropical cyclones have the potential to cause simultaneous high river discharge and storm surge, as exemplified by Cyclone Idai. Storms with prevailing wind directions hitting mountains have also been documented to generate strong sustained winds accompanied with intense rainfall, due to orographic effects (Martius et al., 2016; Svensson and Jones, 2002, 2004). However, the co-occurrence of coastal and river floods can also occur by chance and not be related to an underlying common synoptic weather system. The expected number of co-occurrences happening by chance (i.e. statistical independence) can be determined based on probability theory (Kew et al., 2013; Martius et al., 2016). The impact of a compound flood event is influenced by the magnitude of the river and coastal flood drivers. The presence of a positive and significant statistical dependence between flood drivers indicates a higher probability for the occurrence of extreme combinations of these variables when compared to statistical independence (Diermanse and Geerse, 2012).
A consistent mathematical definition of compound flood events does not exist and multiple statistical methods have been suggested to study this phenomenon (Hao et al., 2018; Tilloy et al., 2019). These methods usually examine the number of joint extremes or the statistical dependence between proxy variables of different flood hazard types such as rainfall and storm surge, river flow and storm surge, and river flow and sea level (Bevacqua et al., 2019; Hendry et al., 2019; Kew et al., 2013; Paprotny et al., 2018a; Sadegh et al., 2018; Svensson and Jones, 2002, 2004; Wahl et al., 2015; Ward et al., 2018; Wu et al., 2018; Zheng et al., 2013). Recent compound flooding studies carried out at the regional to global scale used copula theory to characterise the bivariate joint distribution and assess complex dependence structures, e.g. in the case of upper tail dependence (Bevacqua et al., 2019; Paprotny et al., 2018a; Ward et al., 2018). Possible compound flooding mechanisms are examined by sampling a set of events from the full bivariate time series and then analysing the dependence structure of the latter. Wahl et al. (2015), Moftakhari et al. (2017), and Ward et al. (2018) used conditional sampling to assess the bivariate relationship between a riverine flood driver and a coastal flood driver when one variable was in an extreme state (e.g. by selecting annual maxima or peaks over threshold). Other studies defined compound flood events as pairs based on joint exceedances above a predefined quantile such as the 95th or 97.5th percentile of the respective marginal distribution (Bevacqua et al., 2019; Hawkes, 2008; Kew et al., 2013). However, directly applying such approaches for flood hazard quantification can be difficult due to the challenge of both defining independent and identically distributed events and capturing extremes from both time series (Hawkes, 2008; Hawkes et al., 2008).
A statistically robust analysis of bivariate flood drivers requires an extensive set of high-quality observations. Studies based on observations from gauge data have therefore provided an overview of the compound flood potential globally but are strongly biased towards gauge-rich areas. In the case of Ward et al. (2018), this resulted in a selection of 187 pairs of stations located mainly around the coasts of North America, Europe, Australia, and Japan. Non-stationarities in the observations may be present due to anthropogenic activities, such as water extractions, dam construction, and land-use changes. These factors increase the complexity of the signal and make the attribution of the dependence to synoptic meteorological drivers challenging. One way to address these limitations is by using hydrodynamic models to simulate river discharge and storm surge, and using these simulated time series for the statistical analysis of compound flood potential. Such an approach has been carried out for the European (Bevacqua et al., 2019; Paprotny et al., 2018a) and the Australian coastlines (Wu et al., 2018).
In this paper, we identify potential hotspot regions for compound floods from riverine and coastal floods along the entire global coastline by taking advantage of the extensive spatial and temporal coverage from a global river discharge and a global storm surge model. In doing so, we provide a first statistical assessment of the compound flood potential in areas where observations from discharge and tide gauges are absent or insufficient. We do not limit our analysis to one specific statistical approach but purposefully examine the compound flood potential by analysing both the timing between river discharge and storm surge extremes, and their dependence. We further suggest and apply a new methodology to quantify compound flood hazard that integrates these characteristics while fully capturing both extreme marginal distributions. Finally, we exemplify the critical influence of the dependence structure on the probability of compound discharge and coastal flood events by means of a case study example in Madagascar. Therefore, our global analysis should be considered a first step towards statistically characterising compound flooding from extreme river discharge and storm surge worldwide.
This paper is divided in four parts, as follows. Section 2 introduces the global datasets used and the method developed for this study. Section 3 presents the results and discusses the observed spatial patterns of high (low) compound flood potential based on previous literature. We also emphasise the implication of compound flood events for flood hazard quantification by looking at a selected location in Madagascar and the limitations of our study. The conclusions and outlook for future research are presented in Sect. 4.
We assess the compound flood potential between riverine and coastal flood
drivers using simulated daily river discharge and maximum daily storm surge
as proxy variables, respectively. The latter is a common choice for studying
compound flood hazard analysis in deltas and estuaries (Khanal
et al., 2019; Klerk et al., 2015; Svensson and Jones, 2002; Ward et al., 2018). The research involves the following steps, each of which is described in the following subsections:
selecting global datasets of river discharge and storm surge time series, defining sets of events to analyse compound flooding, quantifying compound flood potential using the defined sets.
We use simulations of instantaneous daily discharge of the river-routing model CaMa-Flood
version 3.62 (Yamazaki et al., 2014) forced by daily average runoff data of the land surface model JULES WRR2 eartH2Observe (Best et al., 2011; Clark et al., 2011; Dutra et al., 2017; Schellekens et al., 2017) available at
Daily river discharge is obtained by routing the mean daily runoff of the
JULES model from the eartH2Observe WRR2 reanalysis data at 0.5
Storm surge, the change in sea level driven by high winds and low atmospheric pressure, is simulated in GTSM with wind speed and atmospheric pressure from ERA-Interim (Muis et al., 2016; Verlaan et al., 2015). The model is a global hydrodynamic model using an unstructured grid with a higher resolution on the shallow continental shelf (up to 3 arcmin) than in deeper parts of the oceans (0.5
We carry out additional validation for both models extending the validation performed by Beck et al. (2017b) and Muis et al. (2016) by looking more specifically at the timing and correlation of discharge and storm surge extremes over a time period of at least 20 years between 1980 and 2014; see Sect. S1 in the Supplement. We calculate the percentage of annual maxima dates correctly predicted, termed “hit rate”, and the Spearman's rank correlation coefficient between observed and simulated annual maxima. For the discharge, we find a relatively high rank correlation globally (median: 0.57; SD: 0.22). Capturing the timing of extreme river discharge is more challenging (median hit rate: 0.21; SD: 0.18), but the hit rate increases close to the coast (see Fig. S1). For the storm surge, we find a higher hit rate (median: 0.34; SD: 0.22) and a lower rank correlation coefficient (median: 0.37; SD: 0.31) than for the discharge. In this case, coastal stations with a high correlation coefficient also capture the timing of storm surge extremes well. As a result, the timing and correlation of extreme storm surge is generally well represented along the European, North American, Japanese, and Australian coasts but not the South African and South American coasts. We further assess how the respective performance of both models can affect the compound flood potential measures defined in Sect. 2.3 (see Sect. S1.3). Due to the selection criteria we impose, this results in 25 pairs of observation stations, which is insufficient for a rigorous comparison at large regional scales. In general, we find that using the simulated discharge and storm surge variables captures the sign of the dependence correctly, but the magnitude of this dependence can largely vary. Moreover, the models tend to correctly identify the locations with the highest number of co-occurring discharge and storm surge annual maxima but overestimate this number. This additional validation showed that the performance of both models vary globally, which as a result can locally bias the compound flood potential. Nevertheless, it provides an acceptable performance on average for the purpose of this study, i.e. to provide a first-cut assessment of the compound flood potential at the global scale.
Finally, each discharge location at the river mouth of coastal catchments
larger than 1000 km
We do not restrict our analysis to one specific set of extreme river
discharge and storm surge events per location but instead define different
sets of events from the paired time series in order to measure the compound
flood potential (presented in Sect. 2.3). In this subsection, we explain
the differences between the sets and illustrate them for an example location
along the coast of Madagascar (Fig. 1c). Figure 1a and b present the paired
time series of simulated daily discharge
To investigate the strength of the dependence between the two variables, we
select the conditional sampling method used in Wahl et al. (2015) and Ward
et al. (2018). We create two sets of events based on the conditional
sampling of the annual maxima of the river discharge
Conversely, the other set is created as follows:
We also examine the co-occurrence of annual maxima by defining another set
of events, the annual maxima pairs of river discharge and of storm surge:
In the extreme case where annual maxima of discharge and surge are always
co-occurring, this means that all the sets of events defined above are
equivalent. In other words,
Maximum daily storm surge
We use the different sets of events constructed from the marginal extremes as defined in Sect. 2.2 to measure the compound flood potential at all the paired locations in three ways. First, we calculate the conditional dependence strength between river discharge and storm surge. Second, we calculate the total number of co-occurring annual maxima from the simulation period globally and analyse the probability of obtaining such a result if discharge and surge were independent. Finally, we calculate the compound flood hazard that corresponds to the probability of observing a co-occurring discharge and storm surge annual maxima above a certain magnitude in a given year. These three approaches are described in the three subsections below.
We characterise the interactions between river discharge and storm surge by
calculating the rank correlation coefficient for the two sets of pairs
(
We analyse the number of co-occurrences of annual maxima of river discharge
and storm surge along the global coastline. To do so, we count the total
number
Probability mass function for the probability of observing
As shown in Fig. 2a, one can expect a higher chance of observing two or more
co-occurrences within the 35 years if both annual maxima are consistently
occurring within the same 90 d season (red curve) compared to if they
occur randomly throughout the whole year (blue curve). This finding can be
summarised as the exceedance probability of obtaining
Finally, we examine the probability
Given the limited temporal coverage of the data, we present the result for a
quantile threshold equivalent to a 5-year discharge magnitude and a 5-year
storm surge magnitude (i.e.
In this section, we present the results for each compound flood potential measure introduced in Sect. 2.3 along the global coastline. We compare the results with respect to existing literature on compound flooding. Relevant meteorological processes likely to lead to the observed regional patterns of high compound flood potential are also discussed.
Figure 3 presents the Spearman's correlation coefficient
We observe clear regional patterns of positive dependence globally. These dependence behaviour patterns are similar to those found by Ward et al. (2018) using observations from river and tide gauges. We obtain more locations for southwestern Japan that exhibit statistically significant dependence when conditional on extreme storm surge (Fig. 3b) than when conditional on extreme discharge (Fig. 3a). We also find positive and statistically significant dependence for locations both on the West and East coasts of the United States (US). However, our results also highlight regions that could not be examined by Ward et al. (2018) due to a lack of gauge observations. Along the South American coastline, we find a cluster of positive dependence along the south coast of Chile. Along the African coastline, the coast of Madagascar is consistently highlighted in both cases, as well as the coast of Morocco. Finally, the coasts of India and large parts of East Asia also show large regions with positive dependence.
Other regions, such as the East Coast of the US, Italy, the United Kingdom (UK) or China, show a more complex dependence behaviour. Riverine flooding in these regions is related to multiple mechanisms, not all of which are related to the mechanisms causing high storm surge. On the East Coast of the US, even though tropical cyclone activity is known to contribute to high storm surge levels and intense precipitation (Villarini and Smith, 2010; Wahl et al., 2015), other river flood generating mechanisms also play a role. For example, high river discharge could also be due to snowmelt or convective storms happening upstream in the catchment (Berghuijs et al., 2016). Similarly, for the eastern coast of the UK, Hendry et al. (2019) found that storms that generate high river discharge are different to the ones that generate high storm surge. This explains why the presence of a statistical significance for discharge conditional on extreme storm surge (Fig. 3a) may be absent when conditionally sampling on extreme discharge (Fig. 3b).
Locations in several regional clusters exhibit a positive statistical dependence in both cases (6 % of all the locations studied) and therefore present the highest potential for compound flooding. Among others, we note the coasts of Madagascar, Portugal, northern Morocco, northern Australia, Vietnam, and Taiwan, which all consistently show a positive and significant dependence between discharge and storm surge.
Spearman's
We examine the total number of co-occurrences between discharge and storm
surge annual maxima obtained from the simulation period and assess their
probability of co-occurrence under the assumption of statistical independence. Figure 4 presents the total
number of co-occurring annual maxima, (
We test the significance of these results globally by comparing the
empirical distribution obtained from the data with the binomial distributions shown in Fig. 2. Figure 5 presents the uncertainty bounds around the empirical distribution using non-parametric bootstrapping with 5000 bootstrap samples and a significance level
Locations where results deviate the most from independence coincide with areas of strong tropical or extratropical cyclone activity. Martius et al. (2016) found that coastal regions affected by frequent tropical cyclones experience the highest number of co-occurring wind and precipitation extremes. Their study highlighted very similar regional patterns compared to the ones presented in Fig. 6 but with less geographical spread. This can be attributed to the fact that they focused on the analysis of climate extremes, whereas we use river discharge and storm surge where these meteorological phenomena are propagated through model chains. In other regions, they identify the interaction of weather systems with regional orographic features that cause compound wind and precipitation extremes. Atmospheric rivers landing on the West Coast of the US have caused recurrent major flood events (Gimeno et al., 2014). Composite analyses of these systems show that they can be accompanied by extreme skew surge (Ridder et al., 2018; Ward et al., 2018). The Iberian Peninsula and the Atlas Mountains contain major orographic features that can block prevailing wind flows and trigger orographic rainfall during low pressure systems, thereby causing high river discharge. These synoptic weather systems were documented to have caused serious flood events in Portugal during windstorms Klaus in 2009, Xynthia in 2010, and Gong in 2013 (Liberato, 2014).
Number of co-occurring annual maxima of discharge and storm surge obtained between 1980 and 2014 using a time window of 3 d. Major rivers are shown in grey.
Comparison of the empirical distribution of the number of co-occurrences from the simulations (in black) with the binomial distributions, assuming that annual maxima can happen randomly within the year (365 d, blue curve) or within the same season (90 d, red curve). Probability mass function
We use Eq. (4) to calculate the probability
Focusing on Europe (inset in Fig. 6), we find the highest compound flood
hazard potential mainly on the western coast, more specifically along the
coasts of Portugal and Ireland, the western coast of the United Kingdom, and the
Straits of Gibraltar. This regional pattern is also observed in the studies
of Bevacqua et al. (2019) and Paprotny et al. (2018a). However, contrary to Bevacqua et al. (2019), we do not find a high compound flood hazard
potential for the coast of France. We attribute this difference to the fact
that we focus on river discharge, whereas their study examined rainfall. In
this area, a cross-correlation analysis on the results from this study (not
shown here) shows a lag between
Probability of a co-occurring annual maxima in a given year
Regions with a high compound flood hazard potential identified in Fig. 6
coincide with regions with a high number of co-occurring annual maxima (Fig. 4) and a strong positive statistical dependence (Fig. 3). This is to be
expected since co-occurring annual maxima (
We assess the sensitivity of the joint return period shown in Fig. 6 to
the selected dependence model by selecting on a location in a region with a
high compound flood hazard potential in Madagascar. The selected paired
location is shown in Fig. 1. Figure 7 presents the probability
Effect of the dependence structure on the joint return period of a discharge and storm surge annual maxima co-occurring in any given year with a magnitude higher than a selected threshold. The latter corresponds to the marginal return period of discharge and of storm surge levels.
While we identify potential compound flooding hotspot regions from extreme discharge and storm surge variables, substantial uncertainties remain as to how this phenomenon will propagate into inland flooding. Flood events in a coastal catchment typically result from the interactions of other drivers not considered in this study, such as local rainfall, wave effects, tidal amplitude, and surge–tide interactions (Saleh et al., 2017; Vousdoukas et al., 2016). Moreover, local characteristics like the bathymetry, catchment properties, and the presence of water control structures further influence the extent to which these interactions contribute to water level extremes at the considered river mouths (Gori et al., 2019; Veldkamp et al., 2018). Specific compound flood event modelling studies (Bilskie and Hagen, 2018; Kumbier et al., 2018) and comprehensive probabilistic simulations at local scales show that multivariate flood drivers result in highly non-linear responses of flood impact variables such as flood depth and flood extent (Couasnon et al., 2018; Serafin et al., 2019). Future assessments of compound flood hazard at the global scale should therefore focus on incorporating those multivariate processes and is left for future work.
In this study, we base our statistical analysis on annual maxima, which
results in 35 data points per paired location. Alternative strategies, such
as a peaks over threshold approach or sampling the
Moreover, the results presented in this study are dependent on the accuracy of the models. The validation performed in this study indicates an acceptable performance on average, albeit with large spatial differences. The timing of the simulated storm surge compared with observations shows a poor performance of the model for stations along the coasts of South America and Africa. Yet, the accuracy of the models in measuring the compound flood potential along large parts of the coastline in South America, Africa, and Asia could not be assessed due to a lack of long-term gauge observations of discharge and sea levels. This was the main motivation for examining different compound flood measures and selecting moderate joint return period conditions, in order to identify regions potentially most exposed to this phenomenon. Nevertheless, important processes for compound flood events may be underestimated or absent in the global models used for this study. Small-scale convective and short-lived processes affecting both wind and precipitation extremes are not fully represented in the weather forcing but may be of critical importance in areas affected by tropical cyclones (Beck et al., 2017a; Martius et al., 2016; Muis et al., 2016). Interactions with ice and snow cover are also currently not resolved at higher latitude, which affects the timing and magnitude of both storm surge heights (Muis et al., 2016) and river discharge (Yamazaki et al., 2011). Therefore, results in northern regions, where we find the lowest compound flood hazard potential, are particularly uncertain and should be interpreted with care. More generally, this emphasises the need for local-scale studies in order to accurately quantify compound flood hazard locally.
Finally, we investigate compound flooding interactions under current climate conditions from hydrometeorological processes only, and neglect anthropogenic changes on the catchment and the climate. Human interventions, such as water extractions, water retention, or flood protection infrastructure can affect the travel time and magnitude of extreme discharges and modify the discharge time series (Allen et al., 2018; Veldkamp et al., 2018). Combined with changes in environmental conditions, e.g. due to sea-level rise and changes in storminess, these additional non-stationary drivers can strongly modulate the multivariate dependence structure between flood drivers and affect compound flood hazard (Moftakhari et al., 2017; Wahl et al., 2015).
This paper provides a global perspective of the compound flood potential from riverine and coastal flood drivers. By selecting time series of flood drivers for both hazard types, we derived a global overview of areas particularly exposed to the co-occurrence of high discharge and storm surge level, and we quantified the strength of the interactions between the two variables. We analysed important characteristics of compound flooding related to the timing and joint dependence between river discharge and storm surge extremes. Regional clusters consistently exhibit a high potential for compound flooding. Hotspot regions such as Madagascar, Portugal, northern Morocco, northern Australia, Vietnam, and Taiwan all show a positive and significant dependence between flood drivers and a large number of co-occurring annual maxima. Using the binomial distribution for different flood season lengths, we showed that the dependence between these variables cannot be explained by random extreme interactions within a season. Instead, we hypothesise that this dependence results from synoptic weather systems and interactions between these weather systems and topography.
Extreme impact events caused by synoptic weather systems, like Hurricane Harvey or Idai, highlight the importance of considering compound flood events for flood protection in coastal communities. Currently, regulatory flood hazard maps, such as those used in the US, often only model flooding due to one flood driver (Federal Emergency Management Agency, 2015; Moftakhari et al., 2019). Because such a methodology discards the interactions between river and coastal floods, it can strongly flaw the representation of flood hazard in deltas and estuaries. This is also the case for current global state-of-the-art flood models, and our study provides a first indication of regions along the global coastline where discharge and storm surge extremes are likely to co-occur. In areas coinciding with rapid economic development, this can have strong implications for emergency responders, reinsurance, and local decision makers.
How compound flood events will affect flood impacts, adaptation strategies, and management operations at local scales is strongly dependent on local conditions and is therefore left for future research. We presented our first insights into how the dependence structure impacts on the probability of hazardous riverine and coastal conditions globally. Such a method could be used to generate stochastic events to explore the impact of unforeseen events within a certain catchment through hydrodynamic and impact model experiments. Similarly, future studies should investigate the importance of synoptic weather conditions with respect to the contributions from local drivers such as estuarine topography, land cover, human interventions, and water management and control in determining the impacts from current and future compound flood events.
The paired daily discharge and storm surge time series at the river mouth locations used for this study are available at:
The supplement related to this article is available online at:
AC wrote the manuscript with assistance and input from all co-authors. DE and SM performed the numerical simulations for the discharge and storm surge data, respectively, and DE compiled it in the dataset used in this study. AC designed the methodology with the help from DE, HCW, and PJW. AC performed the analysis and analysed the results with all co-authors.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Global- and continental-scale risk assessment for natural hazards: methods and practice”. It is a result of the European Geosciences Union General Assembly 2018, Vienna, Austria, 8–13 April 2018.
This research has been supported by the Netherlands Organisation for Scientific Research (NWO) (VIDI; grant no. 016.161.324) and the UK's National Environmental Research Council (grant no. NE/S010262/1).
This paper was edited by James Daniell and reviewed by two anonymous referees.