Scaling precipitation extremes with temperature in the Mediterranean: past climate assessment and projection in anthropogenic scenarios

In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.


Introduction
The complex geography of the Mediterranean region, which features a nearly enclosed sea with high sea surface temperature during summer and fall surrounded by mountains ( Fig. 1), plays a crucial role in steering airflow producing heavy precipitation and catastrophic floods (e.g. Delrieu et al. 2005;Jonkman 2005;Ducrocq et al. 2008;Papagiannaki et al. 2013). In the last 50 years significant increasing trends in daily torrential rains was observed in some western Mediterranean regions, while in the east relatively high interannual variability prevents any trend from being significant (Alpert et al. 2002).
Precipitation extremes and flash floodings are among the most devastating natural hazards in terms of mortality (e.g. Jonkman 2005). Even if flash floods are usually small-scale events, their suddenness and violence account for the high proportion of human losses. According to Jonkman (2005), the European and African continents display the highest mortality rates due to floods or flash floods in the world. Precipitation intensity is projected to increase in most regions under warmer climates, and the increase in precipitation extremes will be larger than that in the mean precipitation (Meehl et al. 2007), with large impacts on society (Oki and Kanae 2006;Pall et al. 2011;Min et al. 2011). In the Mediterranean region, observations (e.g. Quereda Sala et al. 2000;Xoplaki et al. 2003;Shohami et al. 2011) as well as the majority of twenty-first century scenarios show a warming and drying of the Mediterranean. The scenarios show a decrease in average precipitation with a peak signal in summer with either atmosphere-ocean Global Climate Models (GCM) (Giorgi and Bi 2005) or atmosphere Regional Climate Models (Gibelin and Déqué 2003;Déqué et al. 2005;Gao et al. 2006;Ulbrich et al. 2006) and coupled atmosphere-ocean Regional Climate Models (RCM) (Somot et al. 2008). There is, however, no consensus on the evolution of the frequency and intensity of the extreme events over the Mediterranean, even though a decrease of average precipitation and an increase in precipitation variability is expected (Pal et al. 2004;Giorgi 2006;Fig. 1 Target region of the HyMeX/MED-CORDEX regional climate simulations with colors indicating the topography and white circles the locations of the weather stations used in this study 1 3 Lavaysse et al. 2012;Önol et al. 2014) and an increased probability of occurrence of events conducive to floods has been suggested (Gao et al. 2006).
The Clausius-Clapeyron (CC) relation has been used as a predictor of changes in extreme precipitation with global warming. Indeed, the increase in the saturation water vapor pressure associated with warming is described by the CC relation: where e s is the saturation water vapor pressure, T is the absolute atmospheric temperature in Kelvin, L v is the latent heat of vaporization and R v is the gas constant. The CC relation considerably influences the changes in the extreme precipitation intensity under warmer climates (Betts and Harshvardhan 1987;Trenberth et al. 2003;Held and Soden 2006;Pall et al. 2007;Kharin et al. 2007;O'Gorman and Muller 2010;Muller et al. 2011;Romps 2011;Westra et al. 2014). Considering this, observed relations between extreme precipitation and temperature have received considerable attention. In the French Mediterranean region, Drobinski et al. (2016) showed that the scaling of precipitation extreme displays a hook-like shape with a nearly CC-scaling at low temperatures and a negative slope at high temperatures, which they attribute to the arid environment in summer. Deviations from the CC scaling have been observed or simulated for several reasons such as precipitation accumulation duration (Utsumi et al. 2011;Drobinski et al. 2016), use of surface air temperature as a proxy of the tropospheric temperature (Hardwick Jones et al. 2010;Drobinski et al. 2016), change in precipitation type Berg and Haerter 2013), changes in the atmospheric dynamics (O'Gorman and Schneider 2009;Sugiyama et al. 2010;Singleton and Toumi 2013;Muller 2013) which may be related to seasonal variability Drobinski et al. 2016). Hence, the applicability of the CC scaling on the temperature-precipitation extremes relationship is still uncertain. Indeed, there is a need to distinguish two different concepts. On the one hand, there is the long term conception that changes in precipitation extremes in a warming climate scale with the CC relation, such as discussed in Allen and Ingram (2002) and referred as "projected" scaling in the following. On the other hand, there are studies, among which those referenced above, which link precipitation extremes to the daily mean temperature in present-day observations. This relation is referred as "observed" scaling in the following. It is not trivial that these two are connected and that deviations from the thermodynamic CC behavior have similar origins.
This study focuses on the scaling of precipitation extremes with temperature throughout the Mediterranean region, based on quality controlled measurements of (1) ∂e s ∂T = L v e s R v T 2 surface temperature and precipitation at various weather stations around the Mediterranean, fine-scale gridded observational datasets and simulations performed within the context of the HyMeX ) and MED-CORDEX programs (Ruti et al. 2015) for present day and future climate. Our main objective is to investigate the temporal and spatial variability of the scaling and more precisely, to address the following questions: • What is the spatial variability of the "observed" scaling of precipitation extremes with temperature around the Mediterranean basin characterized by coastal areas, plains and high mountain ranges, and under the influence of arid and mid-latitude climates? • What is the uncertainty associated with the simulated "observed" and "projected" scalings of precipitation extremes? • What are the possible skill changes at high resolution relative to low resolution and the possible effects of ocean-atmosphere feedbacks on the simulated "observed" and "projected" scaling of precipitation extremes? • How much and why do the "observed" scaling in present-day and "projected" scaling in anthropogenic scenario differ?
Section 2 details our data (observations and simulations) and analysis methods. Section 3 investigates the ability of the regional climate models to simulate the dependence of precipitation extremes on temperature, its sensitivity to precipitation accumulation duration, model resolution and ocean-atmosphere feedbacks, and its spatial variability within the Mediterranean region. Section 4 analyzes the evolution of the scaling of precipitation extremes in the context of climate change, while Sect. 5 discusses some open research questions needing further investigation.

