Assessment of Future Water Demand and Supply under IPCC Climate Change and Socio ‐ Economic Scenarios, Using a Combination of Models in Ourika Watershed, High Atlas, Morocco

: Climate change will affect the water resources system, on global and regional levels. Over the past thirty years, the High Atlas Mountains in Morocco have experienced severe droughts, which causes a decrease in water supply that affects both agriculture and the urban water system. In this paper, we assess the impact of climate change and socio ‐ economic activities on water supply and demand in the Ourika watershed (High Atlas of Morocco), then we evaluate the efficiency and sustainability of regional adaptation strategies for water supply management. For this, we simulate and analyze the future water situation using the statistical downscaling model (SDSM) and the water assessment and planning tool (WEAP). After the model’s calibration and validation, the precipitation, minimum (Tmin) and maximum (Tmax) temperatures, water demand and unmet water demand were projected for 2100 horizon, using different climate change scenarios. The results revealed that the model’s performance, calibration and validation were found to be satisfactory. The analysis shows that the mean precipitation will decrease by 49.25% and 34.61% by 2100, under A2 and B2 emission scenarios of the Intergovernmental Panel on Climate Change (IPCC). The projected mean Tmax and Tmin will be warmer than the baseline period, with Tmax increasing by 4.2 °C (A2) and 3.6 °C (B2), and Tmin by 3.5 °C (A2) and 2.9 °C (B2) by 2100. The results also show that water demand and the unmet water demand will increase in all scenarios, the pressure on water resources will increase, leading to water scarcity. The results reveal that, under the influence of climate change, future unmet water demand is expected to reach 64 million cubic meters (MCM) by 2100. The results demonstrate that the assessments of the proposed adaptation strategies are effective, but not sufficient to ensure water sustainability for the Ourika watershed.


