Abstract

ESD

Earth System Dynamics

ESD

Earth Syst. Dynam.

2190-4987

Copernicus Publications

Göttingen, Germany

10.5194/esd-16-1029-2025

Concurrent modes of climate variability linked to spatially compounding wind and precipitation extremes in the Northern Hemisphere

Modes of climate variability linked to spatially compounding events

François

Bastien

bastien.francois@knmi.nl

https://orcid.org/0000-0002-9471-026X

Teber

Khalil

https://orcid.org/0009-0006-0075-2052

Brett

Lou

https://orcid.org/0000-0002-1516-3166

Leeding

Richard

https://orcid.org/0000-0002-7717-7678

Gimeno-Sotelo

Luis

https://orcid.org/0000-0001-9305-2325

Domeisen

Daniela I. V.

https://orcid.org/0000-0002-1463-929X

Suarez-Gutierrez

Laura

https://orcid.org/0000-0002-0008-5943

Bevacqua

Emanuele

emanuele.bevacqua@ufz.de

https://orcid.org/0000-0003-0472-5183

1Royal Netherlands Meteorological Institute (KNMI), Research and Development Weather and Climate (RDWK), De Bilt, the Netherlands 2Institute for Earth System Science and Remote Sensing, Leipzig University, Leipzig, 04103, Germany 3Department of Civil and Environmental Engineering, University of Strathclyde, Glasgow, UK 4Department of Earth Sciences, Uppsala University, Uppsala, Sweden 5Departamento de Estatística e Investigação Operacional, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal 6CEAUL – Centro de Estatística e Aplicações, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal 7Centro de Investigación Mariña, Environmental Physics Laboratory, Universidade de Vigo, Ourense, Spain 8Faculty of Geosciences and Environment, Université de Lausanne, Lausanne, Switzerland 9Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland 10Institut Pierre-Simon Laplace, CNRS, Paris, 75252, France 11Department of Compound Environmental Risks, Helmholtz Centre for Environmental Research – UFZ, Leipzig, Germany

These authors contributed equally to this work.

Bastien François (bastien.francois@knmi.nl) and Emanuele Bevacqua (emanuele.bevacqua@ufz.de)

15July2025

16 4 10291051 5July2024 11April2025 3February2025 26July2024

2025

This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/

This article is available from https://esd.copernicus.org/articles/16/1029/2025/esd-16-1029-2025.html

The full text article is available as a PDF file from https://esd.copernicus.org/articles/16/1029/2025/esd-16-1029-2025.pdf

Abstract

Compound wind and precipitation (CWP) extremes often cause severe impacts on human society and ecosystems, such as damage to crops and infrastructure. Spatially compounding events with multiple regions affected by CWP extremes in the same winter can impact the global economy and reinsurance industry; however, our understanding of these events is limited. While climate variability modes such as El Niño Southern Oscillation (ENSO) can influence the frequency of precipitation and wind extremes, their individual and combined effects on spatial co-occurrences of CWP extremes across the Northern Hemisphere have not been systematically examined. Here, by combining reanalysis data and climate model simulations, we investigate how two oceanic and two atmospheric variability modes – ENSO, the Atlantic Multidecadal Variability (AMV), the North Atlantic Oscillation (NAO), and the Pacific North American (PNA) – amplify the wintertime (December–February) frequency of daily CWP extremes and associated spatial co-occurrences across the Northern Hemisphere. We find many hotspot regions where concurrent variability mode anomalies significantly amplify wintertime CWP extreme event frequencies compared to single variability modes. By examining the relationships between frequencies of wintertime CWP extremes across regions, we identify dependencies enabling extreme spatially compounding events, that is winters with many regions experiencing CWP extremes. While ENSO is the most influential variability mode for such extreme spatially compounding events, the occurrence of these events increases further when multiple modes of variability are in anomalous phases. In particular, combinations of modes increase both the number of regions and the population exposed to daily CWP extremes in the same winter. For example, combined ENSO- and NAO+ nearly doubles the number of affected regions compared to neutral conditions on average. Our analysis highlights the importance of considering the interplay between variability modes to improve risk management and adapt to the impacts of spatially compounding CWP extremes.

Deutsche Forschungsgemeinschaft

524780515

HORIZON EUROPE Marie Sklodowska-Curie Actions

101064940

Horizon 2020

101003469

847456

1Introduction

Weather and climate extremes pose severe risks to human health, infrastructure, natural resources, and ecosystem health . Compound weather and climate events, defined as the combination of multiple drivers and/or hazards that contribute to societal or environmental risk, often cause more severe impacts than the respective single hazards . The Intergovernmental Panel on Climate Change (IPCC) Special Report on Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation (SREX) highlighted the importance of studying compound events to improve modeling and risk estimation of weather impacts . Without considering compound extreme weather events, damages from extreme weather events may be either under- or over-estimated e.g.,. A highly studied example of an impactful compound event is hot–dry conditions e.g.,, which are enhanced by land-atmosphere feedback e.g.,. Another example is compound wind and precipitation (CWP) extremes that can cause more damage than high winds or precipitation in isolation e.g.,. For example, the combination of high winds and rainwater can result in severe damage due to the inflow of water through joints or cracks in building structures e.g.,. Furthermore, CWP extremes can lead to a range of other exacerbated impacts, such as increased soil erosion e.g.,, agricultural and forestry losses e.g.,, and damage to energy infrastructure and buildings e.g.,. They can also result in economic impacts through insurance and disaster recovery e.g., and human fatalities e.g.,. Going forward, these CWP extremes are likely to change in frequency and/or intensity due to human-induced climate change .

Alongside co-occurring extremes in the same location, which correspond to multivariate events , another type of compound event involves the simultaneous occurrence of extremes across multiple regions or locations, referred to as spatially compounding extremes . Spatially compounding events are gaining prominence due to the potential to cause widespread impacts on the global food system, disaster management resources, and (re)insurance industries by simultaneously affecting large populations and assets e.g.,. For example, widespread spatially compounding flooding throughout nine European countries in June 2013 placed high pressures on trans-national risk reduction and risk transfer mechanisms and generated EUR 12 billion in losses . In cases where simultaneous extreme hazards across multiple regions occur worldwide, global catastrophe risk pooling could help with reducing the severity of economic shocks . While the characteristics of spatially compounding events such as droughts e.g., or floods e.g., have already been studied, our understanding of spatially compounding CWP extremes and their drivers is limited.

To advance our understanding of CWP extremes and their drivers, this study focuses on the Northern Hemisphere due to the high population density and the severe impacts of CWP extremes in this part of the world e.g.,. Such compound extremes are most frequent in coastal regions of the Pacific and Atlantic oceans e.g., and tend to be more frequent and intense in the winter season e.g.,. Consequently, this study focuses on seasonal frequency of daily CWP extremes during December–February. While multiple drivers can lead to CWP extremes, their influence varies across different geographical regions e.g.,. Typical drivers of CWP extremes include atmospheric rivers e.g., and low-pressure systems e.g.,, including tropical cyclones e.g.,, extratropical cyclones e.g.,, and fronts e.g.,. Accordingly, regions prone to cyclones are particularly exposed to CWP extremes e.g.,. As an example of the potential costs, the extratropical cyclones Anatol, Lothar, and Martin that hit Europe in December 1999 caused damages totaling USD 13.5 billion (adjusted to the value in 2012) and over 150 fatalities .

At the global scale, such drivers are modulated by modes of atmospheric and oceanic variability, which represent variations in atmospheric patterns due to internal climate variability. For example, the North Atlantic Oscillation (NAO) and Pacific North American (PNA) are the leading modes of variability affecting wintertime weather in Europe e.g., and North America e.g.,. During extreme phases of the PNA and NAO, the intensity and location of storms and moisture transport deviate from mean conditions over the Pacific-North American region e.g., and the Euro-Atlantic region e.g.,, respectively. While positive NAO phases intensify westerly winds and shift the North Atlantic storm track toward the northeast, leading to increased storm frequency and intensity over Northern Europe, negative NAO phases weaken the westerlies and amplify storm activity in the Mediterranean region e.g.,. Oceanic modes of variability such as the El Niño Southern Oscillation (ENSO) and the Atlantic Multidecadal Variability (AMV) can further modulate atmospheric modes of variability and CWP extremes. For example, ENSO can influence both the PNA and NAO e.g., as well as weather systems in many regions worldwide e.g.,. The AMV, an alternation of warm and cold sea surface temperatures in the North Atlantic on decadal timescales, has been shown to influence weather in both North America and Europe e.g.,, as well as the long-term variability of the NAO e.g.,. There are complex teleconnections linking the ENSO, NAO, PNA and AMV modes and their influence on weather patterns across the Northern Hemisphere e.g.,. Some of the modes are specific phenomena that are internally driven, such as ENSO, while others are a representation of the dominant weather and climate variability in a particular region, such as the NAO or the PNA.

Understanding the drivers of CWP extremes is crucial for optimizing resource distribution during response efforts, yet research in this area remains limited. Several studies have explored the influence of single or combined mode of variability anomalies on precipitation and wind extremes in isolation, both at regional e.g., and global scales e.g.,. Regarding CWP extremes, some studies have explored the influence of individual variability modes at the regional scale only e.g.,, thus not allowing the assessment of the relation between concurrent climate variability modes and spatial co-occurrences of CWP extremes across several regions. In this study, we use reanalyses data and large ensemble climate model simulations from the CESM General Circulation Model to investigate the effects of ENSO, NAO, PNA, and AMV modes of variability, and their combinations, on the increase in the frequency of December–January–February daily CWP extremes across the Northern Hemisphere. We consider these four modes of variability due to their already-known influence on storm activity and moisture transport in the Northern Hemisphere. As ENSO operates primarily on seasonal timescales e.g.,, seasonal mean indices are considered. Specifically, we (1) analyze how different modes of variability increase wintertime regional frequencies (i.e., seasonal counts) of daily CWP extremes across individual regions of the Northern Hemisphere. We then investigate spatially compounding events with many regions under high CWP extreme frequencies in the same winter by examining the effect of (2) dependencies between CWP extremes across different regions and (3) combinations of modes of variability on such spatially compounding events. Finally, we (4) inspect the atmospheric circulation anomalies associated with these compound events. To the best of the authors' knowledge, the present study is the first to investigate CWP extremes and associated spatially compounding events across the Northern Hemisphere as mediated by multiple large-scale modes of variability. We address this research gap by using large ensemble climate model simulations to provide a robust model-based analysis of the effects of rare concurrent variability mode anomalies on low probability daily CWP extremes e.g.,.

The study is structured as follows: Section describes the data used as well as the methodology to analyze the effects of variability modes on CWP extremes and spatially compounding events. Results are provided in Sect. . Conclusions, discussions, and perspectives for future research are presented in Sect. .

2Data and methods

We examine the influence of four variability modes on CWP extremes across 25 selected regions in the Northern Hemisphere defined in the SREX see Fig. . We chose these regions as they are standard reference in IPCC reports as they encompass areas with relatively homogeneous climatic characteristics . While using these regions does not enable an explicit analysis of dependencies between local-scale CWP extremes and modes of variability, it allows for complementing IPCC assessments.

Figure 1

Presentation of the regions under study which consist of the SREX regions adopted by the Intergovernmental Panel on Climate Change (IPCC) in the Northern Hemisphere. The bottom part shows short and full names of the regions. Regions are clustered in macroareas (see legend). To balance the number of regions between continents, we partitioned the Asian continent into two macroareas (North Asia and Central-South Asia) and included the Arabic Peninsula (ARP) region in the Africa macroarea.

2.1Model and reanalysis data

We employ model simulations from the coupled Community Earth System Model (CESM; spatial resolution of 1.25° by 1.25°), which provides 40 ensemble members . These ensemble members are derived from the CESM model under the same model physics and external forcings, with each of the members starting in 1920 and initialized from slightly different initial states, leading to different evolutions from internal climate variability . Considering multiple members provides a large sample size that allows for assessing the effect of internal climate variability and infrequent combinations of modes of variability on rare compound events e.g.,. Daily total precipitation and daily mean of wind speed are extracted for the historical period (1950–2005). Simulated data are then extended until 2019 using the emission scenario associated with a radiative forcing of +8.5 Wm-2 (RCP 8.5 scenario), resulting in a total of 70×40=2800 years of data. We chose daily averages rather than maxima for wind due to data availability for the CESM simulations.

To evaluate the CESM model, we employed ERA5 reanalysis data (spatial resolution of 0.25°) for the period 1959–2019, from which we also extracted daily means of wind speed and daily total precipitation for consistency.

Seasonal indices for the two oceanic (ENSO and AMV) and two atmospheric (NAO and PNA) variability modes for both CESM and ERA5 data were calculated from monthly data using the National Center for Atmospheric Research (NCAR) data package Climate Variability Diagnostics Package CVDP,. The seasonal indices for the NAO, PNA and AMV were calculated as the mean of the monthly values for December, January and February. For these indices, December was taken from the year n-1, while January and February were taken from year n, and the resulting seasonal index is assigned to year n. For ENSO we proceeded similarly, but used November–January averages to account for lagged effects e.g.,. For each mode, we defined their positive and negative phases when the index was above or below its mean by 1 standard deviation, respectively, otherwise the phase was considered neutral.

2.2Methods 2.2.1CWP extremes

Many techniques have been utilized to characterize CWP extremes, with the selection of a specific method being guided by the research question. For example, the correlation between wind and precipitation has been quantified at daily to seasonal timescales e.g.,. Logistic regression models have been applied to quantify the likelihood of a precipitation extreme occurring given the presence of a wind extreme e.g.,. Alternative approaches include examining tail dependence e.g., and employing impact-focused metrics e.g.,. The most straightforward approaches include counting extreme wind and precipitation co-occurrences above a given percentile e.g., or using extremal dependency measures estimating the probability of one variable being extreme given that the other one is extreme e.g.,.

