Future projections of Mediterranean cyclone characteristics using the Med-CORDEX ensemble of coupled regional climate system models

Here, we analyze future projections of cyclone activity in the Mediterranean region at the end of the twenty-first century based on an ensemble of state-of-the-art fully-coupled Regional Climate System Models (RCSMs) from the Med-CORDEX initiative under the Representative Concentration Pathway (RCP) 8.5. Despite some noticeable biases, all the RCSMs capture spatial patterns and cyclone activity key characteristics in the region and thus all of them can be considered as plausible representations of the future evolution of Mediterranean cyclones. In general, the RCSMs show at the end of the twenty-first century a decrease in the number and an overall weakening of cyclones moving across the Mediterranean. Five out of seven RCSMs simulate also a decrease of the mean size of the systems. Moreover, in agreement with what already observed in CMIP5 projections for the area, the models suggest an increase in the Central part of the Mediterranean region and a decrease in the South-eastern part of the region in the cyclone-related wind speed and precipitation rate. These rather two opposite tendencies observed in the precipitation should compensate and amplify, respectively, the effect of the overall reduction of the frequency of cyclones on the water budget over the Central and South-eastern part of the region. A pronounced inter-model spread among the RCSMs emerges for the projected changes in the cyclone adjusted deepening rate, seasonal cycle occurrence and associated precipitation and wind patterns over some areas of the basin such as Ionian Sea and Iberian Peninsula. The differences observed appear to be determined by the driving Global Circulation Model (GCM) and influenced by the RCSM physics and internal variability. These results point to the importance of (1) better characterizing the range of plausible futures by relying on ensembles of models that explore well the existing diversity of GCMs and RCSMs as well as the climate natural variability and (2) better understanding the driving mechanisms of the future evolution of Mediterranean cyclones properties.