Introduction
Climate change is unequivocal, global and will have severe consequences, even under the most optimistic scenario [1]. Among the most important impacts of climate change are the rise in surface temperature and the reduction in precipitation in some regions [2] and declining water quantity and quality [3,4]. Indeed, water resources are the most vulnerable to climate change [5]. It is estimated that more than a third of the worldʹs population (∼2.4 billion people for 2000) live in countries under water stress and this will rise to two thirds by 2025 [6]. Climate change, population growth, and economic activities are leading to water scarcity, an unpleasant situation in which demand for water increases beyond water supply [7]. Previous studies have considered the Mediterranean region as a global ʺhot spotʺ in terms of climate change, where precipitation is likely to decrease with increasing temperatures [8][9][10]. They are among the countries most threatened by water stress, due to the extreme interannual, spatial, and temporal variability of water resources and the decrease in water inflows expected in the coming decades [11]. When problems related to water quantity and quality occur, it can have repercussions, not only on flora and fauna, but also on the economy, health and well-being of the human population.
Located in one of the driest regions of the world, Morocco stands out as one of the most vulnerable countries to climate change [12]. It has extremely low water resources, a low per capita supply of fresh water, and is considered one of the most stressed countries in terms of water resources, due to rapid population growth. The renewable freshwater per capita is below the United Nations ʺabsolute water scarcityʺ level. In this context of water scarcity, the gap between supply and demand will widen in the years ahead. Morocco is highly dependent on surface water resources, notably for its agricultural sector [9]. The share of agriculture in the national economy is considerable at 8 billion USD, which corresponds to 14.5% of gross domestic product (GDP) for 2014. Water scarcity is a limiting factor in local economic development. The first effects are already visible in the Ourika watershed, where water resources are limited to start with, and are decreasing in both quantity and quality under the combination of decreasing supply and a growing population. Decrease in the river streamflow, agriculture and increased tourism flux are putting even more pressure on the Ourika water resource. An analysis of climate data over the last decades shows an increase in average temperature of 0.76 °C (1990-2001 and 2002-2014), and a decrease in precipitation and the proportion of snow.
Therefore, there is an urgent need to quantify the possible climate change impacts on water resources in this region. Despite the importance of this issue to local population and regional agriculture, to date there are no studies assessing the impact of climate change on water supply and demand using the outputs of global climate models (GCMs).
The water resources of a region depend on both the rainfall amount and temperature. Therefore, for estimating the future impact of climate change, GCMs are used to simulate the climate system. Projections of future climate change are made using a hierarchy of climate models that simulate changes in a response to a set of anthropogenic forcing's scenarios [13]. The results of the GCMs are based on a large-scale grid resolution (>250 km) [14]. Because of this coarse spatial resolution, the output of GCMs needs downscaling, in order to be used to study the environmental and hydrological impacts of climate change at regional and local levels [15,16].
Hydrological models allow one to study the relationship between climate and water resources [17]. Scientists around the world have used various hydrological models for water resource assessments. For example, the soil and water assessment tool (SWAT) watershed model [18], the water resources management model (WRMM) [19], the JAMS J2000 framework [20],the spatial agro hydro salinity model (SAHYSMOD) [21], the modular simulator (ModSim) [22], the hydrologic landsurface processes of mass and energy transfer model (PROMET) [23], and the water assessment and planning tool (WEAP), which has been widely used in different basins of the world in recent decades [24].
To analyze the impact of climate change and different socio-economic scenarios on water supply and demand in the Ourika watershed, and to improve decision-making capacities in water resource management, the WEAP hydrological model was used. The WEAP model has been used in previous studies around the world [25,26], and in Africa [27][28][29]. The objective of this study is to explore the incidence of climate change and socio-economic scenarios on water supply and demand, in the Ourika watershed up to the year 2100, and to analyze the sustainability of the proposed adaptation measures to climate change. In this study, therefore, the statistical downscaling model (SDSM) is used to produce high resolution climate data, using the large-scale climate variables derived from GCMs. The statistically downscaled results are then used to drive the WEAP model. The impact of future climate change and socio-economic scenarios are evaluated by comparing various WEAP responses to different climate projection scenarios and baseline conditions. This study's results are necessary to address future problems threatening a region where 80% of its surface water is used for agriculture. This sector is particularly vulnerable to extensive anthropogenic activities, climate change and its consequent variability in climate and climate extremes, such as droughts and floods. The Ourika watershed is known to face a decreasing trend in rainfall, especially during the rainy season, the availability and quality of water resources, environmental degradation, and population growth.

Study Area
The Ourika watershed is located in the south of Marrakech (Morocco) in a semi-arid to subhumid climate. It is one of the major watersheds in the Tensift basin of the High Atlas mountain (the highest mountains in North Africa). Geographically, the watershed is located between longitudes 7°53' and 7°35' West, and latitudes 31°20 and 31°4' North, and has an area of 507 km 2 ( Figure 1). The size and the height of the mountain make the climate of this region variable [30]. The amount and intensity of precipitation are generally controlled by altitude. The elevation ranges between 974 m and 4001 m above seas level (Jebel n'Tarourt) and about 75% of the watershed area is located between 1600 and 3200 m.a.s.l. Above these high-altitude areas, winter precipitation mainly falls as snow. The average annual precipitation is 541 mm per year at the river outlet and can exceed 700 mm per year at the upstream area of the basin [31]. The average monthly temperatures range from 21.5 to 32 °C for maximums and from 4 to 5.7 °C for minimums. The average annual temperature is 17.6 °C at Aghbalou station for the period of (1970-2014), (Figure 2). The Ourika River originates in the southwestern part of the Atlas Mountains in the North and flows through several villages. Its water is used to meet the daily needs of domestic water for the local residents and tourists, and most importantly to irrigate croplands and orchards (e.g., Ourika Organic Aromatic Gardens, Paradise of saffron). The water fluxes of Ourika river have inter-annual variations with a succession of dry and wet years. Minimum monthly mean flows occur in August. However, in this month, thunderstorms can alter the Ourika river mean summer hydrological regime and cause devastating torrential floods [30]. The maximum flows are recorded in winter and spring respectively. The spring period has the highest flows, due to the high level of precipitation and snowmelt. The months from January to June account for 78% of annual flows ( Figure 3). Geomorphologically, the Ourika watershed is characterized by steep slopes and impermeable and compact rock formations, occupying more than 55% of the total area [30]. They consist mainly of crystalline and fissured soils, semi-permeable soils (Triassic red clays, schist, flysch, granite, and basalt). As for the permeable grounds (limestone, basalt, and alluvium), they represent less than 6%. In terms of land use, the forest dominates much of the catchment area, followed by agriculture and rocky outcrops [30].

