Climate changes can be defined as the Earth or regional global climate negative variations. These changes may be due to intrinsic processes in the Earth, to external influences, or more recently, to human activities  . Current climate changes are not considered as an ordinary change due to its magnitude and speed  .
Since the second Intergovernmental Panel on Climate Change report  has highlighted the lack of information on climate trends and variability, and daily climate extremes  , a number of literature studies of these changes have emerged, both for specific countries      , then synthesize this information between regions and globally    . These studies tended to focus on areas where the daily meteorological observations required for these analyze were already well spaced, controlled and archived  .
In Africa, numerous regional and national studies on recent trends and variability of the monthly climate synthetize that Africa has undergone in the last 50 years, one of the largest variations in rainfall observed, and noticeable disturbance of different climatic parameters      .
In Morocco, the quantification of climate change has become a major issue in terms of the dependence of agricultural production and the local economy on water availability. Any change in the frequency or severity of extreme weather and climate events could have a significant impact on nature and society  .
This research article aims to study the evolution of the observed climate change trends in the two regions Beni Mellal-Khenifra and Darâa-Tafilalet through climatic indices calculated on the basis of observations from the Regional Directorate of Meteorology of Beni Mellal (RDM) and climatological stations.
2. Materials and Methods
2.1. Study Area
Two administrative regions are attached to RDM, Beni Mellal-Khenifra region that is composed by four large natural groups: the mountain ranges, the piedmont or the Dir, the plateau of the Phosphates of Khouribga, and the Tadla Plain  . Darâa-Tafilalt region formed by two large watersheds Darâa and Ziz-Rheris  (Figure 1). Administratively, Beni Mellal-Khenifra region is constitute of five provinces (Azilal, Beni Mellal, Fquih Ben Salah, Khenifra and Khouribga), and 135 communes including 16 municipalities and 119 rural communes, while the Darâa-Tafilalt region is formed by five provinces (Errachidia, Midelt, Tinghir, Ouarzazate and Zagora), and 125 municipalities (16 in urban and 109 in rural areas).
The calculated indices in this research article have been done on the basis of daily temperatures observations (maximum and minimum), and precipitations measurement of 5 stations in the study area. These data were measured over a 46-year period (1970 to 2016). The meteorological stations used in this study (Figure 2) location are summarized in Table 1.
Figure 1. Geographical location of the study area.
Figure 2. Meteorological stations location on a digital terrain model.
Table 1. Used stations information.
Table 2. Martonne aridity index classification.
In this study we used an open source GIS software for cartography and elaborate all the maps with a World Shaded Relief background projected at WGS 84, and some different indexes has benne calculated for estimate climate trends.
2.3.1. Martonne Aridity Index
Aridity index is a number indicating how much more precipitation could be lost by evapotranspirtion if it were available than is actually lost at a givenlocation.
The climate type characterization in Beni Mellal-Khenifra and Darâa-Tafilalt regions has been done using the Martonne aridity index, described as the annual precipitation height (P) and the average annual temperature (T) ratio.
P: Annual rainfall height (in mm)
T: Annual temperature average (in ˚C) 
De Martonne's Index of Aridity has the obvious advantage of showing the transition from one area to another, but it may logically be attacked because of its empirical nature. Nevertheless, it serves as a useful discreet tool in illustrating the slow transition between arid, semiarid, and humid environments.
This index reflects the climatic conditions using the low average annual rainfall and the strong deficit comparison with temperature. The high value of this index due to high rainfall and/or low temperatures led to a wet climate, and the opposite to dry climate as summarized in Table 2.
In order to characterize the types of climate as well as their evolution, the Martonne aridity index was calculated in the observation period for each meteorological station as well as on two different periods 1970-1990, 1991-2010 and a summary period from 1970-2016.
2.3.2. Ombro-Thermic Index
The climate change at the inter-monthly scale characterization has been done using ombro-thermic index, that aims to distinct between dry and wet months that are not easily identifiable, or at least more or less similar.
Wet climate: Io > 0
Dry climate: Io < 0
2.3.3. RClimdex Indices
The used data quality in this work has been controlled by the Rclimdex software to identify any recording errors that may exist on daily data and influence the characterization results. This processing is based on checking and eliminating values that have anomalies, such as:
・ It is not possible to have more than 365 or 366 daily observations per year;
・ The month of February must not have more than 28 or 29 observations whatever the year considered;
・ The missing data is replaced by -99.9 before controlling the data quality;
・ The rain must be between 0 and 200 mm;
・ The maximum temperature must be higher than the minimum temperature.
Table 3. Summary information about the used precipitation indices.
Table 4. Summary information on the used temperature indices.
Based on verified data, the thermal evolution analysis in the Darâa-Tafilalt and Beni-Mellal Khénifra regions has been done using a set of climate indices, calculated annually through RClimdex software for the five stations. These indices allow easy comparison of general climate trends between different regions and climatic zones  . A total of 27 indices are proposed to detect climate change in any study area     . In our study case, we employed only 18 precipitation and temperature indices shown in Table 3 and Table 4. These indices are considered the most representative in climate change trends.
3. Results and Discussion
3.1. Ombro-Thermic Index
A semi-arid climate characterized both Beni-Mellal and Khouribga stations. The two stations located near the Middle Atlas benefit from the generated rain that are coming from Atlantic disturbances by reliefs. The hyper-aridity climate in Errachidia and Ouarzazate stations is mainly due to the location downstream of the Atlas, which constitutes a barrier to the oceanic influence of the Atlantic. Midelt is characterized by an arid desert climate caused by geographical position in Middle Atlas downstream, which are preventing the influence of classical disturbances from the Atlantic.
The Martonne aridity index I (equation (3)) results for the 3 periods showed that there are no changes in climate type for 5 studied stations (Table 5).
To characterize the Ombro-thermic index evolution, four periods (1970-1979, 1980-1989, 1990-1999 and 2000-2016) are used as indicator of general climate trend in the study area.
The ombro-thermal index calculation results (Table 6) showed that the wet period is characterized by six months in Khouribga station and seven months for Beni Mellal station on the totality of the studied period (Table 7).
On the other hand, Table 8 and Table 9 showed that the Errachidia and Ouarzazate stations climates during all the year are dry. Concerning, Midelt station climate (Table 10) is characterized by a short wet period (2 to 3 months) and a long dry period (9 to 10 months).
Table 5. Climate type in studied stations.
Table 6. Ombro-thermic index of Khouribga station.
Table 7. Ombro-thermic index of Beni Mellal station.
Table 8. Ombro-thermic index of Errachidia station.
Table 9. Ombro-thermic index of Ouarzazate station.
Table 10. Ombro-thermic index of Midelt station.
3.2. Indices Cartographies
3.2.1. General Trend of Temperature Indices
The observed indices of climate change are calculated from daily precipitation data and temperatures (maximum and minimum) in Khouribga, Beni-Mellal, Midelt, Errachidia and Ouarzazate stations. The data period used are summarized in Table 11.
Table 11. Reference periods used by RClimdex to calculate climate indices.
We are taking for example the observed trends evolution in Beni Mellal station such as illustrated in Figure 3. The part ((a)-(b)) represent relatively hot days and number of heat wave indices evolution. These two indices have a significant upward trend generally detected on both case. Concerning the rainfall, the part (c) shows a downward trend in the annual cumulative rainfall index, but this trend is not statistically significant.
1) Relatively cool nights (TN10p)
Figure 4 bellow shows a spatial distribution of TN10p. This index is down significantly and suggests −0.024% to −0.325% days/year, which leads to a decrease in the number of cold nights. The highest fall value is located in Beni Mellal city (−0.325% of days/year), Errachidia (−0.32% of days/year), Ouarzazate (−0.295% of days/year), and Midelt (−0.274% of days/year). These trends seem significant in all stations.
Figure 3. Observed trends evolution in Beni Mellal station for: (a) Relatively hot days index; (b) hot sequences durations index; (c) annual totals of precipitations index.
Figure 4. Spatial distribution of the index trends of the fresh nights.
2) Relatively warm days (TX90p)
Figure 5 represents the index evolution trend of the relatively warm TX90p days for the five stations of interest. All trends are positive and statistically significant. The hot days shows a real increase. This observation is valid for all the observations analyzed. This increase varies between 0.568% in days per year (5% days per decade) in Errachidia, Ouarzazate (0.489% of days/year), Khouribga (0.475% of days/year), Midelt (0.397% of days/year) and 0.276% days per year equivalent to 2% day by decade in Beni Mellal.
3) Warm Sequence Duration Indicator (WSDI)/Heat wave
The WSDI index (Figure 6) express the number of hot days in the year with at least six consecutive days or heat wave. Its results are significant at all stations in the study area. The maximum value is observed at Errachidia with 1.301% of days/year, followed by Ouarzazate 1.137% days/year and Khouribga 0.837% days/year. The minimum value is observed at Midelt with 0.581% days/year.
Synthesis: The observed global warming in the two regions (Beni Mellal-Khenifra and Darâa-Tafilalt) is characterized by a significant increase in relatively hot days and for the indicator of successive heat duration in Errachidia and Ouarzazate compared to other stations. The cool nights show a decrease on the whole of stations.
Figure 5. Index evolution trend of the relatively warm.
Figure 6. Warm Sequence Duration index variation in the different stations.
3.2.2. General Trend of Precipitation Indices
1) Maximum of one rainy day (RX1day)
The RX1day index map (Figure 7) expresses the maximum height of the precipitation in one day. This trend is more pronounced at Khouribga station where there has been a decrease in the order of (−0.41 mm) and in Midelt (−0.005 mm). On the other hand, Ouarzazate, Errachidia and Beni Mellal show an increase between 0.177 and 0.052 mm.
2) Very wet days (R95p)
Figure 8 shows that Beni Mellal, Ouarzazate and Errachidia stations have an increasing trend between 0.29 and 0.654 mm. This increase is much accentuated in Beni Mellal with 0.654 mm and in Ouarzazate 0.429 mm. The downward trend is noticeable especially in Midelt with −0.627 mm, and Khouribga −0.299 mm.
3) Annual rainfall totals (PRCPTOT)
Figure 9 represents the PRCPTOT index that is reflecting the total annual precipitations in rainy days. This totals rainfall index is characterized by a downward trend in the Beni Mellal stations −1.343 mm, Khouribga −1.23 mm and Midelt −0.417 mm. While in Errachidia and Ouarzazate stations are characterized by an increase in the order of 1.075 and 0.138 mm.
Figure 7. Maximum height of the precipitation in one day in different stations.
Figure 8. Very wet days index trends.
Figure 9. PRCPTOT index trends of total annual precipitations.
4) Consecutively dry days (CDD)
This index evolution is marked by a very considerable decline in Ouarzazate −0.647 days/year. In the other stations as Errachidia there is a rather low decline is that does not exceed −0.481 days/year. In Beni Mellal and Midelt station they are successively characterized by a rise of 2 days/decade and 3 days/decade. This increase is more pronounced in Khouribga, which has a trend of 7.65 days/decade of consecutive days without rain (Figure 10).
Synthesis: The two regions rainfall (Beni Mellal-Khenifra and Daraa-Tafilalt) is marked by a variable evolution. This recession is reflected in a downward trend in climate change indices for stations and an increase in other stations, particularly for rainfall totals, consecutive wet days and consecutive dry days.
Figure 10. Spatial distribution of consecutive day index trends.
The objective of this work is to provide a summary of extreme and average climate trends observed in the Beni Mellal-Khenifra and Darâa-Tafilalt regions over the last four decades. The indices were calculated based on daily maximum and minimum observations for temperatures and precipitation. Trends were calculated for five weather stations from 1970 to 2016.
For temperature trends, a large number of stations are experiencing significant trends, which would confirm an increase in temperature in the study area and subsequently widening of the hot season and shrinking of the cold season. The results also show a decrease in the number of cold nights and an increase in the number of hot days and indicator of the durations of hot sequences. In addition, precipitation behavior is difficult to analyze given the geographic position and terrain that influences their random distributions.
Precipitation decreases especially in Khouribga and Midelt. This decrease is materialized by a decrease in annual rainfall totals, consecutive wet days, and maximums of a rainy day. An upward trend of consecutively dry days is observed in Beni Mellal, Khouribga and Midelt. Note that these trends are not statistically significant.
In general, the different evolutions of the temperatures and precipitations did not induce a change of the type of climate in the studied regions.
 Bertrand, F. and Larrue, C. (2007) Gestion territoriale du changement climatique -Une analyse à partir des politiques régionales. Programme Gestion et Impacts du Changement Climatique-2. Université UMR CITERES.
 The Intergovernmental Panel on Climate Change (IPCC). http://www.ipcc.ch
 Alexander, L.V., Zhang, X., Peterson, T.C., Caesar, J., Gleason, B., Klein Tank, A.M.G., Haylock, M., Collins, D., Trewin, B. and Rahimzadeh, F. (2006) Global Observed Changes in Daily Climate Extremes of Temperature and Precipitation. Journal of Geophysical Research: Atmospheres, 111.
 Frei, C. and Schär, C. (2001) Detection Probability of Trends in Rare Events: Theory and Application to Heavy Precipitation in the Alpine Region. Journal of Climate, 14, 1568-1584.
 Groisman, P.Y., Knight, R.W., Easterling, D.R., Karl, T.R., Hegerl, G.C. and Razuvaev, V.N. (2005) Trends in Intense Precipitation in the Climate Record. Journal of Climate, 18, 1326-1350.
 Osborn, T.J. and Hulme, M. (2002) Evidence for Trends in Heavy Rainfall Events over the UK. Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, 360, 1313-1325. https://doi.org/10.1098/rsta.2002.1002
 Frich, P., Alexander, L.V., Della-Marta, P.M., Gleason, B., Haylock, M., Tank, A.K. and Peterson, T. (2002) Observed Coherent Changes in Climatic Extremes during the Second Half of the Twentieth Century. Climate Research, 19, 193-212. https://doi.org/10.3354/cr019193
 Kiktev, D., Sexton, D.M., Alexander, L. and Folland, C.K. (2003) Comparison of Modeled and Observed Trends in Indices of Daily Climate Extremes. Journal of Climate, 16, 3560-3571.
 Moberg, A., Sonechkin, D.M., Holmgren, K., Datsenko, N.M. and Karlén, W. (2005) Highly Variable Northern Hemisphere Temperatures Reconstructed from Low- and High-Resolution Proxy Data. Nature, 433, 613-617. https://doi.org/10.1038/nature03265
 Fauchereau, N., Trzaska, S., Rouault, M. and Richard, Y. (2003) Rainfall Variability and Changes in Southern Africa during the 20th Century in the Global Warming Context. Natural Hazards, 29, 139-154. https://doi.org/10.1023/A:1023630924100
 Malhi, Y. and Wright, J. (2004) Spatial Patterns and Recent Trends in the Climate of Tropical Rainforest Regions. Philosophical Transactions of the Royal Society of London, 359, 311-329.
 Sensoy, S., Türkoğlu, N., Akçakaya, A., Ekici, M., Demircan, M., Ulupinar, Y., Atay, H., Tüvan, A. and Demirbaş, H. (2013) Trends in Turkey Climate Indices from 1960 to 2010. 6th Atmospheric Science Symposium, Ankara, 24-26.
 Sohrabi, M.M., Marofi, S. and Ababaei, B. (2009) Investigation of Temperature and Precipitation Indices by Using RClimDex and R Software in Semnan Province. International Conference on Water Resources, 16-18.
 Zhang, X., Aguilar, E., Sensoy, S., Melkonyan, H., Tagiyeva, U., Ahmed, N., Kutaladze, N., Rahimzadeh, F., Taghipour, A. and Hantosh, T.H. (2005) Trends in Middle East Climate Extreme Indices from 1950 to 2003. Journal of Geophysical Research: Atmospheres, 110, D22104.