# Evaluating the dependence structure of compound precipitation and wind speed extremes

Jakob Zscheischler, Philippe Naveau, Olivia Martius, Sebastian Engelke,

Christoph C. Raible

► To cite this version:

Jakob Zscheischler, Philippe Naveau, Olivia Martius, Sebastian Engelke, Christoph C. Raible. Evaluating the dependence structure of compound precipitation and wind speed extremes. Earth System Dynamics, 2021, 12, pp.1-16. 10.5194/esd-12-1-2021. insu-03721921

HAL Id: insu-03721921

https://insu.hal.science/insu-03721921

Submitted on 13 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Evaluating the dependence structure of compound precipitation and wind speed extremes

Jakob Zscheischler<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Philippe Naveau<fnref n="4" />, Olivia Martius<fnref n="1" />,<fnref n="5" />,<fnref n="6" />, Sebastian Engelke<fnref n="7" />, and Christoph C. Raible<fnref n="1" />,<fnref n="2" />

**Correspondence:** Jakob Zscheischler (jakob.zscheischler@ufz.de)

Received: 28 May 2020 – Discussion started: 4 June 2020

Revised: 27 October 2020 – Accepted: 9 November 2020 – Published: 6 January 2021

**Abstract.** Estimating the likelihood of compound climate extremes such as concurrent drought and heatwaves or compound precipitation and wind speed extremes is important for assessing climate risks. Typically, simulations from climate models are used to assess future risks, but it is largely unknown how well the current generation of models represents compound extremes. Here, we introduce a new metric that measures whether the tails of bivariate distributions show a similar dependence structure across different datasets. We analyse compound precipitation and wind extremes in reanalysis data and different high-resolution simulations for central Europe. A state-of-the-art reanalysis dataset (ERA5) is compared to simulations with a weather model (Weather Research and Forecasting – WRF) either driven by observation-based boundary conditions or a global circulation model (Community Earth System Model – CESM) under present-day and future conditions with strong greenhouse gas forcing (Representative Concentration Pathway 8.5 – RCP8.5). Over the historical period, the high-resolution WRF simulations capture precipitation and wind extremes as well as their response to orographic effects more realistically than ERA5. Thus, WRF simulations driven by observation-based boundary conditions are used as a benchmark for evaluating the dependence structure of wind and precipitation extremes. Overall, boundary conditions in WRF appear to be the key factor in explaining differences in the dependence behaviour between strong wind and heavy precipitation between simulations. In comparison, external forcings (RCP8.5) are of second order. Our approach offers new methodological tools to evaluate climate model simulations with respect to compound extremes.

# 1 Introduction

Compound climate extremes such as co-occurring drought and heat or compound precipitation and wind extremes can have a substantial impact on the natural environment and human systems that often exceeds the impact caused by a single extreme (Zscheischler et al., 2014; Raveh-Rubin and Wernli, 2015; Martius et al., 2016; Sippel et al., 2018). Over

recent years a number of compound extremes have been investigated. For instance, several studies have analysed the dependence between storm surge and heavy precipitation (Wahl et al., 2015; Zheng et al., 2013; Bevacqua et al., 2019) or extreme runoff (Ward et al., 2018; Hendry et al., 2019) to estimate the risk of compound flooding in coastal areas. Compound droughts and heatwaves have been studied for dif-

<footnote n="1">¹Oeschger Centre for Climate Change Research, University of Bern, Bern, Switzerland</footnote>

<footnote n="2">²Climate and Environmental Physics, University of Bern, Sidlerstrasse 5, 3012 Bern, Switzerland</footnote>

<footnote n="3">³Department of Computational Hydrosystems,
Helmholtz Centre for Environmental Research – UFZ, Leipzig, Germany</footnote>

<footnote n="4">⁴Laboratoire des Sciences du Climat et de l'Environnement, Gif-sur-Yvette, France</footnote>

<footnote n="5">⁵Institute of Geography, University of Bern, Bern, Switzerland</footnote>

<footnote n="6">⁶Mobiliar Lab for Natural Risks, University of Bern, Bern, Switzerland</footnote>

<footnote n="7">⁷Research Center for Statistics, University of Geneva, Geneva, Switzerland</footnote>
ferent regions and varying temporal scales (Mazdiyasni and AghaKouchak, 2015; Zscheischler and Seneviratne, 2017; Manning et al., 2019; Sutanto et al., 2020; Zscheischler and Fischer, 2020). The occurrence rate of compound precipitation and wind extremes has been estimated for the Mediterranean region (Raveh-Rubin and Wernli, 2015), Europe (De Luca et al., 2020), and at the global scale (Martius et al., 2016). Other studies have investigated the co-occurrence of hot days and hot nights (Wang et al., 2020) or the co-occurrence rate of heavy precipitation and snowmelt to estimate the risk of rain-on-snow events (Musselman et al., 2018; Poschlod et al., 2020). Such a quantification of the occurrence rate of compound extremes is important for assessing the risk of associated impacts today and in the future. Most of the above studies identify compound extremes by thresholding the contributing variables to quantify the occurrence of compound extremes and changes associated with climate change. However, the dependence structure in the tails is only rarely investigated. Due to the rarity of compound extremes, a large number of samples is required to obtain robust estimates, making it difficult to rely solely on observational data (Ridder et al., 2020).

Large ensemble simulations (Deser et al., 2020) offer an opportunity to estimate future changes in the occurrence of compound events without running into data limitations (Poschlod et al., 2020; Champagne et al., 2020). However, such simulations need to be interpreted with care as it is often largely unknown how well the employed models represent observed compound events (Musselman et al., 2018), and differences might be large between models. Climate models are regularly evaluated based on their ability to represent well-known processes in the climate system and predominantly univariate comparisons with key climate variables (Flato et al., 2013), though some multivariate metrics have been explored (Sippel et al., 2017). Yet little is known about the ability of climate models to capture observed occurrence rates of compound extremes (Zscheischler et al., 2018), a challenging task for two primary reasons. First, due to their rarity, a robust quantification of the likelihood of compound extremes requires large amounts of data, thus making it difficult to establish a “ground truth” for many applications. Second, suitable metrics for evaluating multivariate extremes have not been widely tested and applied in a climate context. Such metrics, however, are essential to assess how well models represent compound events, in particular to assess future risks (Zscheischler et al., 2020). When observational data are scarce, process-based model simulations (Couasnon et al., 2020) and reanalysis data (Martius et al., 2016) can be employed to extend or replace purely observational datasets.

To date, model-data comparisons related to compound extremes have been conducted to a very limited extent, often relying on simplifying assumptions and typically confined to precipitation and temperature. For instance, a high likelihood of compound hot and dry summers has been linked to a strongly negative correlation between summer tempera-

ture and precipitation (Zscheischler and Seneviratne, 2017). While there is generally good agreement with respect to this correlation between climate models and observation-based datasets in the Northern Hemisphere, there is strong disagreement in the Southern Hemisphere, for which the models show a much stronger dependence. This finding may suggest that climate models overestimate dependence between summer temperature and precipitation. However, this discrepancy may also be related to the way gridded observation-based datasets are assembled. In particular, for locations without an active measurement station nearby, the mean seasonal cycle is often used to fill gaps in the observational networks (e.g. Mitchell and Jones, 2005). This approach reduces the strength of co-variability between temperature and precipitation in poorly sampled regions, which are mostly in the Southern Hemisphere. Thus, assessing the ability of climate models to represent compound events may reveal underappreciated limitations in gridded observation-based datasets. We are not aware of studies so far that have evaluated the dependence between precipitation and wind speed.