Data Collection
To develop the SDSM-WEAP model, different data sources were used (Table 1), including meteorological stations (observed data), GCMs data from the National Centre for Environmental Prediction (NCEP), Hadley Centre Climate Model-version 3 (HadCM3) for A2 and B2 scenarios, and the digital elevation model (DEM) used for watershed delineation and MODIS land cover datasets. The streamflow measurements from 5 gauging stations and evapotranspiration data were acquired from the Regional Water Agency (Tensift Watershed Agency) and Al Haouz Regional Agricultural Development Office (Regional Office for Agricultural Development of Al Haouz ORMVAH). Supplementary data, including population, livestock, water use rates, water consumption, and irrigation area were acquired from different sources and are subsequently validated by socioeconomic surveys.

Climate Change
SDSM is an open source software that produces high-resolution climate change scenarios. It is a decision support tool [32] for assessing the impact of global climate change at local and regional scale. Two types of daily predictors data required for this study; the 26 predictors of the reanalyzed NCEP data for the period of 1961-2014; and the 26 predictors of HadCM3, for A2 and B2 scenarios, for the period of 1961-2100. These were used to generate future climate projections. A2 and B2 scenarios were used instead of 1. Intergovernmental Panel on Climate Change (IPCC)'s newest emission scenarios, because they are more representative of the trajectory of Moroccan society evolution and the changes associated with climate indicators [33]. The SDSM downscaling process is achieved by computing a statistical relationship between the predictors of NCEP variables and the predictands (local surface variables: precipitation, maximum temperature, and minimum temperatures) [34], using the monthly (SDSM-M) and annual (SDSM-A) sub-models. The NCEP atmospheric predictors and the daily observed local data were used for calibration and validation of the SDSM model.
Daily local climate data are used as predictand variables for downscaling experiments and precipitations and max and min temperatures were measured at the Aghbalou and Oukaimden weather stations. The simulated data were extracted for three periods of 30 years each between 2014 and 2100: the 20 s (2014-2040), 50 s (2041-2070), and 80 s (2071-2100). To ensure the reliability of the data, the SDSM allows quality control and identifies missing or outliers' values before model calibration [35]. Data quality control is followed by screening downscaling predictors [36], to select the most appropriate ones that determine the downscaled climate scenario. Many indicators can be used in this process. In this study, the partial correlation, correlation matrix, and P value were used. The identical approach was used in previous studies [37][38][39].
In general, the calibration of a model is very important; it permits the model to develop the best correspondence between the observed data and simulated data. The data used for calibration in this study covered 30 years (1970-2000), and the period from 2001 to 2014 is used to validate the model. The SDSM calibration calculates the parameters of the multiple regression equations for a set of probable predictors and each predictand [40]. The SDSM outputs and the observed data were compared using several parameters, the coefficient of determination (R 2 ), the explained variance (E%), root mean square error (RMSE), standard error (SE), mean (μ), standard deviation (б), relative error in mean (RE μ), and relative error in standard deviation (RE б) for temperature and precipitation, during the periods of calibration and validation. To improve the quality of outputs from climate models, the correction of biases has become an essential approach in downscaling studies [41]; bias correction is applied to remove or reduce biases in daily time series of downscaled data. In this study, the method developed by [42] was adopted for bias correction. Equations (1) and (2) are used to debias daily Tmax, Tmin, and precipitation data.
where deb are the corrected predictands for (2014-2100). SCEN represents uncorrected downscaled daily data by SDSM for (2014-2100), CONT represents downscaled outputs for the base line period . Obs represent observed values.