Here, to investigate winter season (December–February) CWP extremes at the grid cell level, we derived seasonal counts of daily CWP extremes which were defined as wind and precipitation values simultaneously exceeding high thresholds. This results in one count per season per grid cell, which allows for investigating the effect of seasonally-averaged climate variability modes on the counts. We use the 98th percentile of wind and precipitation over the 1950–2019 period as thresholds for the main analysis based on data from the CESM model. Percentile-based thresholds are frequently used to investigate climate extremes e.g.,. Following and , we chose the 98th percentile, which is a compromise to capture the most extreme events in the CESM simulations while ensuring a sufficiently large sample size for robust statistical analysis. For model evaluation, which involves both the CESM model and ERA5 reanalyses (Figs. S1–S5 in the Supplement), we use the 95th percentiles over the 1959–2019 period – such a lower threshold allows us to present a more robust evaluation. The reason for this is that, given the ERA5's limited period, extremes in the reanalysis data set are more scarce and associated statistics for very extreme events are largely affected by sampling uncertainty . Selecting a slightly lower threshold allows us to reduce this sampling uncertainty and thus improve confidence in assessing the model's ability to simulate extremes e.g.,.

2.2.2Regional and spatially compounding effects of combined variability modes

To investigate the effect of variability modes on wintertime CWP extremes, we use three different metrics. The first metric allows for assessing the effect of modes in individual regions:

Metric 1. Regionally averaged counts of CWP extremes. Wintertime CWP counts are averaged by region over landmasses, weighted by the cosine of latitude to prevent overrepresentation of grid cells closer to the poles. This count-based metric allows us to assess the effects of variability modes on individual regions separately.

As Metric 1 is derived for each region individually, the influence of variability modes on the high regional frequencies of CWP extremes across multiple regions in the same winter (i.e., spatially compounding events) cannot be deduced directly. For example, based on Metric 1, we find that the variability mode phase ENSO+ modulates regionally averaged CWP extreme frequencies for North America and Central Asia. However, a possibility could be that half of the winter seasons with ENSO+ leads to increased CWP extremes for North America only, while the other half affects Central Asia, thus not simultaneously. Examining the dependencies between counts of CWP extremes in different regions can provide preliminary information on spatially compounding events e.g., because regions connected by positive dependencies tend to experience CWP extremes at the same time. Thus, as a first step for investigating spatially compounding events, we analyze dependencies in Metric 1 computed for different regions, so as to provide preliminary information on groups of regions that may be affected by CWP extremes during the same winters.

Then, to examine spatially compounding events, we use two additional metrics. We employ these metrics to investigate the effects of variability modes on regional high frequencies of CWP extremes across multiple regions in the same winter (Metric 2) and on the total population of the Northern Hemisphere exposed to CWP extremes in the same winter (Metric 3).

Metric 2. Total number of affected regions during the same winter. For a given winter, a region is considered as affected by CWP extremes when the regionally averaged count of CWP extremes (Metric 1) is above its 80th percentile derived from the distribution of the 1950–2019 period. Then, similar to , the total number of affected regions during the same winter is counted. Although choosing a higher percentile (>80th percentile) would enable us to focus on cases where regions are more severely affected, it would considerably limit the number of regions reported as affected and, consequently, the statistical robustness of our results.

Metric 3. Total population exposure. To assess the effects of variability modes on the population of the Northern Hemisphere exposed to CWP extremes during the same winter, we calculated for each winter the weighted averages of CWP extremes, wCWP (hereafter referred to as “population-weighted CWP extremes”), using CWP extremes and population counts at the grid cell level as follows: wCWP=∑i=1NgridNCWPi×Popi∑i=1NgridPopi, where NCWPi and Popi are the seasonal counts of CWP events and the population count at grid cell i, respectively, and Ngrid is the total number of grid cells over the Northern Hemisphere. Population weighting is utilized here as a surrogate for the assets at risk that could experience damages due to CWP extremes e.g.,. For this purpose, global population counts for the year 2020 from the Gridded Population of the World (GPW) product have been used .

In addition to enabling quantifying the number of affected regions depending on variability modes, Metric 2 is also used in Fig. b to assess the effect of the dependencies between regions on spatially compounding events. This analysis is performed by comparing the number of affected regions (i.e., Metric 2) from the original dataset with the number obtained after breaking the dependencies via randomly shuffling regional CWP extreme counts using bootstrap in all regions in time .

2.2.3Identifying relevant effects of single and combined variability modes for the regional and spatial metrics

We quantify the effect of a positive (or negative) phase of a single variability mode of interest on wintertime CWP extremes based on the ratio between (i) CWP metrics under the positive (or negative) phase of the mode of interest, while additionally conditioning the other modes in a neutral state, and (ii) CWP metrics under all variability modes in the neutral phase. Hereafter, we referred to this effect as a direct effect. Following , additionally conditioning all other modes in the neutral phase in (i) serves for better isolating the causal effects of the individual variability mode of interest. Specifically, the additional conditioning takes into account confounding effects arising from the considered modes, although some confounding effects may remain. In particular, modes in neutral states still vary within the range of neutral conditions and we do not control for them. In addition, further effects may arise from variability modes not considered in this study. Similarly to the analysis of single variability modes, we quantify combined effects, that is the effect of concurrent variability modes in non-neutral phases based on the ratio between (i) CWP metrics under the concurrent modes of interest while conditioning all other modes to neutral conditions and (ii) CWP metrics under all variability modes in neutral phase. In terms of notation, we refer to concurrent climate variability modes by specifying (positive or negative) phases of NAO, PNA, ENSO, and AMV (in this order), with unspecified modes being in neutral states.

Although the main part of our study is based on the analysis of the metrics defined in Sect. , we provide an overview of the direct (combined) effects of modes on CWP extreme frequencies at the grid cell level in Fig. (Fig. ). Results for the effects of variability modes on regionally averaged metric (Metric 1) and spatial metrics (Metrics 2 and 3) ratios are presented in Figs. , , and . For both regional and spatially compounding cases, we focus on the modes of variability influences leading to an enhancement in the means of the different metrics compared to neutral conditions (hereafter referred to as positive effects) – we focus on the increase, rather than the decrease, as this is of relevance for potential impact to society. This study does not investigate the influence of variability modes on the decrease of the metrics. When results for the different metrics are presented with box plots (Figs. , , and ), the interquartile range and mean of the distribution are displayed.

Figure 2

Direct effect of variability modes on wintertime compound wind and precipitation (CWP) extremes. Ratio of the average wintertime (December–February) CWP extreme frequencies for the (a) positive and (b) negative phases of NAO (while other variability modes are in their neutral phases) compared to neutral conditions (all variability modes being in their neutral phases) based on the CESM model. Corresponding maps are also displayed for (c, d) PNA, (e, f) ENSO, and (g, h) AMV. Numbers in the headers indicate the conditional empirical return period T (in years) for positive and negative phases of the variability modes (while other modes are in their neutral phases). The empirical return period for neutral conditions (i.e., when all modes are in their neutral phases) is T=3 years. Stippling indicates significant differences of mean frequency relative to neutral conditions at the 10 % significance level using permutation tests (two-sided) with Bonferroni correction. The framed regions in blue are those where the direct effects of variability modes significantly increase regionally averaged CWP extreme frequencies compared with neutral conditions (following the methodology defined in Sect. , and reflecting information in Fig. ).

Figure 3

Influence of concurrent variability modes on regional compound wind and precipitation (CWP) extremes. Ratio of the regionally averaged wintertime CWP extreme frequencies for the relevant individual and concurrent variability modes (following the methodology defined in Sect. ) with respect to average frequency for neutral conditions based on the CESM model (orange dots represent associated values of the ratio derived from ERA5, when available, with each dot representing one winter season). Results are displayed for (a) Europe and Africa, (b) North America and Central America, and (c) Northern Asia and Central-South Asia. The conditional empirical return periods (in years; see Methods) for individual and concurrent variability modes are indicated on the top. Gray arrows indicate significant differences in the mean frequency relative to neutral conditions (at the 10 % significance level, with Bonferroni correction). Green boxes indicate distributions for neutral conditions (all variability modes being in their neutral phases) and green horizontal lines indicate ratios equal to one. Blue boxes indicate mode combinations with amplified effects (that is, higher means of regionally averaged CWP extreme frequencies than their underlying mode sub-combinations (methods Sect. )). We remove concurrent variability modes that occur in less than 10 winters and regions without significant effects arrows (EEU, WCE, GIC, CAR, and SCA).

Figure 4

Effects of variability modes on wintertime compound wind and precipitation (CWP) extremes in SREX regions. For different regions (see labels on the left side of the matrix), the ratio of the regionally averaged CWP extreme frequencies for individual and concurrent variability modes (see labels on the bottom side of the matrix) with respect to the average frequency under neutral phases of all modes (Metric 1, see Methods), based on the CESM model. Note that positive effects of concurrent variability modes (ratio above 1) tend to be stronger than negative effects (ratio below 1). The individual and concurrent variability modes with significant effects (at the 10 % level with Bonferroni correction) relative to neutral conditions are indicated with an asterisk. Blue cell borders indicate combinations with higher means of regionally averaged CWP extreme frequencies than their underlying mode sub-combinations (methods in Sect. ). To ensure robustness, we do not display concurrent variability modes that occur in less than 10 winters. Based on the displayed information, the numbers on the top and right margins indicate the total count of regions and mode combinations with significant effects of modes, respectively.

Figure 5

Dependencies between regional wintertime compound wind and precipitation (CWP) extremes. (a) Matrix of Spearman correlations for regionally averaged wintertime (December–February) CWP extreme frequencies between regions based on the CESM model. Stippling indicates significant correlations (correlation values that are outside the 90 % centered confidence interval). A contour is added for regions in the same macroarea. The boxplot shows summaries of the distribution of the correlations shown in the matrix (interquartile range, median and outliers); the gray background shows the confidence interval for no correlation (bootstrap-based 90 % range). (b) Histograms of the number of regions with a high frequency of CWP extremes during the same winter based on the CESM model (orange histogram) and when assuming independence between regional wintertime CWP extreme frequencies (gray histogram; obtained via shuffling the data 1000 times via bootstrap). Vertical lines show the 50 year return levels under dependence and independence. Bootstrap-based confidence intervals for each bin and 50 year return level at 10 % significance level are also shown. (c) Boxplots of the number of affected regions given different numbers of variability modes in anomalous phases. The conditional empirical return periods (in years; see Methods) for the different number of variability modes in anomalous phases are indicated on the right. Mean of the ratios under a different number of modes in anomalous phases that are not significantly different (α=0.10, one-sided permutation test) are indicated with the same letters . The average number of regions affected under neutral conditions is indicated by a vertical line.

Figure 6

Influence of variability modes on spatially compounding wind and precipitation (CWP) extremes. (a) Distributions of the number of regions (left y axis) experiencing a high frequency of CWP extremes during the same winter (December–February) for different individual and concurrent variability modes (x axis) based on the CESM model. The right y axis shows the ratio of the metric to its average under neutral phases of all modes. The combinations presented are selected according to the methodology defined in Sect. . Gray arrows indicate significant differences in the mean with respect to neutral conditions at the 10 % significance level using one-sided permutation tests and Bonferroni correction. Green boxes indicate distributions for neutral conditions (all variability modes being in their neutral phases). Blue boxes indicate combinations exhibiting higher means than their underlying mode sub-combinations (methods Sect. ). (b) The same as panel (a), but for the population-weighted CWP extremes over the Northern Hemisphere. Symbols + and - above the x axis in panels (a) and (b) highlight combinations having a significant effect for a metric containing ENSO+ and ENSO-, respectively.

Figure 7

Compound wind and precipitation (CWP) extremes for concurrent anomalies in variability modes and associated sea level pressure anomalies. (a) Ratio of the average wintertime (December–February) CWP frequencies for the concurrent variability mode NAO-ENSO+ (while other variability modes are in their neutral phases) compared to winter with all variability modes in neutral phases based on the CESM model. (b) Mean sea level pressure anomalies for NAO-ENSO+ (while other variability modes are in their neutral phases) compared to neutral conditions (all variability modes being in their neutral phases) based on the CESM model. Corresponding maps are also displayed for (c, d) PNA+ENSO+, (e, f) ENSO+AMV-, (g, h) PNA-ENSO-, (i, j) ENSO-AMV+, (k, l) NAO+ENSO-AMV+, and (m, n) NAO-PNA+ENSO+. Numbers in the headers indicate the conditional empirical return period T (in years; see Methods) for the different concurrent variability modes stated in the title of the panels, whereas the return period for neutral conditions (all modes are in their neutral phases) is T=3 years. Stippling indicates significant differences in mean frequency relative to neutral conditions at the 10 % significance level using permutation tests (two-sided) with Bonferroni correction. The framed regions are those where concurrent variability modes significantly increase regionally averaged CWP extreme frequencies compared with neutral conditions; by reflecting information in Fig. , the dark framing indicates an amplified effect with respect to underlying mode sub-combinations (methods in Sect. ).