In this study we focus on compound precipitation and wind extremes, which can have severe socio-economic impacts including human fatalities, impaired critical infrastructure, and economical damage (Fink et al., 2009; Lin et al., 2010; Liberato, 2014; Raveh-Rubin and Wernli, 2015; Martius et al., 2016). We investigate differences in the occurrence of compound precipitation and wind extremes for different datasets over a region in central Europe including the Alps. To this end, we introduce a new measure that assesses dissimilarity between the tails of bivariate distributions. We study an experimental design with two factors. The first factor is the type of boundary conditions in a high-resolution regional weather model, either from reanalysis or a global circulation model. The second factor corresponds to the effect of different climate forcing between today and the future under a high-emission scenario. Our object of study under this design is the dependence between heavy precipitation and strong wind in winter over central Europe. In addition, comparisons with a state-of-the-art reanalysis product are implemented.

## 2 Data

We use daily precipitation sums and daily maximum wind speed in the extended winter (November–March) from one reanalysis product and three model simulations over a period of 20 years. The employed reanalysis product is ERA5 (Copernicus Climate Change Service, 2017) from which we use data for the period 1980 to 1999 CE. This reanalysis is generated with an updated numerical weather prediction model and data assimilation system compared to the prior product ERA-Interim (Dee et al., 2011) and integrates additional data sources. The data are available at a resolution of roughly 30 km (spectral resolution of T639), with 137 vertical levels and hourly output.
The three simulations are performed with the Weather Research and Forecasting (WRF) model (Skamarock et al., 2008) which is forced with boundary conditions from (i) ERA-Interim (Dee et al., 2011) (ERAI-WRF), (ii) a period of free-running global climate simulations for the present day (CESM-WRF), and (iii) a period covering the end of the 21st century under Representative Concentration Pathway 8.5 (CESM-WRF-fut, a high-emission scenario). The global climate simulation is performed with the Community Earth System Model (CESM) (Hurrell et al., 2013) for the period 850 to 2100 CE. Details on the setting are described in Lehner et al. (2015) and Raible et al. (2018). In this study we use the period 1980 to 1999 CE as the present day and 2080 to 2099 CE as the future.

The periods of the global simulations and the ERA-Interim period (1980 to 1999 CE) are dynamically downscaled with WRF in version 3.5. WRF is vertically discretized in 40 terrain-following η-coordinate levels. The horizontal resolutions of the four two-way nested domains (Fig. 1) are 54, 18, 6, and 2 km. The innermost domain covers the box [4.75°E, 15.25°E] × [43.25°N, 48.75°N] and is exclusively used in this study. The setup is described in more detail in Gómez-Navarro et al. (2015, 2018) and Messmer et al. (2017, 2020). Important for this a study is that the convection parameterization is disabled for the simulations at 6 and 2 km resolution; at these scales the model is convection-permitting. This is an important step in improving the simulation of precipitation, though some problems still remain (Ban et al., 2014). For adequately simulating wind, the setting of the planetary boundary layer parameterization is key. We use a modified version of the fully non-local scheme developed by Hong and Lim (2020), which specifically treats effects of the unresolved orography (Jimenez and Dudhia, 2012). For the ERAI-WRF simulation we allow analysis nudging of wind, temperature, and humidity above the planetary boundary layer in order to stay close to the large-scale behaviour of the reanalysis data (Gómez-Navarro et al., 2015). For the two simulations driven by CESM, nudging is omitted to allow the regional model to correct potential systematic biases of the CESM (e.g. a zonal atmospheric circulation that is too strong in the mid-latitudes; Bracegirdle et al., 2013). The WRF output is provided in hourly resolution.

We remap the original hourly data to a common regularly spaced grid with 0.25° spatial resolution using conservative remapping and subsequently compute daily precipitation sums and daily wind speed maxima. The 0.25° spatial resolution was chosen as it is closest to the original resolution of the ERA5 reanalysis data. Note, however, that all WRF simulations are run at a much higher convection-resolving resolution. The explicit resolution of convection and a much higher resolution of the topography may result in a more accurate representation of the dependence between precipitation and wind extremes in the simulations than in ERA5. We further note that mean wind speed in ERA5 generally decreases with elevation (Fig. 2a), which is the opposite of the

<figure><img src="image_4.png" /><figcaption>**Figure 1.** The four nested domains used in the dynamical downscaling.</figcaption></figure>

expected behaviour for the response of wind speed to elevation from observations (Graf et al., 2019; Telesca et al., 2020) and what is modelled by WRF (Fig. 2b). The discrepancy in mountainous regions between reanalysis data and observations with respect to wind speed is also evident in other reanalysis datasets such as ERA-Interim (Jones et al., 2017), which is the predecessor of ERA5. In contrast, WRF has been shown to also simulate wind speed reasonably well in mountainous terrain (Stucki et al., 2016). For these reasons – WRF better resolves cloud processes, topography, and wind speed, while ERA5 misrepresents the wind speed gradient with elevation – we use ERAI-WRF as the reference for all analyses.

# 3 A measure for evaluating compound extremes

## 3.1 Measuring tail dependence

The extreme values of a univariate random variable can be analysed with tools from extreme value theory (Embrechts et al., 1997; Coles, 2001; Katz et al., 2002; Naveau et al., 2020). For multivariate random vectors, the dependence between the largest values in the components becomes important (Davison and Huser, 2015; Huser and Wadsworth, 2020; Engelke and Ivanovs, 2021).

We quickly review the concept of bivariate asymptotic tail dependence and independence (Ledford and Tawn, 1997; Poon et al., 2003). Two variables $X_1$ and $X_2$ with cumulative distribution functions $F_1$ and $F_2$, respectively, are asymptotically dependent if

$$ \chi = \lim_{q \to 1} \mathbb{P}(F_1(X_1) > q | F_2(X_2) > q) \in (0, 1] $$
<figure><img src="image_5.png" /><figcaption>Figure 2. Relationship between mean winter wind speed and altitude for ERA5 (a) and the WRF model (ERAI-WRF simulation) (b).</figcaption></figure>

and asymptotically independent otherwise (i.e. if $\chi = 0$). The coefficient $\chi$ is called extremal correlation and represents, after transforming $X_1$ and $X_2$ to the uniform scale, the probability of one variable being extreme given that the other one is extreme. Note that two variables can be dependent at normal levels but asymptotically independent in the extremes, as in the case for a bivariate Gaussian distribution (Sibuya, 1960). To fine tune the rate of decay towards the asymptotically independent case ($\chi = 0$), the residual tail dependence coefficient $\bar{\chi}$ contains additional information (Coles et al., 1999):

$$ \bar{\chi} = \lim_{q \to 1} \frac{\log(\mathbb{P}(F_1(X_1) > q)\mathbb{P}(F_2(X_2) > q))}{\log(\mathbb{P}(F_1(X_1) > q, F_2(X_2) > q))} - 1 \in [-1, 1]. $$

$\bar{\chi}$ is equal to 1 for asymptotically dependent variables, while for asymptotically independent variables $\bar{\chi}$ indicates if $X_1$ and $X_2$ are positively ($\bar{\chi} > 0$) or negatively ($\bar{\chi} < 0$) associated in their extremes. Thus, the pair of coefficients ($\chi, \bar{\chi}$) summarizes the tail dependence structure of $X_1$ and $X_2$.

Because both coefficients $\chi$ and $\bar{\chi}$ are defined as a limit value, a usual way to analyse the behaviour of a bivariate tail dependence structure between two variables is to compute empirical estimates for varying threshold levels $q$ and then visually inspect their behaviour as $q \to 1$. We estimate $\chi$ and