Hydrological Model
The WEAP is a user-friendly tool that takes an integrated approach to water resources planning and policy analysis. The model is based on the water balance concept to reproduce hydrological processes [24]. It can be applied to a single river basin or at the level of a complex basin system [28]. It provides an integrated assessment of hydrology, land use, hydrogeology, climate, water quality, water allocation. WEAP is a conceptual model, which allows the mapping of the physical system [43]. The elements of the natural system (watersheds, rivers, lakes) and the components of the technical system can be represented as a set of interconnected components. The WEAP contains various approaches for modeling water balance processes for a watershed [35]; the "Soil Moisture Model" method was selected for this study. The "Soil Moisture Model" method is considered the most efficient for simulating the processes of the hydrological cycle. The model requires different input variables for hydrological simulation, such as precipitation, temperature, vegetation cover, relative humidity, and latitude, from which water balance components such as infiltration, evapotranspiration, and runoff are estimated [17].

Calibration and Validation of WEAP
Calibration and validation of the models are necessary and used to ensure the quality and accuracy of the results; it involves improving reliability, by modifying the parameters so that the simulated results correspond closely to the measured data [44][45][46]. WEAP includes an interface to a parameter estimation tool (PEST) that allows the user to automate the process of comparing WEAP results with historical observations and to modify the model parameters to improve its accuracy. The PEST tool is a non-linear parameter, considered as a unique calibration tool. PEST executes, as many times as necessary, the control of a model, while changing its parameters until the differences between the simulations of the selected model and a set of observed data are minimized in the direction of weighted least squares [47]. PEST has five methods that allow calibration and validation. The rain-runoff method for model calibration was chosen because all required data are available. The calibration of the WEAP model for the Ourika river was performed using monthly streamflow gauge data at five stations, using the historical data from 2000 to 2008. Data from 2009 to 2014 were used for validation. Statistical parameters, such as Nash-Sutcliffe efficiency (NSE), the coefficient of determination (R 2 ), was taken into account, to evaluate the performance of the calibration results.

Demand Sites and Supply Sources
For modeling purposes, the Ourika watershed was divided into 3 sub-watersheds depending on soil types, drainage network, topography, and differences in climate and land cover. The main demand sites are Setti Fatma, Toubkal, and Oukaimden. These sites represent the counties of the Ourika watershed, where the population and its activities, including agriculture and livestock, are concentrated. Tourism is putting increasing pressure on water resources. The water consumption habits of tourists are very different from those of the local population, and their numbers are increasing annually. Its combination with the water demand of local residents and rural activities leads to the overexploitation of surface water. For agriculture, the irrigated croplands are the major consumers of water. The average water requirement is estimated at 7500 m 3 /ha. Irrigation of agricultural land is still carried out by the traditional method SEGUIAS. It is an open-air irrigation canal, adopted in North Africa, and is often found in oases. Given the climate of the study region, these systems involve significant water losses through evaporation, infiltration, or evacuation. All the sources of supply for existing and planned development in the Ourika watershed are supplied by the surface water resources of the main Ourika River and its tributaries.

