Assessing Future Changes of Climate Extreme Events in the CORDEX-MENA Region Using Regional Climate Model ALADIN-Climate

This study investigates future changes of temperature, precipitation, and associated extreme events in the MENA region using Regional Climate Model ALADIN-Climate over the CORDEX-MENA domain. Model capabilities to reproduce key observed regional climate features are first assessed, including heat waves, drought and high precipitation extremes. Projected changes indicate the intensification of heat waves number, duration and magnitude, and contrasted precipitation changes. A drying is projected in the north-west and moistening in the north-east along the Mediterranean side of the region. Projected regional warming is found at the rate of about 0.2 °C/decade to 0.5 °C/decade over land depending on the scenario. Drought is expected to increase in the northern half of the region independently from the index used, but with a higher rate in the case of the index accounting for both the effect of precipitation and temperature changes. ALADIN-Climate results corroborate previous studies projecting the MENA region to host global hot spots for drought in the late twenty-first century.


Introduction
The Middle East and North Africa (MENA) region extends from the Maghreb in Northwest Africa to the Arabian Peninsula. The MENA region has been proned to climate hazards (Lelieveld et al. 2016;Waha et al. 2017). These include the devastating Syrian Arab Republic's droughts in 1998-2000and 2007(Kelley et al. 2015Flohr et al. 2017), the recurrent droughts registered during 1980-1994 in Morocco (Mansouri 2004;Driouech 2010), the prolonged drought of Greater Horn of Africa (Shongwe et al. 2011;Mwangi et al. 2014;Rowell et al. 2015;Gebremeskel et al. 2019), drought registered in Jordan and Lebanon in 2014 (Dai 2011;Barlow et al. 2016), as well as important changes of climate characteristics are evidenced by several studies focusing on climate trends at regional and local scales (Nasrallah et al. 2004;Zhang et al. 2005;Driouech 2010; AlSarmi and Washington 2011; Dahech and Beltrando 2012;Almazroui et al. 2014;Donat et al. 2014;Filahi et al. 2015;Almazroui 2019a).
The high climate vulnerability of this region is widely recognized and several negative impacts are expected to make the strengthening of resilience and adaptation (Sowers et al. 2011;United Nations Economic and Social Commission for Western Asia (ESCWA) et al. 2017;Waha et al. 2017). Significant and consistent warming trends have been depicted since the middle of the twentieth century across the region with increased frequencies of warm days and warm nights, fewer cold days and cold nights and higher extreme temperature values (Tramblay et al. 2013;Donat et al. 2014;Islam et al. 2015;Merlone et al. 2019;Almazroui 2020). Although less regionally consistent, downward trends of total precipitation amounts have been found in different sub-regions of MENA jointly with an increase in drought (Tanarhte et al. 2012;Driouech et al. 2013;Tramblay et al. 2013;Donat et al. 2014).
Recent scientific and technical advances have explored the issue of the habitability of the Earth under increasing anthropic pressure using General Circulation Models (GCMs), Regional Climate Models (RCMs) and other climate downscaling techniques (IPCC 2018;Déqué et al. 2017;Field et al. 2012). GCMs are the principal tools for investigating the large-scale response of the climate system to various forcings, for making climate predictions on seasonal to decadal time scales and for making projections of future climate over the coming century and beyond (Flato et al. 2013). For example, Almazroui et al. (2020) demonstrated the changes in temperature and precipitation using latest Coupled Model Intercomparison Project (CMIP6) simulations for Africa under the three Shared Socioeconomic Pathways (SSPs). They found that future warning is not uniform over Africa and varies regionally, which demand the necessity of downscaling GCM dataset using RCM. However, GCMs are not suitable for studying extreme phenomena over a region like MENA where climate fluctuations are subject to enhanced spatial heterogeneity, due to their coarse spatial resolution. RCMs and other climate downscaling techniques have been used for ''zooming in" the results of the GCMs on a particular area of the globe. Also, RCMs provide fine-scale resolution and improved physical formulation for mesoscale processes, and the ability to account for local-scale physiographic drivers of regional climate features (Giorgi and Lionello 2008;Rummukainen 2010;Giorgi and Gutowski 2015;Gutowski et al. 2016;Lennard et al. 2018). But, it rapidly appeared that climate change responses from RCM depend both on the emission scenarios, large-scale driving fields from the global model, and the RCM physical formulation Laprise et al. 2013;Llopart et al. 2014;Dosio et al. 2015;Almazroui et al. 2016). The Coordinated Regional Climate Downscaling Experiment (CORDEX) of the World Climate Research Programme provides a novel and standardized approach to assess RCM climate change responses over predefined geographic domains worldwide (Jacob et al. 2014;Giorgi et al. 2014;Llopart et al. 2014;Roudier et al. 2015;Gutowski et al. 2016;Zou and Zhou 2016;Lennard et al. 2018;Mba et al. 2018). The CORDEX program aims at enhancing international cooperation and improving the understanding of relevant regional/local climate phenomena, their variability and changes. Considering the need for a common climate downscaling framework to evaluate future changes and impacts in the MENA region, the CORDEX "Region 13: Middle East North Africa (MENA)" regional domain (CORDEX-MENA) has been established recently under CORDEX. This domain, illustrated in Fig. 1, was initially defined within the framework of the Regional Initiative for the assessment of Climate Change impacts on water resources and socioeconomic vulnerability in the Arab Region (RICCAR; https ://www.escwa .un.org/RICCA R) jointly established by the Arab League and the United Nations Economic and Social Commission for Western Asia (ESCWA). However, much of the CORDEX research in the MENA region has focused on evaluating the RCMs to reproduce historical climate features, with models forced either by large-scale atmospheric conditions from reanalyzes or GCMs. Furthermore, and despite the numerous similarities the region presents, very few climate assessment studies have been performed for the whole region.
The MENA climate is predominantly semi-arid to arid with generally dry and hot summers and mild winters (Donat et al. 2014), and is controlled by oceanic forcing related to heat, humidity and winds coming from surrounding oceans (Almazroui 2019b). The annual total precipitation exhibits significant spatial and temporal variability, and most land areas receive less than 300 mm of annual precipitation (Waha et al. 2017). Over the Arabian Peninsula, spatial variability of precipitation is high, due to the Mediterranean climate's, Indian monsoon and orographic variations (Bucchignani et al. 2015). Precipitation exceeds 900 mm in temperate zones, such as the Atlas Mountains in Morocco, the . Annual mean temperature range from less than 10 °C to more than 30 °C with the highest values located over the desert lands, characterized by very high daily fluctuations and hot summers (reaching 40 °C to more than 50 °C). The temperature is moderated by elevation and by proximity to the sea. The southern part of the region, including southern Mauritania and southern Sudan, benefits from wet tropical climate and receives about 100-600 mm/year with a mean annual temperature of about 18 °C (United Nations Economic and Social Commission for Western Asia (ESCWA) et al. 2017).
In this study, we assess for the first time the climate change responses of the ALADIN-Climate ("Aire Limitée Adaptation dynamique Développement InterNational") RCM over the CORDEX-MENA domain, with a focus on changes of temperature and precipitation extremes. The main objective is to contribute to the improved understanding of future climate extremes in this region. The article is organized as follows. Section 1 describes the key model formulation features and assessment metrics. Section 2 provides a brief evaluation of the model performance in simulating climatological characteristics over the MENA region. Future changes are assessed in Sect. 3 and a conclusion is presented in Sect. 4.