Given the focus on four variability modes, each with three possible phases, there are 81 possible phase combinations, motivating the need for a synthesis of their effects. When presenting the results in Sect. , we focus our analysis on individual and concurrent variability modes having a significant and positive effect on the different metrics using permutation tests (see Sect. for more details on calculating the significance). To support the interpretation of the effects of concurrent variability modes on CWP metrics, in Figs. and we also show the effects of all univariate variability modes that contribute to the concurrent modes with significant effects regardless of whether they exert a significant effect in isolation compared to neutral conditions or not. Furthermore, to ensure robustness in the results, we disregard concurrent variability modes that occur very rarely by only considering concurrent variability modes with an empirical return period, defined as the inverse of their relative frequency, greater than 280 years (i.e., occurring for more than 10 years in our 2800 years dataset). It should be noted here that the empirical return periods for concurrent modes provided in this study are conditional as they are calculated by conditioning all other modes to neutral conditions. Therefore, such conditional return periods, despite providing an indication of the rarity of concurrent modes, should not be considered as absolute return periods as, by construction, they are larger than the unconditional return periods obtained without conditioning the other modes in neutral conditions. To provide an overview of the synthesis of the variability modes, effects of all possible concurrent modes are displayed for regionally averaged metric in Figs. S8–S14 in the Supplement (for a selection of regions for the sake of brevity) and spatial metrics are provided in Fig. S15 of the Supplement.

Finally, to provide insights into possible amplified effects of concurrent variability modes relative to the effects of the modes that contribute to such combination, we identify combinations for which the average metric is higher than that of the underlying mode sub-combinations (note that no test here is performed, so this should be interpreted carefully, see explanations in Sect. ). For example, the concurrent modes NAO-PNA+ENSO+ are deemed to possibly have an amplified effect relative to the underlying mode sub-combinations if the average metric under NAO-PNA+ENSO+ is higher than the underlying mode sub-combinations, which are NAO-, PNA+ and ENSO+ taken in isolation, as well as the bivariate mode sub-combinations NAO-PNA+, NAO-ENSO+ and PNA+ENSO+. The identified combinations with amplified effects are highlighted in Figs. , , , and via dark blue box plots and boxes.

2.2.4Permutation test procedure to assess significance in effects of combined variability modes

As already mentioned in Sect. , the statistical significance of the effects of individual and concurrent variability modes on CWP metrics compared to neutral conditions was assessed using permutation tests e.g.,. Specifically, for a given CWP metric, we test whether the ratio of the average of the metric associated with a given set of phases of interest (e.g., NAO+ENSO-, set as the numerator) to the average of the metric under neutral conditions (set a denominator) is larger than one at significance level α=0.10 based on one-sided tests. Compared to a lower significance level, our chosen level allows the detection of significant effects of modes of variability while reducing false negatives in the context of small sample sizes. To test the hypothesis that the average of the metric under individual and concurrent variability modes is higher than under neutral conditions, we compared the ratio obtained from the original samples with a confidence interval of the ratio obtained from data without an effect of the modes on the CWP metric. Specifically, the latter was obtained via randomly permuting without replacement the samples for both numerator and denominator and re-estimate the ratio from the resampled data. By repeating this procedure, we can then define a confidence interval for the ratio and a critical region for test rejection. If the original ratio is higher than the (1-α)×100th percentile, the average of the metric associated with the set of phases of interest is considered to be significantly higher than that of the neutral conditions. As several tests are carried out for the different concurrent modes, a Bonferroni correction e.g., is applied to control the overall probability of Type I (or false positive) errors. Please note that, while identifying concurrent variability modes with significantly amplified effects relative to the effects of the modes that contribute to such combination is relevant, it has not been done in this study as no statistically consistent procedure including Bonferroni correction has been found for such a statistical problem. For example, to identify the concurrent modes NAO+ENSO+AMV+ as having a significantly amplified effect on the metrics, six permutation tests comparing the effects of NAO+ENSO+AMV+ against those of its six sub-combinations must all reject the equality of effects. Designing such a statistical test procedure while controlling the false discovery rate goes beyond the scope of the present study and is therefore not included in the analysis. In addition, in Figs. and only, we assess whether averages of the CWP extreme frequencies at the grid cell level associated with a given set of phases are significantly different compared to neutral conditions (i.e., not necessarily larger). For these tests performed at the grid cell level, we use two-sided permutation tests. In this case, the critical region for rejection is defined using the (α2)×100th and (1-α2)×100th percentiles.

Regarding the number of samples m for permutation tests, several trials showed that choosing m=100000 allows us to obtain robust results for the three different metrics. Thus, the analysis of the three metrics was carried out using m=100000. However, when applied to the grid cell level for Figs. and , results from permutation tests proved less sensitive to the choice of the number of sample, and m was chosen to be equal to 100.

3Results

Before investigating the effect of individual variability modes and concurrent modes on CWP extremes and associated spatially compounding events in CESM simulations for the 1950–2019 period, we carried out a model evaluation with respect to ERA5 over the 1959–2019 period. Such a model evaluation is mainly performed for the effects of individual variability modes, as a robust assessment of the effects of concurrent variability modes requires a large sample size . For the same reason, the effects of individual variability modes in ERA5 are evaluated without constraining the other modes to be in the neutral phase (which thus differs from the direct effects of variability modes defined in Sect. ).

Results for model evaluation are displayed in Figs. S1–S5. Although the simulated CWP absolute frequencies exhibit some biases with respect to ERA5 over the Northern Hemisphere (Fig. S1), CESM provides an adequate representation of the anomalies in CWP frequencies induced by individual variability modes relative to neutral conditions (Figs. S2–S5). We conclude from this model evaluation that CESM is suitable for further investigating the CWP extremes and their relationships with variability modes. In the following, we assess the effect of modes of variability and their combinations on CWP extremes via the CESM model based on the 98th percentiles and also provide a one-to-one comparison with ERA5 based on the same percentiles.

3.1Direct effects of variability modes on regional CWP extremes

Examining the direct effects of variability modes on CWP extremes in CESM simulations (Fig. ) shows that different phases of NAO, PNA, and ENSO significantly modulate CWP extreme frequencies over multiple regions of the Northern Hemisphere. Within Europe, NAO+ and NAO- significantly increase CWP extremes in Northern Europe (NEU) and the Mediterranean region (MED), respectively (Fig. a and b). NAO+ also significantly influences CWP extremes in the Russian-Arctic (RAR) and Western Siberia (WSB). The low pressure associated with storms under NAO+ is accurately represented in CESM simulations (Fig. S6a in the Supplement), with negative pressure anomalies partially or fully covering the regions experiencing a significant increase in CWP extremes. An exception to such CWP enhancement associated with low pressure is found for Northeastern North America (NEN) during NAO-, where high sea level pressure (SLP) anomalies of ∼+3 hPa are associated with increased CWP extremes. While no regions experience increased CWP extremes under PNA+ (Fig. c), PNA- significantly increases CWP extremes compared with neutral conditions in Eastern Siberia (ESB) and Russian Far East (RFE) within North Asia (Fig. d) and in Northeastern North America (NEN), Northwestern North America (NWN), and Western North America (WNA) in the North America macroarea.

Among the oceanic modes of variability ENSO and AMV, the former has the most relevant and widespread effects. ENSO+ (Fig. e) increases CWP extremes across three regions of Africa: Arabian Peninsula (ARP), Sahara (SAH), and Western Africa (WAF). For ENSO- (Fig. f), CWP extremes are significantly increased for two regions within Central-South Asia: South Asia (SAS) and East Asia (EAS). For the AMV, we find no significant regional effects of AMV+ or AMV- conditions (Fig. g and h), in line with weak sea level pressure anomalies associated with AMV phases, ranging from ∼-1.5 to ∼+1.5 hPa (Fig. S6g and h).

3.2Effects of concurrent anomalies in variability modes on regional CWP extremes

Model simulations (CESM) show that not only individual variability modes can have effects on regional wintertime frequencies of CWP extremes, but also combinations of modes (Fig. ; Fig. S7 in the Supplement shows the same but for absolute frequencies). Figure provides a summary of these effects, including the count of modes combinations that enhance CWP extremes compared to neutral conditions for each region. For example, Northern Europe (NEU) exhibits the highest count of mode combinations that significantly enhance CWP extreme frequencies compared to neutral conditions, with 11 mode combinations that all involve NAO+. The presence of a given mode phase or combination of mode phases that has an effect on CWP extremes in multiple regions suggests the potential for spatially compounding events, which will be examined in Sect. . In the following, we focus on describing the effects of a selection of mode combinations and regions in Figs. –. To maintain clarity and conciseness, we do not discuss all regions and mode combinations in the text, and readers can explore specific regional effects directly in the figures.

For the Mediterranean region (MED) in Europe (Fig. a), five concurrent variability modes containing NAO- have significant positive effects with respect to neutral conditions. Among these five concurrent variability modes, three combinations of modes have an amplified effect relative to the underlying mode sub-combinations (blue boxes). For example, in the MED region, CWP extremes are on average ∼1.5 times more likely under the concurrent mode NAO-ENSO+ than under neutral conditions, while being ∼1.1 and ∼1.2 times more likely than under NAO- and ENSO+ in isolation, respectively. Note that the conditional empirical return period of NAO-ENSO+ (i.e., conditioning all other modes to neutral conditions) is 65 years, while it is equal to 20 and 30 years for NAO- and ENSO+, respectively, which means that the concurrent mode NAO-ENSO+ (with other modes being in neutral conditions) is ∼3 and ∼2 times less likely to occur in any given year than NAO- or ENSO+ in isolation (with other modes being in neutral conditions).

Concurrent modes can also amplify CWP extreme frequencies in other regions. For example, in the Arabian Peninsula (ARP) within Africa, CWP extremes are on average ∼3.5 times more likely than neutral conditions under ENSO+. However, when ENSO+ combines with NAO-PNA+, CWP extremes are on average ∼4 times more likely, while both individual modes (NAO- and PNA+) do not have a significant effect on the mean frequency of CWP extreme frequencies relative to neutral conditions. Interestingly, the combination of NAO-ENSO+ has an amplified effect relative to the underlying mode sub-combinations in all three African regions, increasing approximately the likelihood of CWP extremes frequencies by a factor of ∼4, 1.5, and 2 relative to neutral conditions in the Arabian Peninsula (ARP), Sahara (SAH), and Western Africa (WAF), respectively.

The Central North America (CNA) region illustrates how the effects of individual modes can combine to significantly increase CWP extreme frequencies in a given winter relative to neutral conditions, even if individual modes do not exert significant effects in isolation. For example, despite the fact that both PNA+ and ENSO+ in isolation do not significantly increase CWP extremes (indicated by no gray arrows in Fig. b), combined PNA+ENSO+ increases significantly the likelihood of CWP extremes by a factor of ∼1.5 compared to neutral conditions (indicated by the gray arrow and blue box). Furthermore, within North America, the region with the most combinations of modes that present amplified effects on CWP extreme frequencies is East North America (ENA), with six separate combinations (see blue boxes). Of these combinations, ENSO is always in the positive phase, indicating its potentially large effect on CWP extremes.

In North Asia, concurrent variability modes increase CWP extreme frequencies by ∼1.5 to ∼2 times relative to neutral conditions depending on the combination. In particular, PNA-ENSO- increases CWP extreme frequencies by ∼1.5 times in three of North Asia's four regions (Eastern Siberia ESB, Russian-Arctic RAR, and Russian Far East RFE).

Within Central South Asia, the Tibetan Plateau (TIB) exemplifies again how individual modes without significant effects on CWP extremes (here, AMV- and ENSO+) can combine to significantly increase the likelihood of CWP extremes during winters. Across the five regions of Central South Asia, all concurrent modes that, on average, led to more CWP extremes than their underlying mode sub-combinations (blue boxes) include the variability mode ENSO. Of these concurrent modes, ENSO is always in the positive phase within the inland regions (ECA, TIB, and WCA) while being in the negative phase in the coastal regions (EAS and SAS).

Regarding the comparison of CESM results with those of ERA5, the range of CESM ratios of regionally averaged CWP extreme frequencies with respect to neutral conditions generally covers the values for ERA5 for the different concurrent modes when available (Fig. ).

Results from Fig. and the summary in Fig. cannot be used to conclude whether the effects of concurrent variability modes lead to spatially compounding CWP extremes, that is, to high wintertime frequencies of CWP extreme across multiple regions during the same winter. These figures illustrate the effect of individual and concurrent variability modes on regionally averaged CWP extreme frequencies, which are derived for each region separately. Nevertheless, the number of regions where each mode combination has significant effects in Fig. (see numbers on the top of the matrix) suggests that some mode combinations may potentially lead to spatially compounding CWP extremes. For example, NAO-ENSO+ significantly enhances regional CWP extreme frequencies in eight regions, which means that if these regional effects of NAO-ENSO+ can manifest in the same winter, NAO-ENSO+ would lead to spatially compounding CWP extremes. In the next section, we assess whether individual and concurrent variability modes can lead to concurrent CWP extremes during the same winters across regions.

3.3Effects of concurrent variability modes on spatially compounding CWP extremes

As the first step in the investigation of spatially compounding events, we examine dependencies between counts of CWP extremes in different regions which provides preliminary information on groups of regions that may be affected by CWP extremes during the same winters. Figure a shows Spearman correlations of regionally averaged CWP extreme frequencies (Metric 1) between all pairs of regions. Regions that are geographically close tend to be positively correlated (see correlation values in contoured black boxes) which is in line with spatial autocorrelation. However, some regions that are from different macroareas and therefore more distant can also be correlated, e.g., Central-South Asian regions with African regions or North American regions with North Asian regions. This highlights the potential for underlying effects of modes of variability that can connect distant regions. Notably, Fig. a highlights that most of the pair correlations between CWP extreme counts in different regions are positive.

In line with such dominant positive correlations (Fig. a), we find that dependencies among regions overall enhance the potential for spatially compounding events. Specifically, the distribution of the number of regions affected during the same winters (Metric 2) is different from that obtained by assuming independence between regions (Fig. b), with dependencies elongating both tails of the distribution. For spatially compounding extremes on the right tail of the distribution, this implies a higher 50 year return level for the number of regions affected by high CWP counts in the same winter compared to independence, that is 11 instead of 9 affected regions (see vertical lines). Consistent with findings in , this dependency-driven shift is not observed for the mean of the distribution (not shown).