Scenarios Development
Watershed-based water management is an approach that brings together all water users to develop planned, concerted, and consensual actions in the general interest. This participatory approach brings together people from different fields with sometimes divergent interests. The WEAP allows us to build different current scenarios and to simulate them in the future, to answer various ʺwhat ifʺ questions [48], and to compare them to the current situation. For this study, six scenarios were developed considering socio-economic factors and climate change (population growth, high population growth, A2 and B2 scenarios), to predict the gap between water supply and demand in the future.
We have also developed scenarios to simulate the reliability and sustainability of adaptation strategies that managers plan to implement, to mitigate the impacts of climate change on supply water resources. These management intervention options were identified to meet current and future water demands at a minimal cost. The six scenarios developed are described below: (1) Reference Scenario (RS): The first baseline scenario constructed was the current situation of water supply and demand conditions in the watershed. It serves as a basis for further analysis of the current state and comparison with the other simulated scenarios. The current accounts are chosen to serve, as the base year of the model and the entire information system is entered into the current accounts from which scenarios are created. The scenarios explore possible changes to the system in the years to come after the year of the current accounts. The current situation (2014) is extended in the future (2015-2100). Water demand was assumed to increase with population growth (1.18%) and water irrigation techniques will not change.
(2) High Population Growth Scenario (HPG): Under this scenario population, growth rate is increased from 1.18% to 3.5 %, and all the other parameters are kept as in the reference scenario.
(3) Climate Change Scenario (CC): The third scenario uses the results of downscaled climate change scenarios A2 and B2 with a population growth rate of 1.8% to drive the WEAP for climate change impact assessment.
(4) New Irrigation Techniques Scenario (NIT): Agriculture is the main source of income for the population in the Ourika watershed. This fragile ecosystem is increasingly degraded as a result of a combination of several factors, including climatic disturbances that are causing effects already felt by farmers. We evaluated the impact of drip irrigation as the first adaptation strategy to climate change under A2 and B2 scenarios. This technique allows saving between 25% and 30% of the irrigation water [24,29]. We assume that this technique will cover 50% of agricultural land by 2025 and 100% by 2030.
(5) Improvement of Living Conditions Scenario (ILC): The potential in terms of wastewater discharge at the level of the Ourika watershed is about 3.1 million cubic meters (MCM) in 2014. It would be 4.3 MCM in 2030, with an annual evolution of around 0.08 MCM. This potential is calculated from the populationʹs water consumption, with a return to the sewer system of 80%. Given that there are no wastewater treatment plants, the installation of a water treatment plant is a key adaptation strategy. We evaluated a wastewater treatment plant to save about 30% of domestic water. This strategy is evaluated under both climate change scenarios.
(6) The combination of NIT and ILC adaptation strategies: This scenario simulates the combination of the two assumptions NIT and ILC under A2 and B2 scenarios.

Downscaling Performance
To assess SDSM performance, a number of statistical parameters were used. The E% are 65%, 79%, and 21% for T-max, T-min, and precipitation. During the model calibration process, the R 2 value of the monthly (SDSM-M) and annual (SDSM-A) responses to A2 and B2 scenarios lies between 78.12 and 82.4%, for both T-max and T-min, and between 46% to 60% for precipitation. During the model validation with bias correction, the R 2 value for monthly and annual models was between 88% and 97.3% for both T-max and T-min and between 58% and 74% for precipitation. The bias correction improved almost all statistical indicators. The RMSE decreases to 9.76 mm/month, and RE_б decreases by 5.47%, as shown in Table 2.
These results indicate the strong applicability of SDSM to downscaling Tmax, Tmin, and precipitation under A2 and B2 scenarios after Bias correction was applied in the Ourika watershed. The results are globally satisfactory and comparable to previous studies' results [13,25]. This may be due to the use of a longer calibration period (30 years). A series of observed data of 30 years or older can be used to generate probability distribution functions in SDSM or statistical downscaling procedures [13].

Precipitation Projections
The downscaled precipitation shows a decreasing trend in average annual precipitation for both scenarios over all three-time horizons (Figure 4). Projected mean annual precipitation is shown in table 3. Mean annual precipitation decreased in A2 by 9.98% in the 20s, 33.56% in the 50s and 49.25% in the 80s and the B2 by about 8.21% in the 20s, 23.68% in the 50s and 34.61% in the 80s. The results show that the A2 scenario is more severe than B2 and resulted in more rainfall reduction. These results are consistent with other watershed climate projections studies in Morocco. As an example, we compare our precipitation projections to those of [49] and [50] and found that they present similar trends.