$\bar{\chi}$ with the function taildep from the R package extRemes (Gilleland and Katz, 2016).

To estimate $\chi$ empirically we use a high quantile for which a reasonably large number of data are available. For these reasons we generally estimate $\chi$ at $q = 0.95$. Heavy precipitation events and extreme winds that lead to large damage can be linked through storms or föhn events across neighbouring locations with a lag of several days. To take this aspect into account, we estimate $\chi$ using a local block maxima approach, which is motivated by Ferreira and de Haan (2015). We thus first compute the daily precipitation and wind speed maxima for varying block sizes ranging from 0.25° (approximately 20–30 km) to 1.75° (that is, a maximum of three grid points in any direction, or 100–200 km) and up to 5 d (i.e. maximum of 2 d before and after the day of interest).

We further assess whether estimates of $\chi$ are significantly different from 0. To this end, we bootstrap the data by randomly shuffling the temporal order of one variable to break the dependence structure. The coefficient $\chi$ is then estimated as above. Estimates of $\chi$ are considered significantly different from 0 if they are larger than 95% of the bootstrapped estimates.

## 3.2 Measuring differences in bivariate extremal dependence structures

Classical tail coefficients like $\chi$ are informative summaries to assess the extremal dependence between two univariate random variables, say $X_1$ and $X_2$, but they cannot quantify the difference in extremal dependence between two bivariate random vectors, say $X^{(1)} = (X_1^{(1)}, X_2^{(1)})$ and $X^{(2)} = (X_1^{(2)}, X_2^{(2)})$. For example, $\chi^{(1)}$ can be computed between heavy precipitation and strong winds computed from one dataset, e.g. ERA5, and compared to a $\chi^{(2)}$ for a second dataset, e.g. ERAI-WRF. But it would also be very convenient to have a single number to tell us if the extremal dependences between these two bivariate random vectors are different and, if so, by how much. Recent work by Naveau et al. (2014) showed that the well-known Kullback-Leibler (KL) divergence used in signal processing can be tailored to the framework of extreme value theory. The approach has been applied to cluster climate data according to their bivariate extremal behaviour (Vignotto et al., 2020). However, to our knowledge, multivariate extremal divergence measures have never been applied to the analysis of compound weather and climate events. By complementing tail coefficients, this new tool could shed new light on the joint behaviour of heavy precipitation and strong winds across our different datasets.

The KL divergence is defined on marginals which are normalized to standard Pareto distributions. A risk function ($r: \mathbb{R}^2 \to \mathbb{R}$) is used to describe the extreme region in each of the bivariate distributions. There are different choices for the risk function. Taking the sum or the maximum gives similar results for asymptotically dependent data. In addi-
tion, the minimum also covers asymptotic independence. The sum is defined as $r(\mathbf{x}) = x_1 + x_2$ and the minimum as $r(\mathbf{x}) = \min(x_1, x_2)$, with $\mathbf{x} = (x_1, x_2)$. Hence, we consider those for which the sum (or minimum) of the components exceeds a given high quantile $q_u^{(j)}$ of $r(\mathbf{X}^{(j)})$ corresponding to an exceedance probability $u \in (0, 1)$, $j = 1, 2$ to be extreme points. Varying the threshold $u$ alters the extremal region of interest. For each of the two bivariate distributions, the set $A^{(j)} = \{r(\mathbf{x}) > q_u^{(j)}\}$ and $j = 1, 2$ is partitioned into a fixed number $W$ of disjoint sets $A_1^{(j)}, \dots, A_W^{(j)}$.

For two random samples $(\mathbf{X}_1^{(1)}, \dots, \mathbf{X}_n^{(1)})$ and $(\mathbf{X}_1^{(2)}, \dots, \mathbf{X}_n^{(2)})$ from the distributions $\mathbf{X}^{(1)}$ and $\mathbf{X}^{(2)}$, the empirical proportions of data points belonging to set $A_W^{(j)}$ are computed as