Model and Experimental Design
This study uses version 5 of the limited-area Regional climate Model ALADIN-Climate Radu et al. 2008;Colin et al. 2010;Farda et al. 2010;Herrmann et al. 2011;Nabat et al. 2014;Ruti et al. 2016;Mallet et al. 2019;Spiridonov and Valcheva 2019) with a horizontal resolution of 50 km. The model physical formulation is inherited from ARPEGE-Climat and CNRM-CM5 CGCMs used in the Coupled Model Intercomparison campaigns (Déqué 2010;Voldoire et al. 2013). The dynamical core is bi-spectral with semi-implicit, semi-lagrangian advection scheme, and includes an 11-point wide bi-periodization zone in addition to the more classical 8-point relaxation zone. ALADIN-Climate comprises the Fouquart and Morcrette radiation scheme based on the European Centre for Medium Range Forecasts (ECMWF) model, incorporating effects of greenhouse gases (CO 2 , CH 4 , N 2 O and CFC) and direct effects of aerosols, as well as the first indirect effect of sulphate aerosols. The vertical discretization has 31 layers, with 27 layers from the surface up to the tropopause and a time step of 15 min.
Version 5 of ALADIN-Climate is used in this work, as was the case within the framework of MedCORDEX (Torma et al. 2015;Ruti et al. 2016;Somot et al. 2018;Vaittinada et al. 2016;Fantini et al. 2018). Three simulations have been carried out by the RCM with a horizontal resolution of 50 km (~ 0.44°) driven by lateral boundary conditions from the CMIP5 Global Climate CNRM-CM5 (Salas-Mélia et al. 2005;Voldoire et al. 2013) over the domain shown in Fig. 1. This domain includes the active part (yellow area) where the RCM results are considered usable and the transition zone between the GCM driving boundaries and the RCM (red and green zones). The area considered for analysis (blue rectangle) coincides with the CORDEX-MENA domain. The historical simulation covers the period 1971-2005 and uses observed aerosol field and Greenhouse gas (GHG) concentrations (Tramblay et al. 2013;Farda et al. 2010). Two other simulations were run over the period 2006-2100 based on the Representative Concentration Pathways RCP4.5 and RCP8.5, respectively.
Note that Version 4 of ALADIN-Climate was used for the European ENSEMBLES project where it was inter-compared with several European RCMs at 50-km and 25-km resolution, and was also used for the assessment of future climate change over Morocco (Driouech 2010;Driouech et al. 2010).