The increased count of regions under CWP extremes due to the dependencies among regions can be linked to modes of variability. That is, the high number of regions simultaneously affected by CWP extremes, which is possible due to the dependencies among regional CWP extremes, would also not be possible without the effects of concurrent variability modes (Fig. c). In particular, when all variability modes are neutral, on average around four regions are affected simultaneously by high regionally averaged CWP extremes. Generally, as the number of variability modes in anomalous conditions increases, the likelihood of multiple regions simultaneously experiencing extreme conditions also increases (Fig. 5c; see letters indicating significant differences among distributions, with the exception of the distributions related to three and four modes). Notably, this effect is even more marked when focusing on winter seasons with an extreme number of regions affected (right whiskers of the boxplots). Overall, these results indicate that concurrent anomalies in variability modes are key for spatially compounding events.

Given the dependencies among regions and that combinations of variability modes are essential for spatially compounding CWP extremes, we move to identify which are the relevant mode combinations that enhance the number of regions affected by CWP extremes (Metric 2) and population exposed to CWP extremes (Metric 3). We find that 13 individual and concurrent variability modes significantly increase the number of affected regions compared to neutral conditions (Fig. a). Three bivariate combinations have a significant effect on the total number of affected regions and an amplified effect relative to their underlying mode sub-combinations (see blue boxes for ENSO+AMV-, PNA+ENSO+, and NAO-ENSO+). NAO-ENSO+ nearly doubles the number of regions simultaneously exposed to CWP extremes on average relative to neutral conditions.

In terms of population exposed to CWP extremes, four combinations (ENSO-AMV+, NAO-PNA+ENSO+, NAO+ENSO-AMV+, and PNA-ENSO-) are identified as having a significant effect compared to neutral conditions (Fig. b). Among these four combinations, only one (NAO-PNA+ENSO+) has already been identified as having a significant effect on the number of affected regions in Fig. a. It highlights that the heterogeneous distribution of population density across regions needs to be considered to assess the societal vulnerability to CWP extremes. In particular, variability modes in isolation do not lead to significant effects on the population exposure compared to neutral conditions (Fig. b), indicating the importance of considering combinations of modes to distill the effects of modes of variability on the population affected. Notably, while ENSO+ dominates the influence of modes of variability on the number of affected regions (see “+” sign in Fig. a), ENSO- dominates for the population affected (see “-” sign in Fig. b).

3.4Physical mechanisms underlying spatially compounding CWP extremes

We now move to inspect the physical mechanisms leading to spatial co-occurrences of CWP extremes by analyzing SLP anomalies for concurrent variability modes with significant effects on spatially compounding events (Fig. ). Here, negative SLP anomalies during winter are generally considered indicative of storminess, a key driver of CWP extremes. In the analysis, we restrict to the seven combinations having an amplified effect relative to their underlying mode sub-combinations on Metric 2 and Metric 3 (i.e. combinations with blue boxes in Fig. ): namely, NAO-ENSO+, PNA+ENSO+, and ENSO+AMV- for affected regions (Fig. a) and PNA-ENSO-, ENSO-AMV+, NAO+ENSO-AMV+, and NAO-PNA+ENSO+ for exposed population (Fig. b). These combinations have a reasonably low conditional return period (≤82 years; see Methods), thus are not very unlikely events in the model simulations.

The largest positive effects on the number of affected regions (Fig. a) occurs under concurrent ENSO+ (El Niño phase) and NAO-, that is NAO-ENSO+ (Fig. a and b). Under such a combination of modes of variability NAO-ENSO+, negative SLP anomalies patterns intensify and expand over the North Atlantic Ocean compared to those for modes NAO- and ENSO+ in isolation (Fig. S6b and e), in line with the enhancement of CWP extremes over many regions of North America, North Central America, Europe and Africa.

Although the combination of modes PNA+ENSO+ has smaller significant positive effects on the number of affected regions than that of NAO-ENSO+ (Fig. a), PNA+ENSO+ may be more frequent than NAO-ENSO+ (given their lower conditional return period; see Methods). Such a PNA+ENSO+ combination (Fig. c and d) intensifies the negative SLP anomalies patterns on both sides of North America with respect to PNA+ and ENSO+ modes in isolation (Fig. S6c and e). This is consistent with the amplification of CWP extreme frequencies for regions within North America, Central America, and Africa. In contrast, when ENSO+ combines with AMV- (ENSO+AMV-; Fig. e and f), negative SLP anomalies patterns are quite similar to those of ENSO+ in isolation (Fig. S6e); this suggests that AMV- has a limited amplifying effect on CWP extremes, in line with Fig. a.

Under the resulting combination PNA-ENSO- (Fig. g and h) – which shows the largest positive effects on the population exposure (Fig. b) – negative SLP anomalies are stronger than under PNA- and ENSO- in isolation over the North Asian region (Fig. S6d and f). These intensified low-pressure conditions associated with storminess increase CWP extreme frequencies over a number of densely populated regions in North Asia and East Asia. PNA-ENSO- exemplifies how combined variability modes can lead to spatially compounding events by adding up the effects of individual modes in isolation. Specifically, under PNA-ENSO- the population exposed to extreme CWP extremes under PNA- and ENSO- in isolation partially adds up, resulting in spatially compounding events with a large population affected by CWP extremes.

Although the effects of combined variability modes ENSO-AMV+ (Fig. i and j) are rather limited, when combined with NAO+ (NAO+ENSO-AMV+; Fig. k and l) negative SLP anomalies are intensified over high latitudes, mainly due to NAO+ influence (Fig. S6a) and over South Asia, due to ENSO- (Fig. S6f). These intensified low-pressure systems lead to an amplification of CWP extremes in Northern Europe and South Asia. It is interesting to note that these two regions are not the same as those affected by ENSO-AMV+ (Fig. i and j), suggesting a strong influence of the NAO on regions impacted by CWP extremes.

Similar to what is observed for NAO-ENSO+, NAO-PNA+ENSO+ features a negative SLP anomaly extending across the Southern North Atlantic (Fig. n). However, for some regions, the addition of PNA+ to NAO-ENSO+ tends to weaken the effects of concurrent variability modes on CWP extremes, as some regions exhibit, on average, fewer regional CWP extremes than those obtained under NAO-ENSO+ (see blue boxes in Fig. n and m).

4Discussion

Several results regarding the direct effects of variability modes on CWP extremes in CESM simulations (Sect. ) align with well-established findings in the literature, thereby validating the robustness of CESM modeling. Table S1 in the Supplement provides a concise summary of the agreement between CESM simulations and existing literature. Among these results are the influence of the NAO+ regime favoring CWP extremes in Northern Europe due to stronger than average westerly winds e.g.,, or the influence of PNA- favoring blocking in the Pacific e.g.,, increasing storm frequency in the Northern Pacific. However, several results lack direct support from existing literature. While they may represent novel findings, they should be interpreted cautiously, considering the biases of CESM simulations relative to ERA5. Notably, the direct effect of ENSO+ on CWP extreme frequencies in CESM simulations exhibits some inconsistencies when compared to ERA5 over Northern Africa (Fig. S3). Better understanding and confirming the influence of climate modes on arid regions (e.g., Northern Africa), where CWP events may be less intense than in other areas, can support adaptation and mitigation policies. While CWP extremes can serve as an important source of freshwater e.g.,, they also present a significant flood risk e.g.,.

We analyzed event counts aggregated over winter and at the scale of predefined SREX regions, given that high counts of compound extremes at these scales are expected to have negative effects on society. While the 98th percentile has been used in this study to focus on extremes and is relatively well-established in the literature e.g.,, other higher thresholds could have been chosen to consider more intense extreme events e.g.,. Figures S16–S19 in the Supplement show results from a sensitivity analysis on the influence of variability modes on regional CWP extremes (Metric 1) and spatially compounding events (Metrics 2 and 3) with the 99th and 99.5th percentiles used as thresholds. Although there are some variations in the results compared to those for the 98th percentile, the main conclusions drawn across the different thresholds are broadly consistent for all Metrics. The magnitude of the effects of the combinations are generally consistent across thresholds, and the combinations detected at higher thresholds are generally included among those identified at lower thresholds (Figs. S16–S19). Slight differences may be due to larger sampling uncertainty for higher thresholds limiting the ability to detect significant effects for higher thresholds rather than different physical mechanisms involved for different thresholds. While the sensitivity analyses broadly indicate the robustness of most of our findings, possible relevant differences across thresholds highlight the importance of identifying impact-relevant thresholds, though this task is challenging e.g.,. In addition, the selected SREX regions may not reflect the natural spatial patterns of variation of CWP extremes as they potentially occur at a more localized scale or may span across multiple regions.

We also note that variability modes such as AMV or ENSO can have lagged effects on regional climate extremes e.g.,, an aspect that has been partially taken into account in this study for ENSO. Furthermore, additional variability modes than those considered here e.g., the Indian Ocean Dipole (IOD), may be related to CWP frequencies in some regions, especially in the regions where the different combinations of the variability modes (NAO, ENSO, PNA, and AMV) had no significant effects (in Fig. ). In addition, the IOD is known to co-vary with some of the modes examined here, e.g. with ENSO and the AMV e.g.,. Overall, although our aggregation in time and space may not be optimal for providing a fine-grained analysis of CWP events and additional modes may be relevant in some regions, this study provides a first comprehensive assessment of the interactions between multiple climate variability modes and the frequency of wintertime CWP extremes and associated spatially compounding events across regions of the Northern Hemisphere. In particular, analyzing spatially compound events highlights the potential influence of variability modes that can link distant regions. These long-range relationships modulated by mode combinations can be explored in more detail, for example, with tailored experiments such as nudged atmospheric simulations.

We mainly focused on the increase in the frequency of CWP extremes, however we note that increases in intensity at the same frequency, as well as low frequencies of CWP extremes may also be relevant information for the insurance industry e.g.,. Studying the influence of variability modes on the intensity of the drivers (here, wind and precipitation) could further support mitigation and adaptation strategies e.g.,. Quantifying these effects on wind and precipitation in isolation would allow for a decomposition of the effects of concurrent variability modes, giving an even more complete picture of the processes involved e.g.,. In addition, regarding spatial metrics, considering other relevant measures that combine return period information with the exposed population would allow for focusing on concurrent variability modes that regularly expose the population to CWP extremes.

With the methodology used here, it is not possible to thoroughly examine the causality behind the dependence between different variability modes and between the modes and CWP extremes. It is however clear that there exist important interactions and causal links between oceanic and atmospheric variability modes at different time scales, e.g., between ENSO and the NAO e.g.,, between the AMV and the NAO e.g.,, and between ENSO and the PNA . Therefore, these modes are dependent on each other, a relevant piece of information to take into account when investigating the causal effects of climate variability modes on characteristics of spatially compounding events. Applying statistical methods such as regression techniques e.g., or more advanced approaches such as causal networks e.g., may help to shed light on the complete causal pathway leading to spatially compounding events and better control potential confounding effects. We also note that dependencies among modes influence the return periods of concurrent modes (that is, certain signs of modes will co-occur with others more frequently than others, such as e.g. ENSO+NAO-), with an effect on their associated CWP extreme frequencies. In line with heterogeneous dependencies among different modes, we estimated a wide range of return periods for the different mode combinations (Fig. ).

A natural continuation of this work is the application of the methodology developed in this study to investigate (1) the influence of climate modes on CWP extremes in the Southern Hemisphere, (2) summertime CWP extremes, and (3) changes in CWP extremes over time, as well as their spatial relationships and dependencies with climate variability modes under climate change .

5Conclusions

In this study, we present the first comprehensive assessment of the relation between individual and concurrent climate variability modes and the frequency of wintertime CWP extremes and associated spatially compounding events across regions of the Northern Hemisphere using CESM simulations. We show that simulated concurrent modes are associated with an amplification of CWP extreme frequencies in many individual regions compared to variability modes in isolation. We have identified groups of regions with positive spatial correlations between regional wintertime CWP extremes. These correlations enable extreme spatially compounding events with many regions experiencing CWP extremes during the same winter. We found that such extreme spatially compounding events, which also include a large fraction of the population of the Northern Hemisphere under CWP extremes in the same winter, are possible due to the combinations of multiple variability modes in anomalous phases and their influence on the atmospheric circulation. For example, the combination of ENSO- and NAO+ nearly double the number of affected regions compared to neutral conditions on average. Among the modes, ENSO is found to be the most influential variability responsible for spatially compounding extreme events, which aligns with its known effects on weather and climate extremes worldwide . The high return periods associated with some of the concurrent modes that lead to extreme spatially compounding events make it potentially difficult to factor these findings into long-term planning, e.g., infrastructure development or (re)insurance modeling. However, by identifying the drivers of the most extreme events, the findings raise awareness on the potential for extreme compounding events under certain mode combinations, which could be factored into weekly, sub-seasonal , seasonal , and longer forecasts to better anticipate mitigation action and climate services .

By using two different metrics to characterize spatially compounding events, we highlight that the effects of concurrent variability modes can differ from one spatial metric to another. While combinations with ENSO+ lead to the largest number of affected regions, when weighted by population exposure, combinations with ENSO- lead to higher effects on population. Our analysis thus stresses the importance of considering not only the interplay between variability modes but also a careful choice of metrics, which should be tailored to the ultimate impacts of interest, so as to assess the relevant characteristics of spatially compounding events and improve their risk management and mitigation.

Data availability

The ERA5 reanalysis data can be accessed via the “Climate Data Store” (CDS) web portal 10.24381/cds.adbb2d47 . The CESM Large Ensemble Simulations can be downloaded from the CESM Large Ensemble Community Project website: https://rda.ucar.edu/datasets/d651014/ . Modes of variability for CESM can be accessed via the CVDP-LE Data Repository 10.5065/h7c7-f961 .