$$ \hat{p}_w^{(j)} = \frac{\#\{i : \mathbf{X}_i^{(j)} \in A_w^{(j)}\}}{\#\{i : r(\mathbf{X}_i^{(j)}) > q_u^{(j)}\}}, \quad w = 1, \dots, W. $$

The difference between the extremal behaviours of the two distributions can then be measured as the KL divergence between the two multinomial distributions defined through these proportions, i.e.

$$ d_{12} = D(\mathbf{X}^{(1)}, \mathbf{X}^{(2)}) = \frac{1}{2} \sum_{w=1}^{W} \left( \left( \hat{p}_{w}^{(1)} - \hat{p}_{w}^{(2)} \right) \log \left( \hat{p}_{w}^{(1)} / \hat{p}_{w}^{(2)} \right) \right). \quad (1) $$

Note that this divergence is symmetric and since it is a non-parametric statistic it does not require additional model assumptions. Equation (1) contrasts differences among extremal dependence structures for both asymptotically dependent and asymptotically independent data. The number of partitioning sets $W$ is a free parameter. If it is chosen too high, many sets will be empty, resulting in an undefined KL divergence. If it is too small, only a rough summary is computed but not really an estimate of tail dependence. We chose $W = 5$ in this study based on the simulation study shown in Appendix A. Under suitable assumptions the statistic $d_{12}$ follows a $\chi^2(W-1)$ distribution in the limit as the sample size goes to $\infty$, which allows us to estimate whether distances are significantly different from 0.

The approach is illustrated in Figs. 3 and 4 with $W = 3$. Figure 3 shows daily precipitation sums and maximum wind speed at grid point 9° E, 46.75° N on the original scale (a, d) with margins normalized to the exponential scale (b, e) and to standard Pareto distributions (c, f) for ERAI-WRF (a–c) and CESM-WRF (d–f). The shown grid point reaches the highest tail dependence $\chi$ at $q = 0.95$ in the ERAI-WRF simulation. The colours in all subpanels and the dashed lines in Fig. 3c and f highlight the three disjoint sets $A_1^{(j)}, A_2^{(j)}$, and $A_3^{(j)}$ (see above). At the exponential scale moderate and large extremes can be seen well, whereas at the Pareto scale only very extreme values can be easily visually identified. Figure 4

illustrates $\chi$ (a) and $\bar{\chi}$ (c) for the distributions of the two simulations and the divergence based on Eq. (1) with “sum” (b) and “min” (d) as the risk function, including 95% confidence intervals of the empirical estimates. The estimates of $\chi$ and $\bar{\chi}$ start to diverge somewhat for $q > 0.8$, suggesting different tail behaviour (uncertainty ranges are estimated based on the R function chiplot from the package evd; Stephenson, 2002). This impression is confirmed by the estimates of the KL divergence: for most thresholds $u > 0.8$ and for both choices of the risk function the KL divergence is outside the 95% quantile of the limiting $\chi^2(W-1)$ distribution of the statistic $d_{12}$ under the null hypothesis of equal tail dependence structures. This means that we can conclude that the two distributions have significantly different tail behaviour.

Note that in the bivariate case, a simple approach to quantify the difference in tail dependence would be the difference between $\chi^{(1)}$ and $\chi^{(2)}$. However, for two distributions with the same $\chi$ coefficient but a different dependence structure, it is impossible to distinguish the two cases. In a way, $\chi$ only focuses on the “diagonal”. Furthermore, while in this work we focus on the bivariate case, the KL divergence defined by Eq. (1) could easily be implemented with higher dimensions $d = 3, (4, \dots)$ because it is only based on counting points in different subsets. In contrast, using $\chi$, the number of pairs will increase rapidly with the dimension $d$. In addition, $\chi$ coefficients will only capture pairwise dependencies.

We investigate how well different simulations represent the bivariate tail behaviour of daily precipitation sums and wind speed maxima in winter by comparing ERA5, CESM-WRF, and CESM-WRF-fut against ERAI-WRF with the divergence as defined in Eq. (1) based on the maxima over the spatio-temporal blocks that maximize tail dependence $\chi$ at $q = 0.95$. Using local block maxima ensures that $\chi$, $\bar{\chi}$, and the KL divergence measure very extremal upper tail behaviour. Note, however, that this approach leads to different block sizes depending on the location, which makes a direct comparison in space difficult. For the computation of the KL divergence (Eq. 1) we use $u = 0.9$ and “sum” as the risk function. We further perform a sensitivity test using $u \in \{0.8, 0.85, 0.9, 0.95\}$. Furthermore, the marginals have been transformed into a Pareto scale through ranking. The choice of marginal transformation only has a minor influence on the KL divergence (see Appendix A).

## 4 Results

We first present a simple correlation analysis based on the Spearman's rank correlation coefficient. Daily precipitation sums and maximum wind speed are generally strongly correlated in winter in most areas of the study domain except in the northwest of Italy (Fig. 5). All model simulations show a relatively consistent pattern, whereas ERA5 shows negative correlations at the southern slopes of the Alps along the
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Scatterplots of daily precipitation and wind speed maxima in November–March (1980–1999) for the location with the highest tail dependence $\chi$ ($q = 0.95$) in the ERAI-WRF simulations (**a–c**). CESM-WRF simulations for the same location are shown in (**d–f**). Shown are the original values (**a, d**), values after transformation into exponential marginals (**b, e**), and values after transformation into Pareto marginals (**c, f**). The colours highlight the three separating sets W to compute the KL divergence (see Eq. 1) for a high threshold (see the main text). In (**c** and **f**), the three sets are separated by dashed lines. Note that for the main analyses in the paper we use $W = 5$.</figcaption></figure>

northwestern Italian border (Fig. 5b). Most correlations are significant ($\alpha = 0.05$).

When considering only the dependence in the tails based on $\chi$ and including a spatial and temporal neighbourhood, the spatial patterns look quite different (Fig. 6). The WRF simulations show a highly heterogeneous picture with strong local variations, generally strong dependence over most parts of the Alps and close to the Adriatic coast, and weak dependence otherwise (Fig. 6a, c, d). Overall, ERAI-WRF shows slightly higher tail dependence compared to the WRF simulations driven by CESM. In contrast to the WRF simulations, in ERA5 tail dependence varies rather smoothly in space, with higher values in northeastern Italy and along the eastern border of France (Fig. 6b).

The block sizes that attain the maximum tail dependence $\chi$ for precipitation and wind extremes for each pixel are shown in Fig. 7. On average for 75% of the pixels, the maximum is attained with no temporal lag (blue in Fig. 7). In contrast, there seems to be a shift in space, as maxima tend to co-occur in neighbouring grid points: block sizes with larger than minimal (0.25°) spatial extent occur on average in 60% of all locations (lighter colours in Fig. 7). This means that extremes in daily precipitation sums and wind extremes tend to occur on the same day but potentially not at exactly the same location; they occur with some distance apart. In particular in

the south of the Alps but also in some regions north of the Alps, this distance is 1.75°, or about 100–200 km (very light colours in Fig. 7). The strongest agreement of the dependence patterns is between CESM-WRF and CESM-WRF-fut, which agree for half of the locations in the maximizing block size. In contrast, the agreement is 29% between ERAI-WRF and ERA5 and 39% between ERAI-WRF and CESM-WRF. Note that grid points at the boundaries cannot attain maxima with block sizes larger than one grid point as no data values are available outside the study domain.

The tails between winter daily precipitation sums and wind speed maxima show a significantly different dependence structure between ERAI-WRF and CESM-WRF in 46% of all grid points, mostly in Switzerland and in the north of the study domain but also in many regions in northern Italy (Fig. 8a). The percentage of grid points with significantly different tail behaviour is slightly higher for the comparison of ERAI-WRF and ERA5 (49%), though in this case most of the differences occur in grid points located along a wide diagonal band from the southwest to the northeast through the entire study domain (Fig. 8b). Interestingly, the comparison of ERAI-WRF with CESM-WRF-fut results in only 36% of pixels with significantly different tail behaviour (Fig. 8c). Thus, CESM-WRF-fut agrees better with ERAI-WRF with respect to the tail behaviour than CESM-WRF
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Illustration of the distance metrics between bivariate tails for the location with the highest estimated tail dependence $\chi$ at $q = 0.95$ in ERAI-WRF. (a, c) Tail dependence parameters $\chi$ (a) and $\bar{\chi}$ (c) for daily precipitation sums and daily maximum wind speed for different quantile-based thresholds $q$. Shading highlights the 95% confidence intervals. Grey: ERAI-WRF. Red: CESM-WRF. (b, d) Two different Kullback-Leibler (KL) divergences (Eq. 1 with $W = 5$) for the tails of the bivariate precipitation-wind speed distribution between ERAI-WRF and CESM-WRF (solid lines). Dashed lines highlight the 95% confidence interval of the null hypothesis assuming an equal dependence structure. (b) KL divergence based on the minimum (i.e. $\min(X_1, X_2) > u$). (d) KL divergence based on the sum (i.e. $X_1 + X_2 > u$).</figcaption></figure>

and ERAI-WRF. Finally, only 18% of pixels show significantly different tail behaviour when comparing CESM-WRF and CESM-WRF-fut (Fig. 8d), indicating the pair with the largest number of grid points at which no significant difference in the tail behaviour could be found. The numbers of grid points with significantly different tail behaviour depend somewhat on the threshold $u$ and generally decrease with increasing extremeness (that is, increasing $u$), but the differences between the different pairwise comparisons remain similar (Table 1). In particular, the differences between ERAI-WRF and CESM-WRF and between ERAI-WRF and CESM-WRF-fut are generally larger than the differences between CESM-WRF and CESM-WRF-fut, indicating that the main finding, namely that boundary conditions in WRF appear to be the key factor in explaining differences in the dependence behaviour between wind and precipitation extremes, is robust for different parameter values of the difference measure.

## 5 Discussion

We have introduced a new metric for comparing tail dependence structures between wind and precipitation extremes in reanalysis data and weather model simulations. In our WRF simulations, the type of boundary conditions, either ERAI or CESM, appears to have a stronger effect on the coupling between high wind and heavy precipitation than the change in external forcing (present-day and future) in CESM (Fig. 8). This suggests that the studied dependence structures between the tails of precipitation sums and wind speed maxima in winter are a rather robust feature of the combination of models (boundary conditions plus a high-resolution weather model) and thus also somewhat determined by the boundary conditions. In consequence this also means that here we are probably detecting rather stable dynamical features that are largely independent of strong external forcing such as (much) higher mean temperatures. Because the model setting determines the dependence structure, sampling uncertainties in this dependence, for instance to robustly assess risks under future climate conditions, would require a range of different climate and weather model combinations.
<figure><img src="image_8.png" /></figure>

**Figure 5.** Spearman's rank correlation between daily precipitation sums and maximum wind speed in the extended winter (November–March). (a) ERAI-WRF, (b) ERA5, (c) CESM-WRF, (d) CESM-WRF-fut. Non-significant correlations ($\alpha$ = 0.05) are marked with a cross.

<figure><img src="image_9.png" /></figure>

**Figure 6.** Tail dependence ($\chi$ with $q = 0.95$) between daily precipitation sums and maximum wind speed in the extended winter (November–March). Tail dependence was computed considering block maxima over a maximum range of 5 d temporally and 1.75° spatially. (a) ERAI-WRF, (b) ERA5, (c) CESM-WRF, (d) CESM-WRF-fut. Non-significant values based on bootstraps with the same maximum block size ($\alpha$ = 0.05) are marked with a cross.
<table><caption>**Table 1.** Sensitivity analysis of KL divergence (Eq. 1). Reported is the fraction of grid points with a significantly different ($\alpha$ = 0.05) precipitation-wind speed dependence structure between two datasets for different thresholds $u$ (with $W$ = 5). The case $u$ = 0.90 is shown in Fig. 8.</caption><thead><tr><th></th><th>$u$ = 0.80</th><th>$u$ = 0.85</th><th>$u$ = 0.90</th><th>$u$ = 0.95</th></tr></thead><tbody><tr><td>ERAI-WRF vs. CESM-WRF</td><td>0.61</td><td>0.54</td><td>0.46</td><td>0.32</td></tr><tr><td>ERAI-WRF vs. ERA5</td><td>0.59</td><td>0.53</td><td>0.49</td><td>0.40</td></tr><tr><td>ERAI-WRF vs. CESM-WRF-fut</td><td>0.53</td><td>0.45</td><td>0.36</td><td>0.27</td></tr><tr><td>CESM-WRF vs. CESM-WRF-fut</td><td>0.30</td><td>0.23</td><td>0.18</td><td>0.19</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Blocks for which the maximum tail dependence ($\chi$ with $q$ = 0.95) between daily precipitation sums and maximum wind speed in the extended winter (November–March) is attained (Fig. 6). Block sizes range from 0.25° and 1 d to 1.75° and 5 d. Blue, green, and orange refer to time lags of 1, 3, and 5 d, respectively. Darker shading illustrates higher spatial proximity. The colour bars next to the maps show the number of grid points of that colour in the corresponding map. (a) ERAI-WRF, (b) ERA5, (c) CESM-WRF, (d) CESM-WRF-fut. Grid points with non-significant tail dependence are marked with a cross (see Fig. 6).</figcaption></figure>

The employed block maxima approach (Fig. 6) has the effect that precipitation and wind extremes are considered together even if they might occur some distance apart in either time or space. This is to ensure that extremes in wind and precipitation are considered together if they emerge from the same atmospheric processes (e.g. föhn). At the same time, the block maxima approach can help diagnose why datasets differ in their tail dependence structure of precipitation and wind extremes, for instance if the spatio-temporal blocks for which extremes are attained differ strongly.

Regarding the optimal spatial and temporal lags between wind and precipitation extremes there is generally good agreement that along the southern slopes of the Alps the dependence is maximized for precipitation and wind extremes occurring on the same day and up to 1.75° apart (lightest blue in Fig. 7), which could be related to föhn events that lead to

heavy precipitation north of the mountain range and extreme winds on the southern slopes or vice versa. Indeed, heavy precipitation events on the Alpine southern side are often related to high moisture transport ahead of cold fronts that is associated with moderate winds that are not as strong as potential föhn gusts on the Alpine northern side (Panziera and Germann, 2010).

Most heavy precipitation events in the investigation domain in winter are associated with extratropical cyclones. Within extratropical cyclones, wind speed maxima and precipitation maxima are often linked to fronts and conveyor belts (Parton et al., 2010; Catto and Pfahl, 2013; Pfahl et al., 2014; Pantillon et al., 2020), and this may result in co-located extremes. However, important modulations of both extreme wind and precipitation patterns by the local complex orography are to be expected (Whiteman, 2000; Barry, 2008),
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Locations for which the dependence between the tails of daily precipitation sums and wind speed maxima is significantly different based on the KL divergence according to Eq. (1) with $u = 0.9$ and $W = 5$ (dark grey, with $\alpha = 0.05$). Dependence is assessed for the blocks that attain maximum tail dependence $\chi$ (at $q = 0.95$) (see Fig. 6). Shown are comparisons between (a) ERAI-WRF and CESM-WRF, (b) ERAI-WRF and ERA5, (c) CESM-WRF and CESM-WRF-fut, and (d) ERAI-WRF and CESM-WRF-fut.</figcaption></figure>

and such local föhn effects, channelling effects, flow block-
ing, and many more might be captured by the high-resolution
WRF simulations but not in ERA5.

Overall, ERA5 shows quite a different behaviour for the
Spearman’s rank correlation (Fig. 5) and simple tail depen-
dence χ (Fig. 6) compared to the high-resolution weather
model simulations. Spatial patterns are much smoother, prob-
ably related to the much coarser spatial resolution (30 km
compared to the original 2 km in the WRF simulations).
Furthermore, wind speeds over high mountains are unre-
alistic, as they decrease with height rather than increase
(Fig. 2). These limitations render ERA5 unsuitable as a
benchmark for the tail dependence between precipitation and
wind extremes in the Alpine area with its complex orogra-
phy. Presently, homogenized gridded wind observations of
good quality are not available for this region. Therefore,
driving a well-calibrated high-resolution weather model with
observation-based boundary conditions is currently the best
benchmark to study compound wind and precipitation ex-
tremes.

We would like to note that in our setup ERAI-WRF is
nudged to the driving reanalysis ERA-Interim. The reason for
this is that the simulation should stay close to the large-scale
behaviour of the reanalysis data. As mentioned in the Data
section, we only use wind, temperature, and humidity above
the planetary boundary layer and the nudging is not strong.
Nevertheless, the behaviour of extremes might be changed

due to the modification of the dynamical equations to some
extent, but we think that this effect is minor. Furthermore,
precipitation is not nudged.

Evaluating how well models represent tail dependencies
may help select those models that are fit for purpose (Ma-
raun et al., 2017) regarding the analysis of compound events
(Zscheischler et al., 2020) for a range of different event types
(Ridder et al., 2020). In particular, when the interest lies
in the simulation of impacts, the approach may help decide
when multivariate bias adjustment approaches would need to
be employed (François et al., 2020), as univariate bias adjust-
ment might increase biases in impacts that depend on multi-
ple correlated drivers (Zscheischler et al., 2019).

6 Conclusions

Evaluating the ability of climate models to represent the like-
lihood of compound climate extremes is important for well-
informed climate risk assessments. In this study we investi-
gated differences in the tail behaviour of precipitation and
wind extremes in winter between different weather model
simulations and a reanalysis dataset for a region in central
Europe. Employing a new metric to measure differences in
the tail behaviour of bivariate distributions, we found that
simulations of the same model pair with different external
forcing conditions (climate change conditions) differ less
than simulations for present-day conditions with different
boundary data. Our results further suggest that reanalysis
data are not suitable as a benchmark for the analysis of com-
pound precipitation and wind extremes over complex terrain
such as the Alps. Overall, differences between model sim-
ulations (different boundary conditions and weather and/or
climate models) can be substantial. Our results suggest that
climate impact modelling needs to take uncertainties related
to the simulation of compound extremes into account to pro-
vide robust risk assessments for today and the future.
## Appendix A: Determining W

We simulated $n = 2000$ samples of $X^{(1)}$ and $X^{(2)}$ in the outer power Clayton copula, which is in the domain of attraction of the logistic extreme value distribution. We chose the parameters such that the limiting $\chi$ coefficients are 0.4 and 0.55, which means one model with weaker and one with stronger dependence, respectively. Using the KL divergence for a probability threshold of $u = 0.9$, we compare the samples of $X^{(1)}$ and $X^{(2)}$ for the dependence settings weak-weak, strong-strong, and weak-strong and plot in each case the probability of rejecting the null hypothesis of equal tail dependence structures. Note that the former two cases are in line with the null hypothesis, whereas the latter case does not satisfy the null hypothesis. We conduct the experiment for both known margins and empirically normalized margins, as well as for different numbers of sets $W$ in the KL divergence statistic.

Figures A1 and A2 show the Type I error of rejecting the null hypothesis in the case in which we have the same tail dependence based on 500 repetitions of the simulation based on empirical ranking of the marginals (Fig. A1) and using the true marginals (Fig. A2). For both normalizations the significance level of 5% is in general well attained throughout all numbers of sets. The figures also contain the power of the test when the tail dependence structures are different. After $W = 5$ the power stabilizes and it seems to decrease slightly when the number of sets is chosen too large. We therefore use $W = 5$ throughout the paper. Note that this is only one particular simulation setup, and the results for the optimal number of sets can change depending on sample size and the strength of tail dependence.

<figure><img src="image_12.png" /><figcaption>Figure A1. Simulation study using empirical margins.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Figure A2. Simulation study using true margins.</figcaption></figure>
**Data availability.** ERA5 data are available from the ECMWF website: https://www.ecmwf.int/en/forecasts/datasets/ reanalysis-datasets/era5 (Copernicus Climate Change Ser- vice, 2017). The outputs from the WRF simulations are very large data files and are available from Christoph Raible (christoph.raible@climate.unibe.ch).

**Author contributions.** JZ and PN conceived the idea and study design. PN and SE developed the code for the new metric. CCR pro- vided the model simulations. OM helped with the interpretation of the results. JZ performed all analyses, created all figures, and wrote the first draft. All authors contributed substantially to the writing and revising of the paper.

**Competing interests.** The authors declare that they have no com- peting interests.

**Special issue statement.** This article is part of the special issue “Understanding compound weather and climate events and related impacts (BG/ESD/HESS/NHESS inter-journal SI)”. It is not asso- ciated with a conference.

**Acknowledgements.** This research was supported by a Short- Term Scientific Mission from the European COST Action DAMO- CLES (CA17109). We thank Martina Messmer for creating Fig. 1. The CESM and WRF simulations were performed on the supercom- puting architecture of the Swiss National Supercomputing Centre (CSCS; Lugano, Switzerland). Part of Philippe Naveau’s research was supported by the FRAISE-LEFE-MANU grant and the French Agence National de la Recherche through ANR-Melody and ANR- TRex.

**Financial support.** This research has been supported by the Swiss National Science Foundation (grant nos. 179876, 172745, and 178751), the Helmholtz Initiative and Networking Fund (Young Investigator Group COMPOUNDX, grant agreement VH-NG- 1537), and the European Cooperation in Science and Technology (grant no. CA17109).

**Review statement.** This paper was edited by Gabriele Messori and reviewed by Theophile Caby and two anonymous referees.

**References**

Ban, N., Schmidli, J., and Schaer, C.: Evaluation of the convection- resolving regional climate modeling approach in decade- long simulations, J. Geophys. Res.-Atmos., 119, 889–7907, https://doi.org/10.1002/2014JD021478, 2014.

Barry, R. G.: Mountain weather and climate, Cambridge University Press, Cambridge, UK, 2008.

Bevacqua, E., Maraun, D., Vousdoukas, M. I., Voukouvalas, E., Vrac, M., Mentaschi, L., and Widmann, M.: Higher probability of compound flooding from precipitation and storm surge in Eu- rope under anthropogenic climate change, Science Advances, 5, eaaw5531, https://doi.org/10.1126/sciadv.aaw5531, 2019.

Bracegirdle, T. J., Shuckburgh, E., Sallee, J.-B., Wang, Z., Meijers, A. J. S., Bruneau, N., Phillips, T., and Wilcox, L. J.: Assessment of surface winds over the Atlantic, Indian, and Pacific Ocean sec- tors of the Southern Ocean in CMIP5 models: historical bias, forcing response, and state dependence, J. Geophys. Res-Atmos., 118, 547–562, https://doi.org/10.1002/jgrd.50153, 2013.

Catto, J. L. and Pfahl, S.: The importance of fronts for ex- treme precipitation, J. Geophys. Res-Atmos., 118, 10791–10801, https://doi.org/10.1002/jgrd.50852, 2013.

Champagne, O., Leduc, M., Coulibaly, P., and Arain, M. A.: Winter hydrometeorological extreme events modulated by large-scale atmospheric circulation in southern Ontario, Earth Syst. Dynam., 11, 301–318, https://doi.org/10.5194/esd-11-301-2020, 2020.

Coles, S.: An introduction to statistical modeling of extreme val- ues, Springer, London, UK, https://doi.org/10.1007/978-1-4471- 3675-0, 2001.

Coles, S., Heffernan, J., and Tawn, J.: Dependence measures for extreme value analyses, Extremes, 2, 339–365, 1999.

Copernicus Climate Change Service (C3S): ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate, https: //www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5 (last access: 28 May 2020), Copernicus Climate Change Service, Reading, UK, 2017.

Couasnon, A., Eilander, D., Muis, S., Veldkamp, T. I. E., Haigh, I. D., Wahl, T., Winsemius, H. C., and Ward, P. J.: Measuring compound flood potential from river discharge and storm surge extremes at the global scale, Nat. Hazards Earth Syst. Sci., 20, 489–504, https://doi.org/10.5194/nhess-20-489-2020, 2020.

Davison, A. and Huser, R.: Statistics of Extremes, Annu. Rev. Stat. Appl., 2, 203–235, 2015.

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bid- lot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Källberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011.

De Luca, P., Messori, G., Pons, F. M. E., and Faranda, D.: Dynam- ical systems theory sheds new light on compound climate ex- tremes in Europe and Eastern North America, Q. J. Roy. Meteor. Soc., 146, 1636–1650, https://doi.org/10.1002/qj.3757, 2020.

Deser, C., Lehner, F., Rodgers, K. B., Ault, T., Delworth, T. L., DiNezio, P. N., Fiore, A., Frankignoul, C., Fyfe, J. C., Hor- ton, D. E., Kay, J. E., Knutti, R., Lovenduski, N. S., Marotzke, J., McKinnon, K. A., Minobe, S., Randerson, J., Screen, J. A., Simpson, I. R., and Ting, M.: Insights from Earth system model initial-condition large ensembles and future prospects, Nat. Clim. Change, 10, 1–10, https://doi.org/10.1038/s41558-020-0731-2, 2020.
Embrechts, P., Klüppelberg, C., and Mikosch, T.: Modelling Ex-
tremal Events: for Insurance and Finance, Springer, London, UK,
1997.

Engelke, S. and Ivanovs, J.: Sparse Structures for Multivariate Extremes, Annu. Rev. Stat. Appl., 8, https://doi.org/10.1146/annurev-statistics-040620-041554, 2021.

Ferreira, A. and de Haan, L.: On the block maxima method in ex-
treme value theory: PWM estimators, Ann. Stat., 43, 276–298,
2015.

Fink, A. H., Brücher, T., Ermert, V., Krüger, A., and Pinto, J. G.: The European storm Kyrill in January 2007: synoptic evolution, meteorological impacts and some considerations with respect to climate change, Nat. Hazards Earth Syst. Sci., 9, 405–423, https://doi.org/10.5194/nhess-9-405-2009, 2009.

Flato, G., Marotzke, J., Abiodun, B., Braconnot, P., Chou, S., Collins, W., Cox, P., Driouech, F., Emori, S., Eyring, V., Forest, C., Gleckler, P., Guilyardi, E., Jakob, C., Kattsov, V., Reason, C., and Rummukainen, M.: Evaluation of Climate Models, in: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, USA, 741–866, https://doi.org/10.1017/CBO9781107415324.020, 2013.

François, B., Vrac, M., Cannon, A. J., Robin, Y., and Allard, D.: Multivariate bias corrections of climate simulations: which benefits for which losses?, Earth Syst. Dynam., 11, 537–562, https://doi.org/10.5194/esd-11-537-2020, 2020.

Gilleland, E. and Katz, R. W.: extRemes 2.0: An Extreme
Value Analysis Package in R, J. Stat. Softw., 72, 1–39,
https://doi.org/10.18637/jss.v072.i08, 2016.

Gómez-Navarro, J. J., Raible, C. C., and Dierer, S.: Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain, Geosci. Model Dev., 8, 3349–3363, https://doi.org/10.5194/gmd-8-3349-2015, 2015.

Gómez-Navarro, J. J., Raible, C. C., Bozhinova, D., Martius, O., García Valero, J. A., and Montávez, J. P.: A new region-aware bias-correction method for simulated precipitation in areas of complex orography, Geosci. Model Dev., 11, 2231–2247, https://doi.org/10.5194/gmd-11-2231-2018, 2018.

Graf, M., Scherrer, S. C., Schwierz, C., Begert, M., Martius, O., Raible, C. C., and Brönnimann, S.: Near-surface mean wind in Switzerland: Climatology, climate model evaluation and future scenarios, Int. J. Climatol., 39, 4798–4810, https://doi.org/10.1002/joc.6108, 2019.

Hendry, A., Haigh, I. D., Nicholls, R. J., Winter, H., Neal, R., Wahl, T., Joly-Laugel, A., and Darby, S. E.: Assessing the characteristics and drivers of compound flooding events around the UK coast, Hydrol. Earth Syst. Sci., 23, 3117–3139, https://doi.org/10.5194/hess-23-3117-2019, 2019.

Hong, S. and Lim, J.: The WRF single-moment 6-class micro-
physics scheme (WSM6), Journal of Korean Meteorology So-
ciety, 42, 129–151, 2020.

Hurrell, J. W., Holland, M. M., Gent, P. R., Ghan, S., Kay, J. E., Kushner, P. J., Lamarque, J. F., Large, W. G., Lawrence, D.,

Lindsay, K., Lipscomb, W. H., Long, M. C., Mahowald, N., Marsh, D. R., Neale, R. B., Rasch, P., Vavrus, S., Vertenstein, M., Bader, D., Collins, W. D., Hack, J. J., Kiehl, J., and Marshall, S.: The Community Earth System Model A Framework for Collaborative Research, B. Am. Meteorol. Soc., 94, 1339–1360, https://doi.org/10.1175/BAMS-D-12-00121.1, 2013.

Huser, R. and Wadsworth, J. L.: Advances in Statistical Mod-
eling of Spatial Extremes, WIREs Comput. Stat., e1537,
https://doi.org/10.1002/wics.1537, 2020.

Jimenez, P. A. and Dudhia, J.: Improving the Representation of Re-
solved and Unresolved Topographic Effects on Surface Wind
in the WRF Model, J. Appl. Meteorol. Clim., 51, 300–316,
https://doi.org/10.1175/JAMC-D-11-084.1, 2012.

Jones, P. D., Harpham, C., Troccoli, A., Gschwind, B., Ranchin, T., Wald, L., Goodess, C. M., and Dorling, S.: Using ERA-Interim reanalysis for creating datasets of energy-relevant climate variables, Earth Syst. Sci. Data, 9, 471–495, https://doi.org/10.5194/essd-9-471-2017, 2017.

Katz, R. W., Parlange, M. B., and Naveau, P.: Statistics of extremes in hydrology, Adv. Water Resour., 25, 1287–1304, 2002.

Ledford, A. W. and Tawn, J. A.: Modelling dependence within joint tail regions, J. R. Stat. Soc. Ser. B Stat. Methodol., 59, 475–499, 1997.

Lehner, F., Joos, F., Raible, C. C., Mignot, J., Born, A., Keller, K. M., and Stocker, T. F.: Climate and carbon cycle dynamics in a CESM simulation from 850 to 2100 CE, Earth Syst. Dynam., 6, 411–434, https://doi.org/10.5194/esd-6-411-2015, 2015.

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

Lin, N., Emanuel, K. A., Smith, J. A., and Vanmar-
cke, E.: Risk assessment of hurricane storm surge for
New York City, J. Geophys. Res.-Atmos., 115, D18121,
https://doi.org/10.1029/2009JD013630, 2010.

Manning, C., Widmann, M., Bevacqua, E., Loon, A. F. V., Maraun, D., and Vrac, M.: Increased probability of compound long-duration dry and hot events in Europe during summer (1950–2013), Environ. Res. Lett., 14, 094006, https://doi.org/10.1088/1748-9326/ab23bf, 2019.

Maraun, D., Shepherd, T. G., Widmann, M., Zappa, G., Walton, D., Gutierrez, J. M., Hagemann, S., Richter, I., Soares, P. M. M., Hall, A., and Mearns, L. O.: Towards process-informed bias correction of climate change simulations, Nat. Clim. Change, 7, 764–773, https://doi.org/10.1038/nclimate3418, 2017.

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

Mazdiyasni, O. and AghaKouchak, A.: Substantial increase in concurrent droughts and heatwaves in the United States, P. Natl. Acad. Sci. USA, 112, 11484–11489, https://doi.org/10.1073/pnas.1422945112, 2015.

Messmer, M., Gómez-Navarro, J. J., and Raible, C. C.: Sensitiv-
ity experiments on the response of Vb cyclones to sea surface
temperature and soil moisture changes, Earth Syst. Dynam., 8,
477–493, https://doi.org/10.5194/esd-8-477-2017, 2017.

Messmer, M., Raible, C. C., and Gómez-Navarro, J. J.: Impact of climate change on the climatology of Vb cyclones, Tellus A, 72, 1–18, https://doi.org/10.1080/16000870.2020.1724021, 2020.
Mitchell, T. D. and Jones, P. D.: An improved method of constructing a database of monthly climate observations and associated high-resolution grids, Int. J. Climatol., 25, 693–712, https://doi.org/10.1002/joc.1181, 2005.

Musselman, K., Lehner, F., Ikeda, K., Clark, M., Prein, A., Liu, C., Barlage, M., and Rasmussen, R.: Projected increases and shifts in rain-on-snow flood risk over western North America, Nat. Clim. Change, 8, 808–812, https://doi.org/10.1038/s41558-018-0236-4, 2018.

Naveau, P., Guillou, A., and Rietsch, T.: A non-parametric entropy-based approach to detect changes in climate extremes, J. Roy. Stat. Soc. B, 76, 861–884, 2014.

Naveau, P., Hannart, A., and Ribes, A.: Statistical Methods for Extreme Event Attribution in Climate Science, Annu. Rev. Stat. Appl., 7, 89–110, https://doi.org/10.1146/annurev-statistics-031219-041314, 2020.

Pantillon, F., Adler, B., Corsmeier, U., Knippertz, P., Wieser, A., and Hansen, A.: Formation of Wind Gusts in an Extratropical Cyclone in Light of Doppler Lidar Observations and Large-Eddy Simulations, Mon. Weather Rev., 148, 353–375, https://doi.org/10.1175/Mwr-D-19-0241.1, 2020.

Panziera, L. and Germann, U.: The relation between airflow and orographic precipitation on the southern side of the Alps as revealed by weather radar, Q. J. Roy. Meteor. Soc., 136, 222–238, https://doi.org/10.1002/qj.544, 2010.

Parton, G., Dore, A., and Vaughan, G.: A climatology of mid-tropospheric mesoscale strong wind events as observed by the MST radar, Aberystwyth, Meteorol. Appl., 17, 340–354, https://doi.org/10.1002/met.203, 2010.

Pfahl, S., Madonna, E., Boettcher, M., Joos, H., and Wernli, H.: Warm Conveyor Belts in the ERA-Interim Dataset (1979–2010). Part II: Moisture Origin and Relevance for Precipitation, J. Climate, 27, 27–40, https://doi.org/10.1175/Jcli-D-13-00223.1, 2014.

Poon, S.-H., Rockinger, M., and Tawn, J.: Extreme value dependence in financial markets: Diagnostics, models, and financial implications, Rev. Financ. Stud., 17, 581–610, 2003.

Poschlod, B., Zscheischler, J., Sillmann, J., Wood, R. R., and Ludwig, R.: Climate change effects on hydrometeorological compound events over southern Norway, Weather and Climate Extremes, 28, 100253, https://doi.org/10.1016/j.wace.2020.100253, 2020.

Raible, C. C., Messmer, M., Lehner, F., Stocker, T. F., and Blender, R.: Extratropical cyclone statistics during the last millennium and the 21st century, Clim. Past, 14, 1499–1514, https://doi.org/10.5194/cp-14-1499-2018, 2018.

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, https://doi.org/10.1002/qj.2531, 2015.

Ridder, N., Pitman, A., Westra, S., Ukkola, A., Do, H., Bador, M., Hirsch, A., Evans, J., Luca, A. D., and Zscheischler, J.: Global hotspots for the occurrence of compound events, Nat. Commun., 11, 5956, https://doi.org/10.1038/s41467-020-19639-3, 2020.

Sibuya, M.: Bivariate extreme statistics, Ann. I. Stat. Math., 11, 195–210, 1960.

Sippel, S., Zscheischler, J., Mahecha, M. D., Orth, R., Reichstein, M., Vogel, M., and Seneviratne, S. I.: Refining multi-model projections of temperature extremes by evaluation against land-

atmosphere coupling diagnostics, Earth Syst. Dynam., 8, 387–403, https://doi.org/10.5194/esd-8-387-2017, 2017.

Sippel, S., Reichstein, M., Ma, X., Mahecha, M. D., Lange, H., Flach, M., and Frank, D.: Drought, Heat, and the Carbon Cycle: a Review, Current Climate Change Reports, 4, 266–286, https://doi.org/10.1007/s40641-018-0103-4, 2018.

Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Barker, D. M., Wang, W., and Powers, J. G.: A description of the advanced research WRF version 3, Technical Report, TN-475+STR, National Center for Atmospheric Research, Boulder, CO, USA, 113 pp., 2008.

Stephenson, A. G.: evd, Extreme Value Distributions, R News, 2.0, availabe at: https://CRAN.R-project.org/doc/Rnews/ (last access: 28 May 2020), 2002.

Stucki, P., Dierer, S., Welker, C., Gómez-Navarro, J. J., Raible, C. C., Martius, O., and Brönnimann, S.: Evaluation of down-scaled wind speeds and parameterised gusts for recent and historical windstorms in Switzerland, Tellus A, 68, 31820, https://doi.org/10.3402/tellusa.v68.31820, 2016.

Sutanto, S. J., Vitolo, C., Napoli, C. D., D'Andrea, M., and Lanen, H. A. V.: Heatwaves, droughts, and fires: Exploring compound and cascading dry hazards at the pan-European scale, Environ. Int., 134, 105276, https://doi.org/10.1016/j.envint.2019.105276, 2020.

Telesca, L., Guignard, F., Laib, M., and Kanevski, M.: Analysis of temporal properties of extremes of wind measurements from 132 stations over Switzerland, Renew. Energ., 145, 1091–1103, https://doi.org/10.1016/j.renene.2019.06.089, 2020.

Vignotto, E., Engelke, S., and Zscheischler, J.: Clustering bivariate dependences in the extremes of climate variables, Weather and Climate Extremes, in review, 2020.

Wahl, T., Jain, S., Bender, J., Meyers, S. D., and Luther, M. E.: Increasing risk of compound flooding from storm surge and rainfall for major US cities, Nat. Clim. Change, 5, 1–6, https://doi.org/10.1038/nclimate2736, 2015.

Wang, J., Chen, Y., Tett, S. F., Yan, Z., Zhai, P., Feng, J., and Xia, J.: Anthropogenically-driven increases in the risks of summertime compound hot extremes, Nat. Commun., 11, https://doi.org/10.1038/s41467-019-14233-8, 2020.

Ward, P. J., Couasnon, A., Eilander, D., Haigh, I. D., Hendry, A., Muis, S., Veldkamp, T. I., Winsemius, H. C., and Wahl, T.: Dependence between high sea-level and high river discharge increases flood hazard in global deltas and estuaries, Environ. Res. Lett., 13, 084012, https://doi.org/10.1088/1748-9326/aad400, 2018.

Whiteman, C. D.: Mountain meteorology fundamentals and applications, Oxford University Press, New York, USA, 2000.

Zheng, F., Westra, S., and Sisson, S. A.: Quantifying the dependence between extreme rainfall and storm surge in the coastal zone, J. Hydrol., 505, 172–187, 2013.

Zscheischler, J. and Fischer, E.: The record-breaking compound hot and dry 2018 growing season in Germany, Weather and Climate Extremes, 19, 100270, https://doi.org/10.1007/s00484-020-01951-8, 2020.

Zscheischler, J. and Seneviratne, S. I.: Dependence of drivers affects risks associated with compound events, Science Advances, 3, e1700263, https://doi.org/10.1126/sciadv.1700263, 2017.

Zscheischler, J., Michalak, A. M., Schwalm, C., Mahecha, M. D., Huntzinger, D. N., Reichstein, M., Berthier, G., Ciais, P.,
Cook, R. B., El-Masri, B., Huang, M., Ito, A., Jain, A., King, A., Lei, H., Lu, C., Mao, J., Peng, S., Poulter, B., Ricciuto, D., Shi, X., Tao, B., Tian, H., Viovy, N., Wang, W., Wei, Y., Yang, J., and Zeng, N.: Impact of large-scale climate extremes on biospheric carbon fluxes: An intercomparison based on MsTMIP data, Global Biogeochem. Cy., 28, 585–600, https://doi.org/10.1002/2014GB004826, 2014.

Zscheischler, J., Westra, S., Hurk, B. J., Seneviratne, S. I., Ward, P. J., Pitman, A., Agha Kouchak, A., Bresch, D. N., Leonard, M., Wahl, T., and Zhang, X.: Future climate risk from compound events, Nat. Clim. Change, 8, 469–477, 2018.

Zscheischler, J., Fischer, E. M., and Lange, S.: The effect of univariate bias adjustment on multivariate hazard estimates, Earth Syst. Dynam., 10, 31–43, https://doi.org/10.5194/esd-10-31-2019, 2019.

Zscheischler, J., Martius, O., Westra, S., Bevacqua, E., Raymond, C., Horton, R. M., van den Hurk, B., Agha Kouchak, A., Jézéquel, A., Mahecha, M. D., Maraun, D., Ramos, A. M., Rider, N., Thiery, W., and Vignotto, E.: A typology of compound weather and climate events, Nature Reviews Earth and Environment, 1, 333–347, https://doi.org/10.1038/s43017-020-0060-z, 2020.