Climate Indices
Model capabilities and future climate changes over the region are assessed using a set of sector-specific climate indices derived from daily temperature and precipitation data, and designated to assess simple and complex climate risks (https ://climp act-sci.org/indic es/) in socio-economic sectors such as water resources, agriculture and health. These indices form the standardized set recommended by the World Meteorological Organization (WMO) Commission of Climatology Expert Team on Sector-Specific Climate Indices (ET-SCI), 1 parts of which is based on those defined by the joint Expert Team on Climate Change Detection and Indices (ETCCDI) 2 (Peterson 2005;Zhang et al. 2011;Alexander and Herold 2016). Furthermore, to the annual mean temperature and the total precipitation from wet days, we selected indices relative to extreme temperature and precipitation including heat waves, high precipitation events and drought. These indices were determined by the ET-SCI in representing the acclimatization to heat (EHI(accl.)) and the climatological significance (EHI(sig.). Compared to other heat wave definitions (i.e., based on TX or TN thresholds), the EHF-based heat wave is, from a climatological standpoint, the most appealing as it considers TX and TN within the same index (Alexander and Herold 2016). In addition, the conditions leading up to a given day (i.e., the previous 2 days) are also considered. The given day is also compared to an extreme threshold (via the 90th percentile) to determine its climatological anomaly, and to the mean of the preceding month to determine the anomaly against recent conditions. The EHF is calculated as (Alexander and Herold 2016): where TMi represents the average daily temperature for day i (i.e., TMi-32 is the average daily temperature for day i-32) and TM90i is the 90th percentile of TM over all calendar day i within the user-specified base period, using a 15-day running window. TM is calculated via TM = (TX + TN)∕2.
The EHF is defined from the above two definitions: EHF = EHI(sig.) × max(1,EHI(accl.)). o Heat Wave Duration (HWD, in days) defined as the length in days of the longest heatwave identified by HWN. corresponding to the annual total precipitation from heavy rain days (days > 95th percentile). o Amount of rainfall from extremely wet days (R99p, in mm) corresponding to the annual total precipitation from very heavy rain days (days > 99th percentile).
The drought index SPEI and heat wave indices based on the Excess Heat Factor have not been evaluated before over the MENA region in terms of future changes at high resolution neither in terms of observed trends.
For each index, the statistical significance tests of changes were undertaken, at a confidence level of 95%, on each gridpoint. The Mann-Whitney U test was used, as precipitation has a skewed distribution in many arid regions and the test makes no assumptions about the distribution (James and Washington 2013).
Published in partnership with CECCR at King Abdulaziz University

Model Evaluation
The climate model ability to reproduce the MENA climatic features is assessed by comparing simulated outputs against observed datasets for the historical/reference period. The evaluation of ALADIN-Climate over the MENA domain (described in Sect. 1.1) is conducted using the fifth generation of ECMWF atmospheric reanalysis of the global climate (ERA5) data, recently developed (Hersbach and Dee 2016;Hersbach et al. 2018Hersbach et al. , 2019. ERA5 data cover the Earth on a 30-km fine grid and replaces the ERA-Interim reanalysis which stopped being produced in August 2019. ERA5 improves over ERA-Interim and combines vast amounts of historical observations (including remote sensed and in situ data) into global estimates using advanced modeling and data assimilation systems. For the purpose of this study, ERA5 3 data daily maximum, minimum temperature and precipitation are interpolated onto the 50-km grid of ALADIN-Climate RCM. The climate indices used in the analysis are listed in Sect. 2.2 and calculated for the historical period . Note that a more appropriate evaluation of the RCM requires performing simulation driven by reanalysis (representing "perfect" boundary conditions), free from biases carried over from the GCM. Uncertainties linked to gridded observational datasets may also be significant due to, for example, the density of the meteorological network (Contractor et al. 2015;Prein and Gobiet 2017;Zittis 2018). The absolute mean bias of temperature anomalies averaged over land grid point values is small and does not exceed 0.17 °C. The CNRM5 model exhibits a similar pattern of cold bias indicating that this feature is inherited from the Fig. 2 Climatological mean temperature and hot extreme events simulated by ALADIN-Climate RCM (left panels) and derived from ERA5 reanalysis (middle panels) for the historical period , and their differences (right panels). Annual mean temperature (a, b, c), number of consecutive hot days and nights (TX3TN3) (d, e, f) and number of heat waves (HWN) (g, h, i)

Mean Temperature and Hot Extreme Events
Published in partnership with CECCR at King Abdulaziz University driving GCM (Voldoire et al. 2013;Dosio et al. 2015). Temperature biases have been attributed to systematic errors in non-dynamical processes, while the main features of the large-scale dynamics are well captured in such experiments driven by observed sea-surface temperatures (Cattiaux et al. 2013). Note that CNRM5 is an improved physical formulation of the CNRM3 GCM which also exhibited large cold bias (Voldoire et al. 2013). Cold biases linked to the driving GCMs are also found in simulations undertaken within RICCAR initiative over CORDEX-MENA domain using three different RCMs. (United Nations Economic and Social Commission for Western Asia et al. 2017). RCM outputs have been compared with ERA-Interim Reanalysis over the 25-year period 1980-2004. Other studies have also highlighted cold biases of climate models over the region (Syed et al. 2019;Dike et al. 2015;Ozturk et al. 2018;Zittis et al. 2014). Temperature biases are even higher in some parts of the region although some model simulations were integrated under 'quasi perfect boundary' conditions (driven by reanalysis). Cold biases have been also found in both minimum and maximum temperatures (Zittis et al. 2014), the setup for cloud microphysics seems to have the strongest impact on such biases.
In terms of hot extremes, results show that events with three consecutive hot days and nights (TX3TN3) occur in ERA5 data with a frequency mostly not exceeding 1 time per each 2 years. This means that, hopefully, most of the extreme hot events involving both extreme high maximum and minimum temperatures do not persist generally more than 1 or 2 days in the region. Simulating such kind of events and their evolutions is, however, very useful in the context of global warming due to the damages they can engender and the need for adaptation. For example, the capability for relief and adaptation can have a profound impact on excess mortality and morbidity (Meehl and Tebaldi 2004;Hansen et al. 2008, Ahmadalipour andMoradkhani 2018;Ahmadalipour et al. 2019). These "relatively" low values of TX3TN3 are well reproduced over the entire region exhibiting the model ability to simulate the temporal persistence of hot extreme events (Fig. 2d, e, f). This is further demonstrated by the simulated magnitude of the extended summer (May-September) heat waves (HWM) (not shown), except over the Atlantic coastal zones of Morocco where the HWM reaches its maximum probably due to a high climatological significance (EHI(sig.)) of heat waves (see Sect. 1.2 for calculation of HW). These coastal zones register generally high temperature values during heat waves compared to its 90th percentile which gives a higher EHF (sig.) than in inland zones.
ALADIN-Climate also captures well the numbers of heat waves (HWN) with an absolute mean bias not exceeding 0.4 in approximately all the domain (Fig. 2g, h, i). The length of the longest heat wave (HWD) modeled is as well satisfactory compared with reanalysis despite some overestimation in parts of Saudi Arabia and Egypt and underestimations along the African monsoon regions south of the domain (not shown). However, the biases do not exceed 2 days in nearly all the domain.

Mean Precipitation and High Precipitation Events
ALADIN-Climate reproduces the spatial distribution of annual mean precipitation across the domain and captures the main features including the lowest precipitation amounts in the dry African Sahara and Arabian Desert (precipitation less than 1 mm day −1 ), the rainfall enhancement over complex terrains (such as Zagros Mountains and the Yemen mountains), as well as the precipitation band in the south associated with the Inter-Tropical Convergence Zone (ITCZ) over Africa (Fig. 4a, b, d, g, e, h). Moreover, the humid zone of North Africa associated with the North Atlantic-Mediterranean cyclones is well represented. The simulation exhibits, however, a tendency to slightly overestimate the annual total precipitation amounts over the mountainous areas and along the belt under the influence of African monsoon south of the domain (where a model cold bias is noted). A dry bias is apparent in the most humid part of tropical West Africa located southwest of the domain (Fig. 4c, Zhao et al. (2014) found that CNRM-CM5 performs best for precipitation over North Africa among the models used.
The spatial patterns of total precipitation due to very wet days (R95p) and extremely wet days (R99p) are also well reproduced by ALADIN-Climate. The largest intensities are found around orographic regions and in the tropical wet area resulting from deep convection occurrences (Diedhiou et al. 1998;Sylla et al. 2013;). Consistent with observations, the minima occur in the Sahara and Arabian Desert. The model overestimates the amplitudes of precipitation, although wet biases remain relatively low, over the mountainous areas and along the belt under the influence of the West African Monsoon (WAM). Overestimations of total high precipitation amounts are also found in the case of CORDEX multi-model ensembles compared to GPCP datasets (Huffman et al. 2009). Comparison of CORDEX outputs against TRMM data (Huffman et al. 2007) indicates better estimations, highlighting the potential effect of the observational dataset (Gibba et al. 2018). Several differences have also been found in the magnitudes of the two observational datasets (e.g., GPCP and TRMM), where the TRMM generating larger intensities than the GPCP.    and their respective differences (right panels). Mean annual precipitation (a, b, c), the total precipitation issued from very wet days (d, e, f) and extremely wet days (g, h, i)

3
Published in partnership with CECCR at King Abdulaziz University Focusing on MENA region, Zittis (2018) showed significant discrepancies in various aspects of the five gridded observational datasets used, including the monthly time distribution, interannual variability, mean climatology and precipitation extremes. Significant differences have also been identified between several observational and reanalysis products in the case of exceptional extreme rainfall events for the eastern Mediterranean (Zittis et al. 2017). The CMIP5 GCMs and CORDEX RCMs ensembles, evaluated over Africa, overestimated the frequency of extreme precipitation events and produced little high-intensity precipitations compared to observations (Dosio et al. 2015;Gibba et al. 2018). A previous version of ALADIN-Climate (Version 4) as well as the ARPEGE-Climat model (version 4) have been evaluated in the north-west of the region (Morocco) per comparison to meteorological stations data (Driouech et al. 2009;Driouech 2010). Overall, the model version used in this study shows improved capabilities in reproducing precipitation behaviours across MENA region for both annual total precipitation amounts and high precipitation extremes.

Drought
The model performance in simulating drought across the domain is assessed through the percentage of dry years and the drought indices SPI and SPEI, respectively. For each index, a year is considered as dry if the corresponding index value is below − 1 (McKee et al. 1993;Vicente-Serrano et al. 2010). Our results show that drought occurs in the region with a frequency from 5 to 20% with SPI and from 10 to 25% in the case of SPEI. The comparison, across the region, between the numbers of dry years detected by each of the indices in ERA5 data shows that the SPEI tend to Fig. 5 Percentage of dry years, as depicted through the drought indices SPI (a, b, e) and SPEI (c, d, f), simulated by ALADIN-Climate RCM (a, c) and derived from ERA5 reanalysis (b, d) for the historical period  and their respective differences (e, f) Published in partnership with CECCR at King Abdulaziz University detect more drought events than the precipitation-based index (Fig. 5). This is due to the taking into account by the SPEI of the effect of temperatures that affects the severity of droughts (Vicente-Serrano et al. 2010). The model reproduces well this behavior furthermore to its ability in simulating satisfactorily drought conditions at the annual scale across the domain, independent of the index used (Fig. 5).

Projected Changes
Projected changes over the CORDEX-MENA domain are assessed through contrasting ALADIN-climate results for the future period (2071-2100) against the present-day period  under each of the two RCP scenarios (RCP4.5 and RCP8.5).

Projected Changes of Mean Temperature
and Hot Extremes Figure 6 shows the regional scale climate responses to the increase of greenhouse gases, for annual mean temperature and hot extremes. Results indicate statistically significant changes (at 95% level) consistent with a generalized warming pattern over the domain, with an increase in mean temperature values around 2 °C for RCP4.5 (Fig. 6a) and higher for RCP 8.5 (Fig. 6b). The warming is more moderate along the coastal areas compared to inlands. A warming of at least 4 °C is given with RCP8.5 in North Africa (Eastern Morocco and most Algeria) and the Arabian Peninsula. The Sahel is projected to warm by around 2 °C (RCP4.5) to more than 3 °C warming (RCP8.5). The warming rate is to at least double between the scenario RCP4.5 and RCP8.5, and regional temperature trends over land are estimated to reach 0.2 °C/decade in RCP4.5 and 0.49 °C/decade in RCP8.5 for the period 2020-2100. The number of consecutive hot days and nights (TX3TN3) is projected to increase by about 2-6 for RCP4.5 and 6-10 for RCP8.5 (Fig. 6b). Compared to the historical period corresponding to one event every 5 years (if not less), these changes correspond to frequencies multiplied by at least 10. Changes are more homogeneous across the region and exhibit higher increase over the Arabian Peninsula and west of Iran (+ 4 to + 6 days) compared to the remaining parts (+ 2 to + 4 days).
There is also an intensification of heat waves in terms of number, duration as well as magnitude, generalized to the whole region but with spatial contrast depending on the index and the scenario (Fig. 6c-e). The number of summer (May-September) heat waves increases for RCP4.5 from 3 to about 10 with the highest values (6-10) in the Sahel. The increase is statistically significant at 95% level over the whole region for both scenarios but slightly lower with RCP8.5 over the Arabian Peninsula, Egypt and the Sahel.
The changes are statistically significant and more important with RCP8.5 than RCP4.5 for the duration and the magnitude of heat waves particularly over the Arabian Peninsula and North African countries. Some exception is found with the magnitude of heat waves under RCP4.5 where significant changes do not extend to north of the Arabian Peninsula, Egypt, and part of the Sahel region. The longest heat wave will be more persistent by up to 40 and 80 days for the two scenarios, respectively. The Arabian Peninsula, Egypt and the Horn of Africa reach + 100 days with RCP8.5, corresponding to the twothird of the summer period. The magnitude increases with RCP4.5 by 2 °C 2 -4 °C 2 over southern Morocco and less than 2 °C 2 elsewhere if not insignificant change. Higher emissions result, however, in higher increase extended to nearly the entire region. The heat wave magnitude will intensify the most over Morocco, Western Algeria, Libya, Egypt and Iran (4 °C 2 -5 °C 2 ). It is worth noting that some studies indicated a potential role of soil moisture-linked feedbacks in the increase of extreme temperatures and heat waves in the mid-latitude heat waves (Douville et al. 2016;Vogel et al. 2017).  Figure 7 presents the future changes of annual precipitation as issued from ALADIN-Climate. In contrast to the projected changes of temperature, the response to the GHGs emissions are not uniform across the MENA region. A drying is projected with RCP8.5 in the Mediterranean area, over the western part of the domain including Morocco and north-eastern Algeria, consistent with the previous studies (IPCC 2013). The reduction of total precipitation amounts varies from − 5 to − 20% and exceeds − 40% west of the Atlas Mountains. These changes are significant and consistent with the projected northward shifts of air moisture associated with a stronger North Atlantic Oscillation (NAO) (Waha et al. 2017) and a poleward displacement of the Hadley circulation descending branch (Collins et al. 2013).

Projected Changes of Mean precipitation and High Precipitation Events
On the other hand, precipitation increases in nearly all the remaining parts of the domain. The changes in the Sahara and Arabian Desert correspond to increases exceeding + 40% or even a doubling of the present-day (historical period) few quantities. The exhibited north-south dipole coincides with increases of about + 10% to + 40% over the Sahel and the horn of Africa. With the RCP 4.5, the changes are relatively less intense although they reach a decrease of 5-20% in North-west Africa and increases of about 40-100% southern Arabian Peninsula and Egypt. The northward shift of the ITCZ projected by several previous studies including those summarized in IPCC reports (IPCC 2013) is involved in the increase in precipitation over southern parts of the MENA region (Waha et al. 2017). Based on the Mann-Whitney U test, changes under RCP8.5 are statistically significant at 95% level in most of the region excepting north of the Arabian Desert, south of Algeria, and west of Libya. Significant changes are found for RCP4.5 over the Atlas Mountains, south of the Arabian Peninsula, east of the Sahel and over the Horn of Africa.
The future changes simulated by ALADIN-Climate RCM are in line with projections from previous studies using RCMs or CMIP5 models (Driouech 2010;Collins et al. 2013;Lee and Wang 2014;Gutowski et al. 2016;Waha et al. 2017) including the decrease in total precipitation amounts in the north-west of the region (Morocco) and the increase of precipitation, especially under the intense scenarios, over the Sahara and Southern Arabian Peninsula. Furthermore, using a set of 17 CMIP5 models that includes the CNRM-CM5 (the driving model of ALADIN), Zhao et al. (2014) show also compatible precipitation changes with those issued from ALADIN-Climate. However, it is worth noting that some RCMs project strong decrease in precipitation generalized to all parts of the domain (Ozturk et al. 2018).
Similar behaviors are depicted in the projected changes of very wet days (R95P) with a decrease in precipitation from high events over the band from Morocco and north-west Algeria and more extreme precipitation in the remaining part of the domain (Fig. 7c, d). The highest reductions are given over the northern half of Morocco and reach − 5% to − 40% depending on the scenario and location.
The changes in terms of extremely wet days show with RCP8.5 a large predominance of increases across the region, mostly ranging from + 20% to more than + 80% (Fig. 7e, f). The areas above 25° latitude exhibit contrasted changes with either positive or negative signs depending on the location according to RCP4.5. Compatibly with R95P, the very extremely wet days would be less frequent west of the Atlas Mountains with RCP8.5. In general, the future changes noted in mean precipitations across COR-DEX-MENA domain should be due to those relative to the high precipitation events at least partially. Although statistically not significant at 95% level, the results found here for the very high precipitation events are globally compatible with projected changes issued from other models (Cook et al. 2018). The increase of precipitation at the end of twenty-first century in the similar regions is also projected using the latest CMIP6 simulations ).

Projected Changes of Drought
According to the SPI, more frequent droughts are expected over North Africa, northern parts of the Arabian Peninsula and parts of the Fertile Crescent (Fig. 8a, b). Extend of the increasing area is, paradoxically, reduced with RCP8.5 and the region concerned by a reduction in the percentage of dry years become larger with the changes mostly statically significant. The later includes the southern half of the Arabian Peninsula, most of the Saharan Desert and the Horn of Africa. Both the positive and negative changes vary from 5 to 20% excepting over the northern half of Morocco where it reaches 40% with RCP8.5. Projected increase of drought over Morocco and northern Algeria is statistically significant with RCP8.5 The changes in drought detected using the SPEI (Fig. 8c, d) show more increase in terms of intensity and extent, consistent with the index sensitivity to changes in evaporative demand caused by temperature increase. Temperature rise affects, in fact, the severity of droughts through increased evaporative losses amplifying deficits in soil moisture and runoff, despite fewer precipitation changes. The area above 20° including North Africa and most of the Middle East would experience statistically significant changes corresponding to an increase of 20-40% in the percentage of dry years with the highest values exceeding 40% to a doubling of the present-day ratio. Consistent with the projected increase of mean and extreme precipitation events, less frequent dry years are projected over the remaining southern parts. The significant role of temperature rise has also been found in terms of increased aridity especially in the Mediterranean coastal region due to substantial increase in evapotranspiration (Waha et al. 2017). Land-surface feedbacks tend to increase the drying through contribution to the enhanced warming (Lorenz et al. 2016;Russo et al. 2019).
Our results are in line with various studies that consistently projected the North African countries (notably Morocco, Algeria, and Tunisia), as well as a part of the Middle East area, to become global hot spots for drought by the end of the twentyfirst century (Planton et al. 2008;Collins et al. 2013;Orlowsky and Seneviratne 2013;Sillmann et al. 2013;Prudhomme et al. 2014;Trenberth et al. 2014;Smirnov et al. 2016;Waha et al. 2017. The variation of future drought changes depending on the methodology has also been highlighted (Trenberth et al. 2014) with a larger agreement about the increase of extreme droughts.

Conclusion
This study reports the first application of the ALADIN-Climate version 5 Regional Climate Model over the CORDEX-MENA domain, with a focus on assessing the projected climate change responses for mean variables and extreme events under RCP4.5 and RCP8.5 scenarios. The general validation of ALADIN-Climate simulation for the historical period suggests that the model captures reasonably well several climate extreme indices including the excess heat factor related heatwaves and consecutive extreme hot days and nights (TX3TN3). ALADIN-Climate also reproduces the observed mean regional temperature variations and warming trend, though there are cold biases over mountainous regions, also apparent in the CNRM5 driving GCM. The RCM tends also to overestimate mean precipitation over complex terrains and along the belt under the influence of African Monsoon climate. Although the multi-model ensemble approach is always preferred to climate change response from a single RCM, which may partly reflect the scientific and numerical choices of this particular model, the analysis of ALADIN-Climate responses over the MENA region provides some important insights. Projected regional warming rate ranges from about 0.2 °C/decade to 0.5 °C/decade over land depending on the scenario, and there is an intensification of hot extreme events. Heat waves are projected to intensify in terms of number, duration and magnitude. Precipitation changes are spatially contrasted, with drying over the north-western part. An increase of drought over the northern half of the region is expected, with differences depending on the index and the emission scenario. The effect of warming on the severity of droughts has been detected through the use of SPEI, which is consistent with the finding of previous studies (Cook et al. 2018). Our results are also consistent with several studies that indicate an increase in extreme drought conditions around the Mediterranean, Northern Africa, and the Middle East (Waha et al. 2017;Trenberth et al. 2014).
Due to the wide uncertainty associated with inter-model variability across GCMs, we stress the importance of exploring fully the CORDEX-MENA multi-GCM/multi-RCM ensemble to assess the robustness of climate change responses in terms of climate extremes including drought and heat waves over the MENA region. The effects of increasing spatial resolution on the results should also be investigated. The potential impacts of projected changes in human and natural ecosystems need to be assessed in a rigorous way to inform risk reduction decisions.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.