Methodology
Several methodologies have been used to characterize the scaling of precipitation extremes with temperature. Here we follow the method proposed by Hardwick Jones et al. (2010) and used in Drobinski et al. (2016), and apply it to time series of temperature and precipitation from three different sources: (1) in-situ measurements from surface weather stations across the Mediterranean basin, (2) E-OBS gridded data products and (3) regional climate simulations. At each station for a given precipitation duration (daily or subdaily), the maximum recorded precipitation rate on each rainy period is paired with the daily mean air 1 3 temperature at 2 m above the ground level. Subdaily intensities are also binned using the daily mean 2 m temperature instead of subdaily temperatures, since we are interested in a proxy representing the temperature of the air mass and for consistency with previous studies. For temperature in the E-OBS dataset and the regional climate simulations, we use the output of the grid point closest to the station location. The precipitation-temperature pairs are placed in bins according to temperature, with an equal number of samples in each bin (100 samples) and therefore varying temperature ranges for each bin. A sensitivity analysis was conducted with increasing numbers of samples up to 200 without any significant differences. The value of 100 samples was finally used as the best trade-off between the temperature-range sampling and the accuracy of the precipitation extreme estimate. The mean temperature of the events in each bin is used as the representative temperature for that bin. Within each bin of 100 samples, precipitation intensities are ranked to determine the 99th percentile. An exponential regression is applied to the precipitation values for each percentile (by fitting a least-squared linear regression to the logarithm of precipitation depth) to derive the slope of the temperature-precipitation extremes relationship. When this slope is 0.0596, it is equivalent to a Clausius-Clapeyron-like scaling of 5.96 % • C −1 at 25 °C.