Temperature Projection
The Ourika watershed mean Tmax and Tmin are projected to be warmer than the corresponding baseline (Table 3). Both Tmax and Tmin are projected to increase under A2 and B2 and for all time horizons by about 1.6 °C, 2.7 °C, 4.2 °C and about 1.3 °C, 2.2 °C, 3.6 °C in the 20s, 50s, and 80s, under A2 and B2 scenarios, respectively ( Figure 5). The mean annual Tmin also increased in all two scenarios in all periods; the mean Tmin increase by 2 °C,2.6 °C, 3.5 °C and 1.6 °C, 2.1 °C, 2.9 °C in the 20s, 50s, and 80s, under A2 and B2. Similarly to several other studies ex; [51,52], the projected temperature increase is more severe under A2 than B2. In the Ourika watershed, these changes in precipitation and temperature are expected to reduce the agricultural and domestique water supply and impact water demand significantly.

WEAP Calibration and Validation
Calibration and validation of the WEAP model over the Ourika River were performed using monthly streamflow data from five-gauge stations. Model performance was evaluated with PEST by calculating NSE, R 2 , and r, for all the calibration and validation periods. The results show that the R 2 values ranged from 0.64% to 0.89%, with a mean of 0.76%. The NSE values ranged from 0.61 to 0.85 with a mean of 0.70. The efficiency of Nash-Sutcliffe varies from-∞ to 1 [53]. Effectiveness of 1 means that the modeled flow corresponds perfectly to the observed data [54]. Calibration results showed excellent agreement with the model validation result (Figure 6). This confirms that the WEAP model accurately simulated the streamflow in the study area. Similar results have been reported by [55]. The modelʹs performance is considered satisfactory to simulate long-term changes in water resource availability and to explore its sensitivity to human activities, climate change, and adaptation measures.

References Scenario (RS) and High Population Growth Scenario (HPG)
In Figure 7, we analyze the relationship between water supply and demand, in the reference scenario (RS), and in the high population growth scenario (HPG). The results show that population growth leads to an increase in water demand. The availability of water in a scenario of low population growth would induce insufficient water supply. For the (RS), the results show that, by 2100, the demand for domestic water will have increased. Likewise for the demand for water for agriculture and livestock. Estimates of future water supply and demand conditions indicate that the average water demand coverage in the Ourika watershed will decrease from 93.06% to 61.11% by 2100. The unmet water demand (quantity of water that cannot be physically delivered to the demand site; useful for understanding the extent of the shortage) was 2.5 MCM in 2014 and will be in the order of 14 MCM in 2100; that is, a more than five-fold increase. The reference scenario shows a water shortage in the supply system, while water demand increases. The second scenario (HPG) showed that a high population growth rate will increase unmet demand, as this scenario has a severe effect on water supply and demand. An increased population growth rate from 1.18% (RS) to 3.5% (HPG) induces an increase in water demand for all demand sectors. Consequently, the unmet demand is estimated at 21.8 MCM in 2100. The availability of water (HPG) would result in a much more limited amount of water. Unmet demand in the Ourika watershed is caused largely by population growth, the high per capita water use, as well as declining water inflows and increased human activities (agriculture, tourism, habitat, and deforestation). Hellegers et al. [56] have found that in the Middle East and North African countries, water availability is limited in the face of rapidly increasing water demand. The results obtained are in line with those obtained by [57], who concluded that, at the level of the Upper Indus Basin in Pakistan in 2014, average water demand coverage would decrease from 90% to 75% in 2050.