The supplement related to this article is available online at https://doi.org/10.5194/esd-16-1029-2025-supplement.

Author contributions

EB designed the initial plan of the study with DIVD. EB, LSG, and DIVD supervised the project. EB, LSG, BF, and KT designed the experiments and the statistical analyses with inputs from all co-authors. EB, KT, and LSG provided the code for data pre-processing. BF and KT made all computations and figures. BF, LB, RL, KT, and LGS made the analyses and interpretations with inputs, corrections and additional writing contributions from EB, LSG, and DIVD. BF wrote the first draft of the article with inputs from all co-authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Special issue statement

This article is part of the special issue “Methodological innovations for the analysis and management of compound risk and multi-risk, including climate-related and geophysical hazards (NHESS/ESD/ESSD/GC/HESS inter-journal SI)”. It is not associated with a conference.

Acknowledgements

Emanuele Bevacqua received funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) via the Emmy Noether Programme (grant ID 524780515). Laura Suarez-Gutierrez received funding from the European Union's Horizon Europe Framework Programme under the Marie Skłodowska-Curie grant agreement no. 101064940. The work by Luis Gimeno-Sotelo is partially financed by national funds through FCT – Fundação para a Ciência e a Tecnologia under the project UIDB/00006/2020 (DOI: 10.54499/UIDB/00006/2020) and project UID/00006/2025. This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement no. 101003469. This work emerged from the Training School on Dynamical Modeling of Compound Events organized by the European COST Action DAMOCLES (CA17109). This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement no. 847456 and no. 101003469 (XAIDA)). Support from the Swiss National Science Foundation through project PP00P2_198896 to Daniela I. V. Domeisen is gratefully acknowledged. We thank the NCAR's Climate Analysis Section for producing and making available the Climate Variability Diagnostics Package and the CESM simulations.

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 524780515), the HORIZON EUROPE Marie Skłodowska-Curie Actions (grant no. 101064940), and Horizon 2020 (grant nos. 101003469 and 847456).The article processing charges for this open-access publication were covered by the Helmholtz Centre for Environmental Research – UFZ.

Review statement

This paper was edited by Anais Couasnon and reviewed by four anonymous referees.

References Abeysirigunawardena et al.(2009)Abeysirigunawardena, Gilleland, Bronaugh, and Wong

Abeysirigunawardena, D. S., Gilleland, E., Bronaugh, D., and Wong, P.: Extreme wind regime responses to climate variability and change in the inner south coast of British Columbia, Canada, Atmos. Ocean, 47, 41–62, 10.3137/AO1003.2009, 2009.

Ashok et al.(2003)Ashok, Guan, and Yamagata

Ashok, K., Guan, Z., and Yamagata, T.: A Look at the Relationship between the ENSO and the Indian Ocean Dipole, J. Meteorol. Soc. Jpn. Ser. II, 81, 41–56, 10.2151/jmsj.81.41, 2003.

Berdugo et al.(2020)Berdugo, Delgado-Baquerizo, Soliveres, Hernández-Clemente, Zhao, Gaitán, Gross, Saiz, Maire, Lehmann, Rillig, Solé, and Maestre

Berdugo, M., Delgado-Baquerizo, M., Soliveres, S., Hernández-Clemente, R., Zhao, Y., Gaitán, J. J., Gross, N., Saiz, H., Maire, V., Lehmann, A., Rillig, M. C., Solé, R. V., and Maestre, F. T.: Global ecosystem thresholds driven by aridity, Science, 367, 787–790, 2020.

Bevacqua et al.(2021a)Bevacqua, De Michele, Manning, Couasnon, Ribeiro, Ramos, Vignotto, Bastos, Blesić, Durante, Hillier, Oliveira, Pinto, Ragno, Rivoire, Saunders, van der Wiel, Wu, Zhang, and Zscheischler

Bevacqua, E., De Michele, C., Manning, C., Couasnon, A., Ribeiro, A. F. S., Ramos, A. M., Vignotto, E., Bastos, A., Blesić, S., Durante, F., Hillier, J., Oliveira, S. C., Pinto, J. G., Ragno, E., Rivoire, P., Saunders, K., van der Wiel, K., Wu, W., Zhang, T., and Zscheischler, J.: Guidelines for Studying Diverse Types of Compound Weather and Climate Events, Earths Future, 9, e2021EF002340, 10.1029/2021EF002340, 2021a.

Bevacqua et al.(2021b)Bevacqua, Shepherd, Watson, Sparrow, Wallom, and Mitchell

Bevacqua, E., Shepherd, T. G., Watson, P. A. G., Sparrow, S., Wallom, D., and Mitchell, D.: Larger Spatial Footprint of Wintertime Total Precipitation Extremes in a Warmer Climate, Geophys. Res. Lett., 48, e2020GL091990, 10.1029/2020GL091990, 2021b.

Bevacqua et al.(2022)Bevacqua, Zappa, Lehner, and Zscheischler

Bevacqua, E., Zappa, G., Lehner, F., and Zscheischler, J.: Precipitation trends determine future occurrences of compound hot–dry events, Nat. Clim. Change, 12, 350–355, 10.1038/s41558-022-01309-5, 2022.

Bevacqua et al.(2023)Bevacqua, Suarez-Gutierrez, Jézéquel, Lehner, Vrac, Yiou, and Zscheischler

Bevacqua, E., Suarez-Gutierrez, L., Jézéquel, A., Lehner, F., Vrac, M., Yiou, P., and Zscheischler, J.: Advancing research on compound weather and climate events via large ensemble model simulations, Nat. Commun., 14, 2145, 10.1038/s41467-023-37847-5, 2023.

Bevacqua et al.(2024)Bevacqua, Rakovec, Schumacher, Kumar, Thober, Samaniego, Seneviratne, and Zscheischler

Bevacqua, E., Rakovec, O., Schumacher, D. L., Kumar, R., Thober, S., Samaniego, L., Seneviratne, S. I., and Zscheischler, J.: Direct and lagged climate change effects intensified the 2022 European drought, Nat. Geosci., 17, 1100–1107, 2024.

Blocken and Carmeliet(2004)

Blocken, B. and Carmeliet, J.: A review of wind-driven rain research in building science, J. Wind Eng. Ind. Aerod., 92, 1079–1130, 10.1016/j.jweia.2004.06.003, 2004.

Bloomfield et al.(2023)Bloomfield, J., A., A. L., L. C., F., R., D., A., and P. D.

Bloomfield, H. C., Hillier, J., Griffin, A., Kay, A. L., Shaffrey, L. C., Pianosi, F., James, R., Kumar, D., Champion, A., and Bates, P. D.: Co-occurring wintertime flooding and extreme wind over Europe, from daily to seasonal timescales, Weather Clim. Extremes, 39, 100550, 10.1016/j.wace.2023.100550, 2023.

Bloomfield et al.(2024)Bloomfield, Bates, Shaffrey, Hillier, Champion, Cotterill, Pope, and Kumar

Bloomfield, H. C., Bates, P., Shaffrey, L. C., Hillier, J., Champion, A., Cotterill, D., Pope, J. O., and Kumar, D.: Synoptic conditions conducive for compound wind-flood events in Great Britain in present and future climates, Environ. Res. Lett., 19, 024019, 10.1088/1748-9326/ad1cb7, 2024.

Bonferroni(1936)

Bonferroni, C.: Teoria Statistica Delle Classi e Calcolo Delle Probabilita, Pubblicazioni del R Istituto Superiore di Scienze Economiche e Commerciali di Firenze, Libreria Internazionale Seeber, Florence, Italy, 8, 3–62, 1936.

Bradley(1968)

Bradley, J. V.: Distribution-Free Statistical Tests, Prentice-Hall, Englewood Cliffs, NJ, 283–310, 1968.

Brunner et al.(2021)Brunner, Swain, Gilleland, and Wood

Brunner, M. I., Swain, D. L., Gilleland, E., and Wood, A. W.: Increasing importance of temperature as a contributor to the spatial extent of streamflow drought, Environ. Res. Lett., 16, 024038, 10.1088/1748-9326/abd2f0, 2021.

Calafat et al.(2022)Calafat, Wahl, Tadesse, and Sparrow

Calafat, F. M., Wahl, T., Tadesse, M. G., and Sparrow, S. N.: Trends in Europe storm surge extremes match the rate of sea-level rise, Nature, 603, 841–845, 10.1038/s41586-022-04426-5, 2022.

Camberlin et al.(2001)Camberlin, Janicot, and Poccard

Camberlin, P., Janicot, S., and Poccard, I.: Seasonality and atmospheric dynamics of the teleconnection between African rainfall and tropical sea-surface temperature: Atlantic vs. ENSO, Int. J. Climatol., 21, 973–1005, 10.1002/joc.673, 2001.

Camuffo et al.(2020)Camuffo, Becherini, and della Valle

Camuffo, D., Becherini, F., and della Valle, A.: Relationship between selected percentiles and return periods of extreme events, Acta Geophys., 68, 1201–1211, 10.1007/s11600-020-00452-x, 2020.

Cerveny and Newman(2000)

Cerveny, R. S. and Newman, L. E.: Climatological Relationships between Tropical Cyclones and Rainfall, Mon. Weather Rev., 128, 3329–3336, 10.1175/1520-0493(2000)128<3329:CRBTCA>2.0.CO;2, 2000.

Chongyin and Mingquan(2001)

Chongyin, L. and Mingquan, M.: The Influence of the Indian Ocean Dipole on Atmospheric Circulation and Climate, Adv. Atmos. Sci., 18, 831–843, 10.1007/BF03403506, 2001.

CIESIN(2016)

CIESIN (Columbia University, Center for International Earth Science Information Network): Gridded Population of the World, Version 4 (GPWv4): Population Density, NASA Socioeconomic Data and Applications Center (SEDAC), Palisades, NY, 10.7927/H49C6VHW, 2016.

Ciullo et al.(2023)Ciullo, Strobl, Meiler, Martius, and Bresch

Ciullo, A., Strobl, E., Meiler, S., Martius, O., and Bresch, D. N.: Increasing countries' financial resilience through global catastrophe risk pooling, Nat. Commun., 14, 922, 10.1038/s41467-023-36539-4, 2023.

Coles et al.(1999)Coles, Heffernan, and Tawn

Coles, S., Heffernan, J., and Tawn, J.: Dependence Measures for Extreme Value Analyses, Extremes, 2, 339–365, 10.1023/A:1009963131610, 1999.

Copernicus Climate Change Service, Climate Data Store(2023)

Copernicus Climate Change Service, Climate Data Store: ERA5 hourly data on single levels from 1940 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], 10.24381/cds.adbb2d47 (last access: 19 July 2022), 2023.

Davini et al.(2015)Davini, Hardenberg, and Corti

Davini, P., Hardenberg, J., and Corti, S.: Tropical origin for the impacts of the Atlantic Multidecadal Variability on the Euro-Atlantic climate, Environ. Res. Lett., 10, 094010, 10.1088/1748-9326/10/9/094010, 2015.

DelSole et al.(2017)DelSole, Trenary, Tippett, and Pegion

DelSole, T., Trenary, L., Tippett, M. K., and Pegion, K.: Predictability of Week-3–4 Average Temperature and Precipitation over the Contiguous United States, J. Climate, 30, 3499–3512, 10.1175/JCLI-D-16-0567.1, 2017.

Deser et al.(2017)Deser, Simpson, McKinnon, and Phillips

Deser, C., Simpson, I. R., McKinnon, K. A., and Phillips, A. S.: The Northern Hemisphere Extratropical Atmospheric Circulation Response to ENSO: How Well Do We Know It and How Do We Evaluate Models Accordingly?, J. Climate, 30, 5059–5082, 10.1175/JCLI-D-16-0844.1, 2017.

Domeisen et al.(2019)Domeisen, Garfinkel, and Butler

Domeisen, D. I., Garfinkel, C. I., and Butler, A. H.: The teleconnection of El Niño Southern Oscillation to the stratosphere, Rev. Geophys., 57, 5–47, 2019.

Dowdy and Catto(2017)

Dowdy, A. and Catto, J.: Extreme weather caused by concurrent cyclone, front and thunderstorm occurrences, Sci. Rep.-UK, 7, 40359, 10.1038/srep40359, 2017.

Elsner et al.(2001)Elsner, Bossak, and Niu

Elsner, J. B., Bossak, B. H., and Niu, X.-F.: Secular changes to the ENSO-U. S. hurricane relationship, Geophys. Res. Lett., 28, 4123–4126, 10.1029/2001GL013669, 2001.

Fang et al.(2018)Fang, Chen, Ilvonen, Frank, Pasanen, Holmström, Zhao, Zhang, and Seppä

Fang, K., Chen, D., Ilvonen, L., Frank, D., Pasanen, L., Holmström, L., Zhao, Y., Zhang, P., and Seppä, H.: Time-varying relationships among oceanic and atmospheric modes: A turning point at around 1940, Quatern. Int., 487, 12–25, 10.1016/j.quaint.2017.09.005, 2018.

Fischer et al.(2023)Fischer, Beyerle, Bloin-Wibe, Gessner, Humphrey, Lehner, Pendergrass, Sippel, Zeder, and Knutti

Fischer, E. M., Beyerle, U., Bloin-Wibe, L., Gessner, C., Humphrey, V., Lehner, F., Pendergrass, A. G., Sippel, S., Zeder, J., and Knutti, R.: Storylines for unprecedented heatwaves based on ensemble boosting, Nat. Commun., 14, 4643, 10.1038/s41467-023-40112-4, 2023.

Foulds and Warburton(2007)