Observations: E-OBS dataset
The temperature and precipitation series used in this study come from the E-OBS gridded dataset (Haylock et al. 2008), which was obtained through an interpolation procedure described in Haylock et al. (2008) from the ECA&D recorded series (Klok and Klein Tank 2009). We use here the version 10 of this product (released date: April 2014). The dataset covers Europe from 1950 to 2013 at a daily time-step with a spatial resolution of 0.25°. E-OBS has been widely used for regional climate model evaluation (e.g. Lorenz and Jacob 2010;Nikulin et al. 2011;Samuelsson et al. 2011;Torma et al. 2011;Flaounas et al. 2013;Stéfanon et al. 2014) and has been comprehensively evaluated in the literature (Hofstra et al. 2009(Hofstra et al. , 2010Kysel and Plavcov 2010;Flaounas et al. 2012;Herrera et al. 2012). In the E-OBS version 10 used here substantial improvements were realized, especially by increasing the station density (around 2000 for v1 and v2 vs. more than 4200 stations for temperature and more than 7300 for precipitation in v10). This should lead to some improvements in the drawbacks outlined above, and the E-OBS dataset is the only one that allows an RCM evaluation over the whole Mediterranean region.

Observations: in-situ measurements
Since the E-OBS product only provides daily values of temperature and precipitation we use in-situ measurements at subdaily time sampling for a better insight of the scaling of precipitation extreme with temperature (Utsumi et al. 2011;Drobinski et al. 2016). The in-situ data used in this work were collected at stations in Croatia, France, Israel, Italy, Spain and Greece. The locations of the stations are shown in Fig. 1 and additional details on the data are given in Table 1.

Simulation set-up and models
The list of regional climate models (RCMs) and details about the settings are given in Table 2. The simulations cover the ERA-Interim period 1989-2008 as initially recommended in MED-CORDEX project (Ruti et al. 2015). Some simulations were extended afterwards to cover the updated ERA-Interim period 1979-2008 and some reach 2014 in order to cover the HyMeX observation periods . All models use the 6-hourly ERA-Interim meteorological re-analysis as forcing data at the boundaries. For LMDZ, which is a global model with regional zoom capability, temperature, wind speed and Table 1 Origin of the in-situ data, available time sampling, measurement period used in this work, data provider and references All data provider are regional or national weather services except the Ebro Observatory   Ricard and Royer (1993) and Smith (1990) Noilhan and Planton (1989) and Noilhan and Mahfouf (1996) Ricard and Royer (1993) GLC2000 CNRM ALADIN52/RCM/ MED11 31 Morcrette (1990) Bougeault (1985) Ricard and Royer (1993) and Smith (1990) Noilhan and Planton (1989) and Noilhan and Mahfouf (1996) Ricard and Royer specific humidity are nudged towards the ERA-Interim fields outside the MED-CORDEX domain. It must be noted that the mesh is not regular within the zoom region and the resolution varies between 50 and 30 km. All other limited area models are forced at the boundaries using 3-dimensional re-analyses of wind, humidity, temperature or potential temperature and geopotential height. For CCLM, cloud ice and liquid water are additionally prescribed at the domain boundaries. The IPSL WRF simulations use nudging at all scales within the domain for temperature, wind and humidity above the planetary boundary layer (Salameh et al. 2010;Omrani et al. 2013Omrani et al. , 2015. The other models did not use nudging inside the inner domain. For the downscaling of ERA-Interim reanalyses, 14 simulations were produced from eight models (ALADIN V5.2, Colin et al. 2010;Herrmann et al. 2011;CCLM, Rockel et al. 2008;RegCM V4.3, Giorgi et al. 2012; WRF V3.1.1, Skamarock et al. 2008;Drobinski et al. 2012 for the WRF-NEMO-MED12 coupled system; LMDZ V4, Hourdin et al. 2006 coupled with NEMO-MED8). One of the novel aspects of MED-CORDEX is the use of high spatial resolution and atmosphere/ocean coupling, and one of the objectives of this article is to assess the sensitivity of the models' skill to ocean coupling and resolution. Five simulations were carried out with resolution higher than 50 km (0.11° or 0.18°), and nine with low resolution (0.44°) on the MED-CORDEX domain. Among these two sets, four simulations were made with the same model and similar setups (ALADIN, RegCM4, WRF) at high and low resolutions. Among the 14 simulations, 11 were run in atmosphere-only mode, i.e. at their lower boundary, sea surface temperatures were specified from the ERA-Interim reanalysis. Three runs include an interactive ocean model (using NEMO-MED8 or NEMO-MED12, Lebeaupin Brossier et al. 2011, 2012a, with sea surface temperatures simulated by the coupled ocean model. In the following the simulations are referred to by the name of the modeling group followed by the name and version of the model and the resolution (11 for 0.11°, 18 for 0.18° or 44 for 0.44°). The 3-hourly model outputs are used in the following.
The downscaling of the CMIP5 simulations are performed similarly to what is done for the ERA-Interim reanalyses downscaling. The global/regional model combinations are listed in Table 3. The high end RCP8.5 scenario is considered here, where RCP stands for Representative Concentration Pathways (IPCC 2013). The mean global warming projections for RCP8.5 by the CMIP5 models is 2.0 °C for 2046-2065, with a likely range of 1.4-2.6 °C. It increases up to 3.7 °C with a likely range of 2.6-4.8 °C at the 2081-2100 horizon. Figure 2 shows the 99th percentile of precipitation as a function of daily mean temperature at the surface weather stations in Croatia, France, Israel, Italy, Spain and Greece and from E-OBS at the nearest grid points for all available time sampling. Despite a dry bias of E-OBS with respect to the in-situ measurements, there is a remarkable agreement in the temperature-precipitation extremes relationship, and sub-CC scaling appears to be a common feature across the datasets for high temperatures. The first remarkable feature is that whatever the time sampling, all curves display a hook-shape as shown for France by Drobinski et al. (2016). We can define a temperature break below which the scaling is close to CC, and above which the scaling is sub-CC with a quasi systematic negative slope. In other words, the temperature break corresponds to the temperature at which a sharp variation of the temperature-precipitation extreme slope occurs. The temperature break is around 20 °C in the west (Spain, France) and decreases when going eastwards with a temperature break around 7 °C in Greece. Despite some variability within each country, Fig. 3 shows that the general behaviour going from one country to the other is remarkably reproducible at all available stations. Figure 3 displays the temperature-precipitation extremes relationship at all selected stations in France, Italy, Israel and Greece (for Croatia and Spain, data from only one station is available). In Israel, the slope of the temperature-precipitation relationship is negative for all temperatures (e.g. Fig. 2c). None of the other Mediterranean stations exhibit such a trend for sub-daily timescales. This negative trend is found at several Israeli stations (Fig. 3c). A similar result The crosses indicate the global model used to drive the regional model

Regional models
Global models was found in South-America tropical regions for daily precipitation extremes although sub-daily timescales were not investigated (Maeda et al. 2012). As shown by Utsumi et al. (2011) in Japan and Drobinski et al. (2016) in southern France, the slope of the temperature-precipitation extremes relationship increases as time sampling decreases which can be well explained by the decrease in the duration of the precipitation events at high temperatures, and not by the decrease in the precipitation intensity of individual storms. In Spain, the good quality of the dataset allows the evaluation of the time sampling limit below which the averaging effect is no longer affecting the slope. Since there is no significant differences between the slopes at high temperature for 30, 20 and 10 min time sampling, one can conclude that the time sampling limit is between 30 min and 1 h. In Spain, the sub-CC scaling at high temperatures is therefore real and not the result of time averaging. The fact that the hook-shape is also visible in Croatia for time sampling of 1 h suggests that such behaviour may be robust in this area, even though generalization must be considered with caution. In Spain, super-CC is observed for the smallest timescales between 5 and 10 °C.This is consistent with the observations of Lenderink and Meijgaard (2008) at De Bilt in Netherland, albeit super-CC scaling was found at much warmer temperatures, as also found in France for 3-hourly precipitation (Fig. 2b). This super-CC scaling was interpreted as resulting from a shift in rain type by Haerter and Berg (2009), and Berg and Haerter (2013), though the first reference is a hypothesis, whereas the second only discusses a small number of stations. Berg et al. (2013) and Loriaux et al. (2013) found with smaller statistical effect a super-CC scaling with convective precipitation only, the latter showing an important role of the vertical velocity for super-CC behavior. and from the HyMeX/MED-CORDEX simulations at the nearest grid point. Analysis of the temperature-precipitation extreme curves in Figs. 2 and 4 shows that the temperature break decreases from west to east. It is around 16-20 °C in France, Croatia and Spain, it goes down to around 7-10 °C in Greece and it is no longer detectable in Israel where the slope is always negative. The models display a very similar scaling, with an accurate representation of the slopes and temperature breaks, despite a systematic underestimation of the precipitation extremes which is expected when comparing local precipitation with area averaged precipitation over a grid of about 2500 km 2 (Chen and Knutson 2008). The GUF-CCLM4-8-18-MED44 simulation is the closest to the observed amount of precipitation extreme (as measured by the 99th percentile) while the CNRM-ALADIN52-MED44 simulation always produces much lower extremes than observed. Besides the averaging effect, sub-CC scaling may also be attributed, to aridity which may limit the rainfall amount as suggested by Hardwick Jones et al. (2010) for Central Australia stations. In the northwestern Mediterranean, Drobinski et al. (2016) also attribute the sub-CC scaling to the arid environment with three consequences. First, the 2-m temperature is no longer a proxy of the temperature of the condensed parcel, as precipitation is produced at much higher altitudes due to higher condensation levels. The surface temperature is thus much higher than the temperature of the condensed parcel. Second, the relatively dry air surrounding the precipitating system reduces the precipitation efficiency which thus reduces the slope of the temperature-precipitation extremes relationship. Third, the entrainment of relatively dry unsaturated environmental air dilutes the rising saturated air parcels and causes evaporative cooling, which in turn reduces the buoyancy of the convective parcel and the vertical transport of humidity. These processes reduce the slope of the temperature-precipitation extremes relationship and cause the reduction in temperature break as conditions become drier from Spain to Israel along the coast. Note that the Italian stations do not exhibit a clear temperature break. These stations are located in the mountains at high altitudes (Fig. 1), and as will be shown in Fig. 13, at those altitudes aridity conditions do not occur often and do not limit the rainfall amounts.

Simulation evaluation and analysis
The scalings of precipitation extremes in Italy and Israel display a negative slope above about 10 °C, which is fairly well reproduced by all simulations except CMCC-CCLM4-8-19-MED44 and GUF-CCLM4-8-18-MED44. However, contrary to Israel, the slope of the temperature-precipitation extreme curves in Italy is sensitive to time sampling as CC-scaling is found for hourly precipitation extremes (Fig. 2d). Figure 5 shows the spatial pattern of the slope of the temperature-precipitation extremes relationship for temperatures lower than the temperature break or in the absence of temperature break. The slope at low temperatures is close to the CC-scaling (∼7 % • C −1 ) for most of the northwestern regions. Negative slopes are observed in the south and east of the Mediterranean basin. All models reproduce fairly accurately the limit between the regions with nearly CC-scaling and those with sub-CC scaling, which indicates  Figure 6 displays the spatial pattern of the slope of the temperature-precipitation extremes relationship for temperatures higher than the temperature break. No slope is found in most parts of the eastern and southern Mediterranean, where no temperature break could be estimated (Fig. 5). Negative slope, i.e. sub-CC-scaling, is found everywhere in the Mediterranean in the E-OBS dataset and in the simulations. Sub-CC scaling is found above a temperature mostly ranging between 27 and 3 °C from west to east (Fig. 4), to be compared to a value of ∼20-25 °C in some previous studies (e.g. Utsumi et al. 2011).
As for the daily time sampling, models systematically underestimate the 3-hourly precipitation intensity as shown in Fig. 7. For the Israeli station, the observed temperatureprecipitation extreme curve is too variable so that no conclusion can be drawn. All models still produce a negative slope at all temperatures, except the CNRM-ALADIN52-MED44 simulation which shows a hook shape curve with a temperature break around 10 °C. CMCC-CCLM4-8-19-MED44 and GUF-CCLM4-8-18-MED44 simulations display a slightly positive slope. At all other stations except the French ones, all models produce a flattening of the negative slope at high temperatures. In Italy, the slope retrieved from the in-situ data is positive for temperatures higher than the temperature break, whereas it is still negative in all models. In France, the scaling of precipitation extremes with temperature is better reproduced by all models. Our analysis shows that there is no effect of the horizontal resolution (0.44-0.11) on the scaling of precipitation extremes for daily precipitation (not shown). When 3-hourly precipitation data are used, the sensitivity to the horizontal resolution becomes more significant at some stations, even though the effect is small. Indeed, the ratio of 3-hourly to daily precipitation extremes, which measures the duration of the simulated events (e.g. Utsumi et al.  , appears to be sensitive to the resolution at high temperatures. The ratio equals 1 when the daily averaged precipitation reduced over a period of three hours (i.e. divided by 8) is equal to the 3-hourly accumulated precipitation. This essentially means that it rains throughout the day. Conversely, when the ratio equals 8, it only rains during 3 h during the day. At low temperatures, the ratio is statistically the same between the simulations performed at high and low resolutions, with a value of about 2.5 in good agreement with the observations (not shown). This ratio is quite low, suggesting long-lasting simulated extreme precipitation events probably dominated by non-convective precipitation. The ratio increases with temperature for both simulations and weather station measurements, but with a lower rate for the low-resolution experiments. At the highest temperatures, the difference of ratio between the high and low resolution experiments is generally significant except for UCLM-PROMES. The ratios from the high-resolution simulations are lower, but closer to observations (≃6), than low-resolution models. However, the spread across models increases with temperature (not shown). At high temperatures, the simulated precipitation extremes are mainly of convective nature and their duration is highly sensitive to the convection scheme. The rate at which the ratio increases with temperature is larger for the simulations performed at high than low resolution, and therefore, going from daily (Fig. 4) to 3-hourly (Fig. 7) precipitation extremes, the high resolution runs straighten the temperature-precipitation extreme curve at high temperatures more than the low resolution ones. Regarding ocean-atmosphere feedbacks, no sensitivity to the atmosphere-ocean coupling was found on the scaling of precipitation extremes with temperature (not shown). This is not surprising, since atmosphere-ocean feedbacks affect predominantly the location of the precipitating event rather than the intensity Berthou et al. 2014Berthou et al. , 2015.

Comparison between historical and ERA-Interim simulations
Before investigating climate change effects, we investigate the ability of historical runs to reproduce ERA-Interim results. Figure 8 shows the daily precipitation extremes as a function of the daily mean temperature at the 6 weather stations and in the HyMeX/MED-CORDEX simulations (CNRM-ALADIN52-MED44, ITU-RegCM4-3-BATS-MED44, ITU-RegCM4-3-CLM-MED44, ICTP-RegCM4-3-MED44, LMD-LMDZ4NEMO8-MED44) forced by ERA-Interim reanalysis and by the historical CMIP5 GCM simulations (Table 3). For all stations and precipitation accumulation duration, there is a good agreement between the ERA-Interim and CMIP5 downscaled simulations over the 1979-2005 period (or 1989-2005 for some models; Table 2). Temperature breaks can slightly differ, but these differences are generally not significant and are caused by the different ranges of the temperature bins used to compute the temperature-precipitation extreme relationship. The largest differences between the scalings of precipitation extremes computed from the ERA-Interim and CMIP5 downscaled simulations are found for CNRM-ALADIN52-MED44, where the slope at high temperature is generally steeper in the ERA-Interim downscaled simulations. One can note a general underestimation of the ICTP-RegCM4 precipitation intensity at low temperatures which leads to an overestimation of the slope at these temperatures. Figure 9 shows the daily precipitation extremes as a function of the daily mean temperature at the 6 weather stations and for the GCM-driven HyMeX/MED-CORDEX simulations for the historical period  and the future period 2070-2100. The most striking changes are (1) an increase of the temperature break by about 4-5 °C and therefore a broader temperature range over which the slope is close to the CC-scaling; (2) a change in the value of the sub-CC slope at high temperatures, with steeper or flatter slopes depending on the model. The latter effect is even more pronounced when analyzing the 3-hourly model outputs (not shown). Such evolution of the scaling of precipitation extremes with temperature suggests that a simple CC-scaling may not be used as a predictor.

Projection in the future climate
From a simple CC-scaling argument, if all other aspects of convection remain unchanged with warming (in particular precipitation efficiency or vertical updraft velocities, Drobinski et al. 2016), one would expect a simple shift of the precipitation extremes curve parallel to the CC slope towards higher temperatures (Fig. 10). Indeed, if all other aspects of convection remain constant, the simple thermodynamic expectation is that precipitation extremes should scale with the humidity convergence at cloud base, i.e. P extreme ∝ RHQ sat (Betts and Harshvardhan 1987;O'Gorman and Schneider 2009), where P extreme denotes an extreme rainfall rate, Q sat is the low-tropospheric saturation specific humidity, and RH is the low-tropospheric relative humidity. If RH stays constant in time with warming, and if the temperature increase is uniform in space, then the temperature increase only impacts precipitation extremes through the contribution from Q sat . The latter follows CCscaling (Fig. 10).
The temperature break shifts towards higher temperatures by about 4-5 °C over the whole Mediterranean region and Europe (including Eastern Europe) for all models without much spread across models (not shown). The value of the temperature break shift corresponds to the average temperature change in the region (3-4 °C), which is one of the highest warming rates in the world both in observations (e.g. Quereda Sala et al. 2000;Xoplaki et al. 2003;Solomon et al. 2007) and in the scenario projections (e.g. Giorgi 2006). However, a simple temperature translation by 4-5 °C of the temperature-precipitation extreme curve obtained from the historical runs does not allow the reconstruction of the curve obtained from the downscaled RCP8.5 simulations, since such translation would not reproduce the change of the slope at high temperatures (Fig. 10). Figure 11 shows the warming between historical and future daily mean temperatures for rainy days at the nearest grid point to the weather stations from the HyMeX/ MED-CORDEX simulations at 0.44° resolution. The average shift of the temperature break corresponds approximately to the average temperature difference simulated between 2070-2100 and 1979-2005. In France for instance, ICTP-RegCM4-3-MED44 shifts the temperature break from 20 to 25 °C (Fig. 9d), which also corresponds to the mean temperature change at 25 °C (Fig. 11). Although models produce very similar results for the scaling of precipitation extremes with temperature, the warming relationships in Fig. 11 are characterized by a large spread across models. For instance, in Croatia, all models exhibit temperature warming curves with a "U-shape", but in Greece, the temperature change decreases with temperature for LMD-LMDZ4NEMO8-MED44, it has a "U-shape" for ITU-RegCM4CLM-3-MED44, ITU-RegCM4CLM-3-MED44 and ICTP-RegCM4-3-MED44 and a more irregular shape for CNRM-ALADIN52-MED44.
The fact that the warming T ′ − T is not the same at all temperatures T, where T and T ′ denote current and future climate temperatures respectively, can change the thermodynamic expectation described in Fig. 10. The temperature-precipitation extreme curves obtained from the historical simulations adjusted assuming the CC-scaling with respect to the temperature change are shown in Fig. 12 and are referred as the "adjusted historical" scaling of precipitation extremes. In other words, these "adjusted historical" curves are obtained by shifting the current climate curve along the CC-slope by a shifting amplitude corresponding to the temperature change. They are compared to the curve from the historical simulations and the projected curve from the RCP8.5 simulations. The results are shown at the 6 surface weather stations using the simulations at 0.44° resolution. The results show that the adjusted curves match fairly well the projected curves suggesting that despite a negative slope of the temperature-precipitation extremes relationship at high temperatures, the CC law seems to apply in a warming regional climate. In detail, the best recovery of the projected curve using the CC law are obtained in Spain, France, Italy and Croatia for all models.
The deviation of the "adjusted historical" scaling of precipitation extremes from the scaling in the future climate is larger in Greece, especially at high temperatures for ITU-RegCM4BATS-3-MED44, ITU-RegCM4CLM-3-MED44 and ICTP-RegCM4-3-MED44. In Israel, the "adjusted historical" scaling of precipitation extremes fits well the future scaling of precipitation extremes for ITU-RegCM4BATS-3-MED44, ITU-RegCM4CLM-3-MED44 and LMD-LMDZ4NEMO8-MED44, the latter fitting best the historical scaling of precipitation extremes at cold temperatures. Changes in relative humidity can also impact the thermodynamic expectation described in Fig. 10. Figure 13 shows the relative humidity for the present and future climate from ITU-RegCM4CLM-3-MED44 simulation as a function of the daily mean temperature in the future climate. One single model is shown, but similar results were log Q sat T T T log RH log P extreme T cold climate T' warm climate Fig. 10 Schematic representation of the expected change of the scaling of precipitation extremes with warming (right panel). T and T ′ denote current and future climate temperatures respectively. If all other aspects of convection remain constant, the simple thermodynamic expectation is that precipitation extremes should scale with the humidity convergence at cloud base, i.e. P extreme ∝ RHQ sat . The relative humidity RH is assumed to remain constant at each location (middle panel) but simply shifts to warmer temperatures as each location warms. The temperature increase T ′ − T is assumed to be the same at all locations. In that case precipitation extremes increase with warming only through an increase in the low-tropospheric saturation specific humidity Q sat , which follows CC-scaling (left panel)  Figure 13 also shows the "adjusted historical" relative humidity under the assumption of constant specific humidity. Note that all the relative humidity estimates, even the ones from historical runs, are plotted as a function of the future climate temperature to ease comparison. In addition to the 6 Mediterranean stations, one additional region in Central Europe was added to highlight the differences between the Mediterranean and the continental European climate, in principle less exposed to moisture advection from the oceans. Figure 13 clearly shows that for all Mediterranean stations, relative humidity does not change significantly in the future, even though in Israel, at high temperatures, the future relative humidity is slightly lower than in the present. Figure 13 also shows that the relative humidity is significantly different from the predicted relative humidity assuming constant specific humidity, which means that moisture sources from the oceans compensate the effect of increased temperature on relative humidity. Conversely, more continental regions can display a significant departure of the temperature-relative humidity relationship at high temperatures, as also shown in the work of Lenderink and Attema (2015). Figure 13 shows an example for a large region in Central and Eastern Europe (average over several hundreds of kilometers around Moscow). A decrease of relative humidity can be seen in this region far from large water bodies. At high temperature, the relative humidity in the future tends towards its adjusted value assuming no additional moisture sources (i.e. temperature change effect only). This suggests that at high temperatures in this region, the temperature increase associated with global climate change is not counter-balanced by an increase of humidity.
Why is such behavior not seen at the 6 Mediterranean weather stations is a key question. The Mediterranean Sea is actually a considerable moisture source for the region, and in a warmer future, evaporation is expected to increase with increasing temperature (Somot et al. , 2008Moratiel et al. 2011). Satellite observations show that the integrated water content over the oceans follows the CCscaling with respect to the sea surface temperature (e.g. Stephens 1990) due to energetic constraint (Schneider et al. 2010). Therefore the Mediterranean Sea is the primary source of humidity contributing to maintain the relative humidity constant over parts of the Euro-Mediterranean region. This is consistent with the study of Lebeaupin-Brossier et al. (2015) who showed that the response of precipitation to variations in sea surface temperature affects the continent over a horizontal scale limited by the surrounding mountain ranges.
Despite this finding, we note that the CC correction assuming constant relative humidity does not produce a full match between the adjusted and projected scalings of precipitation extremes. The differences may have various origins which are not accounted for by our simple scaling, such as model uncertainty (e.g. convection and microphysics parameterizations; Singh and O'Gorman 2012), change in precipitation duration (and thus impact on averaging effect), change in precipitation efficiency or vertical moisture transport in a more arid environment (Paltridge et al. 2009;Zhang 2009;Drobinski et al. 2016). The ratio of 3-hourly to daily precipitation extremes in the low-resolution simulations increases with warming (not shown), suggesting shorter extreme precipitation events and thus a steeper negative slope at high temperature for daily precipitation accumulation. However, the significance of these changes is not large.

Conclusion
Expected changes in future extreme precipitation remain a key uncertainty associated with anthropogenic climate change. Extreme precipitation has been proposed to scale with the precipitable water content in the atmosphere. Assuming constant relative humidity, this implies an increase of precipitation extremes at a rate of 7 % • C −1 as indicated by the Clausius-Clapeyron relationship. In this study, the scaling between precipitation extremes and temperature is investigated for the Mediterranean region by assessing against observations over the 1979-2008 period, ensembles of regional climate simulations performed within the frame of the HyMeX and MED-CORDEX programs. Despite quantitative differences in precipitation between the various models, the change in precipitation extremes with respect to temperature is robust and consistent across the models.
This study highlights the spatial variability of the scaling of precipitation extremes with temperature across the Mediterranean basin. In general, this scaling exhibits a hook shape with negative slope at high temperatures and a slope following CC-scaling at low temperatures. The temperature break, i.e. the temperature at which the slope sharply changes, ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. It becomes undetectable further Fig. 13 Relative humidity in the future (blue) and present (red) climate as a function of daily mean future temperature at nearest ITU-RegCM4CLM-3-MED44 grid point from the weather stations from the HyMeX/MED-CORDEX simulations at 0.44° resolution driven by the CMIP5 historical simulations and RCP8.5 scenario performed by Hadgem2-ES. The green curve is the projected relative humidity in the absence of additional moisture source east and south of the Mediterranean, meaning that the slope is always negative in the arid regions of the Mediterranean. The scaling of precipitation extremes is insensitive to ocean-atmosphere coupling and it only slightly depends on resolution at high temperatures for short precipitation accumulation times (where precipitation is largely parameterized). This is due to the precipitation event duration differing at low and high resolutions at short time scales. However, such resolution effect is small and does not affect significantly the scaling of precipitation extremes with temperature.
Our study also investigates the evolution of this scaling in future climate scenario experiments over the 2070-2100 period. The temperature break shifts to higher temperatures by a value which is on average close to the mean regional temperature change due to global warming. The slope of the simulated temperature-precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break. At high temperatures, the curve can be flatter or steeper in the various models compared to present climate, i.e. some models exhibit more intense precipitation extremes in the future. The origin of the slope evolution is in agreement with the Clausius-Clapeyron law suggesting that the relative humidity remains constant in warming regional climate, with simply a distribution shifted to warmer temperature values. This may appear counter-intuitive since the Mediterranean is warming and drying at high rates, however this trend is not only simulated by global and regional climate models but also found in observations. The Mediterranean Sea is a large source of moisture for the region, and as atmospheric water vapor content increases with increasing temperature (following the Clausius-Clapeyron law) the Mediterranean Sea can maintain the relative humidity constant in a warming regional climate at least over a region under its influence. This region is limited by the numerous mountain ridges which surround the Mediterranean Sea, and in fact over continental Europe, at high temperatures, the relative humidity decreases significantly. This has implications for the scaling of precipitation extremes there (weaker extremes expected there at the highest temperatures due to water limitation).
The hook shape of the Mediterranean scaling of precipitation extremes with temperature is robust, and its present shape can be used to approximately predict the evolution of the precipitation extremes in a changing climate. Despite the hook shape deviating significantly from the CC scaling (∼7 % °C −1 ), its evolution appears to follow quite well the CC law. A natural follow-up of this study is (1) a detailed water budget evaluation to identify the humidity pathways that maintain constant relative humidity in parts of the Mediterranean region and (2) a thorough analysis of the climate change induced atmospheric circulation changes that cause precipitation extremes.