Climate Change Scenario (CC)
The WEAP model predicts an increase in future water demand in the Ourika watershed according to the two climate change scenarios A2 and B2. In the (CC) scenario, agricultural and domestic water demand will continue to increase. The results indicate that unmet water demand in 2100 is estimated at 52 MCM and 64 MCM for scenarios B2 and A2 (Figure 8). Future imbalances in water supply and demand are estimated at between 55% and 65% in 2100. A study [29] on the Rheraya watershed in Morocco reports that the average unmet water demand would significantly increase by 70% of the total annual demand by 2100. Another study [58] on water resource trends in the Middle East and North Africa under the impact of climate change shows that demand will increase by 50% and supply will decrease by 12% by 2050, leading to a considerable increase in water scarcity. On the other hand, a study carried in Mali [59] using the WEAP socio-economic and climate change scenarios model reported an increase in water demand from 76 MCM to 224 MCM by 2050. A study conducted in Morocco by [56] suggested that, for the average future climate projection, the water supply and demand imbalances will reach 15.7 billion m 3 per year, or about 60% of total future demand. In the Ourika watershed, potential future climate change may lead to early snowmelt, which could increase water availability at a time of year when demand is relatively low. This could have a negative effect on the dry-months (July-August) supply. The Ourika may experience a temporal opposition of phase between water availability and need. This results in a loss of a significant amount of water every month in the medium and low mountain areas. On the other hand, during summer and autumn, water needs are at a maximum level, while water availability is very low. In addition, increasing temperatures are likely to speed up the evaporation of water and land surfaces and accelerate transpiration. Several studies have shown the possible decrease in water flow in several regions of the world [60,61]. The decrease in water flow in the Ourika watershed is directly related to climate change, followed by anthropogenic over exploitation. Water stress will affect most watersheds in Morocco. In surface-water watershed, already suffering from water stress, population growth, tourism pressure, and climate change exacerbates existing water problems. The WEAP estimated that more than three-quarters of local water resources are consumed for agriculture to ensure local food security. Agriculture is the largest water consuming sector in the world, and the Ourika watershed is not an exception. In 2014 agricultural water demand was about 79%, with an additional 8% for livestock, while domestic water use reached 13%. Our study indicates that climate change is expected to have the strongest impact on the Ourika watershed, where the agriculture is very important for the local economy, and in particular for small farmers.

New Irrigation Technique Scenario (NIT)
This scenario was developed based on the low irrigation efficiency in the study area. The study results show that unmet water demand for agriculture would decrease by 28 MCM to reach 38 and 32 MCM, under scenarios A2 and B2, respectively ( Figure 9). This is due to the widespread use of drip irrigation, which reduces water losses through infiltration and evaporation. This decrease could be the result of farmers awareness and the introduction of new and efficient irrigation techniques, as suggested in the Moroccan Green Plan and the Moroccan National Program of Irrigation water Economy. The SEGUIAS irrigation method seems inadequate to ensure current resource management and mitigate the effects of drought. According to Schilling et al. [9], climate change will lead to a 15-40% drop in primary agricultural production in the 21st century in Morocco, and future droughts are likely to unbalance social stability, by causing severe shocks on food prices. The effects of climate change on water demand are expected to be more important in agriculture than in other sectors [62]. Santillán et al. [63] showed that, in Morocco, as well as in the warmest and driest areas, irrigation may not be enough, and planning strategies would be needed.

Improvement of Living Conditions (ILC)
In this scenario, we have assumed that 30% of wastewater is treated and used for various appropriate purposes such as cleaning, washing, and flushing. The results show that unmet water demand would decrease by 52 and 42 MCM in scenarios A2 and B2, respectively in 2100 ( Figure 10). This decrease could be the result of the design of a sewerage system, water treatment and the re-use of treated water.