Foulds, S. and Warburton, J.: Significance of wind-driven rain (wind-splash) in the erosion of blanket peat, Geomorphology, 83, 183–192, 10.1016/j.geomorph.2006.07.001, 2007.

Gampe et al.(2024)Gampe, Schwingshackl, Böhnisch, Mittermeier, Sandstad, and Wood

Gampe, D., Schwingshackl, C., Böhnisch, A., Mittermeier, M., Sandstad, M., and Wood, R. R.: Applying global warming levels of emergence to highlight the increasing population exposure to temperature and precipitation extremes, Earth Syst. Dynam., 15, 589–605, 10.5194/esd-15-589-2024, 2024.

Gao et al.(2022)Gao, Xu, Wang, Sun, Xie, and Cao

Gao, T., Xu, Y., Wang, H. J., Sun, Q., Xie, L., and Cao, F.: Combined Impacts of Climate Variability Modes on Seasonal Precipitation Extremes Over China, Water Resources Manag., 36, 2411–2431, 10.1007/s11269-022-03150-z, 2022.

Goddard and Gershunov(2020)

Goddard, L. and Gershunov, A.: Impact of El Niño on weather and climate extremes, El Niño Southern Oscillation in a changing climate, 361–375, American Geophysical Union, 10.1002/9781119548164, 2020.

Good(2013)

Good, P.: Permutation Tests: A Practical Guide to Resampling Methods for Testing Hypotheses, Springer Science & Business Media, 10.1007/978-1-4757-2346-5, 2013.

Greeves et al.(2007)Greeves, Pope, Stratton, and Martin

Greeves, C., Pope, V., Stratton, R., and Martin, G.: Representation of Northern Hemisphere winter storm tracks in climate models, Clim. Dynam., 28, 683–702, 10.1007/s00382-006-0205-x, 2007.

Grimm(2011)

Grimm, A.: Interannual Climate Variability in South America: Impacts on Seasonal Precipitation, Extreme Events, and Possible Effects of Climate Change, Stoch. Env. Res. Risk A., 25, 537–554, 10.1007/s00477-010-0420-1, 2011.

Hansen et al.(2019)Hansen, Kruschke, Greatbatch, and Weisheimer

Hansen, F., Kruschke, T., Greatbatch, R. J., and Weisheimer, A.: Factors Influencing the Seasonal Predictability of Northern Hemisphere Severe Winter Storms, Geophys. Res. Lett., 46, 365–373, 10.1029/2018GL079415, 2019.

Hersbach et al.(2020)Hersbach, Bell, Berrisford, Hirahara, Horányi, Muñoz-Sabater, Nicolas, Peubey, Radu, Schepers, Simmons, Soci, Abdalla, Abellan, Balsamo, Bechtold, Biavati, Bidlot, Bonavita, De Chiara, Dahlgren, Dee, Diamantakis, Dragani, Flemming, Forbes, Fuentes, Geer, Haimberger, Healy, Hogan, Hólm, Janisková, Keeley, Laloyaux, Lopez, Lupu, Radnoti, de Rosnay, Rozum, Vamborg, Villaume, and Thépaut

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, 10.1002/qj.3803, 2020.

Hillier et al.(2020)Hillier, Matthews, Wilby, and Murphy

Hillier, J., Matthews, T., Wilby, R., and Murphy, C.: Multi-hazard dependencies can increase or decrease risk, Nat. Clim. Change, 10, 1–4, 10.1038/s41558-020-0832-y, 2020.

Hillier et al.(2015)Hillier, Macdonald, Leckebusch, and Stavrinides

Hillier, J. K., Macdonald, N., Leckebusch, G. C., and Stavrinides, A.: Interactions between apparently “primary” weather-driven hazards and their cost, Environ. Res. Lett., 10, 104003, 10.1088/1748-9326/10/10/104003, 2015.

Hong Lee et al.(2023)Hong Lee, Julien, and Lee

Hong Lee, J., Julien, P. Y., and Lee, S.: Teleconnection of ENSO extreme events and precipitation variability over the United States, J. Hydrol., 619, 129206, 10.1016/j.jhydrol.2023.129206, 2023.

Hu et al.(2023)Hu, Zhang, Jin, Hong, Jiang, and Stuecker

Hu, S., Zhang, W., Jin, F.-F., Hong, L.-C., Jiang, F., and Stuecker, M. F.: Seasonal dependence of the Pacific–North American teleconnection associated with ENSO and its interaction with the annual cycle, J. Climate, 36, 7061–7072, 2023.

Hurrell(1995)

Hurrell, J. W.: Decadal trends in the North Atlantic oscillation: Regional temperatures and precipitation, Science, 269, 676–679, 10.1126/science.269.5224.676, 1995.

Hurrell and Deser(2010)

Hurrell, J. W. and Deser, C.: North Atlantic climate variability: The role of the North Atlantic Oscillation, J. Marine Syst., 79, 231–244, 10.1016/j.jmarsys.2009.11.002, 2010.

Hurrell et al.(2003)Hurrell, Kushnir, Ottersen, and Visbeck

Hurrell, J. W., Kushnir, Y., Ottersen, G., and Visbeck, M.: An Overview of the North Atlantic Oscillation, American Geophysical Union (AGU), 1–35, 10.1029/134GM01, 2003.

IPCC(2012)

IPCC: Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change, edited by: Field, C. B., Barros, V., Stocker, T. F., Qin, D., Dokken, D. J., Ebi, K. L., Mastrandrea, M. D., Mach, K. J., Plattner, G.-K., Allen, S. K., Tignor, M., and Midgley, P. M., Cambridge University Press, Cambridge, UK, and New York, NY, USA, 582 pp., 2012.

IPCC(2021)

IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, in press, 10.1017/9781009157896, 2021.

Iturbide et al.(2020)Iturbide, Gutiérrez, Alves, Bedia, Cerezo-Mota, Cimadevilla, Cofiño, Di Luca, Faria, Gorodetskaya, Hauser, Herrera, Hennessy, Hewitt, Jones, Krakovska, Manzanas, Martínez-Castro, Narisma, Nurhati, Pinto, Seneviratne, van den Hurk, and Vera

Iturbide, M., Gutiérrez, J. M., Alves, L. M., Bedia, J., Cerezo-Mota, R., Cimadevilla, E., Cofiño, A. S., Di Luca, A., Faria, S. H., Gorodetskaya, I. V., Hauser, M., Herrera, S., Hennessy, K., Hewitt, H. T., Jones, R. G., Krakovska, S., Manzanas, R., Martínez-Castro, D., Narisma, G. T., Nurhati, I. S., Pinto, I., Seneviratne, S. I., van den Hurk, B., and Vera, C. S.: An update of IPCC climate reference regions for subcontinental analysis of climate model data: definition and aggregated datasets, Earth Syst. Sci. Data, 12, 2959–2970, 10.5194/essd-12-2959-2020, 2020.

Jeong et al.(2020)Jeong, Cannon, and Morris

Jeong, D. I., Cannon, A. J., and Morris, R. J.: Projected changes to wind loads coinciding with rainfall for building design in Canada based on an ensemble of Canadian regional climate model simulations, Climatic Change, 162, 821–835, 10.1007/s10584-020-02745-y, 2020.

Jiang et al.(2024)Jiang, Tarasova, Yu, and Zscheischler

Jiang, S., Tarasova, L., Yu, G., and Zscheischler, J.: Compounding effects in flood drivers challenge estimates of extreme river floods, Science Advances, 10, eadl4005, 10.1126/sciadv.adl4005, 2024.

Jiménez-Esteve and Domeisen(2019)

Jiménez-Esteve, B. and Domeisen, D.: Nonlinearity in the North Pacific atmospheric response to a linear ENSO forcing, Geophys. Res. Lett., 46, 2271–2281, 2019.

Jiménez-Esteve and Domeisen(2020)

Jiménez-Esteve, B. and Domeisen, D. I. V.: Nonlinearity in the tropospheric pathway of ENSO to the North Atlantic, Weather Clim. Dynam., 1, 225–245, 10.5194/wcd-1-225-2020, 2020.

Jongman et al.(2014)Jongman, Hochrainer-Stigler, Feyen, CJH, Mechler, Bouwer, Pflug, Rojas, and Ward

Jongman, B., Hochrainer-Stigler, S., Feyen, L., CJH, A., Mechler, R., Bouwer, L., Pflug, G., Rojas, R., and Ward, P.: Increasing Stress on Disaster-Risk Finance due to Large Floods, Nat. Clim. Change, 4, 264–268, 10.1038/NCLIMATE2124, 2014.

Kawamura et al.(2004)Kawamura, Suppiah, Collier, and Gordon

Kawamura, R., Suppiah, R., Collier, M. A., and Gordon, H. B.: Lagged relationships between ENSO and the Asian Summer Monsoon in the CSIRO coupled model, Geophys. Res. Lett., 31, L23205, 10.1029/2004GL021411, 2004.

Kay et al.(2015)Kay, Deser, Phillips, Mai, Hannay, Strand, Arblaster, Bates, Danabasoglu, Edwards, Holland, Kushner, Lamarque, Lawrence, Lindsay, Middleton, Munoz, Neale, Oleson, Polvani, and Vertenstein

Kay, J. E., Deser, C., Phillips, A., Mai, A., Hannay, C., Strand, G., Arblaster, J. M., Bates, S. C., Danabasoglu, G., Edwards, J., Holland, M., Kushner, P., Lamarque, J.-F., Lawrence, D., Lindsay, K., Middleton, A., Munoz, E., Neale, R., Oleson, K., Polvani, L., and Vertenstein, M.: The Community Earth System Model (CESM) Large Ensemble Project: A Community Resource for Studying Climate Change in the Presence of Internal Climate Variability, B. Am. Meteorol. Soc., 96, 1333–1349, 10.1175/BAMS-D-13-00255.1, 2015.

Kelder et al.(2022)Kelder, Wanders, van der Wiel, Marjoribanks, Slater, l Wilby, and Prudhomme

Kelder, T., Wanders, N., van der Wiel, K., Marjoribanks, T. I., Slater, L. J., Wilby, R. L., and Prudhomme, C.: Interpreting extreme climate impacts from large ensemble simulations – are they unseen or unrealistic?, Environ. Res. Lett., 17, 044052, 10.1088/1748-9326/ac5cf4, 2022.

Khouakhi et al.(2016)Khouakhi, Villarini, and Vecchi

Khouakhi, A., Villarini, G., and Vecchi, G.: Contribution of Tropical Cyclones to Rainfall at the Global Scale, J. Climate, 30, 359–372, 10.1175/JCLI-D-16-0298.1, 2016.

Kim et al.(2019)Kim, Iizumi, and Nishimori

Kim, W., Iizumi, T., and Nishimori, M.: Global Patterns of Crop Production Losses Associated with Droughts from 1983 to 2009, J. Appl. Meteorol. Clim., 58, 1233–1244, 10.1175/JAMC-D-18-0174.1, 2019.

King et al.(2023)King, Keenlyside, and Li

King, M., Keenlyside, N., and Li, C.: ENSO teleconnections in terms of non-NAO and NAO atmospheric variability, Clim. Dynam., 61, 27172733, 10.1007/s00382-023-06697-8, 2023.

Kirov and Georgieva(2002)

Kirov, B. and Georgieva, K.: Long-term variations and interrelations of ENSO, NAO and solar activity, Phys. Chem. Earth Pt. A/B/C, 27, 441–448, 10.1016/S1474-7065(02)00024-4, 2002.

Klawa and Ulbrich(2003)

Klawa, M. and Ulbrich, U.: A model for the estimation of storm losses and the identification of severe winter storms in Germany, Nat. Hazards Earth Syst. Sci., 3, 725–732, 10.5194/nhess-3-725-2003, 2003.

Knight et al.(2006)Knight, Folland, and Scaife

Knight, J. R., Folland, C. K., and Scaife, A. A.: Climate impacts of the Atlantic Multidecadal Oscillation, Geophys. Res. Lett., 33, L17706, 10.1029/2006GL026242, 2006.

Kossin et al.(2010)Kossin, Camargo, and Sitkowski

Kossin, J. P., Camargo, S. J., and Sitkowski, M.: Climate Modulation of North Atlantic Hurricane Tracks, J. Climate, 23, 3057–3076, 10.1175/2010JCLI3497.1, 2010.

Kretschmer et al.(2021)Kretschmer, Adams, Arribas, Prudden, Robinson, Saggioro, and Shepherd

Kretschmer, M., Adams, S. V., Arribas, A., Prudden, R., Robinson, N., Saggioro, E., and Shepherd, T. G.: Quantifying Causal Pathways of Teleconnections, B. Am. Meteorol. Soc., 102, E2247–E2263, 10.1175/BAMS-D-20-0117.1, 2021.

Kurniadi et al.(2021)Kurniadi, Weller, Min, and Seong

Kurniadi, A., Weller, E., Min, S.-K., and Seong, M.-G.: Independent ENSO and IOD impacts on rainfall extremes over Indonesia, Int. J. Climatol., 41, 3640–3656, 10.1002/joc.7040, 2021.

Lehner(2024)

Lehner, F.: CLIVAR LE project, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory [data set], https://rda.ucar.edu/datasets/d651014/ (last access: 19 July 2022), 2024.

Lenssen et al.(2020)Lenssen, Goddard, and Mason

Lenssen, N. J., Goddard, L., and Mason, S.: Seasonal forecast skill of ENSO teleconnection maps, Weather Forecast., 35, 2387–2406, 2020.

Li et al.(2022)Li, Chen, Messmer, Zhu, Feng, Yin, and Bevacqua

Li, D., Chen, Y., Messmer, M., Zhu, Y., Feng, J., Yin, B., and Bevacqua, E.: Compound Wind and Precipitation Extremes Across the Indo-Pacific: Climatology, Variability, and Drivers, Geophys. Res. Lett., 49, e2022GL098594, 10.1029/2022GL098594, 2022.