Introduction
The Mediterranean region is located at mid-latitudes and includes the continental Southern Europe, Northern Africa and the Middle East, all surrounding the Mediterranean Sea ( Fig. 1). The presence of this relatively large water mass acts as a source of moisture and heat to atmospheric weather systems. Furthermore, it imposes sharp land-sea transitions that strongly affect both atmospheric and oceanic dynamics over the region (see Lionello et al. 2006Lionello et al. , 2012Ulbrich et al. 2012;Flaounas et al. 2021b;Hochman et al. 2021).
Extra-tropical cyclones play a fundamental role in shaping the weather of the Mediterranean region ( Fig. 1, Lionello et al. 2006;Nissen et al. 2010;Ulbrich et al. 2012;Lionello et al. 2012;Flaounas et al. 2013;Reale and Lionello 2013;Lionello et al. 2016 and references inside; Flaounas et al. 2017;Flaounas et al. 2021b). The location and complex orography of the region allow the development of different cyclonic systems with a marked seasonal cycle. These 1 3 include for example: lee cyclones, which occur mostly in winter at the lee of the Alps, Atlas Mountains and Pyrenees, or thermal lows which usually occur in summer (Trigo et al. 1999;Trigo 2006;Lionello et al. 2002;Lionello et al. 2016 and references inside; Flaounas et al. 2018;Flaounas et al. 2021b).
Cyclogenesis in the region has been shown to be driven by the intrusion of upper-tropospheric systems. These result from Rossby wave breaking over the Atlantic Ocean with an active role played by deep convection close to the cyclone center. Deep convection acts as a source of potential vorticity for the cyclone (see Flaounas et al. 2021a, b). Past studies have identified several cyclogenetic areas in the region: the Atlas Mountains, the Gulf of Genoa, the Aegean Sea, the Iberian peninsula, the Black Sea and the Eastern Mediterranean (Alpert et al. 1990;Trigo et al. 1999;Lionello et al. 2006Lionello et al. , 2016Campins et al. 2010;Flaounas et al. 2018Flaounas et al. , 2021bBuzzi et al. 2020). Moreover, approximately 11% of cyclones moving across the region originate in the Atlantic Ocean (Lionello et al. 2016). Mediterranean systems are generally shallower, slower, smaller and less persistent with respect to other mid-latitude cyclones (Lionello et al. 2016).
Cyclone activity in the Mediterranean has important effects on the hydrological cycle and the distribution of the precipitation in the region (Flaounas et al. , 2017. Furthermore, some cyclones are also associated with detrimental extreme events, such as heavy precipitation, windstorms, storm surges and marine storminess (Lionello et al. 2006Nissen et al. 2010;Liberato et al. 2011;Reale and Lionello 2013;Flaounas et al. 2017;Lionello et al. 2019;Lionello et al. 2020;Flaounas et al. 2021b).
Several studies have focused on characterizing Mediterranean cyclone activity in past and present periods using observations, reanalyses, global and/or regional climate models (Trigo et al. 1999;Maheras et al. 2001;Trigo 2006;Lionello et al. 2002Lionello et al. , 2016Flaounas et al. 2018). Others have focused on the future projections of cyclonic activity over the region (Somot 2005;Lionello and Giorgi 2007;Raible et al. 2010;Cavicchia et al. 2014;Nissen et al. 2014;Zappa et al. 2015;Hochman et al. 2018aHochman et al. , b, 2020. For example, Lionello and Giorgi (2007) have demonstrated weaker (stronger) cyclone activity and a decrease (increase) in precipitation over the Eastern (North-Western) Mediterranean for the end of the twenty-first century under SRES emission scenarios A2 and B2. On the other hand, Raible et al. (2010) found in the ECHAM GCM a decrease in cyclone activity over the Western Mediterranean by the end of the century, no significant change in the Eastern Mediterranean and a slight decrease in cyclone intensity for the entire region. Nissen et al. (2014) using a suite of 7 simulations of two coupled GCMs (CMCC-CM and ECHAM5-MPIOM) under different emissions scenarios showed a decrease in the number of cyclones crossing the region and associated windstorms. Zappa et al. (2015) analyzed a suite of 17 CMIP5 models under Representative Concentration Pathway (RCP) 8.5 emission scenario and found an overall reduction of cyclone number in the period 2082-2099 (with respect to 1976-2005) over the Mediterranean region, with a concomitant increase of the associated precipitation over the western part of the Mediterranean basin. Finally, Hochman et al. (2018aHochman et al. ( , 2020 have found a significant decrease in Eastern Mediterranean cyclone frequency (− 35%) and associated daily precipitation (-26%) using 8 CMIP5 models for the end of the twenty-first century under RCP8.5.
There is no general consensus regarding the effect of global warming on cyclone activity over the Mediterranean Fig. 1 The Mediterranean region (blue line) and sub-areas identified for the cyclogenesis processes, AFR (Northern Africa), MED (Mediterranean), ATL (North Atlantic), EEME (Eastern Europe/ Middle East). Original atmospheric domains of the models (Table 1) are shown: ROM (red line), EBUPOM2c (red dashed line), CNRM-RCSM4 (cyan line), LMDZ4NEMOMED8 (green line), COS-MOMED (yellow line). The part of the ERA5 global domain covering the Mediterranean region is shown in blue line basin. Moreover, previous studies were based on a single or an ensemble of low resolution uncoupled/coupled GCMs (ECHAM5 and CMIP5) or single uncoupled regional climate model. This: (i) poses some limitations on the representation of cyclone activity in the region which has been shown to be influenced by the model resolution and the frequency of the atmosphere and ocean coupling (Akhtar et al. 2014;Flaounas et al. 2018;Gaertner et al. 2018), (ii) leads to omit the various sources of uncertainty related for example to the choice of the regional climate models, driving GCM and natural variability. Additionally, all the previous studies have been focused on the analysis of the changes in the spatial distribution of cyclones and related precipitation, without considering changes in the cyclones' features such as size, depth or deepening rate. Since cyclones play a fundamental role in the dynamics of extreme weather in the Mediterranean it is important to consider changes in the features of systems crossing the region and to address the sources of uncertainties affecting the projections in the area by considering different possible climate futures through a multi-model approach.
In recent years the Med-CORDEX initiative (Ruti et al. 2016;Somot et al. 2018; www. medco rdex. eu) has been launched to develop fully-coupled Regional Climate System Models (RCSMs) for the Mediterranean basin under the framework of the Coordinated Regional Downscaling EXperiment (CORDEX, Giorgi et al. 2009). More specifically, the Med-CORDEX RCSMs include a high-frequency coupling of the ocean and atmosphere components (6 min-1 day), both represented at relatively high-resolution (25-50 km for the atmosphere, 6-20 km for the ocean). One of the main aims of this initiative is to use these coupled models to produce coordinated high-resolution ensembles of climate projections for the Mediterranean basin. An initial ensemble has already been generated including seven simulation runs performed with 5 different RCSMs driven by GCM projections from the Coupled Model Intercomparison Project Phase 5 (CMIP5, Taylor et al. 2012) under RCP 8.5 emission scenario, considered as the high end (see Darmaraki et al. 2019 for description and a first use of this ensemble).
In this study we use this ensemble to revisit the still open question concerning the response of Mediterranean cyclone activity to global warming. Our focus is on a variety of cyclone characteristics, such as frequency, intensity, track, associated precipitation yield and wind speed. Particularly, their changes by the end of the twenty-first century are compared to historical reference values. Moreover, we discuss the sources of uncertainties affecting the projections in the region with special attention dedicated to the driving GCM and RCSM setup and internal features.
The paper is organized as follows: Sect. 2 describes the main characteristics of the models, our approach for detecting and tracking Mediterranean cyclones and the procedure adopted to associate precipitation and wind patterns to each system. Section 3 presents an assessment of the ability of RCSMs to reproduce the main observed characteristics of Mediterranean cyclones and related precipitation and wind patterns for present climate conditions, along with an analysis of future projections of cyclone activity in response to increased greenhouse gas concentrations. Section 4 discusses the possible sources of the differences observed among RCSMs projections while Sect. 5 presents our conclusions and the possible future research avenues.

ERA5 reanalysis and Med-CORDEX models
We analyze cyclone activity in the Mediterranean region ( Fig. 1) using a set of 7 RCSM simulations produced as part of the Med-CORDEX initiative (Ruti et al. 2016; www. medco rdex. eu) and the ERA5 reanalysis (Hersbach et al. 2020) as reference dataset. The variables considered here are the mean sea level pressure (MSLP), precipitation and wind fields. ERA5 is used for the sake of consistency, since it is also employed for evaluating cyclone activity. The simulations are based on 5 different coupled RCSMs (ROM, EBUPOM2c, CNRM-RCSM4, LMDZ4NEMOMED8 and COSMOMED) driven by five different GCMs (MPI-ESM-LR, MPI-ESM-MR, CNRM-CM5, IPSL-CM5A-MR and CMCC-CM), being part of CMIP5 (Table 1). One interesting new feature of this ensemble is represented by the relatively high-resolution of the oceanic (eddy-resolving) and atmospheric components as well as the high coupling frequency between the two. The representation of the air-sea exchanges through this high-resolution interface represents a clear advantage for the representation of cyclone dynamics in the region, which is influenced by ocean-atmosphere interactions (Flaounas et al. 2021a). CMIP5 GCM models have been shown to be able to generate a Mediterranean storm track and to reasonably represent the observed precipitation, evaporation, moisture fluxes and wind patterns in the region (e.g. Kelley et al. 2012;Seager et al. 2014;Nissen et al. 2014;Zappa et al. 2015).
All the RCSMs are driven by the respective GCM for a historical period  and for a future period (2006-2100) under emission scenario RCP8.5 (Taylor et al. 2012). The ensemble of simulations we consider consists of the GCM/RCSM combinations shown in Table 1: the LMDZ4NEMOMED8 RCSM (hereafter referred to LMD-MPI, LMD-CNRM and LMD-IPSL) driven by three different GCMs (MPI-ESM-MR, CNRM-CM5, IPSL-CM5A-MR), the CNRM-RCSM4 RCSM driven by the CNRM-CM5 GCM (hereafter CNRM-RCSM4), ROM and EBUPOM2c RCSMs driven by MPI-ESM-LR (hereafter ROM and EBU-POM) and COSMOMED driven by CMCC-CM (hereafter COSMOMED). For more details on each simulation the reader is referred to the bibliography and the information reported in Table 1. Note that an error has been recently reported concerning the CNRM-CM5 GCM files for the historical period for CNRM-RCSM4 and LMD-CNRM (http:// www. umr-cnrm. fr/ cmip5/ spip. php? artic le24). The atmosphere data at the 6-hourly frequency on the model level (6hLev) for the variables ps, ta, ua, va, hus, that were used as atmospheric lateral boundary conditions, do not come from the member r1i1p1 as expected. This simulation has been included in the ensemble given that this error likely has no significant effect on the longterm climate change signal nor on the cyclone analysis performed here (Hochman et al. 2018a;Darmaraki et al. 2019). For additional information on the ocean scheme of the Med-CORDEX RCSMs the reader is referred to Darmaraki et al. (2019) and Soto-Navarro et al. (2020).
Hereafter each simulation analyzed, except if not differently reported, will be identified with the short name of the corresponding RCSM (Table 1, first column).

Cyclone detection and tracking algorithm
Over the last two decades several approaches based on single or multiple automatic tracking algorithms have been applied to characterize cyclones activity in the Mediterranean region (e.g. Trigo et al. 1999;Lionello et al. 2002;Trigo 2006;Nissen et al. 2010;Kouroutzoglou et al. 2010Kouroutzoglou et al. , 2011Reale and Lionello 2013;Lionello et al. 2016;Flaounas et al. 2018). Table 1 Characteristics of the Med-Cordex regional coupled models (RCSMs) All the simulations cover the period 1950-2005 (HIST) and 2006-2100 (RCP8.5). More information about each RCSM and GCM can be retrieved in the bibliography reported in the first column and in Darmaraki et al. (2019). The MEDAR dataset is extensively described in Rixen et al. (2005). PHC is discussed in Steele et al. (2001)  The adoption of different schemes for detecting and tracking cyclones reflects the existence of a wide variety of definitions for a cyclone, which are based on different physical variables such as mean sea level pressure or relative vorticity (e.g. Neu et al. 2013;Hoskins and Hodges 2019a, b).
In this work, we identify cyclones using an automatic scheme specifically developed for the Mediterranean region (Lionello et al. 2002), and later adapted to hemispheric scales Neu et al. 2013). This algorithm has been extensively applied in previous studies (Lionello et al. 2002;Lionello and Giorgi 2007;Reale and Lionello 2013;Neu et al. 2013;Ulbrich et al. 2013;Lionello et al. 2016;Pinto et al. 2016;Flaounas et al. 2018;Lionello et al. 2019;Reale et al. 2019;Reboita et al. 2020;Lionello et al. 2020;Reboita et al. 2021) and it is based on the search of pressure minima in a sequence of MSLP gridded fields. The procedure involves the partition of the MSLP field at time t 0 in a certain number of systems by identifying sets of steepest paths leading to the same MSLP minimum. Each point in the MSLP field is connected to the lowest of the 8 nearest-neighbor points, repeating the procedure until a minimum is reached. All the points crossed by a path leading to the same minimum are assigned to the same cyclone. Inflation of the number of detected cyclones is avoided by merging shallow minima to the large adjacent system. The actual position of the cyclone center is a weighted average of the points around the pressure minimum. The trajectory of each system is built by joining the location of the same cyclone center in successive maps. More specifically a box is associated with each center. This box in the meridional direction is centered at cyclone center with a width given by 1.4 × R where R is the average radius of the system, while in the zonal direction the width is asymmetric, with its western side at a distance of 0.3 × R, and its eastern side at a distance given by 0.8 × R. A cyclone in the map at the time t 0 + 1 is assumed to be the continuation of the cyclone in the map if its center is inside such a box. When no center is identified inside the box, the track is terminated. For a detailed description of the tracking algorithm the reader is referred to Lionello et al. (2002) and Reale and Lionello (2013). The final outcome of this automatic algorithm is a list of cyclone tracks, their position as a function of time and the temporal evolution of variables such as: MSLP minimum, Laplacian gradient, depth and radius. In particular the radius of a generic storm (R) is defined as: where x 2 i represents the distance of a generic point from the center of the cyclone, with N representing the total number of points associated to a cyclone. (1) Recent studies based on IMILAST dataset (e.g. Neu et al. 2013;Ulbrich et al. 2013;Pinto et al. 2016;Lionello et al. 2016;Reale et al. 2019) have shown that the choice of a specific tracking scheme could influence the representation of the characteristics of the systems in that region, suggesting that the use of different tracking algorithms on the same original dataset could be a better approach for describing the life cycle of extratropical cyclones. In this study however, we focus on the ability of the ensemble members to reproduce the main characteristics of cyclones in the Mediterranean, rather than on the sensitivity of track climatology to the different cyclone definitions. The sensitivity of our climatology and climate change signal observed with respect to the definition of cyclone is out of the scope of the present manuscript and will be eventually discussed in a future study.
The different RCSMs simulations were completed over a different spatial domain and with different resolution, ranging from 25 to 50 km in their atmospheric component ( Fig. 1 and Table 1). In order to provide a fair comparison between ERA5 and each RCSM (e.g. Kouroutzoglou et al. 2010Kouroutzoglou et al. , 2011Flaounas et al. 2018) all the 6-hourly model outputs were interpolated onto the ERA5 grid (0.25°, approximately 30 km) and on the region shown in blue line in Fig. 1. Moreover, we identified four sub-areas encompassing the main sources of the cyclones crossing the region: North Atlantic (ATL), Northern Africa (AFR), Mediterranean region (MED) and Eastern Europe/Middle East (EEME; Lionello et al. 2016).
The use of ERA5 represents an advantage with respect to the use of former global reanalysis due to its finer resolution, more suitable to describe small scale systems such as the Mediterranean cyclones (Lionello et al. 2016). On the other hand, the use of high-resolution data as input in the tracking algorithm can lead to an inflation in the number of spurious and weak systems which, however, are unlikely to have a strong impact on the climate dynamics and extremes of the Mediterranean region Flaounas et al. 2018). In order to filter out all these systems, we consider here cyclones that remain for > = 24 h in the region and whose maximum intensity (measured by the MSLP) is reached within the Mediterranean area (MED, Fig. 1). Second, we retain only cyclones with a lifetime longer than 24 h originating above 20° N , with a maximum adjusted deepening rate (Trigo 2006) lower than zero  and the difference between the minimum and maximum of MSLP along their life cycle larger than 12 hPa (Lionello et al. 2016).
Finally, the period October-March (hereafter ONDJFM) has been chosen to analyze cyclone-related precipitation and wind patterns as it is the period of the year where the influence of cyclone activity on these weather variables has been extensively studied and characterized (see Toreti et al. 1 3 2010; Donat et al. 2011;Lionello et al. 2006;Ulbrich et al. 2012;Reale and Lionello 2013;Nissen et al. 2014;Zappa et al. 2015). Although different techniques exist to associate precipitation and wind patterns to a cyclone Nissen et al. 2014;Scoccimarro et al. 2014;Villarini et al. 2014;Zappa et al. 2015;Zhang et al. 2020) here we adopted a procedure similar to that of Zappa et al. (2015), i.e.: precipitation and wind speed values in each grid point are attributed to a cyclone if it occurs within a predefined radius from its center. Based on the size distribution of storms analyzed in Sect. 3.1 and following previous analysis (Zappa et al. 2015), the radius has been chosen equal to ten horizontal degrees. Precipitation and wind speed observed beyond this distance are considered as "background" precipitation and wind and not associated with the cyclone activity itself (Zappa et al. 2015). In the case of precipitation, to better focus on possible changes in the intensity of a rainy event on a short time scale all the values are shown in mm/h.
The climate change signal is evaluated considering the difference between a 25 years period at the end of the twenty-first century in the scenario simulation (2075-2099) and a 25 years period at the end of the historical simulation . Both periods partially encompass the IPCC-AR5 end of the 21st century and baseline time slices respectively (IPCC 2014) and other time slices adopted in previous studies (e.g. Zappa et al. 2015). The statistical significance of the observed differences between the two periods were assessed by means of a Mann-Whitney test at the 95% significance level. Following Lionello et al. (2016) the Mediterranean region has been divided in cells of 1.5° × 1.5°. Then, the number of cyclones crossing each cell is counted only once, even if it stays or crosses the cell for more than one time step. The 1.5° criterion has been chosen to better visualize the spatial distribution of cyclone tracks in the region and to reduce the noise due to the high resolution of original data. Moreover, in order to assess the level of agreement among RCSMs, we use the normalized standard deviation, which is the ratio between the standard deviation of the multi-model mean and the multi-model mean (std multi-model-mean /multi-model mean; Lionello et al. 2016;Reale et al. 2019) and is shown in the panel (j). A small value of normalized standard deviation corresponds to a high level of agreement among the RCSMs (Lionello et al. 2016).

Evaluation of present cyclone activity in the Mediterranean region
The relatively low number of tracks observed in ERA5 (a) with respect to what found in previous multi-tracking/multi-model studies using ERA-Interim as reference (e.g. Lionello et al. 2016;Flaounas et al. 2018) is a direct consequence of the strict criteria adopted for identifying cyclones crossing the region (see Sect. 2.2). The spatial patterns of cyclones tracks are qualitatively reproduced in the Mediterranean region by the multi-model mean (b) and each RCSM (c-i) compared to ERA5 (a). Relative maxima are captured in the Gulf of Genoa, Western Mediterranean, North Africa and in proximity to Cyprus in agreement with previous studies (e.g. Alpert et al. 1990;Trigo et al. 1999;Lionello et al. 2006Lionello et al. , 2016Campins et al. 2010;Ulbrich et al. 2012;Flaounas et al. 2018). The spatial correlation coefficient between each RCSM and ERA5 is high, ranging from 0.73 (CNRM-RCSM4) to 0.89 (LMD-CNRM). The level of agreement among RCSMs is relatively high in the Western and Central Mediterranean region, while it is lower around Cyprus and in general at the boundaries of the Mediterranean Sea (Fig. 2j). Figure 3 shows the relative differences between the multi-model mean (a) and each RCSM (b-h) with ERA5 (Fig. 2a). Most of the RCSMs and the multi-model ensemble underestimate more than 50% the maxima observed in ERA5 over the Gulf of Genoa and the Iberian Peninsula. Of the same magnitude is the overestimation of the signal over Northern Africa and in the Cyprus area. LMD-CNRM (Fig. 3e) and LMD-MPI (Fig. 3f) show higher values of the density of tracks in the western and central part of the basin. On the other hand, LMD-IPSL (Fig. 3h) has an opposite behavior with lower values over most of the two regions. Besides, EBU-POM in Fig. 3c (ROM in Fig. 3b) underestimates (overestimates) the maximum observed in the area around Cyprus. Most of the differences observed are not statistically significant except for a limited number of areas and RCSMs such as in the Gulf of Genoa in the case of EBU-POM (Fig. 3c), CNRM-RCSM4 (Fig. 3d), LMD-IPSL (Fig. 3g) and COSMOMED (Fig. 3h),over the Ionian Sea in the case of EBU-POM (Fig. 3c) and LMD-IPSL (Fig. 3g), and around Cyprus in the case of the multimodel mean (Fig. 3a) and ROM (Fig. 3b).
In all the RCSMs, the majority of cyclones crossing the region (more than 50%) are formed within the Mediterranean basin (Fig. 4a). RCSMs also agree that most of the remaining systems crossing the region originate either over the Atlantic or over Northern Africa, while the contribution from the Middle East and Eastern Europe areas is quite small. However, RCSMs tend to overestimate the percentage of cyclones originating in North Africa (up to 15% in CNRM-RCSM4 with respect to ERA5). At the same time ROM, EBU-POM, LMD-MPI, LMD-CNRM  Despite some relevant overestimation/underestimation observed in each single month (for example − 14% observed in December in LMD-IPSL), the seasonal cycle of cyclone frequency is captured by all RCSMs (Fig. 4b) (Fig. 4d). In the LMD, COSMOMED and ROM simulations, cyclones tend to be deeper than those in ERA5 by about 10 hPa. On the other hand, in the EBU-POM simulation, cyclones are shallower than those in ERA5. Moreover, both LMD and COSMOMED simulations are also characterized by cyclones with higher ADR than ERA5. The differences observed, with the only exception of ROM, are statistically significant.
Finally, concerning the average radius (Fig. 4e) most RCSMs are in agreement with ERA5 in showing that the most likely value for this variable falls in the interval (1100; 1300) km. ROM and COSMOMED tend to simulate wider  Precipitation patterns associated with cyclones are qualitatively well captured by each RCSM in comparison to ERA5. Some relative maxima are observed along the Algerian, Eastern Adriatic, Ionian, Anatolian and Eastern Mediterranean coastlines (e.g. Zappa et al. 2015). Moreover, other precipitation maxima are observed over the Alps and the Pyrenees. The level of agreement among RCSMs is relatively higher over the Mediterranean Sea while it tends to decrease over land.
However, most of the simulations overestimate (more than 50%) precipitation over mountain peaks, in agreement with previous studies (Fig. 6; Fantini et al. 2018;Reale et al. 2020). Another common feature among RCSMs is the underestimation of the precipitation over part of the Adriatic Sea, the Balkan area and Aegean Sea. Moreover, ROM (Fig. 6b) and EBU-POM (Fig. 6c) display underestimations of about 50% precipitation over Western and Central part of the Mediterranean region. CNRM-RCSM4 (Fig. 6d), LMD-CNRM (Fig. 6e), LMD-MPI (Fig. 6f) and COSMOMED (Fig. 6h) simulates values of precipitation almost double with respect to ERA5 (Fig. 5a) over the Iberian Peninsula. Drier conditions are observed in this area in the other ensemble members.
Direction and speed of the winds associated with cyclones are also qualitatively well captured by the multimodel mean (Fig. 7b) and by each RCSM (Fig. 7 c, i) in comparison to ERA5 (Fig. 7a). Multi-model-mean and each RCSM simulation captures the northerly direction of the winds associated with cyclone activity over the Aegean Sea, the westerly winds flowing along the Northern Africa coastlines, the northerly and easterly flows over the gulf of Lions and Northern Adriatic Sea, which correspond to the location of Mistral and Bora wind patterns respectively. Most RCSMs capture the spatial pattern of wind magnitude in the region with some relative maxima observed in the gulf of Lions, Northern African coastlines, Northern Adriatic, Southern Ionian and over the Aegean Sea. As in the case of precipitation, the spread measured by the normalized standard deviation is low over the Mediterranean Sea, while it increases over the land with some maxima. These findings are spatially coherent with the presence of some topographic features such as the Alps and the Pyrenees (Fig. 7j). However, RCSMs tend also to overestimate (the values simulated are almost double with respect to those observed in ERA5) the intensity of wind flowing through the strait of Otranto (Fig. 8). At the same time, they tend to underestimate (~ 50%) the intensity of wind flowing through the Aegean Sea (Fig. 8). Moreover, some differences are also observed among RCSMs concerning the maximum associated with the Mistral wind in the Gulf of Lions. Finally, CNRM-RCSM4 underestimates the intensity of wind flowing over the Southern Ionian and Levantine basin (Fig. 8d).
In summary, although some noticeable quantitative biases exist in the ensemble of Med-CORDEX RCSMs simulations, the models qualitatively capture key cyclone characteristics, in particular, the spatial patterns and there are no extreme outliers. We therefore conclude that all simulations can be considered as plausible representations of Mediterranean cyclones and we therefore use them all for the analysis of the response of cyclone activity to increased greenhouse gases concentrations in the next section. .We define a climate change signal as "robust" in our ensemble if it fulfills two criteria: (i) more than 5 out of 7 models agree on the sign of the change and (ii) the normalized standard deviation of changes observed is (in absolute value) lower than 1 (Akperov et al. 2019). The ensemble average shows a general decrease in the number of cyclones crossing the region of up to ~ 50% (Fig. 9a). The observed decrease is robust over the Tyrrhenian Sea, Southern Italy, part of the Gulf of Genoa, over part of Northern Africa, the Eastern part of the Ionian Sea, the Anatolian Peninsula and in the area between Crete and Cyprus. Despite this overall decrease, an increase in some areas is also evident. ROM (Fig. 9b), EBU-POM (Fig. 9c), CNRM-RCSM4 (Fig. 9d) and LMD-MPI (Fig. 9f) show an increase in cyclone track density over most of Western Europe and the Iberian Peninsula. At the same time, an increase in density of tracks over the Ionian Sea and most of the Levantine basin is projected by CNRM-RCSM4 (Fig. 9d) and only over part of the Ionian Sea in LMD-CNRM (Fig. 9e). EBU-POM shows a positive signal located over the Middle-East (Fig. 9c). In general, the observed statistically significant differences between the two periods in each RCSM are mostly negative and located over the Mediterranean Sea (Fig. 9b-h). Among the RCSMs the only statistically significant signal encompassing most of the Mediterranean storm track is observed in LMD-IPSL (Fig. 9g). The level of agreement among RCSMs (not shown) is low over France, Northern Africa, Gulf of Lions, over the Ionian Sea and in the Levantine basin, around the Cyprus area and over the Middle-East. Figure 10a and Fig. 10b show the changes in the area of the origin and the seasonality of cyclones crossing the region, respectively. Most RCSMs project a decrease of the number of cyclogenesis processes in the Mediterranean A general agreement on the sign of the changes in the minimum lifetime MSLP is observed (Fig. 10c). All RCSMs project an overall decrease in the intensity of cyclones moving through the region. Moreover, the signal observed is statistically significant in 5 out of 7 RCSMs (ROM, LMD-CNRM, LMD-IPSL, LMD-MPI and COSMOMED). The models show a relatively low agreement in the case of changes in maximum life-time ADR (Fig. 10d). On the other hand, in the case of mean radius the level of agreement is relatively higher (Fig. 10e). For example, RCSMs laterally forced by CNRM-CM5 (CNRM-RCSM4 and LMD-CNRM) project a shift towards significantly lower values of maximum ADR and higher values of the radius in the range [1300; 1500] km. COSMOMED and LMD-IPSL simulates a significant increase in the frequency of systems in the interval [− 6; − 3] hPa/6 h and [− 9; − 6] hPa/6 h respectively, and a decrease of the size of the systems crossing the region. Finally, EBU-POM simulates a significant increase in frequency of the systems in the interval [− 9; − 6] and To summarize, except for the decrease in the number of cyclone tracks over most of the Western part of Mediterranean region and a general tendency of the systems crossing the region to become shallower, there is no general agreement among RCSMs about the future changes in the cyclone activity in the Eastern basin, genesis areas, seasonality and in some key variables such as adjusted deepening rate. For the mean size five out of seven RCSMs project a statistically significant decrease. Figure 11 shows the distribution of changes in annual cyclone-related hourly-precipitation. Here, we consider the extended winter season ONDJFM for 2075-2099 and the changes (in %) with respect to the same season in 1981-2005. In the multi-model-mean a robust climate change signal is observed only in limited areas of the Mediterranean region (Fig. 11a). It is possible to observe an increase over part of the Alps, the central part of Italy and the Adriatic Sea (which coincides with the decrease of cyclone track density shown in Fig. 9a), and the Balkan area. Moreover, a decrease is observed over limited areas of the southern Iberian Peninsula, northern Africa and Levantine basin. Also, in each single RCSM the variations observed are hardly statistically significant and mostly limited to some small areas. LMD-IPSL (Fig. 11g) projects a decrease over most of the eastern part of the Mediterranean Sea (which also coincides with the decrease of cyclone track density observed in Fig. 9g). On the other hand, ROM (Fig. 11b), LMD-CNRM (Fig. 11e), LMD-MPI (Fig. 11f) and COSMOMED (Fig. 11h) project an overall increase over Italy (spatially coherent with the decrease of cyclone track density in these models shown in Fig. 9 b, e, f, g, h). CNRM-RCSM4 (Fig. 11d) and LMD-CNRM (Fig. 11e) simulates another increase over part of the Ionian Sea, which corresponds spatially to a previously observed increase of density of tracks (Fig. 9d, e). Moreover, CNRM-RCSM4 (Fig. 11d) also shows a decrease of cyclone-related precipitation, which is coherent with a decrease/increase of density track over the Tyrrhenian/ Levantine basin (Fig. 9d).

Cyclone activity in the Mediterranean region at the end of twenty-first century
RCSMs agree in showing a decrease in the wind speed over part of Western France and around the Pyrenees areas, part of the Ionian Sea, Gulf of Gabes area, Northern Africa and around Cyprus (Fig. 12a). An increase of the wind speed is observed in small areas around Sardinia, Tyrrhenian Sea and Eastern Europe (Fig. 12a). LMD-IPSL (Fig. 12g) projects a significant decrease of cyclonerelated wind speed over most of the eastern Mediterranean Sea (which also coincides with the decrease of both cyclone track density and precipitation observed in Fig. 9g and Fig. 11g). On the other hand, ROM, LMD-CNRM and LMD-MPI show a statistically significant increase of the wind speed on the Tyrrehian and part of Italy, which is spatially coherent with the variation of the precipitation previously discussed (Fig. 12b, e, f, respectively).

Discussion with respect to GCM and RCSM characteristics
In this section we discuss the sources of uncertainty in the characteristics (Sect. 3.1) and in the sign of future changes (Sect. 3.2) of the cyclone activity in the Mediterranean region.
In a recent study, Flaounas et al. (2021a) showed that the dynamics and intensity of the Mediterranean cyclones is influenced by both baroclinic forcing and diabatic processes. In fact, the more a cyclone is baroclinically driven, the less diabatic processes play a role to its development and intensity and vice versa. As a result, physical parameterization might play an important role in cyclone intensity differences among the models only for a fraction of the systems. As a token of this fraction, Galanaki et al. (2016) showed that only about a third of intense cyclones may produce deep convection in their center. For the rest of the cyclones where baroclinic forcing is important, the dynamics and intensity of the cyclone will depend on the driving GCM and more precisely on the dynamical characteristics of the intruding troughs that trigger cyclogenesis.
This explains the differences in the density of tracks between LMD-MPI (Fig. 2g) and LMD-IPSL (Fig. 2h) or the strong negative bias observed in the cyclone frequency (− 14%) in December in LMD-IPSL (Fig. 4b).
The IPSL GCM is, in fact, characterized by anticyclonic biases, which results in an overestimation of the blocking in the Northern Atlantic and reduced cyclone activity in the Mediterranean region (Tuel et al. 2020).
Some additional factors can explain the spread among RCSMs and biases observed in Sect. 3.1. One of these is represented by the regional model physics. In fact, the physics of the model and the resulting SST biases and surface temperature patterns (Darmaraki et al. 2019;Soto-Navarro et al. 2020) influences the air-sea heat exchanges and therefore partially the cyclogenesis processes together with the vertical stability of the air column and thus the convection and related precipitation. For example, ROM is characterized by an overall cold bias with respect to satellite SST of approximately 1.1 °C (Darmaraki et al. 2019) and this explains the underestimation of the precipitation observed over most of the Mediterranean Sea (Fig. 8b). The other factors are represented by a group of RCSM internal features such as the ratio between the size of the relaxation zone at the boundary and the lateral coupling frequency (Li et al. 2021), the RCSM internal variability (Sanchez-Gomez  the observed density of tracks and related precipitation. The so-called Regional Climate Model (RCM) internal variability could represent an additional mechanism driving the eastward increase of the spread among RCSMs observed in Fig. 2. In fact Sanchez-Gomez and  have shown that in the Mediterranean region, using the same boundary conditions and RCM numerical setting but slightly different initial conditions, the spread affecting the cyclone activity among different runs done with the same RCM is higher in the Eastern part with respect to the Western part of the region. This means that the eastward increase of spread among RCSMs density of tracks, mainly around the Cyprus area, could be driven by the internal variability of each member of the ensemble, further amplifying the differences eventually related to the factors discussed above. Finally, the different size of the domains could provide a clue of why for example ROM and COSMOMED, whose domains are the largest (Fig. 1), are characterized by higher percentages of Atlantic cyclones compared, for example, to LMD which has the smallest domain ( Fig. 1) and the highest percentage of cyclones formed in the Mediterranean region (Sect. 3.1).
The importance of the driving GCM can also explain the differences between LMD-MPI ( Fig. 11f and Fig. 12f) and LMD-IPSL ( Fig. 11g and Fig. 12g) in the projected changes for the precipitation and wind speed. The former projects a decrease of the density of tracks and an increase of hourly precipitation and wind speed over Central Italy, the latter a decrease of the cyclone track density and associated precipitation and wind speed over the Eastern basin. Further evidence of the importance of the driving GCM is represented by the comparison of the projections of CNRM-RCSM4 and LMD-CNRM. Both tend to show at the end of the twenty-first century an increase in the density of tracks over the Ionian and Levantine basin (Fig. 9d-e) together with an increase in the precipitation over the Ionian Sea ( Fig. 11d-e) and a decrease(increase) of the wind speed over the northern(southern) Ionian (Fig. 12d-e). On the other hand, the two simulations differ in the terms of climate change signal sign and significance in the Levantine basin, with CNRM-RCSM4 showing a decrease in the precipitation and wind speed over the area (Fig. 11d and Fig. 12d) and LMD-CNRM projecting an increase in the precipitation over Cyprus and an increase of the wind speed along southern Anatolian coastlines ( Fig. 11e and Fig. 12e).This further supports the importance of internal physics of each RCSM in modulating the climate change signal imposed by the driving GCM.

Summary and conclusions
Cyclones are an important component of the atmospheric circulation in the Mediterranean region, influencing the energy balance and water budget of the regional climate system. They also largely determine the timing and intensity of extreme events in the area (Nissen et al. 2010;Reale and Lionello 2013;Flaounas et al. 2016Flaounas et al. , 2017Flaounas et al. , 2021bLionello et al. 2019Lionello et al. , 2020Hochman et al. 2021).
Here, we use a set of 7 state-of-the-art fully-coupled Regional Climate System Models (RCSMs) simulations from the Med-CORDEX initiative (Ruti et al. 2016;Darmaraki et al. 2019) to analyze cyclone activity and associated precipitation and wind fields over the Mediterranean region for present day  and future projection (2075-2099; RCP8.5).
The RCSMs used here are able to capture the main characteristics of cyclone activity in the Mediterranean region in present-day climate conditions and related precipitation and wind patterns, despite some quantitative differences such as an overestimation (underestimation) of number of tracks over Northern Africa and Cyprus area (Gulf of Genoa and Iberian Peninsula), of the precipitation over the mountain chains (over part of the Adriatic Sea, the Balkan area and Aegean Sea) and of the wind flowing through the Otranto strait (over the Aegean Sea). However, no dramatic biases with respect to ERA5 reanalysis have been found. We have therefore considered these stateof-art RCSMs as plausible representations of cyclone activity in the Mediterranean region in present and future scenarios.
Our RCSMs ensemble projects for the end of the twentyfirst century a robust (according to the definition introduced in Sect. 3.2) decrease of the number and intensity of cyclone crossing the central part of Italy, Tyrrhenian Sea, part of the Anatolian Peninsula, Balkan area and part of Northern Africa and an overall weakening of the systems crossing the region. We also found a robust increase in the cyclone-related precipitation and wind intensity in the central part of the Mediterranean region. Conversely, we observed a change of opposite sign in both variables in the south eastern part of the region. Both signals are spatially coincident with the decrease of the number of tracks observed before. Our ensemble projects an overall decrease in total accumulated precipitation in winter over most of the Mediterranean region (not shown here). This means that, using the decomposition technique adopted in Zappa et al. 2015, the decrease of the accumulated precipitation in winter over most of the Mediterranean will be driven by the decrease in the number of cyclones crossing the area (Fig. 9) and will be only partially compensated by the increase in the intensity of the rainy events associated with each cyclone (Fig. 11). In the Eastern part of the region the drier conditions observed in winter will be driven by both decrease in the number of cyclones and intensity of each rainy event. These findings are in agreement with the findings of Zappa et al. (2015), which used 17 CMIP5 models to show that the northern part of the Mediterranean region (40 u N-45° N) is projected to have a decrease of the number of Mediterranean cyclones compensated by an increase in the amount of precipitation generated by each system and an overall decrease of cyclone-related precipitation in the Levantine area, which further amplifies the decrease of precipitation due to the reduction in the number of cyclones crossing that region (Hochman et al. 2020). Two of them (CMCC-CM and IPSL-CM5A-MR) are used as driving GCMs in this study. Additionally, the storyline approach adopted in Zappa et al. (2017) further supports the cyclone-related precipitation response to climate change observed in Fig. 11 and the importance of the driving GCM in explaining the signal observed. According to the classification discussed in Zappa et al. (2017), MPI-ESM-LR, MPI-ESM-MR, IPSL-CM5-MR, CMCC-CM are associated with high tropical amplification and weak polar stratospheric vortex. Conversely CNRM-CM5 is associated with low tropical amplification and strong stratospheric vortex. In the first case the storyline involves an overall reduction in the precipitation over the Eastern Mediterranean, as is shown Fig. 11b, c, f, g, h, linked to a persistent anticyclonic pattern insisting over the area. In the second case the storyline shows a large reduction of the precipitation over the Iberian peninsula as also shown in this case in Fig. 11d, e. It is relevant to note that, unfortunately, only two out of the four Zappa's storylines have been explored in Med-CORDEX using RCSMs, leading to unexplored possible futures.
Differences among the RCSMs emerge in the sign and magnitude of the changes concerning the cyclone's area of origin, their seasonality, deepening rate and associated precipitation and wind speed. The differences observed have been assessed with respect to different plausible sources, such as the driving GCM, the internal variability and physics of each RCSM and the domain size of the atmospheric module. The magnitude of the projected changes (also with RCSMs driven by the same GCM) varies considerably across models, showing the importance of the model physics representation in modulating the projected changes driven by GCM, especially in winter (e.g. Akperov et al. 2019).
Our results indicate that there are still important uncertainties in RCSM-based projections of changes in cyclone characteristics over the Mediterranean region. This uncertainty needs to be fully characterized through the use of large model ensembles sharing a consistent simulation protocol able to cover the whole climate change spread. Our results are based on the largest available, but still relatively small set of RCSMs. Another possible approach to better characterize the uncertainty in the future evolution of the Mediterranean cyclone activity would be a thorough and process-relevant sub-sampling of the whole a-priori ensemble to determine self-consistent, region-relevant and plausible future pathways of the regional climate system for a given phenomenon following the aforementioned storyline approach (Zappa and Shepherd 2017;Shepherd 2019).
Future work will include a larger set of simulations currently being produced under the Med-CORDEX framework to better characterize the uncertainty, but also to better understand the driving mechanisms of the projected changes in cyclone characteristics.