Combination of the Two Adaptation Strategies (NIT + ILC)
The combination of the two previous adaptation strategies will probably have a positive impact on the management and conservation of water resources in the Ourika watershed. The unmet demand will remain present all along the century. The results show that unmet water demand will drop to 34 and 26 MCM in scenarios A2 and B2, respectively. This scenario's results show that unmet water demand would decrease by about 50% and 44% in scenarios B2 and A2, respectively ( Figure  11).
Increasing resilience to climate change, rational and rigorous management of water resources in the Ourika watershed, is of paramount importance for the sustainable development of the region, to meet the needs of the population, to develop the potential of agricultural land and to avoid water scarcity to be an obstacle to local socio-economic progress. To mitigate the impacts of climate change and strengthen the adaptive management capacity of the Ourika watershed, several management intervention options were simulated to meet water demands with minimal possible cost and environmental impact. The results obtained should be taken into account in planning water supply. Schilling et al. [9] mentioned that Moroccoʹs adaptive capacities are weaker than in Europe, but greater than the capacity of the African continent, which makes these countries even more vulnerable to climate change. Therefore, the combination of the two strategies does not ensure a complete supply in the Ourika watershed. Other studies showed similar results; Amin et al. [57] estimated that the installation of dams as an adaptation strategy would reduce unmet demand in 2025 to about 60% in the study basin; Birhanu et al. [64] indicated that no water management strategy would cover 100% of the water demand by 2025 in the city of Addis Ababa, also, the results of Guermazi et al. [65] showed that it is not enough to use a single optimal cropping strategy to avoid aquifer depletion by 2050, in the Regueb watershed in Tunisia. The planning and maintenance of water supply systems are another important aspect to increase water system resilience to climate change. Moreover, [66] confirms the usefulness of the planning water supply system based on climate change scenarios to ensure an optimal resource distribution in response to the expected reductions in water availability while satisfying current demands. Vulnerability measures for water distribution systems, as developed in [67], can inform about the capacity of the water distribution system to meet users' needs and help to carry out proper management and maintenance of the system to ensure adequate response to natural and made-man damages.
In the Ourika watershed, the scenario of the adaptation strategy for irrigation water management is more promising for water conservation in the future. It provided better results, in the face of pressures exacerbated by climate change and population growth. This was also proven by [49], who ranked irrigation first among the 14 adaptation measures studied. The above findings indicate that, by implementing appropriate management strategies, it may be possible to improve water resource management, the irrigation deficit, and save water consumption. We consider that the results also demonstrate that the assessments of the proposed adaptation strategies are effective, but not sufficient to ensure water sustainability for the Ourika watershed.

Conclusions
This study combines two models, the SDSM for future climate projection and the WEAP to simulate water supply and demand in the Ourika watershed, to assess the impact of climate change and socio-economic scenarios. Like the Mediterranean basins, the SDSM results show that the Ourika watershed will face an increase in temperature and a decrease in precipitation. The developed WEAP model shows that, in the future, climatic variables and socio-economic factors will be the two main drivers of water scarcity in the basin. The results indicate that climate change and population growth exacerbate the difference between water demand and unmet demand. The unmet water demand will be 64 MCM in 2100, under the A2 climate change scenario. All scenarios reveal water shortages. Even though the proposed adaptation strategies showed a positive impact on the long-term preservation of water resources, they are not sufficient to ensure water sustainability for the Ourika watershed. Thus, we conclude that better management policies are needed in the watershed, to ensure the sustainability of water resources in the long term.
In regard to the used data and modeling approach, as with all modeling studies, there are uncertainties in the GCM results for precipitation and temperature. Although parameter calibration ensures model applicability, the conceptual nature of the SDSM tends to overestimate temperature values and underestimate precipitation values. These sources of uncertainty are taken into account in the validation and analysis of the model output. Among this modeling approach limitations, we highlight the small number of climate stations in the study area. We also suggest that future studies in the Ourika watershed use a set of GCMs other than NCEP with the most recent RCP scenarios (representative concentration pathways). These studies could provide a better understanding of modeling and simulation uncertainties. However, the high-resolution climatological data presented in this study are shown to be useful for future climate change and hydrological studies in this region.