Li et al.(2024)Li, Zscheischler, Chen, Yin, Feng, Freund, Qi, Zhu, and Bevacqua

Li, D., Zscheischler, J., Chen, Y., Yin, B., Feng, J., Freund, M., Qi, J., Zhu, Y., and Bevacqua, E.: Intensification and Poleward Shift of Compound Wind and Precipitation Extremes in a Warmer Climate, Geophys. Res. Lett., 51, e2024GL110135, 10.1029/2024GL110135, 2024.

Li et al.(2011)Li, Zhai, and Cai

Li, W., Zhai, P., and Cai, J.: Research on the Relationship of ENSO and the Frequency of Extreme Precipitation Events in China, Advances in Climate Change Research, 2, 101–107, 10.3724/SP.J.1248.2011.00101, 2011.

Li et al.(2017)Li, Liu, Mauzerall, Li, Fan, Horowitz, He, Yi, and Tao

Li, Z., Liu, J., Mauzerall, D. L., Li, X., Fan, S., Horowitz, L. W., He, C., Yi, K., and Tao, S.: A potential large and persistent black carbon forcing over Northern Pacific inferred from satellite observations, Sci. Rep.-UK, 7, 43429, 10.1038/srep43429, 2017.

Liberato(2014)

Liberato, M. L.: The 19 January 2013 windstorm over the North Atlantic: large-scale dynamics and impacts on Iberia, Weather. Clim. Extremes, 5–6, 16–28, 10.1016/j.wace.2014.06.002, 2014.

Liu et al.(2013)Liu, Chen, Chen, Lian, and Wu

Liu, B., Chen, J., Chen, X., Lian, Y., and Wu, L.: Uncertainty in determining extreme precipitation thresholds, J. Hydrol., 503, 233–245, 10.1016/j.jhydrol.2013.09.002, 2013.

Liu et al.(2024)Liu, Li, Chen, Feng, Qi, and Yin

Liu, H., Li, D., Chen, Q., Feng, J., Qi, J., and Yin, B.: The multiscale variability of global extreme wind and wave events and their relationships with climate modes, Ocean Eng., 307, 118239, 10.1016/j.oceaneng.2024.118239, 2024.

Lodise et al.(2022)Lodise, Merrifield, Collins, Rogowski, Behrens, and Terrill

Lodise, J., Merrifield, S., Collins, C., Rogowski, P., Behrens, J., and Terrill, E.: Global Climatology of Extratropical Cyclones From a New Tracking Approach and Associated Wave Heights From Satellite Radar Altimeter, J. Geophys. Res.-Oceans, 127, e2022JC018925, 10.1029/2022JC018925, 2022.

Luca et al.(2017)Luca, Hillier, Wilby, Quinn, and Harrigan

Luca, P. D., Hillier, J. K., Wilby, R. L., Quinn, N. W., and Harrigan, S.: Extreme multi-basin flooding linked with extra-tropical cyclones, Environ. Res. Lett., 12, 114009, 10.1088/1748-9326/aa868e, 2017.

Maher et al.(2021)Maher, Milinski, and Ludwig

Maher, N., Milinski, S., and Ludwig, R.: Large ensemble climate model simulations: introduction, overview, and future prospects for utilising multiple types of large ensemble, Earth Syst. Dynam., 12, 401–418, 10.5194/esd-12-401-2021, 2021.

Manning et al.(2018)Manning, Widmann, Bevacqua, Loon, Maraun, and Vrac

Manning, C., Widmann, M., Bevacqua, E., Loon, A. F. V., Maraun, D., and Vrac, M.: Soil Moisture Drought in Europe: A Compound Event of Precipitation and Potential Evapotranspiration on Multiple Time Scales, J. Hydrometeorol., 19, 1255–1271, 10.1175/JHM-D-18-0017.1, 2018.

Manning et al.(2024)Manning, Kendon, Fowler, Catto, Chan, and Sansom

Manning, C., Kendon, E. J., Fowler, H. J., Catto, J. L., Chan, S. C., and Sansom, P. G.: Compound wind and rainfall extremes: Drivers and future changes over the UK and Ireland, Weather Clim. Extremes, 44, 100673, 10.1016/j.wace.2024.100673, 2024.

Maraun(2016)

Maraun, D.: Bias Correcting Climate Change Simulations – A Critical Review, Curr. Clim. Chang. Reports, 2, 211–220, 10.1007/s40641-016-0050-x, 2016.

Martius et al.(2016)Martius, Pfahl, and Chevalier

Martius, O., Pfahl, S., and Chevalier, C.: A global quantification of compound precipitation and wind extremes, Geophys. Res. Lett., 43, 7709–7717, 10.1002/2016GL070017, 2016.

Matthews et al.(2014)Matthews, Murphy, Wilby, and Harrigan

Matthews, T., Murphy, C., Wilby, R. L., and Harrigan, S.: Stormiest winter on record for Ireland and UK, Nat. Clim. Change, 4, 738–740, 10.1038/nclimate2336, 2014.

Messmer and Simmonds(2021)

Messmer, M. and Simmonds, I.: Global analysis of cyclone-induced compound precipitation and wind extreme events, Weather Clim. Extremes, 32, 100324, 10.1016/j.wace.2021.100324, 2021.

Mirrahimi et al.(2015)Mirrahimi, Lim, and Surat

Mirrahimi, S., Lim, C. H., and Surat, S.: Review of method to estimation of wind-driven rain on building façade, Advances in Environmental Biology, 9, 18–23, 2015.

Müller and Roeckner(2006)

Müller, W. A. and Roeckner, E.: ENSO impact on midlatitude circulation patterns in future climate change projections, Geophys. Res. Lett., 33, L05711, 10.1029/2005GL025032, 2006.

Nowack et al.(2020)Nowack, Runge, Eyring, and Haigh

Nowack, P., Runge, J., Eyring, V., and Haigh, J. D.: Causal networks for climate model evaluation and constrained projections, Nat. Commun., 11, 1415, 10.1038/s41467-020-15195-y, 2020.

O'Reilly et al.(2024)O'Reilly, Drouard, Ayarzagüena, Ambaum, and Methven

O'Reilly, C. H., Drouard, M., Ayarzagüena, B., Ambaum, M. H. P., and Methven, J.: The role of storm-track dynamics in the intraseasonal variability of the winter ENSO teleconnection to the North Atlantic, Q. J. Roy. Meteor. Soc., 150, 2069–2086, 10.1002/qj.4691, 2024.

Osman et al.(2023)Osman, Domeisen, Robertson, and Weisheimer

Osman, M., Domeisen, D., Robertson, A. W., and Weisheimer, A.: Sub-seasonal to decadal predictions in support of climate services, Climate Services, 30, 100397, 10.1016/j.cliser.2023.100397, 2023.

Owen et al.(2021)Owen, Catto, Stephenson, and Dunstone

Owen, L. E., Catto, J. L., Stephenson, D. B., and Dunstone, N. J.: Compound precipitation and wind extremes over Europe and their relationship to extratropical cyclones, Weather Clim. Extremes, 33, 100342, 10.1016/j.wace.2021.100342, 2021.

Pearl(2013)

Pearl, J.: Linear Models: A Useful “Microscope” for Causal Analysis, Journal of Causal Inference, 1, 155–170, 10.1515/jci-2013-0003, 2013.

Phillips et al.(2014)Phillips, Deser, and Fasullo

Phillips, A., Deser, C., and Fasullo, J.: Evaluating Modes of Variability in Climate Models, Eos T. Am. Geophys. Un., 95, 453–455, 10.1002/2014EO490002, 2014.

Phillips et al.(2020)

Phillips, A., Deser, C., Fasullo, J. T., Schneider, D. P., and Simpson, I. R.: Assessing Climate Variability and Change in Model Large Ensembles: A User's Guide to the “Climate Variability Diagnostics Package for Large Ensembles” Version 1.0, NSF NCAR [data set], 10.5065/h7c7-f961, 2020.

Pilorz et al.(2023)Pilorz, Laskowski, Surowiecki, and Łupikasza

Pilorz, W., Laskowski, I., Surowiecki, A., and Łupikasza, E.: Fatalities related to sudden meteorological events across Central Europe from 2010 to 2020, Int. J. Disast. Risk Re., 88, 103622, 10.1016/j.ijdrr.2023.103622, 2023.

Pinto et al.(2009)Pinto, Zacharias, Fink, Leckebusch, and Ulbrich

Pinto, J. G., Zacharias, S., Fink, A. H., Leckebusch, G. C., and Ulbrich, U.: Factors contributing to the development of extreme North Atlantic cyclones and their relationship with the NAO, Clim. Dynam., 32, 711–737, 10.1007/s00382-008-0396-4, 2009.

Priestley et al.(2023)Priestley, Stephenson, Scaife, Bannister, Allen, and Wilkie

Priestley, M. D. K., Stephenson, D. B., Scaife, A. A., Bannister, D., Allen, C. J. T., and Wilkie, D.: Return levels of extreme European windstorms, their dependency on the North Atlantic Oscillation, and potential future risks, Nat. Hazards Earth Syst. Sci., 23, 3845–3861, 10.5194/nhess-23-3845-2023, 2023.

Qian et al.(2023)Qian, Ye, Bevacqua, and Zscheischler

Qian, C., Ye, Y., Bevacqua, E., and Zscheischler, J.: Human influences on spatially compounding flooding and heatwave events in China and future increasing risks, Weather Clim. Extremes, 42, 100616, 10.1016/j.wace.2023.100616, 2023.

Qiu et al.(2014)Qiu, Cai, Guo, and Ng

Qiu, Y., Cai, W., Guo, X., and Ng, B.: The asymmetric influence of the positive and negative IOD events on China's rainfall, Sci. Rep.-UK, 4, 4943, 10.1038/srep04943, 2014.

Raible(2007)

Raible, C. C.: On the relation between extremes of midlatitude cyclones and the atmospheric circulation using ERA40, Geophys. Res. Lett., 34, L07703, 10.1029/2006GL029084, 2007.

Ralph et al.(2006)Ralph, Neiman, Wick, Gutman, Dettinger, Cayan, and White

Ralph, F. M., Neiman, P. J., Wick, G. A., Gutman, S. I., Dettinger, M. D., Cayan, D. R., and White, A. B.: Flooding on California's Russian River: Role of atmospheric rivers, Geophys. Res. Lett., 33, 13, 10.1029/2006GL026689, 2006.

Rappaport(2000)

Rappaport, E.: Loss of Life in the United States Associated with Recent Atlantic Tropical Cyclones, B. Am. Meteorol. Soc., 81, L13801, 10.1175/1520-0477(2000)081<2065:LOLITU>2.3.CO;2, 2000.

Rasmijn et al.(2018)Rasmijn, Schrier, Bintanja, Barkmeijer, Sterl, and Hazeleger

Rasmijn, L., Schrier, G., Bintanja, R., Barkmeijer, J., Sterl, A., and Hazeleger, W.: Future equivalent of 2010 Russian heatwave intensified by weakening soil moisture constraints, Nat. Clim. Change, 8, 1–5, 10.1038/s41558-018-0114-0, 2018.

Raveh-Rubin and Catto(2019)

Raveh-Rubin, S. and Catto, J. L.: Climatology and dynamics of the link between dry intrusions and cold fronts during winter, Part II: Front-centred perspective, Clim. Dynam., 53, 1893–1909, 10.1007/s00382-019-04793-2, 2019.

Raveh-Rubin and Wernli(2015)

Raveh-Rubin, S. and Wernli, H.: Large-scale wind and precipitation extremes in the Mediterranean: a climatological analysis for 1979–2012, Q. J. Roy. Meteor. Soc., 141, 2404–2417, 10.1002/qj.2531, 2015.

Raymond et al.(2022)Raymond, Suarez-Gutierrez, Kornhuber, Pascolini-Campbell, Sillmann, and Waliser

Raymond, C., Suarez-Gutierrez, L., Kornhuber, K., Pascolini-Campbell, M., Sillmann, J., and Waliser, D. E.: Increasing spatiotemporal proximity of heat and precipitation extremes in a warming world quantified by a large model ensemble, Environ. Res. Lett., 17, 035005, 10.1088/1748-9326/ac5712, 2022.

Renwick and Wallace(1996)

Renwick, J. A. and Wallace, J. M.: Relationships between North Pacific Wintertime Blocking, El Niño, and the PNA Pattern, Mon. Weather Rev., 124, 2071–2076, 10.1175/1520-0493(1996)124<2071:RBNPWB>2.0.CO;2, 1996.

Ridder et al.(2022)Ridder, Ukkola, Pitman, and Perkins-Kirkpatrick

Ridder, N., Ukkola, A., Pitman, A., and Perkins-Kirkpatrick, S.: Increased occurrence of high impact compound events under climate change, npj Clim. Atmos. Sci., 5, 2065–2074, 10.1038/s41612-021-00224-4, 2022.

Roberts et al.(2014)Roberts, Champion, Dawkins, Hodges, Shaffrey, Stephenson, Stringer, Thornton, and Youngman

Roberts, J. F., Champion, A. J., Dawkins, L. C., Hodges, K. I., Shaffrey, L. C., Stephenson, D. B., Stringer, M. A., Thornton, H. E., and Youngman, B. D.: The XWS open access catalogue of extreme European windstorms from 1979 to 2012, Nat. Hazards Earth Syst. Sci., 14, 2487–2501, 10.5194/nhess-14-2487-2014, 2014.

Ruiz-Barradas et al.(2000)Ruiz-Barradas, Carton, and Nigam

Ruiz-Barradas, A., Carton, J. A., and Nigam, S.: Structure of Interannual-to-Decadal Climate Variability in the Tropical Atlantic Sector, J. Climate, 13, 3285–3297, 10.1175/1520-0442(2000)013<3285:SOITDC>2.0.CO;2, 2000.

Schär et al.(2016)Schär, Ban, Fischer, Rajczak, Schmidli, Frei, Giorgi, Karl, Kendon, Tank, O'Gorman, Sillmann, Zhang, and Zwiers

Schär, C., Ban, N., Fischer, E. M., Rajczak, J., Schmidli, J., Frei, C., Giorgi, F., Karl, T. R., Kendon, E. J., Tank, A. M. G. K., O'Gorman, P. A., Sillmann, J., Zhang, X., and Zwiers, F. W.: Percentile indices for assessing changes in heavy precipitation events, Climatic Change, 137, 201–216, 10.1007/s10584-016-1669-2, 2016.

Schmidt et al.(2001)Schmidt, Lipp, Rose, and Luther

Schmidt, N., Lipp, E. K., Rose, J. B., and Luther, M. E.: ENSO Influences on Seasonal Rainfall and River Discharge in Florida, J. Climate, 14, 615–628, 10.1175/1520-0442(2001)014<0615:EIOSRA>2.0.CO;2, 2001.

Sedgwick(2014)

Sedgwick, P.: Multiple hypothesis testing and Bonferroni's correction, BMJ, 349, 3, 10.1136/bmj.g6284, 2014.

Seneviratne et al.(2012)Seneviratne, Nicholls, Easterling, Goodess, Kanae, Kossin, Luo, Marengo, McInnes, Rahimi, Reichstein, Sorteberg, Vera, and Zhang

Seneviratne, S., Nicholls, N., Easterling, D., Goodess, C., Kanae, S., Kossin, J., Luo, Y., Marengo, J., McInnes, K., Rahimi, M., Reichstein, M., Sorteberg, A., Vera, C., and Zhang, X.: Changes in climate extremes and their impacts on the natural physical environment, in: Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation, edited by: Field, C. B., Barros, V., Stocker, T. F., Qin, D., Dokken, D. J., Ebi, K. L., Mastrandrea, M. D., Mach, K. J., Plattner, G.-K., Allen, S. K., Tignor, M., and Midgley, P. M., Cambridge University Press, 109–230, 10.1017/CBO9781139177245, 2012.

Singh et al.(2021)Singh, Ashfaq, Skinner, Anderson, and Singh

Singh, J., Ashfaq, M., Skinner, C. B., Anderson, W. B., and Singh, D.: Amplified risk of spatially compounding droughts during co-occurrences of modes of natural ocean variability, npj Clim. Atmos. Sci., 4, 7, 10.1038/s41612-021-00161-2, 2021.

Stuecker et al.(2017)Stuecker, Timmermann, Jin, Chikamoto, Zhang, Wittenberg, Widiasih, and Zhao

Stuecker, M. F., Timmermann, A., Jin, F.-F., Chikamoto, Y., Zhang, W., Wittenberg, A. T., Widiasih, E., and Zhao, S.: Revisiting ENSO/Indian Ocean dipole phase relationships, Geophys. Res. Lett., 44, 2481–2492, 2017.

Toniazzo and Scaife(2006)

Toniazzo, T. and Scaife, A. A.: The influence of ENSO on winter North Atlantic climate, Geophys. Res. Lett., 33, L24704, 10.1029/2006GL027881, 2006.

Trascasa-Castro et al.(2021)Trascasa-Castro, Ruprich-Robert, Castruccio, and Maycock

Trascasa-Castro, P., Ruprich-Robert, Y., Castruccio, F., and Maycock, A. C.: Warm Phase of AMV Damps ENSO Through Weakened Thermocline Feedback, Geophys. Res. Lett., 48, e2021GL096149, 10.1029/2021GL096149, 2021.

van den Hurk et al.(2023)van den Hurk, White, Ramos, Ward, Martius, Olbert, Roscoe, Goulart, and Zscheischler

van den Hurk, B. J., White, C. J., Ramos, A. M., Ward, P. J., Martius, O., Olbert, I., Roscoe, K., Goulart, H. M., and Zscheischler, J.: Consideration of compound drivers and impacts in the disaster risk reduction cycle, iScience, 26, 106030, 10.1016/j.isci.2023.106030, 2023.

van der Wiel et al.(2019)van der Wiel, Bloomfield, Lee, Stoop, Blackport, Screen, and Selten

van der Wiel, K., Bloomfield, H. C., Lee, R. W., Stoop, L. P., Blackport, R., Screen, J. A., and Selten, F. M.: The influence of weather regimes on European renewable energy production and demand, Environ. Res. Lett., 14, 094010, 10.1088/1748-9326/ab38d3, 2019.

van Oldenborgh et al.(2005)van Oldenborgh, Philip, and Collins

van Oldenborgh, G. J., Philip, S. Y., and Collins, M.: El Niño in a changing climate: a multi-model study, Ocean Sci., 1, 81–95, 10.5194/os-1-81-2005, 2005.

Van Stan et al.(2011)Van Stan, Siegert, Levia, and Scheick

Van Stan, J., Siegert, C., Levia, D., and Scheick, C.: Effects of wind-driven rainfall on stemflow generation between codominant tree species with differing crown characteristics, Agr. Forest Meteorol., 151, 1277, 10.1016/j.agrformet.2011.05.008, 2011.

Vignotto et al.(2021)Vignotto, Engelke, and Zscheischler

Vignotto, E., Engelke, S., and Zscheischler, J.: Clustering bivariate dependencies of compound precipitation and wind extremes over Great Britain and Ireland, Weather Clim. Extremes, 32, 100318, 10.1016/j.wace.2021.100318, 2021.

Wahl et al.(2015)Wahl, Jain, Bender, Meyers, and Luther

Wahl, T., Jain, S., Bender, J., Meyers, S., and Luther, M.: Increasing risk of compound flooding from storm surge and rainfall for major US cities, Nat. Clim. Change, 5, 1093–1097, 10.1038/nclimate2736, 2015.

Waliser and Guan(2017)

Waliser, D. and Guan, B.: Extreme winds and precipitation during landfall of atmospheric rivers, Nat. Geosci., 10, 179–183, 10.1038/ngeo2894, 2017.

Wallace and Gutzler(1981)

Wallace, J. M. and Gutzler, D. S.: Teleconnections in the Geopotential Height Field during the Northern Hemisphere Winter, Mon. Weather Rev., 109, 784–812, 10.1175/1520-0493(1981)109<0784:TITGHF>2.0.CO;2, 1981.

Wang(2019)

Wang, C.: Three-ocean interactions and climate variability: a review and perspective, Clim. Dynam., 53, 5119–5136, 10.1007/s00382-019-04930-x, 2019.

Wang et al.(2023)Wang, Chen, Tett, Stone, Nie, Feng, Yan, Zhai, and Ge

Wang, J., Chen, Y., Tett, S. F. B., Stone, D., Nie, J., Feng, J., Yan, Z., Zhai, P., and Ge, Q.: Storyline attribution of human influence on a record-breaking spatially compounding flood-heat event, Science Advances, 9, eadi2714, 10.1126/sciadv.adi2714, 2023.

Whan and Zwiers(2017)

Whan, K. and Zwiers, F.: The impact of ENSO and the NAO on extreme winter precipitation in North America in observations and regional climate models, Clim. Dynam., 48, 1401–1411, 10.1007/s00382-016-3148-x, 2017.

White et al.(2022)White, Domeisen, Acharya, Adefisan, Anderson, Aura, Balogun, Bertram, Bluhm, Brayshaw et al.

White, C. J., Domeisen, D. I., Acharya, N., Adefisan, E. A., Anderson, M. L., Aura, S., Balogun, A. A., Bertram, D., Bluhm, S., Brayshaw, D. J., Browell, J., Büeler, D., Charlton-Perez, A., Chourio, X., Christel, I., Coelho, C. A. S., DeFlorio, M. J., Delle Monache, L., Di Giuseppe, F.,García-Solórzano, A. M., Gibson, P. B., Goddard, L., González Romero, C., Graham, R. J., Graham, R. M., Grams, C. M., Halford, A., Huang, W. T. K., Jensen, K., Kilavi, M., Lawal, K. A., Lee, R. W., MacLeod, D., Manrique-Suñén, A., Martins, E. S. P. R., Maxwell, C. J., Merryfield, W. J., Muñoz, Á. G., Olaniyan, E., Otieno, G., Oyedepo, J. A., Palma, L., Pechlivanidis, I. G., Pons, D., Ralph, F. M., Reis, D. S., Remenyi, T. A., Risbey, J. S., Robertson, D. J. C., Robertson, A. W., Smith, S., Soret, A., Sun, T., Todd, M. C., Tozer, C. R., Vasconcelos, F. C., Vigo, I., Waliser, D. E., Wetterhall, F., and Wilson, R. G.: Advances in the application and utility of subseasonal-to-seasonal predictions, B. Am. Meteorol. Soc., 103, E1448–E1472, 2022.

Williams(1978)

Williams, G. P.: Bank-full discharge of rivers, Water Resour. Res., 14, 1141–1154, 10.1029/WR014i006p01141, 1978.

Xie et al.(2020)Xie, Sun, and Gao

Xie, N., Sun, Y., and Gao, M.: The Influence of Five Teleconnection Patterns on Wintertime Extratropical Cyclones over Northwest Pacific, Atmosphere-Basel, 11, 1248, 10.3390/atmos11111248, 2020.

Xing et al.(2022)Xing, Yu, Wei, Zhang, Ma, Yi, Ju, Wang, Laux, and Kunstmann

Xing, Z., Yu, Z., Wei, J., Zhang, X., Ma, M., Yi, P., Ju, Q., Wang, J., Laux, P., and Kunstmann, H.: Lagged influence of ENSO regimes on droughts over the Poyang Lake basin, China, Atmos. Res., 275, 106218, 10.1016/j.atmosres.2022.106218, 2022.

Xue et al.(2022)Xue, Luo, Zhang, and Yamagata

Xue, J., Luo, J.-J., Zhang, W., and Yamagata, T.: ENSO–IOD Inter-Basin Connection Is Controlled by the Atlantic Multidecadal Oscillation, Geophys. Res. Lett., 49, e2022GL101571, 10.1029/2022GL101571, 2022.

Yeh et al.(2018)Yeh, Cai, Min, McPhaden, Dommenget, Dewitte, Collins, Ashok, An, Yim, and Kug

Yeh, S.-W., Cai, W., Min, S.-K., McPhaden, M. J., Dommenget, D., Dewitte, B., Collins, M., Ashok, K., An, S.-I., Yim, B.-Y., and Kug, J.-S.: ENSO Atmospheric Teleconnections and Their Response to Greenhouse Gas Forcing, Rev. Geophys., 56, 185–206, 10.1002/2017RG000568, 2018.

Yin et al.(2023)Yin, Gao, Chen, Yu, Wilby, Wright, Ge, Bricker, Gong, and Guan

Yin, J., Gao, Y., Chen, R., Yu, D., Wilby, R., Wright, N., Ge, Y., Bricker, J., Gong, H., and Guan, M.: Flash floods: why are more of them devastating the world's driest regions?, Nature, 615, 212–215, 2023.

Zhang et al.(2011)Zhang, Alexander, Hegerl, Jones, Tank, Peterson, Trewin, and Zwiers

Zhang, X., Alexander, L., Hegerl, G. C., Jones, P., Tank, A. K., Peterson, T. C., Trewin, B., and Zwiers, F. W.: Indices for monitoring changes in extremes based on daily temperature and precipitation data, WIREs Climate Change, 2, 851–870, 10.1002/wcc.147, 2011.

Zhang et al.(2021)Zhang, Sun, and Chen

Zhang, Y., Sun, X., and Chen, C.: Characteristics of concurrent precipitation and wind speed extremes in China, Weather Clim. Extremes, 32, 100322, 10.1016/j.wace.2021.100322, 2021.

Zhao(2015)

Zhao, H.: On the distinct interannual variability of tropical cyclone activity over the eastern North Pacific, Atmósfera, 28, 161–178, 10.20937/ATM.2015.28.03.02, 2015.

Zscheischler and Lehner(2021)

Zscheischler, J. and Lehner, F.: Attributing compound events to anthropogenic climate change, B. Am. Meteorol. Soc., 103, E936–E953, 10.1175/BAMS-D-21-0116.1, 2021.

Zscheischler and Seneviratne(2017)

Zscheischler, J. and Seneviratne, S.: Dependence of drivers affects risks associated with compound events, Sci. Adv., 3, e1700263, 10.1126/sciadv.1700263, 2017.

Zscheischler et al.(2018)Zscheischler, Westra, Hurk, Seneviratne, Ward, Pitman, AghaKouchak, Bresch, Leonard, Wahl, and Zhang

Zscheischler, J., Westra, S., Hurk, B., Seneviratne, S., Ward, P., Pitman, A., AghaKouchak, A., Bresch, D., Leonard, M., Wahl, T., and Zhang, X.: Future climate risk from compound events, Nat. Clim. Change, 8, 469–477, 10.1038/s41558-018-0156-3, 2018.

Zscheischler et al.(2020)Zscheischler, Martius, Westra, Bevacqua, Raymond, Horton, Hurk, AghaKouchak, Jézéquel, Mahecha, Maraun, Ramos, Ridder, Thiery, and Vignotto

Zscheischler, J., Martius, O., Westra, S., Bevacqua, E., Raymond, C., Horton, R., Hurk, B., AghaKouchak, A., Jézéquel, A., Mahecha, M., Maraun, D., Ramos, A., Ridder, N., Thiery, W., and Vignotto, E.: A typology of compound weather and climate events, Nat. Rev. Earth Environ., 1, 1–5, 2020.