minimale exigée de 70 % (Hair et al., 1998). Cette dernière préoccupation est aussi bien confirmée par le tracé d'effondrement, qui est illustré dans la Figure 17.

<figure><img src="image_19.png" /><figcaption>Figure 17 : Tracé d'effondrement de valeurs propres en axes de l'ACP sur les variables périodiques de tendances des précipitations au périmètre de Tadla (Période de 1958 à 2018).</figcaption></figure>

Le « coude de Cattell » (Hair et al., 1998) selon son allure sur graphique, coïncide effectivement en coordonnée à l'abscisse pour la composante F2 de l'ACP envisagée (Figure 17).

### d) Eléments caractéristiques de l'ACP sur les variables traitées

Les valeurs des composantes de l'ACP sur les variables analysées, leurs qualités de représentations et les contributions respectives sont consignées au tableau ci-dessous.

<table><caption>Tableau 59 : Valeurs de composantes, qualités de représentation et contributions en ACP pour les variables périodiques des précipitations au périmètre de Tadla (Période de 1958 à 2018).</caption><thead><tr><th rowspan="2">Variables de tendances périodiques de précipitations</th><th colspan="2">Composantes</th><th colspan="3">Qualités de représentation</th><th colspan="2">Contributions (%)</th></tr><tr><th>F1</th><th>F2</th><th>Cos<sup>2</sup> (1)</th><th>Cos<sup>2</sup> (2)</th><th>Extraction</th><th>F1</th><th>F2</th></tr></thead><tbody><tr><td></td><td>(1)</td><td>(2)</td><td>(3)=(1)<sup>2</sup></td><td>(4)=(2)<sup>2</sup></td><td>(3) + (4)</td><td>(5)=[100*(1)<sup>2</sup>/(Valeur propre F1=3,714)]</td><td>(6)=[100*(2)<sup>2</sup>/(Valeur propre F2=1,141)]</td></tr><tr><td>Septembre</td><td>0,397</td><td>-0,678</td><td>0,157</td><td>0,460</td><td>0,617</td><td>4,238</td><td>40,279</td></tr><tr><td>Décembre</td><td>0,765</td><td>0,077</td><td>0,585</td><td>0,006</td><td>0,591</td><td>15,742</td><td>0,520</td></tr><tr><td>Avril</td><td>0,860</td><td>0,262</td><td>0,739</td><td>0,069</td><td>0,808</td><td>19,911</td><td>6,006</td></tr><tr><td>Juin</td><td>-0,273</td><td>0,730</td><td>0,074</td><td>0,533</td><td>0,608</td><td>2,002</td><td>46,757</td></tr><tr><td>MAM</td><td>0,845</td><td>0,265</td><td>0,714</td><td>0,070</td><td>0,784</td><td>19,223</td><td>6,168</td></tr><tr><td>SemH</td><td>0,810</td><td>-0,052</td><td>0,655</td><td>0,003</td><td>0,658</td><td>17,647</td><td>0,236</td></tr><tr><td>SemS</td><td>0,888</td><td>0,002</td><td>0,788</td><td>5E-06</td><td>0,788</td><td>21,227</td><td>5E-04</td></tr></tbody></table>
Les points proéminents des analyses des résultats aboutis sont les suivants (Tableau 59) :

i) Qualités des représentations :

*Les qualités totales d'extraction des variables (somme de cos²) varient dans la fourchette de 0,59 à 0,81 et indiquent alors que le maximum d'information est bien expliqué.

*Sur l'axe F1, la qualité d'extraction (cos²) des variables est importante pour le groupe de 5 variables : Décembre, Avril, MAM, SemH et SemS, ayant des valeurs comprises entre 0,59 et 0,79, mais cette qualité demeure très faible pour les deux variables de Septembre et Juin, à valeurs inférieures à 0,20.

*Sur l'axe F2, la qualité d'extraction (cos²) est élevée pour les deux variables de Septembre et de Juin, qui est de l'ordre de la valeur de 0,50. Comparées à l'autre groupe en 5 variables restantes (Décembre, Avril, MAM, SemH et SemS) qui présentent de faible niveau en cette valeur de moins de 0,07.

ii) Contributions des variables traitées :

*En axe F1, les fortes contributions à l'extraction d'information qui concernent les 5 variables : Décembre, Avril, MAM, SemH et SemS, sont dans la gamme du taux de 16 à 21 % par variable ; contre la faible contribution unitaire de moins 4 % pour les deux variables de Septembre et Juin.

*En axe F2, la contribution des deux variables de Septembre et de Juin, est d'une forte qualité d'extraction avec le pourcentage de 40 à 47 %, comparées à l'autre groupe de 5 variables (Décembre, Avril, MAM, SemH et SemS) ayant une qualité d'extraction plus réduite et qui ne dépasse pas le taux unitaire de 6%.

e. Représentations graphiques sur axes factoriels de l'ACP

Le graphique du « cercle des corrélations » qui caractérise notre ACP pour les variables périodiques des tendances des précipitations périodiques pour le périmètre de Tadla (zone d'étude) est montré à la Figure 18.
<figure><img src="image_20.png" /><figcaption>Figure 18 : Cercle des corrélations de l'ACP pour les variables périodiques de tendances de précipitations au périmètre de Tadla (Période de 1958 à 2018)</figcaption></figure>

L'exploration graphique de l'ACP sur les axes F1-F2 des tendances périodiques pluviométriques, nous renseigne largement par les éléments suivants (Figure 18) :

- L'importance des représentations des variables est matérialisée par la proximité à la circonférence du « cercle de corrélations » pour les 3 variables mensuelles (Septembre, Avril et Juin), pour la seule variable trimestrielle Mars-Avril-Mai (MAM) et pour les deux variables du semestre sec (SemS) et du semestre humide (SemH).

- L'ampleur des représentations est relativement réduite pour la seule variable restante à savoir pour la période mensuelle de Décembre.

- Les délimitations distinctes sont plus marquantes sur le graphique pour deux groupes homogènes de variables, dont le premier groupe (I) concerne les 5 variables : Décembre, Avril, MAM, SemH et SemS ; le second groupe (II) est consacré aux 2 variables se rapportant aux tendances des précipitations mensuelles de Septembre et de Juin.
## V-4-Spatialisation des résultats de l'ACP pour les tendances des précipitations périodiques

Le graphique du « *tracé double* » de l'ACP menée sur les variables des tendances des précipitations périodiques pour le périmètre de Tadla, a été conçu en considérant les dites variables comme les « *variables de l'ACP* » et les stations des mesures pluviométriques de la zone d'étude en tant que les « *individus de l'ACP* ». L'illustration de ce « *tracé double* » est portée à la Figure 19.

<figure><img src="image_21.png" /><figcaption>**Figure 19 :** Tracé-double de l'ACP sur les tendances des variables des précipitations périodiques au périmètre de Tadla (Période d'observations de 1958 à 2018).</figcaption></figure>

A l'issue de la configuration spatiale sous le support du « *tracé double* » de l'ACP (Figure 19), les groupements des variables des tendances des précipitations périodiques sont alors affectés aux sites géographiques des stations de mesures pluviométriques de la zone étudiée. Les faits sont récapitulés aux Tableaux 60 et 61.

<caption>**Tableau 60 :** 1er groupement géographique des sites des stations pluviométriques en fonction des variables des tendances des précipitations périodiques, selon l'ACP au périmètre de Tadla</caption>

(Période d'observations de 1958 à 2018)

### *Groupement 1 (ACP)*

*Repéré par les 13 sites pluviomètres :*

-Stations de la zone des Béni Amir : S3 ; S5 ; S7 ; S8 ; S10 ; SAH

-Stations de la zone des Béni Moussa : S20 ; S21 ; S23 ; S25 ; S27 ; S35 ; SOG

*Variables périodiques des précipitations :*

+2 mensuelles : Décembre et d'Avril

+1 trimestrielle : Mars-Avril-Mai, et

+ 2 semestrielles : sèche et humide

En l'occurrence, la régionalisation des tendances périodiques des précipitations est marquée par deux groupements spatiaux homogènes à savoir: le groupement G I, qui est géo-référencé par 13
sites de stations pluviomètres et caractéristiques des 5 variables (mois de Décembre et d'Avril, trimestre printanier (Mars-Avril-Mai) et les 2 semestres) (Tableau 60). Par ailleurs, l'autre groupement G II qui est jalonné par 14 sites de mesures pluviométriques, se rapporte aux deux variables mensuelles de Septembre et de Juin (Tableau 61).

<caption>Tableau 61 : 2ème groupement géographique des sites des stations pluviométriques en fonction des variables des tendances des précipitations périodiques, selon l'ACP au périmètre du Tadla</caption>

(Période d'observations de 1958 à 2018)

Groupement 2 (ACP)

*Repéré par 14 sites pluviomètres :

-Stations de la zone des Béni Amir : S1 ; S4 ; S11 ; S12

-Stations de la zone des Béni Moussa : S26 ; S28 ; S29 ; S30 ; S31 ; S32 ; S33 ; S34 ; S36 ; S38

*Variables périodiques des précipitations :

+2 mensuelles : Septembre et Juin

La disposition géographique de ces groupements spatiaux d'ordre homogène pour les tendances périodiques des pluies est montrée sur la carte de la Figure 20.

<figure><img src="image_22.png" /><figcaption>Figure 20 : Régionalisation de tendances périodiques de précipitations au périmètre du Tadla</figcaption></figure>

(Période d'observations de 1958 à 2018).

Il est constaté visiblement sur l'espace physique du périmètre de Tadla que les tendances des pluies périodiques au sein d'une même année agricole, se cantonnent en de deux groupements identifiés avec leur délimitation géographique correspondant et à savoir les suivants (Figure 20) :
- 1er groupement géographique :

+Se rapporte à cinq variables de tendances des pluviométriques périodiques: Décembre ; Avril ; trimestre Mars-Avril-Mai ; semestres sec et humide) :

+Représente spatialement la partie Est du périmètre de Tadla, y compris une portion de la zone du piémont (Dir) au limitrophe de la chaine montagneuse du Moyen Atlas marocain.

- 2ème groupement géographique Gg2 :

+Concerne deux variables de tendances des pluviométriques périodiques : Septembre et Juin) :

+Est spécifique à la partie Ouest du périmètre de Tadla, aux voisinages à l'Est de la plaine du Haouz et au Sud de la plaine du Chaouia.
# Chapitre III : VARIABILITES THERMIQUE ET DE LA DEMANDE CLIMATIQUE EN EVAPORATION AU PERIMETRE DE TADLA

## I-Descriptifs thermique et de la demande climatique en évaporation au périmètre de Tadla

### I-1- Caractéristique de la température au périmètre du Tadla

La température annuelle moyenne au périmètre de Tadla est aux alentours de 19°C (Tableau 39C). La zone se caractérise par la température moyenne des trois mois les plus froids de l'année de 4,7°C contre 35,5°C pour les trois mois les plus chauds, soit un écart thermique de l'ordre de 30 °C (Tableau 62).

<caption>Tableau 62 : Températures moyennes annuelles des mois les plus froids et des mois les plus chauds de l'année par stations au périmètre de Tadla</caption>

(Période d'observations de 1970 à 2018)

<table><thead><tr><th rowspan="2">Stations des mesures</th><th colspan="4">Trois mois les plus froids de l'année</th><th colspan="4">Trois mois les plus chauds de l'année</th><th rowspan="2">Moyenne annuelle</th></tr><tr><th>Décembre</th><th>Janvier</th><th>Février</th><th>Moyenne</th><th>Juin</th><th>Juillet</th><th>Août</th><th>Moyenne</th></tr></thead><tbody><tr><td>AS</td><td>5,2</td><td>4,1</td><td>5,5</td><td>4,9</td><td>31,9</td><td>36,4</td><td>36</td><td>34,8</td><td>18,9</td></tr><tr><td>FBS</td><td>5,1</td><td>4</td><td>5,7</td><td>4,9</td><td>32,9</td><td>37,9</td><td>37,2</td><td>36,0</td><td>19,7</td></tr><tr><td>OG</td><td>4,6</td><td>3,4</td><td>5</td><td>4,3</td><td>32,4</td><td>37,2</td><td>36,7</td><td>35,4</td><td>18,9</td></tr><tr><td>MBA</td><td>4,9</td><td>3,6</td><td>5,2</td><td>4,6</td><td>32,6</td><td>37,4</td><td>36,9</td><td>35,6</td><td>19,2</td></tr><tr><td>Moyenne</td><td>5,0</td><td>3,8</td><td>5,3</td><td>4,7</td><td>32,4</td><td>37,2</td><td>36,7</td><td>35,5</td><td>19,2</td></tr><tr><td>CV (%)</td><td>5%</td><td>7%</td><td>5%</td><td>5%</td><td>1%</td><td>1%</td><td>1%</td><td>1%</td><td>2%</td></tr></tbody></table>

Il est observé que l'évolution de tendance de la température annuelle est dans le sens de la hausse sur la période de cinq dernières décennies (entre 1970 et 2018), comme illustré au graphique de la Figure 21.

<figure><img src="image_23.png" /><figcaption>Figures 21 : Evolution de température moyenne annuelle au périmètre de Tadla (Région Béni Mellal-Khénifera, Maroc)</figcaption></figure>

(Période d'observations de 1970 à 2018)

Ceci trouve explication dans la tendance à l'augmentation de la moyenne annuelle des maxima de température en l'air ambiant du périmètre de Tadla, représenté sur graphique de la Figure 22.
<figure><img src="image_24.png" /><figcaption>**Figures 22: Evolution de température maximale moyenne annuelle**
**au périmètre de Tadla (Région Béni Mellal-Khénifera, Maroc)**
(Période d'observations de 1970 à 2018)</figcaption></figure>

De même pour, la moyenne annuelle des minimas de température qui a évolué dans la tendance à l'élévation sur le périmètre de Tadla, représenté graphiquement en Figure 23.

<figure><img src="image_25.png" /><figcaption>**Figures 23: Evolution de température minimale moyenne annuelle**
**au périmètre de Tadla (Région Béni Mellal-Khénifera, Maroc)**
(Période d'observations de 1970 à 2018)</figcaption></figure>

A l'examen des températures extrêmes journalière enregistrées au périmètre de Tadla, il est retenu en principe que le niveau minimum absolu de l'année est atteint à la moyenne de -1,3 +C sur un intervalle périodique restreint le plus souvent entre le 30 Décembre et le 25 Janvier.
Quant à la valeur maximale absolue de température de l'année, elle passe à la limite supérieure moyenne de 45,4 °C sur l'intervalle saisonnier plus fréquent entre le 10 Juillet et le 12 Août en chaque année (Tableau 63).

<caption>Tableau 63 : Températures extrêmes moyennes de l'année par stations au périmètre de Tadla</caption>

(Période d'observations de 2014 à 2018)

<table><thead><tr><th rowspan="2">Stations des mesures</th><th colspan="2">Température maximale absolue</th><th colspan="2">Température minimale absolue</th></tr><tr><th>Valeur (°C)</th><th>Intervalle de fréquence d'enregistrement</th><th>Valeur (°C)</th><th>Intervalle de fréquence d'enregistrement</th></tr></thead><tbody><tr><td>AS</td><td>45,2</td><td>17 Juillet au 03 Août</td><td>-0,6</td><td>30 Décembre au 25 Février</td></tr><tr><td>FBS</td><td>45,7</td><td>10 Juillet au 05 Août</td><td>-0,8</td><td>01 Janvier au 11 Février</td></tr><tr><td>OG</td><td>45,3</td><td>15 Juillet au 07 Août</td><td>-2,2</td><td>30 Décembre au 06 Février</td></tr><tr><td>MBA</td><td>45,4</td><td>23 Juin au 12 Août</td><td>-1,5</td><td>02 Janvier au 06 Février</td></tr><tr><td>Moyenne</td><td>45,4</td><td>10 Juillet au 12 Août</td><td>-1,3</td><td>30 Décembre au 25 Février</td></tr><tr><td>CV (%)</td><td>0,4%</td><td>-</td><td>-49,4%</td><td>-</td></tr></tbody></table>

## I-2- Caractéristique de l'évapo-transpiration de référence (ETo) au périmètre de Tadla

L'ETo du périmètre de Tadla est d'une moyenne annuelle de 3,6 mm/jour, soit l'équivalent de 1 314 mm par an ou rapporté à la surface est pour valeur de 13 140 m3 à l'hectare (Tableau 64).

<caption>Tableau 64 : Evapotranspiration de référence (ETo) annuelle par stations du périmètre de Tadla</caption>

(Période d'observations de 1970 à 2018)

<table><thead><tr><th>Stations des mesures</th><th>Moyenne annuelle de l'ETo<br>(en mm/jour)</th></tr></thead><tbody><tr><td>AS</td><td>3,1</td></tr><tr><td>FBS</td><td>3,6</td></tr><tr><td>OG</td><td>3,5</td></tr><tr><td>MBA</td><td>3,9</td></tr><tr><td><strong>Moyenne</strong></td><td><strong>3,6</strong></td></tr><tr><td>CV (%)</td><td>7%</td></tr></tbody></table>

Egalement, la tendance de la demande climatique en évaporation exercée sur le végétal et qui est simulée par l'évapotranspiration de référence, a été mesurée au périmètre de Tadla depuis le début des années 1970 et illustré à la Figure 24.
<figure><img src="image_26.png" /><figcaption>Figures 24 : Evolution de l'évapo-transpiration de référence moyenne annuelle au périmètre du Tadla (Région Béni Mellal-Khénifera, Maroc)
(Période d'observations de 1970 à 2018)</figcaption></figure>

L'évapotranspiration de référence au périmètre de Tadla évolue dans un sens de l'accroissement tendanciel, comme visualisé graphiquement sur la Figure 24.

## II-Résultats des tests de tendance thermiques et de la demande climatique en évaporation

### II-1-Tests préliminaires de vérification de la normalité des distributions des séries de mesures

Les résultats des tests préliminaires statistiques sur les séries de mesures d'indicateurs agro-climatiques montrent que les distributions suivent la loi normale pour la majorité des situations (Taux de validation : 70%). Le détail des données de ces tests est rapporté dans le Tableau 65

<table><caption>Tableau 65 : Résultats des tests de la normalité précédant les tests d'homogénéité sur les séries agro-climatiques au périmètre du Tadla.
(Période d'observations de 1970 à 2018)</caption><thead><tr><th>Indicateur ou<br>indice climatique</th><th>Station<br>métrologique</th><th>Distribution</th></tr></thead><tbody><tr><td rowspan="4">Température minimale</td><td>OG</td><td>Normale</td></tr><tr><td>AS</td><td>Normale</td></tr><tr><td>MBA</td><td>Normale</td></tr><tr><td>FBS</td><td>Normale</td></tr><tr><td rowspan="4">Température maximale</td><td>OG</td><td>Non normale</td></tr><tr><td>AS</td><td>Non normale</td></tr><tr><td>MBA</td><td>Normale</td></tr><tr><td>FBS</td><td>Normale</td></tr><tr><td rowspan="4">Température moyenne</td><td>OG</td><td>Normale</td></tr><tr><td>AS</td><td>Normale</td></tr><tr><td>MBA</td><td>Normale</td></tr><tr><td>FBS</td><td>Normale</td></tr><tr><td rowspan="4">Evapo-transpiration</td><td>OG</td><td>Non normale</td></tr><tr><td>AS</td><td>Non normale</td></tr><tr><td>MBA</td><td>Non normale</td></tr><tr><td>FBS</td><td>Non normale</td></tr><tr><td rowspan="4">Indice de De Maritonne</td><td>OG</td><td>Normale</td></tr><tr><td>AS</td><td>Normale</td></tr><tr><td>MBA</td><td>Normale</td></tr><tr><td>FBS</td><td>Normale</td></tr></tbody></table>
## II-2-Tests de tendance de variabilité et d'homogénéité thermique et de la demande climatique en évaporation

Il est présenté ci-après les résultats de tendance de la variabilité thermique et de la demande climatique en évaporation selon la méthode Mann-Kendall et aussi concernant les tests en quatre méthodes menées séparément au sujet de l'homogénéité des données ; ce en quête d'identifier d'éventuelle rupture sur les séries chronologiques périodiques étudiées.

### II-2-1- Test de Mann-Kendall sur les températures et sur l'évapo-transpiration

Les résultats de ce test sont consignés dans le Tableau 66.

<caption>Tableau 66 : Résultats des tests statistiques de tendance par la méthode de Mann-Kendall sur les indicateurs agro-climatiques au périmètre du Tadla</caption>

(Période observations de 1970 à 2018).

<table><thead><tr><th rowspan="2">Indicateur ou indice climatique</th><th rowspan="2">Station Météorologique</th><th colspan="3">Paramètres du test de Mann-Kendall</th></tr><tr><th>U Mann-Kendall</th><th>Hypothèse Ho: Absence de tendance</th><th>Niveau de signification du test</th></tr></thead><tbody><tr><td rowspan="4">Température minimale</td><td>OG</td><td>4,9133</td><td>Non</td><td>Hautement significatif</td></tr><tr><td>AS</td><td>2,4998</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>BAM</td><td>4,7064</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>FBS</td><td>3,1031</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td rowspan="4">Température maximale</td><td>OG</td><td>5,1719</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>AS</td><td>5,5684</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>BAM</td><td>5,4822</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>FBS</td><td>5,6374</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td rowspan="4">Température moyenne</td><td>OG</td><td>6,1890</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>AS</td><td>5,890</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>BAM</td><td>6,4994</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>FBS</td><td>6,7752</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td rowspan="4">Evapo-Transpiration</td><td>OG</td><td>3,6548</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>AS</td><td>6,1546</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>BAM</td><td>3,9824</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>FBS</td><td>5,3443</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td rowspan="4">Indice de De Maritonne</td><td>OG</td><td>3,9996</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>AS</td><td>3,1376</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>BAM</td><td>4,2065</td><td>Non</td><td>Très hautement significatif</td></tr><tr><td>FBS</td><td>2,4480</td><td>Non</td><td>Très hautement significatif</td></tr></tbody></table>

La totalité des tests par la méthode Mann-Kendall, met en évidence l'existence en la quasi majorité des cas de variabilités à des niveaux très hautement significatifs pour les tendances des paramètres traités. Ceci approuve que des modifications effectives sont survenues au cours du temps sur l'évolution sur de différentes séries climatiques en question.

Faits qui ouvrent une large voie pour mener en relais d'autres tests encore plus poussés et notamment sur homogénéité temporelle de ces séries périodiques agro-climatiques.
## II-2-2- Tests d'homogénéité sur la température et sur l'évapo-transpiration

Les résultats des quatre tests réalisés sur l'homogénéité des tendances des paramètres agro-climatiques, qui ont concernés le test de Pettitt, le test de Buishand, le test de Lee et Heghinian et la méthode de segmentation de Hubert sont reportés et analysés ci-après.

### a) Test de rupture de Pettitt

Les données abouties par ce test sont récapitulées dans le tableau suivant.

<caption>Tableau 67 : Tests de détection de rupture de Pettitt sur les chronologies agro-climatiques au périmètre du Tadla</caption>

(Période d'observations de 1970 à 2018)

<table><thead><tr><th rowspan="2">Indicateur ou indice climatique</th><th rowspan="2">Station météorologique</th><th rowspan="2">U<sub>Pettitt</sub></th><th colspan="4">Paramètres du test de Pettitt</th><th rowspan="2">Dates de rupture</th></tr><tr><th>Hypothèse Ho: Absence de tendance</th><th>Probabilité dépassement</th><th>Niveau de signification du test</th></tr></thead><tbody><tr><td rowspan="4">Température minimale</td><td>OG</td><td>-426</td><td>Non</td><td>2,30E-04</td><td>Très hautement significatif</td><td>1989</td></tr><tr><td>AS</td><td>-</td><td>Oui</td><td>-</td><td>Non significatif</td><td>-</td></tr><tr><td>BAM</td><td>-402</td><td>Non</td><td>6,21E-04</td><td>Très hautement significatif</td><td>1989</td></tr><tr><td>FBS</td><td>-255</td><td>Non</td><td>7,76E-10</td><td>Significatif</td><td>1985</td></tr><tr><td rowspan="4">Température maximale</td><td>OG</td><td>-</td><td>Oui</td><td>-</td><td>Non significatif</td><td>-</td></tr><tr><td>AS</td><td>-570</td><td>Non</td><td>8,39E-08</td><td>Très hautement significatif</td><td>1996</td></tr><tr><td>BAM</td><td>-499</td><td>Non</td><td>7,88E-06</td><td>Très hautement significatif</td><td>1997</td></tr><tr><td>FBS</td><td>-564</td><td>Non</td><td>2,49E-07</td><td>Très hautement significatif</td><td>1997</td></tr><tr><td rowspan="4">Température moyenne</td><td>OG</td><td>-378</td><td>Non</td><td>1,36E-06</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>AS</td><td>-506</td><td>Non</td><td>5,54E-06</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>BAM</td><td>-531</td><td>Non</td><td>1,52E-06</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>FBS</td><td>-571</td><td>Non</td><td>1,68E-07</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td rowspan="4">Evapo-transpiration (ETo)</td><td>OG</td><td>-281</td><td>Non</td><td>3,86E-02</td><td>Très hautement significatif</td><td>1983</td></tr><tr><td>AS</td><td>-345</td><td>Non</td><td>5,22E-03</td><td>Très hautement significatif</td><td>2000</td></tr><tr><td>BAM</td><td>-</td><td>Oui</td><td>-</td><td>Non significatif</td><td>-</td></tr><tr><td>FBS</td><td>-324</td><td>Non</td><td>1,05E-02</td><td>Hautement significatif</td><td>1983</td></tr><tr><td rowspan="4">Indice de De Maritonne</td><td>OG</td><td>343</td><td>Non</td><td>5,50E-03</td><td>Très hautement significatif</td><td>1991</td></tr><tr><td>AS</td><td>286</td><td>Non</td><td>3,35E-02</td><td>Hautement significatif</td><td>1998</td></tr><tr><td>BAM</td><td>364</td><td>Non</td><td>2,66E-03</td><td>Très hautement significatif</td><td>1991</td></tr><tr><td>FBS</td><td>-</td><td>Oui</td><td>-</td><td>Non significatif</td><td>-</td></tr></tbody></table>

Il est identifié par ce test des ruptures sur les séries chronologiques à propos de tous les paramètres étudiés. En effet, il a été repéré pour les températures minimales, maximales et moyennes, pour l'ETo et pour l'Indice d'aridité de De Maritonne des ruptures chronologiques circonscrites aux intervalles respectifs suivants : 1985-1989; 1996-1997; {1994}; 1983-2000; 1991-1998 (Tableau 67).

Il est à noter aussi que ce test s'est accompli avec la prédominance du niveau « très hautement significatif » pour les cas envisagés.

### b) Test de rupture de Buishand

Les résultats de ce test sont synthétisés dans le tableau ci-après.
<caption>Tableau 68 : Tests de la détection de rupture de Buishand sur les séries temporelles
agro-climatiques au périmètre du Tadla</caption>

(Période d'observations de 1970 à 2018)

<table><thead><tr><th rowspan="2">Indicateur ou indice Climatique</th><th rowspan="2">Station Météorologique</th><th rowspan="2">U<sub>Max</sub></th><th colspan="3">Paramètres du test de de Buishand</th></tr><tr><th>Hypothèse Ho: Absence de tendance</th><th>Niveau signification en « ellipse du Bois »</th><th>Dates de rupture</th></tr></thead><tbody><tr><td rowspan="4">Température minimale</td><td>OG</td><td>-9,6735</td><td>Non</td><td>Très hautement significatif</td><td>1988</td></tr><tr><td>AS</td><td>-4,5184</td><td>Non</td><td>Significatif</td><td>1977</td></tr><tr><td>BAM</td><td>-9,2776</td><td>Non</td><td>Très hautement significatif</td><td>1989</td></tr><tr><td>FBS</td><td>-5,700</td><td>Non</td><td>Très hautement significatif</td><td>1985</td></tr><tr><td rowspan="4">Température maximale</td><td>OG</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td>AS</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td>BAM</td><td>-16,7184</td><td>Non</td><td>Très hautement significatif</td><td>1999</td></tr><tr><td>FBS</td><td>-0,1072</td><td>Non</td><td>Très hautement significatif</td><td>1997</td></tr><tr><td rowspan="4">Température moyenne</td><td>OG</td><td>-13,1122</td><td>Non</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>AS</td><td>-11,8612</td><td>Non</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>BAM</td><td>-13,9367</td><td>Non</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td>FBS</td><td>-15,8837</td><td>Non</td><td>Très hautement significatif</td><td>1994</td></tr><tr><td rowspan="4">Evapo-transpiration</td><td>OG</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td>AS</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td>BAM</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td>FBS</td><td>-</td><td>Oui</td><td>Non significatif</td><td>-</td></tr><tr><td rowspan="4">Indice de De Maritonne</td><td>OG</td><td>38,9163</td><td>Non</td><td>Très hautement significatif</td><td>1991</td></tr><tr><td>AS</td><td>29,6265</td><td>Non</td><td>Très hautement significatif</td><td>1979</td></tr><tr><td>BAM</td><td>41,8347</td><td>Non</td><td>Très hautement significatif</td><td>1979</td></tr><tr><td>FBS</td><td>29,3571</td><td>Non</td><td>Très hautement significatif</td><td>1979</td></tr></tbody></table>

En l'occurrence, le test de Buishand repère aussi des ruptures sur les séries temporelles pour tous les paramètres mise en épreuve, à savoir pour les températures minimales, maximales et moyennes et pour l'Indice de De Maritonne, sur les intervalles propices respectifs : 1977-1989; 1997-1999 ; {1994} ; {1994} ; 1979-1991.

Il est à signaler qu'en présence de ruptures en question, le niveau du test se qualifie de très hautement significatif ; avec l'exception pour le paramètre de l'ETo, dont le test n'a décelé aucune rupture sur sa série pluri-annuelle de données (Tableau 68).

### c) Test de rupture de Lee et Heghinian

Les résultats essentiels du test sont assemblés dans le tableau ci-dessous.
<caption>Tableau 69 : Tests de détection de rupture de Lee et Heghinian sur séries chronologiques agro-climatiques au périmètre du Tadla</caption>

(Période d'observations de 1970 à 2018)

<table><thead><tr><th rowspan="2">Indicateur ou<br>indice climatique</th><th rowspan="2">Station<br>météorologique</th><th colspan="2">Paramètres et résultats du test de de Buishand</th></tr><tr><th>Probabilité associée</th><th>Dates de<br>rupture</th></tr></thead><tbody><tr><td rowspan="4">Température minimale</td><td>OG</td><td>0,2358</td><td>1989</td></tr><tr><td>AS</td><td>0,2546</td><td>1977</td></tr><tr><td>BAM</td><td>0,1692</td><td>1989</td></tr><tr><td>FBS</td><td>0,7046</td><td>1997</td></tr><tr><td rowspan="4">Température maximale</td><td>OG</td><td>-</td><td>-</td></tr><tr><td>AS</td><td>-</td><td>-</td></tr><tr><td>BAM</td><td>0,4617</td><td>1999</td></tr><tr><td>FBS</td><td>0,7046</td><td>1997</td></tr><tr><td rowspan="4">Température moyenne</td><td>OG</td><td>0,3832</td><td>1997</td></tr><tr><td>AS</td><td>0,3648</td><td>1994</td></tr><tr><td>BAM</td><td>0,3421</td><td>1994</td></tr><tr><td>FBS</td><td>0,4704</td><td>1994</td></tr><tr><td rowspan="4">Evapo-<br>Transpiration (ETo)</td><td>OG</td><td>-</td><td>-</td></tr><tr><td>AS</td><td>-</td><td>-</td></tr><tr><td>BAM</td><td>-</td><td>-</td></tr><tr><td>FBS</td><td>-</td><td>-</td></tr><tr><td rowspan="4">Indice de De<br>Maritonne</td><td>OG</td><td>0,1540</td><td>1978</td></tr><tr><td>AS</td><td>0,1220</td><td>1979</td></tr><tr><td>BAM</td><td>0,2831</td><td>1979</td></tr><tr><td>FBS</td><td>0,1442</td><td>1979</td></tr></tbody></table>

Il a été constaté également par ce test des clivages dans la chronologie des séries périodiques pour les indicateurs de températures minimales, maximales et moyennes, respectivement aux gammes annuelles suivantes: 1977-1997 ; 1997-1999 ; 1994-1997 ; 1978-1979. Concernant l'ETo, le test n'a révélé par contre aucune telle rupture (Tableau 69).

### d) Segmentation de Hubert

Les principaux résultats des calculs par méthode sont consignés dans le tableau suivant.
<table><caption>Tableau 70 : Segmentation de Hubert appliquée sur les séries temporelles agro-climatiques au périmètre du Tadla
(Période d'observations de 1970 à 2018)</caption><thead><tr><th rowspan="2">Indicateur ou indice Climatique</th><th rowspan="2">Station météorologique</th><th rowspan="2">Dates de rupture</th><th colspan="6">Paramètres et résultats de la segmentation de Hubert</th></tr><tr><th>Nombre rupture (NP)</th><th>Avant rupture Valeur moyenne (1)</th><th>Après rupture Valeur moyenne (21)</th><th>Ecart entre après et avant rupture [ (2) - (1) ]</th></tr></thead><tbody><tr><td rowspan="3">Température minimale (en mm)</td><td>OG</td><td>1977<br>2000</td><td>2</td><td>10,38<br>11,17</td><td>0,39<br>0,52</td><td>11,17<br>11,81</td><td>0,52<br>0,64</td><td>0,79<br>0,64</td></tr><tr><td>AS</td><td>1977</td><td>1</td><td>10,93</td><td>0,37</td><td>11,60</td><td>0,58</td><td>0,67</td></tr><tr><td>BAM</td><td>1989<br>1977</td><td>1<br>1</td><td>11,13<br>11,53</td><td>0,55<br>0,42</td><td>11,91<br>12,33</td><td>0,57<br>0,59</td><td>0,78<br>0,80</td></tr><tr><td rowspan="3">Température maximale (en mm)</td><td>OG</td><td>1993<br>2009</td><td>2</td><td>25,95<br>26,84</td><td>0,54<br>0,81</td><td>26,84<br>27,76</td><td>0,81<br>0,69</td><td>0,89<br>0,92</td></tr><tr><td>AS</td><td>1996<br>1999</td><td>1<br>1</td><td>25,11<br>26,17</td><td>0,55<br>0,56</td><td>27,48<br>27,61</td><td>0,63<br>0,74</td><td>1,97<br>1,44</td></tr><tr><td>FBS</td><td>1997</td><td>1</td><td>26,55</td><td>0,58</td><td>28,28</td><td>0,61</td><td>1,73</td></tr><tr><td rowspan="3">Température moyenne (en mm)</td><td>OG</td><td>1978<br>1997</td><td>2</td><td>18,06<br>18,66</td><td>0,31<br>0,48</td><td>18,66<br>19,54</td><td>0,48<br>0,44</td><td>0,60<br>0,88</td></tr><tr><td>AS</td><td>1978<br>1997</td><td>2</td><td>18,12<br>18,69</td><td>0,32<br>0,50</td><td>18,69<br>19,48</td><td>0,50<br>0,41</td><td>0,57<br>0,79</td></tr><tr><td>BAM</td><td>1978<br>2000</td><td>2</td><td>18,19<br>18,84</td><td>0,30<br>0,49</td><td>18,84<br>19,81</td><td>0,49<br>0,49</td><td>0,65<br>0,97</td></tr><tr><td rowspan="3">Evapo-Transpiration (en mm/jour)</td><td>FBS</td><td>1978<br>1997</td><td>2</td><td>18,59<br>19,21</td><td>0,34<br>0,50</td><td>19,21<br>20,31</td><td>0,50<br>0,41</td><td>0,62<br>1,10</td></tr><tr><td>OG</td><td>1973</td><td>1</td><td>2,93</td><td>0,05</td><td>3,62</td><td>0,14</td><td>0,69</td></tr><tr><td>AS</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Indice de De Maritonne</td><td>FBS</td><td>1973<br>2012</td><td>2</td><td>2,90<br>3,61</td><td>0,00<br>0,29</td><td>3,61<br>4,30</td><td>0,29<br>0,61</td><td>0,71<br>0,69</td></tr><tr><td>OG</td><td>1971</td><td>1</td><td>17,200</td><td>2,828</td><td>10,422</td><td>2,792</td><td>6,778</td></tr><tr><td>AS</td><td>1991<br>2013</td><td>2</td><td>16,200<br>9,633</td><td>2,263<br>2,764</td><td>9,633<br>5,300</td><td>2,764<br>2,978</td><td>6,567<br>4,333</td></tr><tr><td rowspan="3"></td><td>FBS</td><td>1974<br>2010</td><td>2</td><td>15,400<br>8,717</td><td>3,201<br>2,696</td><td>8,717<br>5,231</td><td>2,696<br>2,404</td><td>6,683<br>3,486</td></tr><tr><td></td><td>1979</td><td>1</td><td>12,150</td><td>2,701</td><td>8,462</td><td>3,140</td><td>3,688</td></tr></tbody></table>

Il est repéré par cette méthode de fréquentes segmentations statistiques sur les séries chronologiques agro-climatiques. En effet, il y a au moins deux ruptures successives pour chaque paramètre étudié, avec le caractère répétitif de l'espacement dans le temps de 24 à 33 ans. Ces coupures temporaires sont exprimées aux bornes des intervalles des indicateurs des températures minimales, maximales et moyennes, de l'ETo et des valeurs de l'Indice de De Maritonne respectivement en 1977-2000, 1993-2009, 1978-2000, 1973-2012 et en 1971-2013 (Tableau 70).

## II-2-3-Determination des dates optimales de ruptures des séries chronologiques thermique et de demande climatique en évaporation

A l'issue des résultats des tests précités d'homogénéité, sont déterminées graphiquement les dates optimales des ruptures des séries chronologiques correspondant aux paramètres ainsi : températures minimale en 1989, température maximale en 1997, température moyenne en 1995 et ETo en 1978; comme sont illustrées aux graphes de la Figures 25. Pour l'Indice de De Maritonne, via de méthodologie similaire, l'année de rupture optimale est en 1988.
<figure><img src="image_27.png" /><figcaption>Graphique A : Température maximale</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>Graphique B : Température minimale</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>Graphique C : Température moyenne</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>Graphique D : ETo</figcaption></figure>

<caption>Figures 25 : Graphiques des ruptures temporelles pour paramètres : Température maximale (graphique A), Température minimale (graphique B), Température moyenne (graphique C), ETo (graphique D) au périmètre de Tadla</caption>

(Période d'observations de 1970 à 2018).

### III-Cractérisation de tendance thermique et de la demande climatique en évaporation pour le périmètre de Tadla

A la lumière des précédents tests d'homogénéité réalisés et après avoir pris connaissance des dates butoirs des ruptures chronologiques sur les séries de données de paramètres agro-climatiques étudiés, il est présenté ci-après leurs évolutions sous forme de graphiques à histogrammes et sur des cartes SIG en des thématiques spatio-temporelles.
### III-1-Tendances thermiques au périmètre de Tadla

Les analyses des tendances des paramètres thermiques de températures minimale, maximale et moyenne, sont relatées plus bas.

#### III-1-1-Températures minimales

La température minimale moyenne du périmètre de Tadla a évolué de part et d'autre de la date de la rupture en année 1989 de la valeur de 11,3 à 11,9°C ; soit un réchauffement du milieu de 0,6°C (4%), constat et qui perdure en période de 1990 à 2018 (Figure 26).

<figure><img src="image_31.png" /><figcaption>Figures 26 : Tendance de la température minimale au périmètre du Tadla en avant et en après la rupture chronologique</figcaption></figure>

(Période d'observations de 1970 à 2018)

De point de vue tendance spatio-temporelle à l'échelle de ce périmètre, il se profile en après la date de rupture, un contraste thermique en sa cotation minimale, qui est relativement plus chaud sur les terrains de la zone bour du Nord du périmètre, comparé à la zone Sud en piedmont à proximité de la chaine des montagnes du Moyen Atlas (voir les cartes de la Figure 27).
<figure><img src="image_32.png" /><figcaption>Carte A : Température minimale en situation d'avant la rupture chronologique</figcaption></figure>

<figure><img src="image_33.png" /><figcaption>Carte B : Température minimale en situation en après la rupture chronologique</figcaption></figure>

Figure 27: Cartes de tendance spatio-temporelle de la température minimale en avant (Carte A) et en après rupture chronologique (Carte B) au périmètre de Tadla (Période d'observations de 1970 à 2018).
### III-1-2-Températures maximales

La température maximale du périmètre est passée après la rupture de l'année 1997 de 26,1 à 27,6 °C ; soit un important relatif réchauffement du milieu ambiant qui est évalué à 1,5 °C (6%). Le fait plausible qui se maintient entre sur dernière période de 1998 et 2018 (Figure 28).

<figure><img src="image_34.png" /><figcaption>Figures 28 : Tendances de la température maximale au périmètre du Tadla en avant et en après la rupture chronologique
(Période d'observations de 1970 à 2018).</figcaption></figure>

Selon la vision spatio-temporelle à l'échelle du périmètre de Tadla, la tendance après la rupture en question (1997) est constatée d'une manière presque homogène dans cet espace. Toutefois une nuance relative est signalée pour les trois zones suivantes : au Centre-Nord (partie de la zone bour), au Sud-Ouest (partie de la zone irriguée) et au Sud-Est (grande partie de la zone Dir), avec respectivement à caractère relativement chaud, moins chaud et faiblement chaud (voir les cartes de la Figure 29).
<figure><img src="image_35.png" /><figcaption>Carte A : Température maximale en situation d'avant la rupture chronologique</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>Carte B : Température maximale en situation en après la rupture chronologique</figcaption></figure>

Figure 29: Cartes de tendance spatio-temporelle de la température maximale en avant (Carte A) et en après rupture chronologique (Carte B) au périmètre de Tadla (Période d'observations de 1970 à 2018).
### III-1-3-Températures moyennes

La température moyenne du périmètre a sursauté après sa rupture d'évolution en année de 1995 de 18,7 à 19,7 °C ; mettant aussi en évidence le réchauffement moyen du milieu de 1 °C (5%). La circonstance qui a sévie durant les deux dernières décennies (1996-2018) (Figure 30).

1970

<figure><img src="image_37.png" /><figcaption>Figures 30 : Tendances de la température moyenne au périmètre du Tadla en avant et en après rupture chronologique</figcaption></figure>

(Période observation de 1970 à 2018)

La cartographie de la température moyenne révèle que la tendance du réchauffement du milieu touche l'ensemble du territoire de ce périmètre, mais en de différentes ampleurs. L'accroissement thermique est relativement plus accentué en zone Sud-Ouest, laquelle abrite une partie relativement importante de la zone potentielle irriguée de ce périmètre (voir les cartes de la Figure 31).
<figure><img src="image_38.png" /><figcaption>Carte A : Température moyenne en situation d'avant rupture chronologique</figcaption></figure>

<figure><img src="image_39.png" /><figcaption>Carte B : Température moyenne en situation en après rupture chronologique</figcaption></figure>

Figure 31: Cartes de tendance spatio-temporelle de la température moyenne en avant (Carte A) et en après rupture chronologique (Carte B) au périmètre de Tadla (Période d'observations de 1970 à 2018).
Dans l'objectif de développer d'éventuelles analyses thermiques implicites en relation avec la pluviométrie mensuelle dans le périmètre de Tadla, selon une cadence mensuelle, des courbes de régression des précipitations mensuelles en fonction du temps ont été élaborées pour la période en après la rupture chronologique de la température moyenne (année de scission chronologique : 1988) (Figures 32 et 33).

<figure><img src="image_40.png" /><figcaption>Figure 32 : Courbes de régression de relation « Année-Pluviométrie mensuelle » au périmètre de Tadla définie par mois de Septembre à Février
(Période : 1989-2018 ; après rupture chronologique d'ETo).</figcaption></figure>
Il est à signaler que les courbes des tendances en question suivent le meilleur modèle sélectionné, qui est de type polynomial, avec des coefficients de la corrélation en des niveaux variant entre la moyenne et la bonne acceptabilité (Figures 32 et 33).

<figure><img src="image_41.png" /><figcaption>Figure 33 : Courbes de régression de relation « Année-Pluviométrie mensuelle » au périmètre du Tadla définie par mois de Mars à Août
(Période : 1989-2018 ; après rupture chronologique d'ETo).</figcaption></figure>
### III-2-Tendance en demande climatique en évaporation au périmètre du Tadla

#### III-2-1-Evolution de tendance d'évapo-transpiration de référence au périmètre de Tadla

Le paramètre d'évapo-transpiration de référence, l'indicateur renseignant sur la demande climatique en évaporation d'un sol enherbé, a subi précocement sa rupture chronologique aux séries des mesures en l'année 1978, coïncidant avec la date de la rupture pluviométrique. L'ETo évolue après cette rupture de 3,2 à 3,6 mm/jour ; soit pour un écart positif de 0,4 mm/jour (13%) (Figure 34).

<figure><img src="image_42.png" /><figcaption>Figure 34 : Tendances de l'évapo-transpiration de référence au périmètre du Tadla en avant et en après rupture chronologique
(Période d'observations de 1970 à 2018)</figcaption></figure>

#### III-2-2-Cartographie de tendance de l'évapo-transpiration de référence au périmètre de Tadla

La carte SIG de l'évapo-transpiration de référence illustre la hausse relative en valeur moyenne de ce paramètre auprès de la zone Sud-Ouest, suivi de la large zone du Centre avec des valeurs intermédiaires et enfin en zone d'extrême Est du périmètre à des valeurs relativement les plus basses. Après rupture chronologique, on assiste à la persistance des mêmes niveaux de la répartition spatiale, ceci en dépit de l'accroissement constaté de l'ETo (voir cartes de la Figure 35).
<figure><img src="image_43.png" /><figcaption>Carte A : Evapo-transpiration en situation d'avant rupture chronologique</figcaption></figure>

<figure><img src="image_44.png" /><figcaption>Carte B : Evapo-transpiration en situation en après rupture chronologique</figcaption></figure>

Figure 35: Cartes de tendance spatio-temporelle de l'évapo-transpiration de référence en avant (Carte A) et en après rupture chronologique (Carte B) au périmètre de Tadla (Période d'observations de 1970 à 2018).

### III-3-Relation entre la température et l'évapo-transpiration au périmètre du Tadla

L'expression de l'évapo-transpiration de référence en fonction par la température à travers des courbes-abaques par mois, en se référant aux données de la station météorologique d'Ouled Gnaou sur le périmètre de Tadla en après-rupture chronologique thermique moyenne (en 1989) est donnée aux graphiques des (Figures 36 et 37).
<figure><img src="image_45.png" /></figure>

<figure><img src="image_46.png" /></figure>

<figure><img src="image_47.png" /></figure>

<figure><img src="image_48.png" /></figure>

<figure><img src="image_49.png" /></figure>

<figure><img src="image_50.png" /><figcaption>Figure 36 : Courbes de régression de relation Température-Evapotranspiration au périmètre de Tadla pour les mois de Septembre à Février</figcaption></figure>

(Période : 1989-2018 ; après rupture chronologique d'ETo).

Il est à noter que la relation entre de l'évapo-transpiration de référence avec la température est selon la méthode régressive comporte des globalement des coefficients de corrélation dans les niveaux requis (Figures 36 et 37).
<figure><img src="image_51.png" /><figcaption>Figure 37 : Courbes de régression de relation Température-Evapotranspiration au périmètre de Tadla pour les mois de Mars à Août</figcaption></figure>

(Période : 1989-2018 ; après rupture chronologique d'ETo).

## IV-Bilan hydrique au sol au périmètre de Tadla

Le bilan hydrique au sol au périmètre de Tadla ayant été déterminé selon une cadence mensuelle par la différence entre l'ETo et les précipitations efficaces, a été doublement établi par mois et d'une manière cumulative au cours de l'année et pour les deux situations d'avant et en après-rupture chronologique des paramètres en question.
Les résultats afférents des analyses sont présentés et commentés ci-après avec illustration par des graphiques.

## IV-1-Bilan hydrique mensuel au sol au périmètre de Tadla

### IV-1-Evaluation du bilan hydrique au sol mensuel global au périmètre de Tadla en avant et en après la rupture chronologie de l'évapo-transpiration de référence.

En tenant compte seulement des précipitations interceptées au périmètre de Tadla, le bilan hydrique au sol qui concerne la végétation du milieu à caractère pluvieux, est constamment négatif pour tous les mois de la campagne. Ceci avait lieu à différentes amplitudes des mois et qui ne dépassent guère la barre mensuelle de -150 mm/mois avant la rupture chronologique de l'année 1978.

Après la rupture chronologique (précipitations et ETo), le déficit hydrique au sol se creuse davantage à partir du mois de Mars, avec des dépassements supérieurs par rapport au seuil de l'avant rupture de - 150 mm/mois, notamment pour les mois de la saison estivale en la phase de l'après rupture (voir graphes de la Figure 38).

### IV-1-2-Bilan hydrique cumulatif annuel caractéristique du périmètre de Tadla

En raisonnant selon le bilan hydrique au sol cumulatif annuel sur base mensuelle, la disparité entre les situations d'avant et en après la rupture chronologique de l'année 1978 (précipitations et ETo), la situation devient plus délicate dès le mois de Février en continuant à s'accentuer au fil des mois restant de la campagne agricole (Figure 39).
<figure><img src="image_52.png" /><figcaption>Graphe A : Bilans hydriques au sol mensuels avant-rupture chronologique de la pluie et de l'ETo</figcaption></figure>

<figure><img src="image_53.png" /><figcaption>Graphe B : Bilans hydriques au sol mensuels après-rupture chronologique de la pluie et de l'ETo</figcaption></figure>

Figure 38 : Tendance des bilans hydriques au sol mensuels le périmètre de Tadla en avant-rupture (graphe A) et en après-rupture (graphe B) chronologique de la pluie et de l'ETo (Période : 1970-2018).

A l'achèvement de la campagne agricole, le bilan hydrique au sol cumule une valeur déficitaire de - 805 mm/an (équivalent à -8 050 m3 d'eau par hectare et par an), en période d'avant-rupture chronologique (en l'année 1978), contre un volume aussi déficitaire mais plus important et qui
est évalué à - 1 093 mm/an (équivaut à -10 930 m3/ha/an) en période d'après-rupture ; soit alors
un taux d'accroissement toujours dans le sens de la négativité de 35,78 % (Figure 39).

<figure><img src="image_54.png" /><figcaption>Figure 39 : Tendance de la variabilité du bilan hydrique au sol cumulatif d'une manière comparative entre l'avant et l'après rupture chronologique des paramètres agro-climatiques au périmètre de Tadla

(Période d'observations de 1970 à 2018).</figcaption></figure>
# Chapitre IV : VARIABILITE EN NIVEAU D'ARIDITE AU PERIMETRE DU TADLA

## I-Diagrammes de Gaussen

L'élaboration des diagrammes de Gaussen pour les stations météorologiques retenues pour le périmètre du Tadla avait pour assise les trois périodes interposées entre les années successives des ruptures chronologiques. La première période est de 1970-1978, spécifique seulement pour le paramètre des précipitations (rupture en 1978). La deuxième période est de 1979-1995 en intermédiaire à propos du couple de paramètres précipitations-température moyenne. La troisième période est de 1996-2018 et spécifie uniquement à la température moyenne (rupture en 1995). Les résultats sont récapitulés au Tableau 71.

Ces aspects sont analysés ci-après et illustrés aux Figures 43 et 44 pour les différents diagrammes établis pour les périodes desdites ruptures et en fonction des stations de mesures.

Les durées des périodes sèches en fonction des trois périodes des ruptures chronologiques précitées, extraites à partir des Diagrammes de Gaussen réalisés sont rapportées au Tableau 71.

<table><caption>Tableau 71 : Tendances de la variabilité de la période sèche sur le diagramme de Gaussen au périmètre de Tadla
(Période de 1970 à 2018)</caption><thead><tr><td>Intervalle annuel séparant les ruptures chronologiques des indicateurs agro-climatiques</td><td>Durée de la période sèche (en nombre de jours)</td></tr></thead><tbody><tr><td>1970-1978</td><td>221 (+/- 43)</td></tr><tr><td>1979-1995</td><td>192 (+/- 15)</td></tr><tr><td>1996-2018</td><td>235 (+/- 28)</td></tr></tbody></table>

Il se profile de ces résultats qu'en passant de la dernière date de la rupture chronologique des indicateurs agro-climatiques, coïncidant avec celle de la température moyenne du périmètre (année de 1995), que la période sèche a significativement évolué de **192 à 235 jours**, soit un allongement de **43 jours** (Tableau 71 ; Figures 40 et 41).
<figure><img src="image_55.png" /><figcaption>Graphique A : Période de 1970 à 1978</figcaption></figure>

<figure><img src="image_56.png" /><figcaption>Graphique B : Période de 1979 à 1995</figcaption></figure>

<figure><img src="image_57.png" /><figcaption>Graphique C : Période de 1996 à 2018
STATION D'OULED GNOU</figcaption></figure>

<figure><img src="image_58.png" /><figcaption>Graphique D : Période de 1970 à 1978</figcaption></figure>

<figure><img src="image_59.png" /><figcaption>Graphique E : Période de 1979 à 1995</figcaption></figure>

<figure><img src="image_60.png" /><figcaption>Graphique F : Période de 1996 à 2018
STATION AHL SOUSS</figcaption></figure>

Figure 40 : Diagrammes de Gaussen au périmètre de Tadla aux niveaux des stations agro-climatiques Ouled Gnou (graphes A à C) et Ahl Souss (graphes D à F) (Période d'observations de 1970 à 2018).
En effet, il se traduit fidèlement et d'une manière conjointe le réchauffement atmosphérique tendanciel avec la tendance de la réduction graduelle des précipitations, surtout en fin de la saison printanière pour le périmètre de Tadla (Figures 40 et 41).

<figure><img src="image_61.png" /><figcaption>Graphique A : Période de 1970 à 1978</figcaption></figure>

<figure><img src="image_62.png" /><figcaption>Graphique B : Période de 1979 à 1995</figcaption></figure>

<figure><img src="image_63.png" /><figcaption>Graphique C : Période de 1996 à 2018</figcaption></figure>

STATION MABROUK

<figure><img src="image_64.png" /><figcaption>Graphique D : Période de 1970 à 1978</figcaption></figure>

<figure><img src="image_65.png" /><figcaption>Graphique E : Période de 1979 à 1995</figcaption></figure>

<figure><img src="image_66.png" /><figcaption>Graphique E : Période de 1996 à 2018</figcaption></figure>

STATION FKIH BEN SALAH

<caption>Figure 41 : Diagrammes de Gaussen au périmètre de Tadla aux niveaux des stations agro-climatiques Mabrouk (graphes A à C) et Fkih Ben Salah (graphes D à F) (Période d'observations de 1970 à 2018).</caption>
## II-Indice de De Maritonne

Les résultats numériques de l'évolution de l'indice d'aridité de De Maritonne pour le périmètre de Tadla sont restitués au Tableau 72.

<table><caption>Tableau 72 : Tendance de la variabilité de l'indice d'aridité de De Maritonne au périmètre de Tadla
(Période de 1970 à 2018).</caption><thead><tr><th>Intervalle annuel séparant les ruptures chronologiques des indicateurs agro-climatiques</th><th>Indice moyen de Maritonne I.M. (en mm/°C)</th></tr></thead><tbody><tr><td>1970-1988</td><td>11 (+/- 1)</td></tr><tr><td>1989-2018</td><td>9 (+/- 1)</td></tr></tbody></table>

La valeur moyenne de l'indice d'aridité de De Maritonne pour le périmètre de Tadla passe entre l'avant à l'après-rupture chronologique (année de 1988) de 11,2 à 8,6 mm/°C (Tableau 72). Ceci est à l'origine de l'accentuation d'aptitude à l'aridité du milieu lors de l'intervalle mensuel de Janvier à Mars, comme est illustré dans la Figure 42.

<figure><img src="image_67.png" /><figcaption>Figure 42: Tendance de la variabilité de l'indice mensuel d'aridité de De Maritonne au périmètre de Tadla
(Période d'observations de 1970 à 2018).</figcaption></figure>

De point de vue spatio-temporelle, après le clivage dans la chronologie de tendance de l'indice de De Maritonne, repérée en l'année 1988, l'aridité a fortement augmenté à différents endroits du périmètre du Tadla, avec l'accélération dans sa zone Sud-Est, dont les cartes thématiques sont sur la Figure 43.
<figure><img src="image_68.png" /><figcaption>Carte A : Niveaux d'aridités en situation d'avant rupture chronologique</figcaption></figure>

<figure><img src="image_69.png" /><figcaption>Carte B : Niveaux d'aridités en situation en situation en après rupture chronologique</figcaption></figure>

Figure 43 : Cartes d'évolution spatio-temporelle de l'aridité au périmètre du Tadla entre l'avant (Carte A) et l'après (Carte B) de la rupture chronologique sur paramètres agro-climatique

(Période d'observations de 1970 à 2018).
### III-Quotient d'Emberger

Les résultats des calculs de l'évolution du Quotient d'Emberger (Q.E.) pour le périmètre de Tadla, sont données au Tableau 73.

<table><caption>Tableau 73 : Tendances de la variabilité de l'indice mensuel d'aridité de De Maritonne au périmètre de Tadla
(Période de 1970 à 2018).</caption><thead><tr><th>Intervalle annuel séparant les ruptures chronologiques des indicateurs agro-climatiques</th><th>Quotient d'Emberger (Q.E.) en mm/ (°C) x(°C)</th></tr></thead><tbody><tr><td>1970-1978</td><td>28 (+/- 5)</td></tr><tr><td>1979-1995</td><td>21 (+/- 4)</td></tr><tr><td>1996-2018</td><td>19 (+/- 3)</td></tr></tbody></table>

Le périmètre du Tadla se caractérisait en avant la rupture chronologique en question avec un Q.E de 28 mm/°C². Ce quotient s'est replié en fin de compte après deux scissions temporaires pour un écart négatif de 11 unités à l'au-delà de l'année de fatidique de 1995 (Tableau 73).
# Chapitre V : CHANGEMENTS CLIMATIQUES ET LEURS IMPACTS SUR LES BESOINS EN EAU DE L'IRRIGATION DES CULTURES AU PERIMETRE DE TADLA

Au préalable des projections d'évolution des besoins en eau d'irrigation des cultures du périmètre du Tadla sous les impacts éventuels des changements climatiques prévisionnels, il est objet de déterminer leurs exigences actuelles en eau d'irrigation selon les systèmes d'irrigation usuels adoptés. A ce propos, il est présenté les résultats des enquêtes agricoles réalisées, suivi des évaluations par les calculs et à travers les expérimentations pour la détermination des besoins nets et bruts en eau d'irrigation des cultures. Ceux-ci sont relayés par les extrapolations au périmètre de Tadla pour les besoins en eau des assolements culturaux actuels et projetés à long terme par systèmes d'irrigation adoptés, tout en tenant compte de l'éventualité des changements climatiques.

## I-Résultats des enquêtes agricoles auprès des agriculteurs sur l'usage comparatif des deux systèmes d'irrigation gravitaire et localisée.

A l'issue des deux enquêtes complémentaires de caractérisation du processus d'irrigation des cultures menées au périmètre de Tadla, auprès d'agriculteurs superviseurs directs de leurs exploitations agricoles, des renseignements ont été tirés et sont récapitulés au Tableau 74.

<caption>Tableau 74 : Principaux résultats des enquêtes sur l'opération de l'irrigation auprès d'exploitations agricoles du périmètre de Tadla avec la comparaison des systèmes d'irrigation gravitaire et localisée</caption>

(Campagnes agricoles des enquêtes de 2015/2016 à 2017/2018)

<table><thead><tr><th rowspan="2">Cultures</th><th>Volume eau<br>irrigation</th><th>Rendement</th><th>Efficience<br>agronomique<br>d'eau (EUE)</th><th>Valeur de<br>production</th><th>Marge<br>nette</th><th>Valeur<br>ajoutée</th><th>Valorisation<br>eau irrigation</th></tr><tr><th>m3 / ha</th><th>T / ha</th><th>Kg / m3 eau</th><th>Dhs / ha</th><th>Dhs / ha</th><th>Dhs / ha</th><th>Dhs / m3 eau</th></tr></thead><tbody><tr><td><strong>Céréales d'automne</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>6 720</td><td>5,40</td><td>0,80</td><td>19 900</td><td>5 586</td><td>14 440</td><td>2,15</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>3 450</td><td>6,58</td><td>1,91</td><td>23 030</td><td>6 618</td><td>15 490</td><td>4,49</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-49%</strong></td><td><strong>22%</strong></td><td><strong>137%</strong></td><td><strong>16%</strong></td><td><strong>18%</strong></td><td><strong>7%</strong></td><td><strong>109%</strong></td></tr><tr><td><strong>Betterave à Sucre</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>8 950</td><td>66,00</td><td>7,37</td><td>31 189</td><td>10 229</td><td>24 049</td><td>2,69</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>5 170</td><td>86,08</td><td>16,65</td><td>36 765</td><td>12 535</td><td>27 845</td><td>5,39</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-42%</strong></td><td><strong>30%</strong></td><td><strong>126%</strong></td><td><strong>18%</strong></td><td><strong>23%</strong></td><td><strong>16%</strong></td><td><strong>100%</strong></td></tr><tr><td><strong>Luzerne</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>11 520</td><td>55,00</td><td>4,77</td><td>32 140</td><td>19 040</td><td>30 200</td><td>2,62</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>8 940</td><td>65,00</td><td>7,27</td><td>37 984</td><td>23 614</td><td>34 602</td><td>3,87</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-22%</strong></td><td><strong>18%</strong></td><td><strong>52%</strong></td><td><strong>18%</strong></td><td><strong>24%</strong></td><td><strong>15%</strong></td><td><strong>48%</strong></td></tr><tr><td><strong>Mais Ensilage</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>7 000</td><td>45,00</td><td>6,43</td><td>18 279</td><td>4 639</td><td>6 759</td><td>0,97</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>5 370</td><td>55,39</td><td>10,31</td><td>22 500</td><td>8 204</td><td>10 524</td><td>1,96</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-23%</strong></td><td><strong>23%</strong></td><td><strong>60%</strong></td><td><strong>23%</strong></td><td><strong>77%</strong></td><td><strong>56%</strong></td><td><strong>103%</strong></td></tr><tr><td><strong>Agrumes</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>12 000</td><td>27,00</td><td>2,25</td><td>53 859</td><td>26 344</td><td>49 139</td><td>4,09</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>8 000</td><td>42,00</td><td>5,25</td><td>84 450</td><td>45 038</td><td>79 205</td><td>9,90</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-33%</strong></td><td><strong>56%</strong></td><td><strong>133%</strong></td><td><strong>57%</strong></td><td><strong>71%</strong></td><td><strong>61%</strong></td><td><strong>142%</strong></td></tr><tr><td><strong>Oliviers</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>7 680</td><td>4,20</td><td>0,55</td><td>18 840</td><td>8 606</td><td>16 804</td><td>2,19</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>5 450</td><td>8,00</td><td>1,47</td><td>40 340</td><td>13 638</td><td>37 599</td><td>6,90</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-29%</strong></td><td><strong>90%</strong></td><td><strong>168%</strong></td><td><strong>114%</strong></td><td><strong>58%</strong></td><td><strong>124%</strong></td><td><strong>215%</strong></td></tr><tr><td><strong>Oignon</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>6 000</td><td>57,50</td><td>9,58</td><td>28 750</td><td>9 032</td><td>24 150</td><td>4,03</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>5 800</td><td>83,57</td><td>14,41</td><td>33 428</td><td>11 806</td><td>27 998</td><td>4,83</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-3%</strong></td><td><strong>45%</strong></td><td><strong>50%</strong></td><td><strong>16%</strong></td><td><strong>31%</strong></td><td><strong>16%</strong></td><td><strong>20%</strong></td></tr><tr><td><strong>Melon</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Syst. Gravitaire (1)</td><td>6 720</td><td>32,50</td><td>4,84</td><td>37 000</td><td>16 000</td><td>32 850</td><td>4,89</td></tr><tr><td>Syst. Irr. Localisée (2)</td><td>4 830</td><td>52,70</td><td>10,91</td><td>78 750</td><td>42 010</td><td>56 510</td><td>11,70</td></tr><tr><td><strong>Ecrat : (2)/(1)</strong></td><td><strong>-28%</strong></td><td><strong>62%</strong></td><td><strong>126%</strong></td><td><strong>113%</strong></td><td><strong>163%</strong></td><td><strong>72%</strong></td><td><strong>139%</strong></td></tr></tbody></table>
Les analyses des données de ces enquêtes sur milieu réel de l'agriculture du périmètre de Tadla,
sont données ci-après, modulées en fonction des aspects directeurs phyto-techniques et agro-
économiques.

* **Consommation en eau d'irrigation par les cultures**

-Les cultures en irrigué bénéficient d'apports d'eau par les agriculteurs sous le système d'irrigation gravitaire, qui oscillent entre un minimum pour le maraichage de 6 000 m3/ha en oignon et 6 720 m3/ha en melon à niveau maximum pour les cultures pérennes d'agrumes à 12 000 m3/ha) et pour la luzerne à 11 520 m3/ha.

-Le système d'irrigation localisée contribue effectivement à économiser l'eau d'irrigation dans un taux variant entre 3% à 49 % selon les cultures.

* Productivités des cultures

Le recours au système d'irrigation localisée par les agriculteurs leur a permis d'améliorer les rendements de leurs cultures, avec un taux de gain allant de 18 % pour la luzerne à 90 % pour les oliviers. Aussi sur le plan productif, les plantations d'oliviers se trouvent désormais mieux valorisées avec la nouvelle technique d'irrigation plus efficiente.

• **Efficiences agronomique de l’utilisation d’eau d’irrigation**

-Les efficiences agronomiques de l'utilisation de l'eau d'irrigation (EUE) au niveau des exploitations agricoles du périmètre de Tadla montrent des disparités pour les cultures. En effet, l'EUE passe de 0,80 Kg/m3 d'eau en cultures des céréales d'automne à 9,58 Kg/m3 pour la culture maraichère de référence d'oignon d'hivers.

-A l'exception des cultures céréalières, les autres cultures affichent des valeurs d'EUE plus élevées par rapport aux résultats des expérimentations, ce compte-tenu que les agriculteurs utilisent massivement de l'eau du pompage pour le complément d'irrigation en plus de l'eau du barrage fourni à une stricte mesure. Cette attitude de l'agriculteur est dans le sens de maximaliser leurs rendements de campagne surtout pour la culture industrielle de la betterave à sucre, contractée avec l'unique usine sucrière locale.

-L'utilisation de l'irrigation localisée chez les agriculteurs s'est montrée plus bénéfique dans l'accroissement de l'EUE des cultures.

* Valeurs des productions agricoles

En milieu irrigué, les recettes de vente des productions des cultures sont évaluées dans la fourchette de e 18 840 à 37 000 Dhs/ha sous l'irrigation gravitaire. Les accroissements financiers
productifs après la reconversion au système d'irrigation localisée sont à des taux fluctuant entre 16% pour la céréaliculture et au multiple de 1,14 fois pour l'olivier et à 1,13 fois la culture maraichère de référence du melon.

* **Marge nette agricole**

Le bénéfice net réalisé à l'échelle de l'exploitation agricole en irrigué évolue du minimum de 4 639 Dhs/ha pour la culture du maïs ensilage, ayant un court cycle productif de 4 à 5 mois,) à la valeur au maximum de 26 344 Dhs/ha pour les plantations agrumicoles. Le système d'irrigation localisée, intervient dans l'accroissement de la marge nette agricole dans une gamme du taux de 16% pour le melon et de 18 % pour les céréales à élévations maximales en taux de 71% pour les agrumes et de 77% pour le maïs d'ensilage.

* **Valeur ajoutée agricole**

Les cultures conduites sous l'irrigation gravitaire dégagent des valeurs ajoutées dans un intervalle de 6 759 Dhs/ha pour le maïs d'ensilage à 49 139 Dhs/ha en agrumiculture. Par ailleurs, l'irrigation localisée s'implique à augmenter davantage le taux de gain de 7% pour la culture de blés et à le multiplier par facteur de 1,24 fois pour les vergers des oliviers.

* **Valorisation économique de l'eau d'irrigation**

L'eau d'irrigation se valorise économiquement sous le système gravitaire en son plus bas niveau de 0,87 Dhs par le mètre cube d'eau en présence du maïs ensilage en passant à des niveaux plus élevés pour les cultures du melon, des agrumes et du maraichage (oignon) respectivement aux valeurs de 4,89 ; 4,09 et 4,03 Dhs/m3. La valorisation économique de l'eau est améliorée pour toutes les cultures en intégrant l'irrigation localisée, en réalisant un surplus par rapport à la technique du gravitaire de 20% pour le maraichage (oignon) jusqu'au multiple de 2,15 fois pour la spéculation oléicole.

## II-Evaluation des besoins nets en eau actuels des cultures du périmètre de Tadla.

Les données des besoins nets en eau qui ont été déterminés pour les principales cultures du périmètre de Tadla par systèmes d'irrigation gravitaire et localisée, sont présentées d'une manière détaillée sur des canevas des calculs et les récapitulatifs au Tableau 75.
<caption>Tableau 75: Besoins nets en eau d'irrigation des principales cultures du périmètre de Tadla</caption>

(Campagnes de référence : 2014/2015 à 20117/2018)

<table><thead><tr><th>Culture</th><th>Système d'irrigation</th><th>Campagne de référence</th><th>Rang annuel</th><th>Besoins nets en eau d'irrigation (mm)</th><th>Ecarts entre système d'irrigation localisée (2)/gravitaire (1)</th></tr></thead><tbody><tr><td rowspan="2">Céréales d'automne (Blés)</td><td>Gravitaire (1)</td><td>2016/2017</td><td>A3</td><td>445,30</td><td rowspan="2">-12%</td></tr><tr><td>Localisée (2)</td><td>2016/2017</td><td>A3</td><td>391,20</td></tr><tr><td rowspan="4">Betterva à sucre</td><td rowspan="2">Gravitaire</td><td>2014/2015</td><td>A1</td><td>575,83</td><td rowspan="4">-8%</td></tr><tr><td>2015/2016</td><td>A2</td><td>514,34</td></tr><tr><td rowspan="2">Localisée</td><td>2014/2015</td><td>A1</td><td>491,82</td></tr><tr><td>2015/2016</td><td>A2</td><td>466,38</td></tr><tr><td rowspan="2">Luzerne</td><td>Gravitaire</td><td>2016/2017</td><td>A3</td><td>845,09</td><td rowspan="2">-6%</td></tr><tr><td>Localisée</td><td>2016/2017</td><td>A3</td><td>791,67</td></tr><tr><td rowspan="2">Maïs ensilage</td><td>Gravitaire</td><td>2014/2015</td><td>A1</td><td>491,82</td><td rowspan="2">-7%</td></tr><tr><td>Localisée</td><td>2014/2015</td><td>A1</td><td>458,44</td></tr><tr><td rowspan="2">Agrumes</td><td>Gravitaire</td><td>2017/2018</td><td>A4</td><td>701,07</td><td rowspan="2">-7%</td></tr><tr><td>Localisée</td><td>2017/2018</td><td>A4</td><td>652,79</td></tr><tr><td rowspan="2">Oliviers</td><td>Gravitaire</td><td>2017/2018</td><td>A4</td><td>504,21</td><td rowspan="2">-5%</td></tr><tr><td>Localisée</td><td>2017/2018</td><td>A4</td><td>479,65</td></tr><tr><td rowspan="2">Maraichage (Pomme de Terre)</td><td>Gravitaire</td><td>2016/2017</td><td>A3</td><td>560,48</td><td rowspan="2">-12%</td></tr><tr><td>Localisée</td><td>2016/2017</td><td>A3</td><td>494,80</td></tr></tbody></table>

Les besoins annuels nets d'irrigation des cultures se situent dans la gamme de 445 à 845 mm par cycle productif cultural en cas du système d'irrigation gravitaire et de 391 à 792 mm/cycle productif en cas du système d'irrigation localisée. Les réductions d'eau correspondant sont donc de 5 à 12 % en passant du premier au deuxième système d'irrigation (Tableau 75).

Il est alors approuvé que le système d'irrigation localisée confère une relative meilleure performance par le rationnement dans le temps des besoins nets en eau pour les cultures à l'échelle parcellaire. Ceci étant avant l'exécution de l'opération de distribution d'eau au niveau de la parcelle agricole. Sachant que cette distribution intègre un complément d'eau pour la compensation de l'écart à l'origine de l'efficience de la technique d'irrigation adoptée.

En fin du compte, le raisonnement déjà initié pour les besoins nets en eau des cultures sera parachevé ci-après en intégrant les besoins bruts en eau d'irrigation des spéculations agricoles envisagées.
### III-Evaluation des besoins bruts en eau actuels de l'irrigation et les performances agronomiques liées à l'eau pour les cultures du périmètre de Tadla.

Il est objet d'exploiter les résultats des expérimentations des essais culturaux réalisés durant la période de 2014/2015 à 2017/2018 au niveau de la station expérimentale de référence nationale SEHA-Ouled Gnou, située au sein du périmètre irrigué du Tadla, dont les fiches des protocoles expérimentaux sont portées en annexe et les illustrations *in situ* sont sur les Figures 44-46.

En effet, ces essais ont permis d'aboutir à des caractérisations agronomiques spécifiques à l'utilisation de l'eau d'irrigation pour les principales cultures pratiquées dans le périmètre de Tadla. L'essentiel des résultats expérimentaux est évoqué ci-après pour chaque culture étudiée.

#### III-1- Caractéristiques agronomiques en relation avec l'eau pour les céréales d'automne

Les principaux résultats des essais en céréales d'automne sont sur le Tableau 76 avec illustration à la Figure 44.

<table><caption>Tableau 76 : Résultats des essais culturaux sur céréales d'automne dans la Station SEHA-Ouled Gnaou au périmètre de Tadla</caption><thead><tr><th rowspan="2">Objectif de l'expérimentation</th><th rowspan="2">Référence de l'essai</th><th rowspan="2">Campagne agricole de</th><th colspan="3">Résultats de l'essai</th></tr><tr><th>Volume eau irrigation (en m3/ha)</th><th>Rendement (en Qx/ha)</th><th>Efficience agronomique d'utilisation eau irrigation EUE (en Kg/m3 eau)</th></tr></thead><tbody><tr><td rowspan="2">Comportement variétal</td><td>1<sup>er</sup> Essai cultural</td><td>2016/2017</td><td>6 720</td><td>62,21 ± 1,01.</td><td>0,93 ± 0,01.</td></tr><tr><td>2<sup>ème</sup> Essai cultural</td><td>2017/2018</td><td>1 920</td><td>67,17 ± 2,33</td><td>3,50 ± 0,12.</td></tr></tbody></table>

Les céréales d'automne (blés) ont une productivité similaire pendant deux campagnes successives, soit de 62 à 67 quintaux à l'hectare, mais avec des exigences en eau d'irrigation très variables en raison de l'état très pluvieux de la campagne 2017/2018, ce qui s'est répercuté sur l'EUE fluctuant donc entre 1,0 à 3,5 Kg en grain de blé par mètre cube d'eau fourni par l'irrigation (Tableau 76).

#### III-2 Caractéristiques agronomiques en relation avec l'eau pour la betterave à Sucre

La culture de la betterave à sucre a été soumise au périmètre du Tadla à deux types d'expérimentation à savoir les suivantes (Figure 44):

* Comparaison de trois systèmes d'irrigation localisée, aspersion et gravitaire avec réalisation des essais culturaux en deux versions pendant les campagnes de 2014/2015 et de 2015/2016.
* Modalités du pilotage des irrigations de la culture de la betterave à sucre avec le système d'irrigation localisée, moyennant des tests de réduction d'eau par rapport aux exigences. Les essais culturaux ont été menés pendant les campagnes de 2016/2017 et de 2017/2018. Les aboutissements de ces expérimentations sont formulés ci-après

### III-2-1 Comparaison des trois techniques d'irrigation sur la betterave à sucre

Les résultats des travaux de comparaison des trois systèmes d'irrigation appliqués à la betterave à sucre sont présentés dans le Tableau 77.

<table><caption>Tableau 77 : Résultats des essais culturaux des techniques d'irrigation sur la betterave à sucre à la Station SEHA-Ouled Gnaou au périmètre de Tadla (Ière version)</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th rowspan="2">Système d'irrigation</th><th rowspan="2">Volume d'irrigation (m3/ha)</th><th colspan="2">Produit: Racine</th><th colspan="2">Produit: Sucre</th></tr><tr><th>Rendement (T/ha)</th><th>Efficience agronomique (Kg/m3 eau)</th><th>Rendement (T/ha)</th><th>Efficience agronomique (Kg/m3 eau)</th></tr></thead><tbody><tr><td rowspan="9">2014/2015</td><td rowspan="3">Indicateurs agronomiques</td><td>Localisée</td><td>5 400 a</td><td>64,92 a</td><td>12,02 a</td><td>11,89 a</td><td>2,20 a</td></tr><tr><td>Aspersion</td><td>6 110 b</td><td>61,05 ab</td><td>9,99 b</td><td>11,01 ab</td><td>1,80 ab</td></tr><tr><td>Gravitaire</td><td>8 950 c</td><td>59,48 b</td><td>6,65 c</td><td>10,26 b</td><td>1,15 b</td></tr><tr><td rowspan="2">Ecarts par rapport au système d'irrigation gravitaire</td><td>Localisée</td><td>-39,66%</td><td>9,15%</td><td>80,75%</td><td>15,89%</td><td>91,30%</td></tr><tr><td>Aspersion</td><td>-31,73%</td><td>2,64%</td><td>50,23%</td><td>7,31%</td><td>56,52%</td></tr><tr><td rowspan="9">2015/2016</td><td rowspan="2">Indicateurs agronomiques</td><td>Localisée</td><td>5 170 a</td><td>75,00 a</td><td>14,51 a</td><td>13,17 a</td><td>2,55 a</td></tr><tr><td>Aspersion</td><td>5 957 b</td><td>74,25 a</td><td>12,46 b</td><td>12,3 a</td><td>2,06 b</td></tr><tr><td rowspan="3">Ecarts entre les systèmes d'irrigation localisée et d'aspersion</td><td></td><td>-13,21%</td><td>1,01%</td><td>16,39%</td><td>7,07%</td><td>23,37%</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Les analyses en guise de ces essais culturaux comparatifs sur la betterave à sucre à propos des techniques comparatives d'irrigation sont données ci-dessous.

#### a) Comparaison des trois techniques de systèmes d'irrigation : Gravitaire-Aspersion-Localisée

La culture de la betterave à sucre exige un volume d'eau d'irrigation de 5 400 à 8 950 m3/ha selon le système d'irrigation. Le système d'irrigation localisée économise effectivement de l'eau d'irrigation par rapport au gravitaire à 40% et reste presque similaire par rapport au système d'irrigation en aspersion.

Quant à l'efficience agronomique de l'utilisation de l'eau par la culture betteravière, elle varie de 6,7 à 12,0 kg/m3 d'eau en production racinaire et de 1,2 à 2,2 kg/m3 en production du sucre. En notant que comparé au gravitaire, le système du localisé réalise de meilleure performance aussi
bien qu'en produit racine (81 % de gain) qu'en produit sucre (+91 % de gain), suivi par la
technique de l'aspersion avec des accroissements respectifs de 50 % et de 57% (Tableau 77).

b) Comparaison de deux techniques d'irrigation : Aspersion-Localisée

Dans cet essai à caractère répétitif pour la deuxième campagne par rapport à l'essai précédent, en se limitant aux deux systèmes d'irrigation à savoir goutte à goutte et l'aspersion, lesquels ont déjà attestés leur efficience améliorée par rapport au système d'irrigation gravitaire. La culture de la betterave exprime un besoin en volume d'eau d'irrigation variant entre 5 170 et 5 957 m3/ha selon la technique d'irrigation.

Le système d'irrigation localisée approuve encore son aptitude à l'économie d'eau en un taux de gain de 13%. Quant à son efficience agronomique d'utilisation de l'eau d'irrigation, elle est aussi plus avantageuse en production racinaire (+16 %) et en production du sucre (+23 %). Ceux-ci ont eu lieu malgré que les rendements (racine et sucre) étaient presque rapprochés pour les deux systèmes d'irrigation comparées (Tableau 77).
<figure><img src="image_70.png" /><figcaption>1er Essai cultural des céréales d'automne
(campagne 2016-2017)</figcaption></figure>

<figure><img src="image_71.png" /><figcaption>2ème Essai cultural des céréales d'automne
(campagne 2017-2018)</figcaption></figure>

<figure><img src="image_72.png" /><figcaption>Système d'irrigation localisée</figcaption></figure>

<figure><img src="image_73.png" /><figcaption>Système d'irrigation par aspersion</figcaption></figure>

Essai cultural comparatif des systèmes d'irrigation sur betterave à sucre (2014-2015)

<figure><img src="image_74.png" /><figcaption>Système d'irrigation localisée</figcaption></figure>

<figure><img src="image_75.png" /><figcaption>Système d'irrigation par aspersion</figcaption></figure>

Essai cultural comparatif des systèmes d'irrigation sur betterave à sucre (2015-2016)

<caption>**Figure 44 : Illustrations sur essais culturaux de céréaliculture et de la betterave à sucre**
*(expérimentations 1) réalisés dans la station SEHA-Ouled Gnaou (Périmètre de Tadla ; période : 2014-2018)*</caption>
### III-2-2 Effet de la réduction d'eau pour la culture de la betterave à sucre

Les résultats des expérimentations sur l'impact de la réduction d'eau au cours du cycle productif de la betterave à sucre sont présentés au Tableau 78 et en des images sur les lieux à la Figure 45.

<table><caption>Tableau 78 : Résultats des essais de la réduction d'eau sur la production de la betterave à sucre à la Station SEHA-Ouled Gnaou dans le périmètre de Tadla (2ème version)</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Régime hydrique</th><th rowspan="2">Pluie (mm)</th><th rowspan="2">Volume eau irrigation (m3/ha)</th><th colspan="2">Produit: Racine</th><th colspan="2">Produit: Sucre</th></tr><tr><th>Rendement (T/ha)</th><th>Efficience agronomique d'utilisation d'eau (kg/m3)</th><th>Rendement (T/ha)</th><th>Efficience agronomique d'utilisation d'eau (kg/m3)</th></tr></thead><tbody><tr><td rowspan="4">2016/2017</td><td>100% ETM</td><td>195,5</td><td>5 350,0</td><td>68,07 b</td><td>12,72 a</td><td>13,28 a</td><td>2,48 a</td></tr><tr><td>90% ETM</td><td>195,5</td><td>4 890,0</td><td>60,02 ab</td><td>12,27 a</td><td>14,59 a</td><td>2,98 a</td></tr><tr><td>80% ETM</td><td>195,5</td><td>4 430,0</td><td>57,27 a</td><td>12,93 a</td><td>10,92 a</td><td>2,47 a</td></tr><tr><td>70% ETM</td><td>195,5</td><td>3 970,0</td><td>42,53 a</td><td>10,71 a</td><td>8,64 a</td><td>2,18 a</td></tr><tr><td rowspan="4">2017/2018</td><td>100% ETM</td><td>378,6</td><td>4 124,3</td><td>59,48 a</td><td>14,42 a</td><td>2,68 a</td><td>0,65 a</td></tr><tr><td>90% ETM</td><td>378,9</td><td>3 769,3</td><td>60,66 a</td><td>16,09 ab</td><td>3,10 a</td><td>0,82 b</td></tr><tr><td>80% ETM</td><td>378,6</td><td>3 414,4</td><td>61,44 a</td><td>18,0 ab</td><td>2,46 a</td><td>0,72 ab</td></tr><tr><td>70% ETM</td><td>378,6</td><td>3 026,2</td><td>61,39 a</td><td>20,29 b</td><td>3,26 a</td><td>1,08 ab</td></tr></tbody></table>

A l'issue d'essais entrepris sur le système d'irrigation gravitaire, il a été constaté que la réduction de l'eau pour la betterave sucrière est susceptible d'abaisser le rendement racinaire, vu qu'une diminution de l'apport d'eau nécessaire de 30 % de l'ETM, induit la régression du tonnage produit de 68,1 à 42,5 T/ha, alors que le rendement en sucre reste indifférent à l'approvisionnement en eau.

L'efficience de l'utilisation d'eau par la culture betteravière (EUE), est significativement la meilleure sous le régime hydrique à 70% ETM pour le produit racine (20,3 Kg/m3) et pour le produit sucre. L'EUE reste irrégulière vis-à-vis la réduction de l'eau d'irrigation (Tableau 78).

### III-3 Caractéristiques agronomiques en relation avec l'eau pour la Luzerne

Les résultats saillants des essais culturaux réalisés sur la culture fourragère la luzerne, sont présentés dans le Tableau 79 et sur Figure 45.

<caption>Tableau 79 : Résultats des essais culturaux sur luzerne</caption>

(Station SEHA-Ouled Gnaou, périmètre de Tadla)

<table><thead><tr><th rowspan="3">Objectif de l'expérimentation</th><th rowspan="3">Campagne agricole</th><th colspan="4">Résultats de l'essai</th></tr><tr><th rowspan="2">Volume eau irrigation (en m3/ha)</th><th colspan="2">Rendement (en T/ha)</th><th colspan="2">Efficience agronomique d'utilisation eau irrigation EUE (en Kg/m3 eau)</th></tr><tr><th>Matière fraîche</th><th>Matière sèche</th><th>Matière fraîche</th><th>Matière sèche</th></tr></thead><tbody><tr><td>Comportement variétal</td><td>2016/2017</td><td>7 680</td><td>35,13 ± 0,29</td><td>8,48 ± 0,06</td><td>4,57 ± 0,04</td><td>1,10 ± 0,01</td></tr></tbody></table>
<figure><img src="image_76.png" /><figcaption>Essai de réduction hydrique sur betterave à sucre sous irrigation localisée (2016-2017)</figcaption></figure>

<figure><img src="image_77.png" /><figcaption>Essai cultural de la luzerne (2016-2017)</figcaption></figure>

<figure><img src="image_78.png" /><figcaption>1-Régime hydrique à 100 % ETM</figcaption></figure>

<figure><img src="image_79.png" /><figcaption>2-Régime hydrique à 90 % ETM</figcaption></figure>

<figure><img src="image_80.png" /><figcaption>3-Régime hydrique à 80 % ETM</figcaption></figure>

<figure><img src="image_81.png" /><figcaption>4-Régime hydrique à 70 % ETM</figcaption></figure>

Essai de réduction hydrique sur betterave à sucre sous irrigation localisée (2017-2018)

**Figure 45** : Illustrations sur essais culturaux de luzerne et de la betterave à sucre (expérimentations 2) réalisés dans la station SEHA-Ouled Gnaou (Périmètre de Tadla ; période : 2016-2018)

La jeune luzernière installée d'âge seulement de 16 mois, a été soumise à 4 coupes durant la campagne d'essai et a exigé un volume d'eau d'irrigation sous le système d'irrigation gravitaire
de 7 680 m3/ha pour réaliser un rendement en frais de l'ordre de 35 T/ ha, correspondant à 8,5 T/ha en matière sèche et correspondant aux EUE respectives de 4,6 et 1,1 Kg/m3.

### III-4 Caractéristiques agronomiques en relation avec l'eau pour la culture du maïs

Les résultats de deux essais culturaux sur le maïs à des fins de la production d'ensilage pour l'affourragement des animaux d'élevages (ovins et bovins), sont récapitulés dans le Tableau 80 et conjointement en des photos sur la Figure 46.

<table><caption>Tableau 80 : Résultats des essais culturaux sur le Maïs ensilage réalisés à Station SEHA-Ouled Gnaou, périmètre du Tadla</caption><thead><tr><th rowspan="2">Campagne</th><th rowspan="2">Objectif de l'expérimentation</th><th rowspan="2">Variétés/ Régime hydrique</th><th rowspan="2">Pluie (en mm)</th><th rowspan="2">Volume eau irrigation (m3/ha)</th><th colspan="2">Produit: Grain</th><th colspan="2">Produit: Ensilage</th></tr><tr><th>Rendement (T/ha)</th><th>Efficience agronomique d'utilisation d'eau (kg/m3)</th><th>Rendement (Qx/ha)</th><th>Efficience agronomique d'utilisation d'eau (kg/m3)</th></tr></thead><tbody><tr><td rowspan="5">2016/2017</td><td rowspan="4">Impact de réduction de l'eau</td><td>100% ETM</td><td>30,4</td><td>5 460</td><td>19,6 a</td><td>3,60 a</td><td>105 a</td><td>1,92 a</td></tr><tr><td>125% ETM</td><td>30,4</td><td>6 590</td><td>22,4 b</td><td>3,40 a</td><td>118 a</td><td>1,79 a</td></tr><tr><td>75% ETM</td><td>30,4</td><td>4 110</td><td>14,8 c</td><td>3,60 a</td><td>69 b</td><td>1,68 a</td></tr><tr><td>50% ETM</td><td>30,4</td><td>3 170</td><td>11,0 d</td><td>3,47 a</td><td>47 c</td><td>1,48 a</td></tr><tr><td rowspan="2">Inscription au Catalogue National</td><td>Précoces</td><td>25</td><td>1 056</td><td>55,26 a</td><td>52,33 a</td><td>21,39 a</td><td>20,26 a</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>Tardives</td><td>25</td><td>1 152</td><td>75,22 b</td><td>65,30 b</td><td>26,25 b</td><td>22,79 b</td></tr></tbody></table>

La variation de l'apport d'eau à la culture du maïs par rapport à ses exigences normales sous conduite du système d'irrigation localisée, avec des volumes d'irrigation allant de 3 170 à 6 590 m3/ha s'est accompagnée dans le même sens par la fluctuation des niveaux des rendements en grain et en ensilage. Par ailleurs, l'efficience EUE a été maintenue constante pour les deux cas de la production en grain et d'ensilage du maïs (Tableau 80).

Il est aussi constaté que les variétés tardives du maïs ont les meilleures performances productives et elles valorisent agronomiquement mieux l'eau par apport aux variétés précoces (Tableau 80).
<figure><img src="image_82.png" /><figcaption>Epis récoltés à la maturité en culture
du maïs d'ensilage conduite sous le système d'irrigation localisée</figcaption></figure>

<figure><img src="image_83.png" /><figcaption>Culture du maïs sous d'irrigation localisée</figcaption></figure>

<figure><img src="image_84.png" /><figcaption>Disposition du système d'irrigation gravitaire</figcaption></figure>

## Essai de réduction d'eau pour culture du Maïs ensilage-grain (2017-2018)

<figure><img src="image_85.png" /><figcaption>1er Essai cultural de la Pomme de Terre
(2014-2015)</figcaption></figure>

<figure><img src="image_86.png" /><figcaption>2ème Essai cultural de la Pomme de Terre
(2017-2018)</figcaption></figure>

## Essais culturaux en culture de la Pomme de Terre

<caption>**Figure 46 :** Illustrations sur essais culturaux du maïs et maraichage (expérimentations 3) réalisés dans la station SEHA-Ouled Gnaou (Périmètre de Tadla ; période : 2016-2018)</caption>
### III-5 Caractéristiques agronomiques en relation avec l'eau pour le maraichage

Les résultats des essais culturaux de la spéculation maraichère modèle la Pomme de Terre sont résumés au Tableau 81. Aussi des photos illustrant les expérimentation sur terrain sont sur la Figure 46.

<table><caption>Tableau 81 : Résultats d'essais culturaux sur la Pomme de Terre accomplis au niveau de la Station SEHA-Ouled Gnaou au périmètre du Tadla</caption><thead><tr><th rowspan="2">Objectif de l'expérimentation</th><th rowspan="2">Référence de l'essai</th><th rowspan="2">Campagne agricole</th><th colspan="3">Résultats de l'essai</th></tr><tr><th>Volume eau irrigation (en m3/ha)</th><th>Rendement (en T/ha)</th><th>Efficience d'utilisation d'eau EUE (en Kg/m3)</th></tr></thead><tbody><tr><td rowspan="2">Inscription au Catalogue Officiel National des variétés agricoles</td><td>1er Essai</td><td>2014/2015</td><td>6 720</td><td>33,51 ± 0,55</td><td>4,55 ± 0,08</td></tr><tr><td>2ème Essai</td><td>2016/2017</td><td>7 680</td><td>35,82 ± 0,13.</td><td>4,66 ± 0,02.</td></tr></tbody></table>

L'exigence en eau d'irrigation par la Pomme de Terre en culture d'hivers et sous irrigation gravitaire est évaluée de 6 720 à 7680 m3/ha, en réalisant un rendement de 33,5 à 35,8 T/ha et une efficience d'utilisation d'eau qui est de l'ordre de 4,6 Kg/m3 d'eau (Tableau 81).

### III-6-Ordres de grandeurs des efficiences agronomiques d'utilisation de l'eau par les cultures dans le périmètre du Tadla déterminées expérimentalement

Le récapitulatif des résultats des essais culturaux précités concernant les efficiences agronomiques d'utilisation d'eau d'irrigation (EUE) est présenté dans le Tableau 82.

<table><caption>Tableau 82 : EUE des cultures menées dans le périmètre de Tadla
(Période : campagnes agricoles de 2015 à 2018)</caption><thead><tr><th>Cultures soumises aux expérimentations</th><th>Système d'irrigation</th><th>Efficience EUE (en kg/m3)</th></tr></thead><tbody><tr><td>Céréales d'automne</td><td>Gravitaire</td><td>0,93 à 3,5</td></tr><tr><td rowspan="3">Betterave à sucre</td><td>Gravitaire</td><td>0,65 à 2,48</td></tr><tr><td>Localisée</td><td>2,20 à 2,55</td></tr><tr><td>Aspersion</td><td>1,80 à 2,06</td></tr><tr><td>Luzerne</td><td>Gravitaire</td><td>1,10</td></tr><tr><td rowspan="2">Mais ensilage</td><td rowspan="2">Gravitaire</td><td>Précoce (frais) : 32,77 à 65,74</td></tr><tr><td>Tardif (frais): 44,01 à 93,48</td></tr><tr><td>Maraichage (Pomme de Terre)</td><td>Gravitaire</td><td>4,55 à 4,66</td></tr></tbody></table>

L'efficience EUE pour une même culture avec même système d'irrigation reste plus dépendante aux conditions du déroulement de l'essai cultural, dont il est cité essentiellement pour les causes-à-effets météorologiques de chaque campagne agricole. Il est à noter, aussi, que cette efficience s'améliore davantage avec le système d'irrigation localisée sur culture de la betterave sucrière
(rendement sucre) et augmente avec le caractère tardif des variétés exploitées en culture du maïs destiné à l'ensilage (Tableau 86).

## IV-Détermination de l'efficience d'application de l'eau à la parcelle pour les systèmes d'irrigation gravitaire et localisée au périmètre de Tadla

En se référant aux données déterminées ultérieurement concernant les besoins nets en eau d'irrigation pour les principales cultures conduites dans le périmètre de Tadla (Tableau 75) et en exploitant les résultats des expérimentations réalisées (Tableaux 76-81), il est évalué les efficiences comparatives des deux systèmes d'irrigation gravitaire et localisée, dont résultats sont cités au Tableau 83.

<table><caption>Tableau 83 : Besoins en eau bruts et efficiences des systèmes d'irrigation gravitaire et localisée pour les principales cultures conduites dans le périmètre de Tadla, en se référant aux enquêtes et expérimentations
(Campagnes de référence de 2014 à 2018)</caption><thead><tr><th rowspan="2">Culture</th><th rowspan="2">Système d'irrigation</th><th rowspan="2">Campagne de référence</th><th rowspan="2">Rang annuel</th><th colspan="2">Besoins bruts en eau d'irrigation</th><th colspan="2">Efficience du système d'irrigation</th></tr><tr><th>Volume d'eau (m3/ha)</th><th>Ecarts entre système d'irrigation Localisée (2)/Gravitaire (1) (en %)</th><th>Rapport entre besoins nets et besoins bruts d'eau d'irrigation (%)</th><th>Ecarts entre système d'irrigation Localisée (2)/Gravitaire (1) (en nombre de points)</th></tr></thead><tbody><tr><td rowspan="2">Céréales d'automne (Blés)</td><td>Gravitaire (1)</td><td>2016/2017</td><td>A3</td><td>6 720</td><td rowspan="2">-37%</td><td>66,25%</td><td rowspan="2">26</td></tr><tr><td>Localisée (2)</td><td>2016/2017</td><td>A3</td><td>4 240</td><td>92,26%</td></tr><tr><td rowspan="4">Betterva à sucre</td><td rowspan="2">Gravitaire</td><td>2014/2015</td><td>A1</td><td>8 950</td><td rowspan="4">-43%</td><td>64,34%</td><td rowspan="4">32</td></tr><tr><td>2015/2016</td><td>A2</td><td>9 640</td><td>53,35%</td></tr><tr><td rowspan="2">Localisée</td><td>2014/2015</td><td>A1</td><td>5 400</td><td>91,08%</td></tr><tr><td>2015/2016</td><td>A2</td><td>5 170</td><td>90,21%</td></tr><tr><td rowspan="2">Luzerne</td><td>Gravitaire</td><td>2016/2017</td><td>A3</td><td>11 480</td><td rowspan="2">-27%</td><td>73,61%</td><td rowspan="2">21</td></tr><tr><td>Localisée</td><td>2016/2017</td><td>A3</td><td>8 390</td><td>94,36%</td></tr><tr><td rowspan="2">Maïs ensilage</td><td>Gravitaire</td><td>2014/2015</td><td>A1</td><td>7 680</td><td rowspan="2">-36%</td><td>64,04%</td><td rowspan="2">30</td></tr><tr><td>Localisée</td><td>2014/2015</td><td>A1</td><td>4 878</td><td>93,98%</td></tr><tr><td rowspan="2">Agrumes</td><td>Gravitaire</td><td>2017/2018</td><td>A4</td><td>10 000</td><td rowspan="2">-28%</td><td>70,11%</td><td rowspan="2">20</td></tr><tr><td>Localisée</td><td>2017/2018</td><td>A4</td><td>7 230</td><td>90,29%</td></tr><tr><td rowspan="2">Oliviers</td><td>Gravitaire</td><td>2017/2018</td><td>A4</td><td>6 800</td><td rowspan="2">-24%</td><td>74,15%</td><td rowspan="2">19</td></tr><tr><td>Localisée</td><td>2017/2018</td><td>A4</td><td>5 158</td><td>92,99%</td></tr><tr><td rowspan="2">Maraichage (Pomme de Terre)</td><td>Gravitaire</td><td>2016/2017</td><td>A3</td><td>7 680</td><td rowspan="2">-33%</td><td>72,98%</td><td rowspan="2">23</td></tr><tr><td>Localisée</td><td>2016/2017</td><td>A3</td><td>5 130</td><td>96,45%</td></tr></tbody></table>

Il est retenu que les cultures les plus exigeantes pour les besoins bruts en eau de l'irrigation à l'échelle parcellaire, quel que soit le système d'irrigation adopté sont : la luzerne (8 390 à 11 480 m3/ha), suivie par les agrumes (7 230 à 10 000 m3/ha) et ensuite la betterave à sucre (5 170 à 9 640 m3/ha). Les cultures les moins nécessiteuses en eau sont les céréales d'automne (4 240 à 6720 m3/ha) et les plantations d'oliviers (5 158 à 6 800 m3/ha) (Tableau 83).
Il est à signaler qu'après calcul de l'efficience d'application de l'eau à la parcelle (rapport entre les besoins nets et les besoins bruts en eau de l'irrigation), il est démontré que la technique d'irrigation localisée réalise de meilleure efficience avec un taux allant de 90,3 % à 96,5 % par rapport à l'irrigation gravitaire dont le taux ne dépasse pas au maximum la barre de 74 % (Tableau 83).

# V-Occupation du sol par les cultures et l'intensification culturale en relation avec la source d'eau d'irrigation au périmètre de Tadla

## V-1-Occupation du sol par les cultures dans la globalité de la zone irriguée au périmètre de Tadla

La zone irriguée du périmètre de Tadla est un support d'occupations du sol par diverses cultures sur un total de superficie annuelle qui s'élève à 101 922 ha en campagne de référence de 2017/2018. La grande part de l'assolement (82,2 %) est prédominée essentiellement par les cultures des céréales d'automne, la luzerne et les oliviers en des contributions respectives de 22,8 % ; 26,8% et 21 ?6 % ; alors que les 17,8 % restant sont répartis entre les cultures d'agrumes (9,8%), la betterave à sucre (7,2%) ; le maraichage (0,7%) et le maïs ensilage (0,7%). Il est rapporté pour ce périmètre, qu'il est à son juste niveau requis de l'exploitation des assolements culturaux, sachant que le taux de l'intensification cultural (T.I.C) actuel est à sa valeur équilibrée de 1,04.

Par ailleurs, il est noté que les besoins bruts en eau d'irrigation pour l'essentiel des cultures irriguées afférentes au réseau hydro-agricole ont été évalués à environ 846 millions de m3 par an, dont la quasi dominance à 93 % des superficies incombe à la culture fourragère de la luzerne (36%), les plantations arboricoles fruitières (32%) et aussi aux céréales d'automne (36%) (Tableau 83).

Les prélèvements d'eau pour l'ensemble de la zone du réseau hydro-agricole (gravitaire et goutte à goutte), à partir des barrages locaux sont en campagne de 2017/2018 pour un volume annuel de 622 millions de m3 accordé par l'Agence du Bassin Hydraulique d'Oum Rbia (ORMVA du Tadla, 2020 b). Ce volume d'eau qui arrive à peine à couvrir les besoins bruts en eau des cultures de la zone réseau pour seulement les 73% des exigences. Le reste des besoins en eau non satisfaits (27%) est apporté par les agriculteurs pour combler le déficit, via le pompage sur place à partir de la nappe d'eau souterraine d'eau ou à défaut sont délaissées passivement les cultures aux faits réels des déficits hydriques inéluctables (Tableau 83).
Le volume annuel d'eau nécessaire à l'hectare assolé dans le périmètre du Tadla est estimé à 8 297 m3/ha, contre la dotation en eau du barrage pour volume de 6 003 m3/ha/an, l'écart à restituer donc est de 2 214 m3/ha/an, qui constitue alors le volume potentiel unitaire de pompage rapporté à la l'unité de la surface à exploiter en zone irriguée du réseau hydro-agricole (Tableau 84).

<table><caption>Tableau 84 : Assolements culturaux et intensification culturale en relation avec la source d'eau d'irrigation dans le périmètre du Tadla (Campagne de référence : 2017/2018).</caption><thead><tr><th rowspan="3">Cultures</th><th rowspan="2">Superficie (en ha)</th><th rowspan="2">Part de la culture dans l'assolement</th><th rowspan="2">Besoins bruts totaux eau irrigation par culture</th><th colspan="2">Eau de barrage</th><th colspan="2">Eau de pompage de nappe souterraine</th><th rowspan="2">Volume total des besoins en eau d'irrigation à l'hectare assolé</th><th rowspan="2">Part maximale de contribution du pompage au volume eau irrigation exigée</th></tr><tr><th>Dotation totale au périmètre</th><th>Dotation eau à l'hectare assolé</th><th>Potentiel total de complément eau exigée</th><th>Prélevement eau par hectare assolé</th></tr><tr><th>(en ha)</th><th>(en %)</th><th>(en 1000 m3)</th><th>(en 1000 m3)</th><th>(en m3/ ha)</th><th>(en 1000 m3)</th><th>(en m3/ ha)</th><th>(en m3/ ha)</th><th></th></tr></thead><tbody><tr><td>Céréales d'automne</td><td>34 516</td><td>16,93%</td><td>209 374</td><td rowspan="7">620 000</td><td rowspan="7">6 083</td><td rowspan="7">225 646</td><td rowspan="7">2 214</td><td rowspan="7">8 297</td><td rowspan="7">27%</td></tr><tr><td>Betterave à Sucre</td><td>7 325</td><td>3,59%</td><td>60 021</td></tr><tr><td>Luzerne</td><td>27 323</td><td>13,40%</td><td>300 660</td></tr><tr><td>Mais Ensilage</td><td>90</td><td>0,04%</td><td>601</td></tr><tr><td>Agrumes</td><td>9 966</td><td>4,89%</td><td>111 619</td></tr><tr><td>Oliviers</td><td>22 031</td><td>10,81%</td><td>159 372</td></tr><tr><td>Maraichage</td><td>671</td><td>0,33%</td><td>3 999</td></tr><tr><td>TOTAL</td><td>101 922 (1)</td><td>100,00%</td><td>845 646</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Superficie zone réseau hydro-agricole équipée</td><td colspan="9">98 300 (2)</td></tr><tr><td>Taux d'intensification cultutal : T.I.C. [(1)/(2)]</td><td colspan="9">1,04</td></tr></tbody></table>

## V-2-Evolution comparative de l'occupation du sol entre avant et après projet de la reconversion collective à l'irrigation localisée selon les méthodes télédétection spatiale et statistique terrestre

Sur l'espace pilote au niveau des réseaux M4-M7 relevant du CDA 536 en zone d'action de l'ORMVA du Tadla, d'une superficie agricole utile (SAU) de 3 662 ha, l'évolution de l'occupation du sol par les cultures entre avant et après projet de la reconversion collective à l'irrigation localisée, a été approché d'une manière comparative par la méthode de la télédétection spatiale avec indice NDVI et par la méthode statistique terrestre, dont les résultats sont exposés et analysés ci-après.
En effet, les occupations du sol de la zone précitée pour les deux campagnes agricoles références de 2013/2014 et de 2017/2018 aux dates respectives du 17/07/2014 et du 15/07/2018, avec appui de cartes SIG suite au traitement d'images satellitaires Landsat_OLI8 avec indice NDVI sont présentées aux Figures 47 et 48.

<figure><img src="image_87.png" /><figcaption>Figure 47 : Carte d'occupations du sol en mi-Juillet 2014 en zone des réseaux M4-M7 du périmètre de Tadla, en avant-projet de la reconversion collective à l'irrigation localisée.</figcaption></figure>

En passant de l'avant (année: 2014) à l'après (année: 2018) lancement de la reconversion collective à l'irrigation localisée en cette zone, le taux de l'occupation du sol par les cultures en date de mi-juillet a évolué de 49,75 à 59,80, % (Figures 47 et 48).

<figure><img src="image_88.png" /><figcaption>Figure 48 : Carte d'occupations du sol en mi-Juillet 2018 en zone des réseaux M4-M7 du périmètre de Tadla, en après-projet de la reconversion collective à l'irrigation localisée.</figcaption></figure>

Les analyses comparatives entre les deux méthodes de la télédétection spatiale et par les statistiques terrestres concernant l'évaluation de l'occupation du sol par les cultures sont rapportées dans le Tableau 85.
<table><caption>Tableau 85 : Occupations du sol entre avant/ après projet de reconversion collective à l'irrigation localisée par méthodes comparatives de la télédétection spatiale et des statistiques terrestres aux campagnes 2013/2014 et 20117/2018 sur la zone modèle des réseaux M4-M7 (périmètre de Tadla)</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="3">Campagne 2013/2014</th><th colspan="3">Campagne 2017/2018</th></tr><tr><th>Septembre-<br>Avril (1)</th><th>Mai-Août (2)</th><th>Total emblavé<br>[(1)+(2)]</th><th>Septembre-<br>Avril (1)</th><th>Mai-Août (2)</th><th>Total emblavé<br>[(1)+(2)]</th></tr></thead><tbody><tr><td rowspan="10">Statistiques/CDA</td><td>Emblavements ( en ha)</td><td>3 391</td><td>3 338</td><td>3 932</td><td>3 939</td><td>3 861</td><td>4 713</td></tr><tr><td>Céréles d'automne</td><td>891</td><td>891</td><td>891</td><td>1 138</td><td>1 138</td><td>1 138</td></tr><tr><td>Betterave</td><td>363</td><td>363</td><td>363</td><td>290</td><td>290</td><td>290</td></tr><tr><td>Luzerne</td><td>689</td><td>689</td><td>689</td><td>816</td><td>816</td><td>816</td></tr><tr><td>Mais</td><td>0</td><td>79</td><td>79</td><td>0</td><td>148</td><td>148</td></tr><tr><td>Plantations</td><td>844</td><td>854</td><td>844</td><td>842</td><td>842</td><td>842</td></tr><tr><td>Agrumes</td><td>339</td><td>348</td><td>339</td><td>336</td><td>336</td><td>336</td></tr><tr><td>Oliviers</td><td>506</td><td>506</td><td>506</td><td>507</td><td>507</td><td>507</td></tr><tr><td>Maraichage</td><td>604</td><td>461</td><td>1 065</td><td>852</td><td>626</td><td>1 478</td></tr><tr><td>SAU</td><td></td><td>3 662</td><td>3 662</td><td>3 662</td><td>3 662</td><td>3 662</td><td>3 662</td></tr><tr><td>TIC (CDA)</td><td></td><td></td><td></td><td>1,07</td><td></td><td></td><td>1,29</td></tr><tr><td rowspan="5">Télédétection<br>spatiale par indice<br>NDVI</td><td>Emblavements ( en ha)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Occupation du sol par NDVI</td><td></td><td>2 119</td><td></td><td></td><td>2 547</td><td></td></tr><tr><td>Betterave à sucre restante</td><td></td><td>73</td><td></td><td></td><td>58</td><td></td></tr><tr><td>Nouveaux emblavements entre<br>Mai-Août</td><td></td><td>503</td><td></td><td></td><td>831</td><td></td></tr><tr><td>Assolement total/campagne</td><td></td><td></td><td>3 894</td><td></td><td></td><td>4 769</td></tr><tr><td></td><td>TIC (NDVI)</td><td></td><td></td><td>1,06</td><td></td><td></td><td>1,30</td></tr><tr><td rowspan="3">Ecart entre<br>statistiques CDA et<br>Télédétection<br>spatiale</td><td>Statistiques CDA<br>(Mi-Juillet) (3)</td><td></td><td>2 157</td><td></td><td></td><td>2 491</td><td></td></tr><tr><td>Télédetection (Mi-Juillet) (4)</td><td></td><td>2 119</td><td></td><td></td><td>2 547</td><td></td></tr><tr><td>Ecart [(4)-(3)] en ha</td><td></td><td>-38</td><td></td><td></td><td>56</td><td></td></tr><tr><td></td><td>Ecart [(4)-(3)]/(3)] en%</td><td></td><td>-1,75%</td><td></td><td></td><td>2,25%</td><td></td></tr></tbody></table>

Il est constaté que l'occupation du sol est impacté par la technique de l'irrigation localisée, plus efficiente, en substitut au système d'irrigation gravitaire, sur l'amélioration du taux de l'intensification cultural, qui est passé de 1,06 à 1,30. La méthode de la télédétection spatiale par indice NDVI s'est avérée pratique et plus précise, vu que l'écart relatif avec les données des statistiques terrestres n'a pas dépassé le faible taux de 2 % soit au sens d'augmentation ou au sens de diminution de la superficie mesurée (Tableau 85).
# VI-Projections alternatives à long terme des besoins en eau des cultures au périmètre de Tadla

## VI-1-Modalités des projections à long terme des précipitations et de l'évapo-transpiration dans le périmètre du Tadla.

Une vision futuriste en agriculture serait difficile pour de longue durée, notamment pour les prédicitions des besoins en eau des cultures qui sont tributaires aux choix des assolements culturaux, des aléas pluviométriques, du contexte technologique en évolution et aussi aux tendances du marché pour la commercialisation des productions agricoles. A cet effet et en considérant aussi l'implication effective des changements climatiques pour les décennies ultérieures, il serait judicieux de tabler sur une durée minimale d'une trentaine à une quarantaine d'années prochaines raisonnable pour le secteur agricole du périmètre de Tadla.

Ainsi, l'horizon de l'an 2050 est fixé pour mener les investigations nécessaires pour établir les projections des précipitations et en évapo-transpiration de référence, moyennant le simulateur fiable et usuel pour l'espace marocain à savoir le système MOSAICC (FAO et INRA-Maroc, 2015). A travers les manipulations répétitives du dernier système, moyennant ses seize alternatives possibles, la sélection la plus favorable s'est portée sur la combinaison prévisionnelle la plus optimiste qui associe le scénario RCP 4.5 avec le modèle CANESM.

Les résultats essentiels de la démarche prévisionnelle agro-climatique sont les suivants :

<table><tr><td>-Période de référence initiale</td><td>: 2010-2020.</td></tr><tr><td>-Période fixée des projections</td><td>: 2050-2060.</td></tr><tr><td>-Réchauffement thermique prétendu</td><td>: Ecart positif de 1 à 3°C</td></tr><tr><td>-Variabilité prévisionnelle des précipitations Annuelles</td><td>: Réduction de 11 mm sur le cumul pluviométrique annuelle.</td></tr><tr><td>-Variabilité prévisionnelle de l'évapo-transpiration de référence</td><td>: 1 mm/jour sur la moyenne ETo journalière annuelle.</td></tr></table>

## VI-2-Evaluations des projections à l'horizon 2050 des besoins nets en eau d'irrigation pour les cultures du périmètre du Tadla.

Les résultats précités des projections agro-climatiques ont été exploités dans les évaluations prévisionnelles des besoins nets en eau d'irrigation pour les principales cultures projetées à long terme dans le périmètre de Tadla. Le détail des calculs est présenté aux tableaux en annexes et le récapitulatif des prévisions en question sont données au Tableau 86.
<table><caption>Tableau 86 : Situation actuelle et projections à l'horizon 2050 des besoins en eau d'irrigation pour principales cultures du périmètre de Tadla (Base de calculs simulatifs des campagnes de 2014/2015 à 201720/18)</caption><thead><tr><th rowspan="3">Cultures en irrigué</th><th colspan="4">Situation actuelle</th><th colspan="4">Situation projetée</th><th colspan="2">Ecarts à l'ETM</th><th colspan="2">Ecarts sur besoins nets eau irrigation</th></tr><tr><th colspan="4">Période de référence: 2015-2018</th><th colspan="4">Horizon de: 2050</th><th colspan="2">[(2) / (1)]</th><th colspan="2">[(3) / (4)]</th></tr><tr><th>ETo (mm)</th><th>ETM (mm) (1)</th><th>Pluie (mm)</th><th>Besoins nets eau irrigation (mm) (3)</th><th>ETo (mm)</th><th>ETM (mm) (2)</th><th>Pluie (mm)</th><th>Besoins nets eau irrigation (mm) (4)</th><th>Volume (mm)</th><th>%</th><th>Volume (mm)</th><th>%</th></tr></thead><tbody><tr><td><strong>Céréales d'automne</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Blés</td><td>666,88</td><td>560,72</td><td>239,90</td><td>445,30</td><td>849,86</td><td>714,57</td><td>230,78</td><td>617,79</td><td>153,85</td><td>27,44%</td><td>172,49</td><td>38,73%</td></tr><tr><td><strong>Cultures industrielles</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Betterave à sucre</td><td>672,34</td><td>677,40</td><td>179,60</td><td>575,83</td><td>856,82</td><td>863,27</td><td>172,77</td><td>734,88</td><td>185,87</td><td>27,44%</td><td>159,05</td><td>27,62%</td></tr><tr><td><strong>Cultures fourragères</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Luzerne</td><td>1 330,25</td><td>959,19</td><td>289,20</td><td>845,09</td><td>1 695,25</td><td>1 222,37</td><td>278,20</td><td>1 105,27</td><td>260,18</td><td>30,79%</td><td>260,18</td><td>30,79%</td></tr><tr><td>Mais ensilage</td><td>503,83</td><td>544,14</td><td>89,10</td><td>491,82</td><td>642,07</td><td>693,44</td><td>85,71</td><td>642,10</td><td>149,30</td><td>27,44%</td><td>150,28</td><td>30,56%</td></tr><tr><td><strong>Ariculture fruitière</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Agrumes</td><td>1 238,84</td><td>902,31</td><td>378,60</td><td>701,07</td><td>1 581,49</td><td>1 151,94</td><td>364,20</td><td>888,49</td><td>249,63</td><td>27,67%</td><td>187,43</td><td>26,73%</td></tr><tr><td>Olivier</td><td>1 238,84</td><td>664,56</td><td>378,60</td><td>504,21</td><td>1 581,49</td><td>848,40</td><td>364,20</td><td>715,49</td><td>183,85</td><td>27,66%</td><td>211,28</td><td>41,90%</td></tr><tr><td><strong>Maraichage</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pomme de Terre d'hivers</td><td>518,26</td><td>596,44</td><td>78,90</td><td>560,48</td><td>660,46</td><td>760,10</td><td>75,90</td><td>714,46</td><td>163,65</td><td>27,44%</td><td>153,98</td><td>27,47%</td></tr></tbody></table>

A l'horizon 2050, les besoins nets annuels en eau pour les principales cultures, en conduites sous le système d'irrigation gravitaire, seraient dans la gamme de **618 à 1 105 mm/an** (équivalents de 6 180 à 11 050 m3/ha/an) respectivement en des valeurs extremes pour les céréales d'automne et pour la luzerne. Les taux d'accroissement projetés de ces besoins hydriques par rapport à la situation de référence (période : 2014-2018) seraient dans un intervalle de **27 à 39 %** en des cas de limites inférieure et supérieure respectives des cultures du maraichage et des plantations agrumicoles. Les explications résident aux changements climatiques, qui s'impacteraient probablement en des augmentations effectives de l'ETM des cultures, qui ont été estimées pour taux ascendants simulés de **27 à 31 %** en fonction des cultures considérées (Tableau 86).

## VI-3-Projections à l'horizon 2050 des exigences en volumes d'eau d'irrigation par hectare assolé dans le périmètre irrigué de Tadla.

Compte-tenu des données antérieures abouties, les projections des besoins en eau d'irrigation des cultures du périmètre de Tadla, sont selon trois situations, à savoir la première situation de départ en adoptant les campagnes référentielles de 2014 à 2018, suivi du 1er scenario des projection à l'horizon de l'an 2050 avec le maintien des cultures en conduite exclusivement sous le système
d'irrigation gravitaire sans sa reconversion et achevé avec le 2ème scenario en des projections à l'an 2050 avec spécificité de la reconversion totale du système d'irrigation gravitaire par le système d'irrigation localisée (plus efficient).

Les données détaillées de ces situations alternatives sont portées dans les Tableaux 87-89 et les analyses synthétiques sont les suivantes :

* Situation de départ :

Les besoins bruts en eau d'irrigation pour les cultures en zone du réseau hydro-agricole avec le système d'irrigation gravitaire sont évalués à la moyenne de 8 356 m3 par hectare assolé, dont la contribution avec les eaux des barrages hydrauliques desservant le périmètre est pour volume de 6 083 m3/ha équipé, soit le taux de la couverture à 72,8% de ces besoins en eau.

* Situation projetée du 1er scenario à horizon de l'an 2050 :

Les exigences en eau d'irrigation seraient pour un volume de 11 283 m3/ha assolé avec l'accroissement de 35% par rapport à la situation de départ.

* Situation projetée du 2ème scenario à horizon de l'an 2050 :

A ce sujet, il est posé l'hypothèse qui consiste à l'achèvement total à l'horizon de l'an 2050 du processus déjà lancé au périmètre de Tadla depuis l'année 2015 en vue de la reconversion du système d'irrigation gravitaire en système d'irrigation localisée. A l'horizon de l'an 2050, il est prétendu le volume d'eau nécessaire à l'irrigation des cultures pour l'hectare assolé d'une moyenne de 7 831 m3, soit une diminution de 6 % par rapport à la situation actuelle de départ.

Toutefois, il persisterait encore le déficit à couvrir par la dotation en eau de barrage, mais avec un volume plus réduit d'une moyenne de 1 748 m3/ha, ce qui représente le taux de 22% du total dû en eau, contre un taux de 27 % exprimé par le système d'irrigation gravitaire en sa situation de départ (année : 2018). Ceci serait bien attendu en tenant compte des incidences à long-terme des changements climatiques au niveau du périmètre de Tadla (zone d'étude).
<table><caption>Tableau 87 : Situation de départ des besoins en eau d'irrigation pour les principales cultures conduites sous le système d'irrigation gravitaire en période de 2014-2018 sur la globalité du périmètre de Tadla</caption><thead><tr><th>Cultures</th><th>Superficie par culture<br>(en ha)</th><th>Besoins unitaires nets en eau irrigation<br>(en m3/ha)</th><th>Effice du système d'irrigation<br>en %</th><th>Besoins unitaires bruts en eau irrigation<br>(en m3/ha)</th><th>Besoins totaux bruts en eau irrigation par culture<br>(en Millions m3)</th><th>Besoins totaux bruts en eau irrigation<br>(en Millions m3)</th><th>Volume total des besoins en eau irrigation à l'hectare assolé<br>(en m3/ha)</th><th>Volume total d'eau irrigation accordée<br>(en Millions m3)</th><th>Volume accordé en eau à l'hectare assolé<br>(en m3/ha)</th><th>Volume déficit en eau à l'hectare assolé (pompage)<br>(en m3/ha)</th><th>Volume total du déficit en eau (pompage)<br>(en Millions m3)</th></tr></thead><tbody><tr><td>Céréales d'automne</td><td>33 416</td><td>4 453</td><td>66%</td><td>6 722</td><td>224,6</td><td rowspan="6">735,0</td><td rowspan="6">8 355,6</td><td rowspan="6">535</td><td rowspan="6">6 083</td><td rowspan="6">2 273</td><td rowspan="6">199,9</td></tr><tr><td>Betterave à Sucre</td><td>6 575</td><td>5 758</td><td>64%</td><td>8 950</td><td>58,8</td></tr><tr><td>Luzerne</td><td>26 123</td><td>8 451</td><td>74%</td><td>11 481</td><td>299,9</td></tr><tr><td>Mais Ensilage</td><td>50</td><td>4 918</td><td>64%</td><td>7 680</td><td>0,4</td></tr><tr><td>Agrumes</td><td>876</td><td>7 011</td><td>70%</td><td>10 000</td><td>8,8</td></tr><tr><td>Oliviers</td><td>20 703</td><td>5 042</td><td>74%</td><td>6 800</td><td>140,8</td></tr><tr><td>Maraichage</td><td>221</td><td>5 605</td><td>72%</td><td>7 736</td><td>1,7</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<caption>Tableau 88 : Projections à long terme des besoins en eau d'irrigation pour les principales cultures sur la globalité du périmètre de Tadla en comparaison avec la situation actuelle de référence aux campagnes de 2014/2015 à 2017/2018</caption>

Situation projetée au scenario 1 : Cultures maintenues sous le système d'irrigation gravitaire en période de 2040-2050

<table><thead><tr><th>Cultures</th><th>Superficie par culture</th><th>Besoins unitaires nets en eau irrigation</th><th>Efficace du système d'irrigation</th><th>Besoins unitaires bruts en eau irrigation</th><th>Besoins totaux bruts en eau irrigation par culture</th><th>Besoins totaux bruts en eau irrigation</th><th>Volume total des besoins en eau irrigation à l'hectare assolé</th><th>Volume total d'eau à accorder</th><th>Volume à accorder d'eau à l'hectare assolé</th><th>Volume déficit en eau prévu à l'hectare assolé (pompage)</th><th>Volume total prévisionnel du déficit en eau (pompage)</th></tr><tr><th></th><th>(en ha)</th><th>(en m3/ha)</th><th>en %</th><th>(en m3/ha)</th><th>(en Millions m3)</th><th>(en Millions m3)</th><th>(en m3/ha)</th><th>(en Millions m3)</th><th>(en m3/ha)</th><th>(en m3/ha)</th><th>(en Millions m3)</th></tr></thead><tbody><tr><td>Céréales d'automne</td><td>33 416</td><td>6 178</td><td>66%</td><td>9 325</td><td>311,6</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Betterave à Sucre</td><td>6 575</td><td>7 349</td><td>64%</td><td>11 422</td><td>75,1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Luzerne</td><td>26 123</td><td>11 053</td><td>74%</td><td>15 015</td><td>392,2</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mais Ensilage</td><td>50</td><td>6 421</td><td>64%</td><td>10 027</td><td>0,5</td><td>992,5</td><td>11 283,0</td><td>535</td><td>6 083</td><td>5 200</td><td>457,4</td></tr><tr><td>Agrumes</td><td>876</td><td>8 885</td><td>70%</td><td>12 673</td><td>11,1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Oliviers</td><td>20 703</td><td>7 155</td><td>74%</td><td>9 649</td><td>199,8</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Maraichage</td><td>221</td><td>7 145</td><td>72%</td><td>9 861</td><td>2,2</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<caption>Tableau 89 : Projections à long-terme des besoins en eau d'irrigation pour les principales cultures sur la globalité du périmètre de Tadla en comparaison avec la situation actuelle de référence aux campagnes de 2014/2015 à 2017/2018.</caption>

Situation projetée au scenario 2 : Equipements des parcelles des cultures en l'après reconversion totale du système d'irrigation gravitaire en système d'irrigation localisée à l'horizon de l'an 2050

<table><thead><tr><th>Cultures</th><th>Superficie par culture</th><th>Besoins unitaires nets en eau irrigation</th><th>Efficace du système d'irrigation</th><th>Besoins unitaires bruts en eau irrigation</th><th>Besoins totaux bruts en eau irrigation par culture</th><th>Besoins totaux bruts en eau irrigation</th><th>Volume total des besoins en eau irrigation à l'hectare assolé</th><th>Volume total d'eau à accorder</th><th>Volume à accorder d'eau à l'hectare assolé</th><th>Volume déficit en eau prévu à l'hectare assolé (pompage)</th><th>Volume total prévisionnel du déficit en eau (pompage)</th></tr><tr><th></th><th>(en ha)</th><th>(en m3/ha)</th><th>en %</th><th>(en m3/ha)</th><th>(en Millions m3)</th><th>(en Millions m3)</th><th>(en m3/ha)</th><th>en Millions m3</th><th>(en m3/ha)</th><th>(en m3/ha)</th><th>(en Millions m3)</th></tr></thead><tbody><tr><td>Céréales d'automne</td><td>33 416</td><td>5 423</td><td>92%</td><td>5 878</td><td>196,4</td><td rowspan="7">688,8</td><td rowspan="7">7 830,9</td><td rowspan="7">535</td><td rowspan="7">6 083</td><td rowspan="7">1 748</td><td rowspan="7">153,7</td></tr><tr><td>Betterave à Sucre</td><td>6 575</td><td>5 952</td><td>91%</td><td>6 535</td><td>43,0</td></tr><tr><td>Luzerne</td><td>26 123</td><td>10 354</td><td>94%</td><td>10 973</td><td>286,7</td></tr><tr><td>Mais Ensilage</td><td>50</td><td>5 985</td><td>94%</td><td>6 369</td><td>0,3</td></tr><tr><td>Agrumes</td><td>876</td><td>8 273</td><td>92%</td><td>8 964</td><td>7,9</td></tr><tr><td>Oliviers</td><td>20 703</td><td>6 806</td><td>92%</td><td>7 399</td><td>153,2</td></tr><tr><td>Maraichage</td><td>221</td><td>6 307</td><td>96%</td><td>6 539</td><td>1,4</td></tr></tbody></table>
DISCUSSIONS

Les résultats constatés des fréquences respectives des précipitations annuelles en des années sèches, normales et humides de : 33%; 35% et 32%, obtenus par des méthodes d'intervalle de confiance et d'indice IPS sont trouvés relativement comparables par rapport aux résultats antérieus traitant la même thématique, couvrant soit partiellement ou soit totalement la même zone d'étude, selon Jouilil et al. (2013) sur la période de 1975 à 2009 et selon Krimissa et al. (2017) sur la période d'observations de 1982 à 2015.

Le constat de la tendance du régime pluviométrique en baisse, décelé par le test de Mann-Kendall et par le test de Sen's, se confirme par les travaux similaires sur large étendue de la zone du bassin hydraulique d'Oum Er Rbia par Jouilil et al. (2013) et sur le grand territoire du Maroc septentrional par Sebbar et al. (2011), y compris notre zone d'étude.

La précision de l'étude de tendance des précipitations par méthode de pondération avec moyenne mobile au filtre passe-bas de Hanning d'ordre 2, ayant identifié une période pluviométrique transitoire de 1980 à 1987 avec des alternances irrégulières d'années sèches, normales et déficitaires, est bien confirmée par Krimissa et al. (2017) sur l'intervalle-temps de 1982 à 2006, concernant la tendance des pluies dans la zone Est du bassin hydraulique Oum Er Rbia, pour de ratios comparatifs des années sèches par rapport aux années à pluviométrie normale.

Le taux de réduction de la pluviométrie constaté par notre étude de 36% est plus élevé par rapport à celui rapporté par Sebbar et al. (2011) de 8 à 15 % sur la même zone d'étude, qui s'explique par les décalages sur la période d'observations et aussi pour les étendues différentes de ses périodes identifiées d'avant-rupture de 1935 à 1977 (43 années) et en après-rupture de 1978-2004 (27 années). Par contre en notre propre cas, la période d'observations de 1958 à 2018 s'est éclatée en deux périodes d'avant et d'après rupture de séries chronologiques des pluies de 1958 à 1978 (20 années) et de 1979 à 2018 (40 années).

A propos des ruptures chronologiques pluviométriques détectées qui sont centrées sur l'année de 1978, suite à nos tests d'homogénéité, elles se ne confirment pas avec celles de Krimissa et al. (2017) portant sur une période observée plus récente de 1983 à 2015 sur l'Est du bassin hydraulique Oum Er Rbia. Nos résultats sont comparables avec ceux de Sebbar et al. (2011) à propos de la rupture en régime des pluies sur la décennie des années 1980. Les mêmes tendances à la baisse de la pluviométrie ont été déjà rapportées au Nord-Ouest algérien entre fin des années 1960 et début des années 1970 (Taibi et al., 2013 ; Ghenim et al., 2010). Egalement les cas en Afrique Occidentale, les ruptures ont été a priori identifiées au haut bassin versant du Sénégal
(Faty et al., 2017) et aux zones du Centre et du Sud de la Côte d'Ivoire (Fossou et al., 2015 ; Ahoussi et al., 2013).

Notre interpolation spatiale de la répartition pluviométrique dans la zone d'étude, avec ses quatre sous-zones homogènes, concorde effectivement avec la carte des précipitations annuelles du Maroc établie par El Ajhar et al. (2018) pour la période de 1960-2015, avec un semblable positionnement du périmètre de Tadla (zone d'étude) entre les isohyètes de 200 à 300 mm/an.

A propos des périodicités intra-annuelles des précipitations, nos constats de tendance de la diminution des pluies par mois entre Novembre et Juin, concordent avec ceux de Hanchane (2013) pour les régions continentales marocaines à climat semi-aride (y compris la zone d'étude) et précisément pour les deux mois de Mars et de Juin. En outre, des travaux de recherche analogues de Norrant et Douguédroit (2004) ont précédemment décelé de telle tendance aux mois d'Octobre, Novembre et Mars au Sud français sur le large bassin méditerranéen du fleuve de Rhône. A citer aussi des similitudes observées de la diminution progressive des pluies mensuelles au Nord du territoire algérien (Hebal et Remini, 2012).

Quant à notre seul constat de la tendance significative d'accroissement pluvieux au mois de Septembre au début de l'année agricole, il n'a été nulle part signalée sur le territoire marocain ou aux niveaux de pays méditerranéen de proximité selon de la bibliographique en la matière disponible et consultée.

Concernant les saisons de l'année, sachant que nos travaux ont montré une réduction significative des pluies aux trois trimestres: Décembre-Janvier-Février, Mars-Avril-Mai et Juin-Juillet-Août, le constat se confirme avec résultats de Driouech (2010) sur la globalité du territoire marocain pour uniquement la saison d'hivers et pour les quatre saisons prises individuellement d'une manière restreinte à l'extrême Nord marocain en région de Tanger et à en plaine du Saïs (région de Meknès-Fès). La même tendance a été relevée par Sebbar et al., (2011) en zone du Centre du Maroc pour l'hivers et le printemps, mais à l'inverse avec une légère augmentation des pluies en saison d'automne. Les faits constatés sont identiquement relevés sur l'espace algérien pour le Centre et Nord-Ouest rapportés par Hebal et Remini (2012) et au Nord selon Touhami (2017).

Considérant la périodicité du semestre, dont nos travaux affichent le repli pluviométrique, est conforme à ce qui a été détecté pour l'hivers étendu (Octobre-Mars) suite à l'étude climatique de Driouech (2010) sur la globalité du territoire marocain.

Au pas annuel des précipitations qui est en cadence de la régression selon notre étude sur le périmètre de Tadla, a été aussi rapportée sur le territoire marocain presque similaire par Driouech (2010). Par ailleurs, cette tendence est qualifiée aux zones du Centre du Maroc. par Sebbar et al.
(2011) d'une manière irrégulière, voir avec des contrastes en amplitude (augmentation ou diminution). Aussi, la tendance à la baisse des cumuls annuels des pluies est soulevée en de diverses zones marocaines semi-arides à caractère continental par Hanchane (2013) au Nord algérien par Hebal et Remini (2012) et par Taibi et al. (2013) au Centre-Nord algérien par (Touhami, 2017) et également en Afrique occidentale auprès de la vaste région agricole d'Abidjan-Agboville en Côte d'Ivoire par Ahoussi (2013).

Par ailleurs, l'approche menée de la régionalisation de tendance périodiques des précipitations par la méthode de l'ACP, ayant abouti à deux zones géographiques pluvieuses homogènes, trouve sa confirmation concernant Sebbar et al. (2011), qui ont identifié cinq zones climatiques au centre du Maroc en une tendance négative des précipitations de l'hiver et du printemps, où notre zone d'étude s'y intègre complètement. D'une manière analogue, pour la zone du bassin versant du Tensift au voisinage de la zone d'étude, comportant quatre sous-zones à contrastes pluvieux annuels plus marquants, rapportés par Salama (2011).

Aussi, il est à noter pour la région du Nord algérien un gradient de dessèchement pluviométrique en allant d'Est en Ouest (Hebal et Remini, 2012), mais dans le sens inversé par rapport à notre constat au présent travail de recherche, avec un gradient de baisse de la pluviométrie annuelle dans la direction Ouest-Est. En signalant que l'espace algérien se subdivise sur lui-même en six régions pluviométriques homogènes (Taibi et al., 2013). Comparativement, au sujet de la pluie mensuelle, l'espace du Sud-Est méditerranéen français comporte trois régions homogènes (Norrant et Douguédroit, 2004), contre deux sous-zones qu'on a identifiées sur notre zone d'étude.

A propos de la périodicité des pluies du semestre, les résultats de nos travaux se confirment avec ceux aboutis en région algérienne semi-aride du bassin versant de la Macta concernant la baisse pluviométrique graduelle des précipitations en se déplaçant du Nord vers le Sud (Chibane et al., 2015). D'une manière conforme, nos constats rejoignent les cas soulevés en régions du Centre et du Nord-Ouest algériens présentant eux aussi trois premières saisons de l'année agricole avec la baisse pluvieuse et des sous délimitation en cinq sous-espaces géographiques homogènes pour la tendance pluviométrique saisonnière (Touhami, 2017).

Il y a une concordance dans les résultats trouvés par rapport aux données sur l'Afrique Occidentale à propos de la régionalisation cartographiée de la tendance pluviométrique par la méthode d'ACP. Ainsi, trois grandes régions ont été auparavant jalonnées pour les saisons sèches de l'année sur les territoires de divers pays de cette grande zone africaine (Khodja, 1998). Le cas modèle de la Côte d'Ivoire avec ses trois zones de la région d'Abidjan-Agboville qui se caractérisent par le déclin annuel progressif des précipitations (Ahoussi, 2013). Au niveau de ce
dernier pays, avec une cadence de pluie mensuelle vulnérable à la baisse sur l'étendue bassin versant du N'Zi, il est récemment identifié de même deux zones homogènes pluviométriques (Koudou, 2015).

Après la rupture chronologique thermique repérée en campagne 1994/1995 au périmètre du Tadla, qui marque le passage de la température annuelle moyenne de 18,7 à 19,7 °C en 2017/2018 en enregistrant une anomalie du réchauffement de 1 °C ; il y a coïncidence avec la même date et pour le même accroissement de la température à l'échelle mondiale (Lenssen et al., 2019).

L'évapo-transpiration de référence qui est l'indicateur adopté pour le rapprochement effectif au comportement du végétal en realtion avec ses exigences en eau (Doorebos et Pruitt, 1977; Doorebos et Kassam, 1979; Alien et al. 1998), ayant évolué dans le périmètre du Tadla après rupture chronologique en 1977/1978, avec un taux de progression de 12,5 % ; est trouvé plus proche du taux de 11,4% constaté en zone du Centre tunisien souis le même climat semi-aride (Hajri, 2006).

Le bilan hydrique (entre ETo et les précipitations efficaces) qui évolue en zone d'étude après la rupture chronologique de 1977/1978 pour sa valeur déficitaire de 805 à 1 093 mm/an, est comparable dans le même sens dev l'évolution en région méditerranéenne de la Sardaigne en Italie durant la période de 1951 à 2010, où le déficit d'eau s'etat aussi accentué de 899 à 2 544 mm/an (Pulina, 2012). En se référant à la classification bioclimatique avec ses six étages pour la zone méditerranéenne (E.M.) basée sur le bilan hydrique (Gammar et Hlaoui, 2012), le périmètre de Tadla qui fait contraste entre l'avant et l'après rupture de 7 à 10 mois en non satisfaction de point de vue hydrique, évolue alors dans une translation de la zone d'étude du 3ème Etage E.M. Semi-aride (sous-étage supérieur) au 5ème Etage E.M. Aride.

A l'issue des diagrammes ombro-thermiques de Gaussen obtenus après la rupture sur la tendance de la température moyenne en année de 1994/1995 et ayant montré que l'extension de la période sèche passait de 6 mois à 8 mois, implique selon la classification usuelle afférente (Gaussen et Bagnouls, 1957), que notre zone d'étude bascule actuellement du climat thermo-méditerranéen vers le climat xérothermo-méditerranée, dans un sens d'accroissement de l'aridité du milieu agro-écologique.

L'indice d'aridité de De Maritonne (I.M.) au périmètre de Tadla a régressé après la rupture tendancielle en 1987/1988 de 11,2 à 8,6 points. En se référant à la classification climatique selon cet indice (De Parcevaux et al, 1990), la zone d'étude qui a été répertoriée auparavant au-dessus de la limite inférieur de la classe bioclimatique semi-aride (I.M. : de 10 à 20 unités) se déplace vers la limite supérieure de la classe bioclimatique aride (I.M : de 5 à 10 unités).
Avant la rupture thermique de 1994/1995, le périmètre de Tadla avait un quotient d'Emberger de l'ordre de 28 unités et qui se positionne après rupture temporelle à la valeur 38 unités. Par conséquent et avec la référence à la classification bioclimatique tablée sur ce quotient (De Parcevaux et al, 1990), le périmètre (zone d'étude) régresse depuis la limite supérieure vers le centre de la classe aride, sachant l'intervalle en Q.E pour cette classe est de 12 à 30 ; approuvant encore une accentuation du phénomène de l'aridité depuis plus de deux décennies.

L'efficience agronomique de l'utilisation d'eau (EUE) par les spéculations de nos essais culturaux en station expérimentale, relevée pour les céréales d'automne dans le périmètre de Tadla à des niveaux de 0,96 à 3,50 Kg/m3 d'eau, s'avère plus améliorée par rapport aux données des enquêtes sur ce même périmètre de Bouaziz et Belabbes (2002), pour des valeurs de l'EUE de 0,70 à 1,80 Kg/m3. Quant à la betterave à sucre, les EUE constatées respectivement pour les systèmes d'irrigation gravitaire, localisée et aspersion sont de 0,65 à 2,48 ; de 2,20 à 2,55 et de 1,80 à 2,06 Kg/m3 et qui sont comparables aux résultats des essais culturaux réalisées sur le même site en campagnes de 2011/2012 à 2013/2014 avec des EUE moyennes correspondantes de 1,16 à 1,69 ; de 2,28 à 3,26 et de 1,82 à 2,48 Kg/m3 (ORMVA du Tadla, 2019).

Pour la culture de la luzerne, son EUE trouvée en expérimentation sur le système d'irrigation gravitaire de 1,10 Kg/m3 d'eau, se compare aux alentours des données en la même zone par Bouzzama et al. (2015) en une EUE de 1,44 à 1,46 Kg/m3. Au sujet de la culture du maïs pour l'ensilage, en tenant compte de teneurs en matière sèche sur les rendements en frais de 30 à 35 % selon les normes en vigueur (Ait Houssa et al., 2008), l'EUE de nos essais culturaux qui varie entre 32,77 et 83,48 Kg/m3, s'avère plus importante par rapport aux résultats des expérimentations de Bouzzama et al. (2013) dans la même zone en 2008/2009 et en 2009/2010 avec EUE déterminée sur matière sèche de la productivité culturale de 2,85 à 2,99 Kg/m3. L'explication de cet écart réside à la faiblesse des volumes d'eau d'irrigation utilisés dans nos expérimentations (de 1 052 à 1 056 m3/ha) contre les volumes d'eau afférents aux auteurs précités (de 5 350 à 6 190 m3/ha).

La valorisation économique de l'eau d'irrigation constatée dans notre étude est plus améliorée en cas d'adoption du système d'irrigation localisée par rapport au système d'irrigation gravitaire ; avec des données pour ce dernier système plus proches de celles rapportées par Lionboui (2018) pour la culture du maïs fourrager (1,3 Dhs/m3), pour les plantations des agrumes (3,5 Dhs/m3) et pour les vergers oliviers en irrigué (1,4 Dhs/m3).

A l'issue de nos projections de l'ETo au niveau du périmètre de Tadla, où la moyenne annuelle passerait de 3,6 à 4,6 mm/jour à l'horizon de l'an de 2050, soit à un écart journalier de 1
mm/jour, correspondant au taux d'accroissement de 27%, obtenu selon le Modèle MOSAICC (FAO et INRA-Maroc, 2015) ; il y a une convergence avec les projections du projet « Explore » à l'horizon de l'an de 2070 qui projette pour la zone modèle des bassins de Rhône-Méditerranée et de la Corse, à l'Ouest-méditerranéenne, une augmentation de l'ETo de l'ordre de 20% à l'horizon de 2046 à 2065 (MEDDE-France, 2012).

Les besoins en eau d'irrigation pour les cultures du périmètre de Tadla conduites sous le système d'irrigation gravitaire, qui sont selon nos projections à l'horizon 2050 de 8 356 à 11 283 m3/ha assolé, soit un écart supplémentaire exigé en eau de 2 927 m3/ha qui reste dans les grandeurs des projections pour la zone du bassin versant du Rhône au Sud de la France sur rive Nord-Méditerranéenne pour de tels besoins qui augmenteraient à l'horizon 2070-2100 en un écart moyen annuel dans la gamme de 730 à 2 560 m3/ha (Dayon, 2015).

Par ailleurs, les projections des besoins en eau pour l'irrigation des cultures du périmètre de Tadla en conduites sous système d'irrigation localisée à l'horizon 2050 et avec l'alternative d'avoir reconverti la totalité de la superficie actuelle en cette nouvelle technique d'irrigation en substitut à l'actuel mode moins efficient de l'irrigation gravitaire, en prédiction de volume moyen annuel de 7 831 m3/ha, soit une diminution relative de 6 % par rapport aux besoins initiaux de 2018. Notre cas est plus optimiste par rapport à la référence à l'échelle mondiale, où il est prétendu à l'horizon 2050 la réduction de 17 % des ressources en eau à allouer à l'irrigation en agriculture, qui serait à l'origine probable de la forte concurrence des autres secteurs économiques prépondérants, en premier lieu les activités industrielles et les productions énergiques, en étroite relation avec l'augmentation soutenue à long terme de la démographique humaine du globe terrestre (OCDE, 2012).
CONCLUSION GENERALE ET PERSPECTIVES

La caractérisation des changements climatiques en relation avec l'eau en agriculture en milieu à
bioclimat semi-aride vulnérable à la raréfaction hydrique, a été appréhendée sur le modèle du
périmètre de Tadla. Sachat que ce dernier espace, est un support d'excellence pour une
agriculture potentielle irriguée en mode intensif, cependant il est assujettie à des contraintes
récurrentes des disponibilités en eau pour l'irrigation.

En effet, ce périmètre de référence, situé en plein centre continental marocain, a fait l'objet d'un
système intégré d'étude de la tendance agro-climatique se rapportant aux : variabilité des
distributions des données enregistrées, analyses d'homogénéité chronologiques des séries des
mesures, analyses multivariées de la pluviométrie, évolution des valeurs des indices climatiques
usuels et aussi exploitation des cartographies spatio-temporelles des paramètres afférents au
climat. Ceux-ci ont permis de cerner la démarche scientifique sur le facteur eau dans son aspect
climatique, en se focalisant sur les indicateurs annuels et saisonniers de la pluviométrie, sur la
température ambiante, sur l'évapo-transpiration de référence ainsi que sur les indices
conventionnels pour le suivi de l'aridité d'un espace agro-écologique.

En l'occurrence et à travers l'exécution des analyses précitées, il a été concrétisé pratiquement la
convenance et les inter-complémentarités des diverses méthodes adoptées. Il s'agit en effet des
tests de rupture chronologique, non paramétriques et fiables (Test de Pettitt, test de Lee et
Heghinian, test Buishand, segmentation d'Hubert). Aussi ont été impliquées les méthodes
statistiques tendancielles sur les séries temporelles (Test de Mann-Kendall couplé au test de
Sen's), la méthode des intervalles de confiance sur l'évolution des pluies. A cet effet, Il a été
question d'emploi des indices pluviométriques à usage universel, tel le cas pour l'Indice
Standardisé des Précipitations (ISP) en ses deux versions, l'une de base et l'autre plus finalisée
selon la méthode de la pondération pluviométrique en des moyennes mobiles avec le filtre passe-
bas de Hannig d'ordre 2.

L'achevement du processus méthodologique se trouve contraint d'être complété par l'étude de
la saisonnalité intra-annuelle des pluies par la méthode de l'ACP. En vue de pousser davantage
les analyses des résultats des méthodes sus-indiquées, il s'est aussi imposé d'ordre opportun de
faire des comparaisons entre les situations de l'avant et en après-rupture des séries agro-
climatiques, concernant les évolutions spatiale et dans le temps du phénomène de l'aridité,
moyennant l'indice de De Maritonne, l'indice d'Emberger et les Diagrammes Ombro-thermiques de Gaussen.

Les résultats aboutis ont montré qu'au niveau du périmètre de Tadla, il y a une catégorisation des pluviométries annuelles en des fréquences équitables entre les années déficitaires (33% des cas), les années normales (35%) et les années excédentaires (32%). Il est essentiellement retenu pour cette zone semi-aride modèle que les précipitations sont dans une perpétuelle tendance à la baisse. En effet, il est montré l'existence de deux périodes distinctes, séparées par l'année de clivage chronologique de 1978. Le taux de réduction pluvieux après cette rupture est évalué au taux moyen de 36 %. Plus encore, la configuration géographique de la variabilité spatio-temporelle pluviométrique a mis en évidence un gradient d'accroissement des précipitations dans le sens Ouest-Est pour le périmètre de Tadla. En l'occurrence, après la rupture précitée, un front de dessèchement pluvieux est identifié et qui touche d'une manière uniforme le territoire régional observé, avec pour constat plausible le déplacement des isohyètes selon une ampleur déterminée de $\Delta=135$ mm/an.

Quant à la saisonnalité des précipitations, sa variabilité qui a été élucidée au périmètre de Tadla, a révélé que pour la majorité des moins de l'année, au nombre de 9 sur 12, affiche des fluctuations significatives pour les pluies mensuelles. Ce qui a été constaté dans le sens de la hausse pour l'unique mois de Septembre, alors que pour les autres huit mois, pris au singulier, de Novembre à Juin, le sens est inversé à la basse pluviométrique. Concernant, les trois mois restants d'Octobre (en mi-saison automnale) et de Juillet-août (en fin saison estivale), le niveau des précipitations persistent à l'état constant entre la fin dés années cinquante du dernier siècle jusqu'à la fin de l'année de 2018. De point de vue spatio-temporelle, la saisonnalité des précipitations du périmètre étudié est spécifié en deux groupes homogènes, dont le premier groupe est le représentant de sa partie Ouest, identifiable aux deux variables-mois de Septembre et de Juin, le deuxième groupe est la caractéristique de sa partie Est pour la large combinaison suivante en : des deux variables pluvieuses mensuelles de Décembre et d'Avril, de la variable trimestre « Mars-Avril-Mai » et aussi en intégrant des deux variables-semestre en humide (Septembre-Février) et en sec (Mars-Août).

A l'instar des précipitations, des ruptures chronologiques ont été déterminées pour les autres paramètres agro-climatiques à savoir : la température minimale en 1989, la température maximale en 1997, la température moyenne en 1995, l'évapo-transpiration de référence en 1978 et l'indice d'aridité de De Maritonne en 1988. Ainsi, Il a été abouti d'une manière plus avancée à caractériser les manifestations des changements climatiques, à l'issue des ruptures
chronologiques sus-mentionnées, en indiquant essentiellement le réchauffement de l'air de 1°C. En conséquenc, il est constaté une accentuation du déficit du bilan hydrique au sol au périmètre de Tadla, plus propice entre les mois de Mars à Août, en impactant alors la régression du cumul annuel de ce blan de la valeur de - 805 à la valeur de - 1 093 mm/an. Il est, en effet, observé un allongement anormal de la période sèche annuelle de 43 jours supplémentaires et aussi l'augmentation relative de l'aridité du périmètre, sachant que l'indice de De Maritonne a évolué à la baisse de 11,2 à 8,6 unités et le Quotient d'Emberger est de même passé de 27,8 à 18,7 unités.

En conséquence, les relations entre les changements climatiques et le facteur eau sont effectivement argumentés par les résultats précités. Dans cette optique, les besoins en eau d'irrigation pour les principales cultures conduites sous l'actuel système d'irrigation gravitaire moins efficient et plus dominant au périmètre de Tadla, subiraient à long terme (à horizon de l'an 2050), si aucune mesure structurante n'est prise, une augmentation excessive srait probablement en volume éxigé en eau de 8 356 à 11 283 m3 rapporté à l'hectare assolé, correspondant à un important écart vers la hausse au taux de 35%.

Il se profile alors la crise latente à satisfaire des exigences de plus en plus accrues en eau d'irrigation des cultures, ce qui se conjugue avec la tendance de la raréfaction de ressources en eaux mobilisables sur le milieu agricole, fortement tributaires à un régime pluviométrique régional en sa phase de dérèglement défavorable. La situation est sans doute susceptible d'empirer davantage avec le maintien du système d'irrigation classique gravitaire, plus gaspillant d'eau. Il serait donc opportun de concenter ultérieurement les efforts pour achever à terme le projet lancé de la reconversion de l'ancien système d'irrigation classique par le système d'irrigation localisée, plu efficient et plus pratique dans son plilotage à l'échelle parcellaire. Le processus de cette reconversion engagée, qui est déjà lancé depuis quelques années dans le périmètre de Tadla, vu son économie effective éventuel en eau d'irrigation, est considéré alors d'un rôle prometteur dans la rationalisation des besoins en eau d'irrigation pour les prochaines décennies. Ceci serait accéssible malgré les revers des changements climatiques attendus, argumenté par la projection établie qui prétend une fois le projet entrepris serait mené jusqu'au bout, il y aurait une réduction d'environ 6% des besoins en eau d'irrigation des cultures prioritaires à l'horizon 2050, en partant de la situation actuelle retenue comme référence de départ.

Les résultats du présent travail de recherche, à caractère ambitieux, constituent alors une assise consistente et diversifiée de la caractérisation bioclimatique de la région semi-aride, qui serait
plus utile à l'appui aux études de conception de programmes éventuels de la résilience aux changements climatiques, qui prônent le développement durable et la preservation de la biodiversité écologique.

Il s'agit en outre d'une contribution aux reflexions sur les perspectives des utilisations plus rationnelles et efficientes des ressources en eau disponibles dans des espaces plus vulnérables aux pénuries des eaux surtout d'origine pluviale, dans le souci majeur d'exploiter durablement la diversité des cultures praticables couramment sur les sols des périmètres agricoles potentiels similaires, de vocation d'agriculture irriguée et intensive, mais sous les risques énormes à encourir par les changements climatiques, notamment pour l'aspect stratégique de l'eau.

En guise de la clôture du document de la présente étude, il est recommandé de mener ultérieurement des travaux de recherche de continuité des investigations et des analyses déjà entreprises, surtout au sujet de la variabilité des précipitations journalières extrêmes et leurs impacts éventuels sur les systèmes culturaux plus axées sur les zones bioclimatiques vulnérables. Il serait aussi question d'entreprendre des activités scientifiques plus poussées sur les incidences actuelles et aussi à long terme du réchauffement atmosphérique sur les aspects physiologiques clés des végétaux, aussi bien pour les cultures usuellement pratiquées que pour les écosystèmes naturels, avec focalisation sur les facteurs de l'eau et de la nutrition minérale végétale.

Les finalités escomptées pour les études ultérieures recommandées, seraient dans la perspective de prodiguer des résiliences d'ordre pratique contre les méfaits du climat en phase de la profonde mutation. Les préoccupations majeures futures porteraient donc sur les éventualités nécessaires de la refonte du système de l'exploitation des spéculations stratégiques des productions végétales et sylvicoles sur les deux milieux pluvial et irrigué. Les attentes pour de ces travaux escomptés, seraient en vue de baliser les interventions de la pérennisation des périmètres agricoles potentiels, tels les cas vulnérables sur le territoire marocain et pour d'autres inombrables cas à l'au-delà des frontières. Ce serait en réponse pour la sauvegarde urgente d'espaces agricoles, de biotopes naturels, de domaines agro-forestiers et d'aires pastorales, les plus fragiles aux agressions changements climatiques.
REFERENCES BIBLIOGRAPHIQUES

ABDOU A, THIERRY L et ABOU A. (2008). Signification et usage de l'indice pluviométrique au Sahel. Revue Sécheresse 19: 22735.

ABIDINE O. M. M., EL ABOUDI A., INEJIH. C. A., O. SOULE O. A. (2017). Télédétection et SIG en appui à la caractérisation du climat en Mauritanie: Cas du Parc National du Diawling (PND) et sa zone périphérique. Communication au 30ème colloque de l'Association Internationale de Climatologie, du 03 au 07/07/2017 sous le thème : climat, ville et environnement. Université de Sfax (FLSH). 2017. pp : 423-428.

ABOUATALLAH A., SALGHI R., EL FADL A., AFFI N., GHNIZAR Y., ZARROUK A. and HAMMOUTI B. (2012). Impact assessment of dripper's distribution around the tree on soil moisture, roots and fruits growth of Citrus. *Der Pharma Chemica* 4, n° 5 (2012): 1969-1981.

ADEME (Agence de la transition écologique) (2020). un établissement public à caractère industriel et commercial (EPIC) français créé en 1991, site web : www.ademe.fr/

AGBANOU T., ABDOULAYE D., SABI OROU BOGO G. A., PAEGELOW M. et TENTE B. (2018). Variabilité pluviométrique et son impact sur le couvert végétal dans le secteur Natitingou-Noukombé au Nord-Ouest du Bénin, *Afrique Science : Revue Internationale des Sciences et Technologie* 14, n°3 (2018), pp : 182-191.

AHOUSSI K. E., KOFFI Y.B., KOUASSI A. M., SORO G., SORO N. et BIÉMI J. (2013). Étude de la variabilité hydro-climatique et de ses conséquences sur les ressources en eau du Sud forestier et agricole de la Côte d'Ivoire: cas de la région d'Abidjan-Agboville. Intern. *J. Pure & Ap. Bio-science Research.* 1(6): 3050 (2013).

AIT HOUSSA A., MOUTID S, BELBASRI M., HSAYNI M. et LOULTITI M. R. (2008). Productivité et rentabilité du maïs ensilage conduit en goutte à goutte aux sables de Larache. Bulletin mensuel d'information et de liaison : *Transfert de technologie en agriculture*, n°169, Octobre 2008, IAV Hassan II (Rabat).

AIT HOUSSA A., DRISSI S., ASFERS A., OUBAKI L., HADDOUCHE Z., AMLAL F. et CHRAIBI H. (2020). Éléments d'aide au contrôle de gestion et au management de l'agrumiculture au Maroc. *Revue Marocaine des Sciences Agronomiques et Vétérinaires* 8, n° 1 (2020).

AIT YACINE Z. (1995). Irrigation d'appoint du blé (*Triticum aestivum*) : Gestion de l'eau et techniques culturales adaptées. DESA Faculté des Sciences- Université Chouaib Doukkali El Jadida 1995.

AKOGNONGBE A., ABDOULAYE D., VISSIN E. W. et BOKO M. (2014). Dynamique de l'occupation du sol dans le bassin versant de l'Oueme à l'exutoire de Bétérou (Bénin). *Afrique Science : Revue Internationale des Sciences et Technologie* 10, n° 2 (2014), pp 228-242.

ALLEN R. G., PEREIRA L. S., RAES D. and SMITH M. (1998). Crop evapo-transpiration-Guidelines for computing crop water requirements- Ir. & drainage paper 56. FAO, Rome (Italy) 300, n° 9 (1998): D05-109.

ASSANI A. A. (1999). Analyse de la variabilité temporelle des précipitations (1916-1996) à Lubumbashi (Congo-Kinshasa) en relation avec certains indicateurs de la circulation atmosphérique (oscillation australe) et océanique (El Niño/La Niña). Science et changements planétaires. *Sécheresse* 10, no 4 (2000): 245-252.

BALAH B. et AMARCHI H. (2016). Variabilité des séries pluviométriques du bassin versant de la Seybouse du Nord-Est algérien. *Rev. Sciences et de la Technologie* 383, n° 3473 (2016) : 1-12.

BARONETTI A., GONZÁLEZHIDALGO J.C., VICENTESERRANO S.M., ACQUAOTTA F. et FRATIANNI S. (2019). Les épisodes de sècheresse sur la plaine du fleuve Po en Italie. Communication au 32ème colloque Internationale de l'AIC sous thème : « Le changement climatique, la variabilité et les risques climatique » ; Thessaloniki-Grèce, 29 mai au 1 juin 2019 ; Département de Météorologie et de Climatologie-Collège de Géologie Université d'Aristote de Thessaloniki ; Grèce, pp : 157-162.

BEDOUM A., BOUKA B. C., ALLADOUM M., ADOUMI I. et BAOHOUTOUL L. (2014). Impact de la variabilité pluviométrique et de la sécheresse au Sud du Tchad : Effet du changement climatique. *Revue Ivoirienne des Sciences et Technologie*, volume 23 (2014): 13-20.

BEKKAR Y., KUPER M., HAMMANI A., DIONNET M., et EL IAMANI A. (2007). Reconversion vers système d'irrigation localisée au Maroc. *Revue Hommes, terre et eaux* 137 (2007): 7-20.

BELGHITI M. (2008). Le programme national d'économie et de valorisation de l'eau en irrigation marocaine. Communication de l'Administration du Génie Rural au ministère de l'agriculture du Maroc, en journée mondiale de l'alimentation, organisée par la FAO à Rabat, date du 14/11/2008.

BENHADJ I., DUCHEMIN B., MAISONGRANDE P., SIMONNEAUX V., KHABBA S. et CHEHBOUNI A. (2012). « Automatic un-mixing of MODIS multi-temporal data for inter-annual monitoring of land use at a regional scale (Tensift, Morocco) ». *International journal of remote sensing* 33, N° 5 (2012): 1325-1348.

BERDAI H., KARROU M., CHATI M.T., BOUTFIRASS M. and BEKAOUI A. (2011). Irrigation water management in Morocco. Paper at the symposium on: "Improving water and land productivities in rainfed systems". Community Based Optimization of the Management of Scarce Water Resources in Agriculture in CWANA. Report, 2011, edited by ICARDA, Aleppo, (Syria), pp: 11-36.

BERRADA A. (2009). Assessment of drip irrigation in Morocco with particular emphasis on the plain of Tadla research grant report. US Fulbright Scholar, November 2008 April 2009 Research Scientist/Manager
Southwestern Colorado Res. Center Agr. Exper. Station & Department of Soil and Crop Sciences Colorado State University. Technical bulletin (Colorado Agricultural Experiment Station: LTB); LTB091, 2009.

BODIAN A., DACOSTA H. et DEZETTER A. (2011). « Caractérisation spatio-temporelle du régime pluviométrique du haut bassin du fleuve Sénégal dans un contexte de variabilité climatique ». *PhysioGéo. Géographie physique et environnement*, N° volume 5 (2011) : 107-124.

BOIS P. (1971). Une méthode de contrôle de séries chronologiques utilisées en climatologie et en hydrologie. Publication du Laboratoire de Meca. Flu. Univ. Grenoble I, mai 1971, 33p.

BOUAZIZ A. et BELABBES K. (2002). Efficience productive d'eau en irrigué au Maroc. *H.T.E.*,124 : 5772.

BOUAZZAMA B. et BAHRI A. (2009). Détermination des coefficients culturaux de la variété d'oranger « Maroc late » au périmètre irrigué Tadla ». Actes du Symposium International : « Agriculture durable en région méditerranéenne (AGDUMED) », pp : 14-16, 2009.

BOUAZZAMA B., DIMITRI D., BOUAZIZ A., RUELLE P. and MAILHOL J.C. (2012). Effect of water stress on growth, water consumption and yield of silage maize under flood irrigation of Tadla (Morocco), 2012. *J. Biotechnologie, Agronomie, Société et Environnement*, 2012, 16 (4), 468-477.

BOUAZZAMA B., BOUAZIZ A., XANTHOULIS D. et BAHRI A. (2015). Effet du déficit hydrique sur la croissance, le rendement et l'efficience d'utilisation de l'eau chez la luzerne (Medicago sativa L.) au Tadla. *Revue Marocaine des Sciences Agronomiques et Vétérinaires* 3, n° 2 (2015): 16-26.

BOUNDI Ah. (1989). Contribution à la gestion optimale de l'irrigation localisée d'un verger d'agrumes établis dans la région du Souss (Maroc). Mémoire de fin d'étude pour l'obtention du diplôme d'Ingénieur Horticole. Institut Agronomique et Vétérinaire Hassan II-Campus Agadir.

BOUNDI Ah. (1995). Besoins en eau, économie eau d'irrigation et diversification variétale de la culture de l'oignon (*Allium cepa* L.) dans la région du Tadla. Mémoire de 3ème cycle pour l'obtention du diplôme d'Ingénieur Agronome (option : Horticulture). Institut Agronomique et Vétérinaire Hassan II-Campus Agadir.

BOUTHIBA, A. et DEBAEKE P. (2019). Besoins en eau de différentes variétés de blé dur en conditions semi-arides. In *Symp. Intern. : « Agriculture durable en région méditer. (AGDUMED) », Rabat, Maroc*, 14-16, 2009.

BOYER J.F. (1998). Khronostat statistical time series analyses software. Montpellier : UMR 5569 Hydro-sciences, IRD-Maison des sciences de l'eau.

BRAUMER J. S. (1997). Nonparametric estimation of slope : Sen's method in environnemental pollution, Environnement sampling & monitoring primer.

B.R.G.M. (Bureau de Recherches Géologiques et Minières-Maroc), 1993. Étude du système multicouche de la plaine de Tadla. Description d'aquifères Maroc, Rapport R35-610 4S/Eau92, 129 p.

BUISHAND, T. A. (1982). Some methods for testing the homogeneity of rainfall records. Hydrol. (58): 11-27.

BUISHAND, T. A. (1984). Tests for detecting a shift in the mean of hydrological time series. J. Hydrol. (73) : 51-69.

Centre Canadien de Cartographie et d'Observation de la Terre (CCCOT) (2020). Ministère canadien des ressources naturelles. Adresse ministériel : 5700, 4e Avenue Ouest Québec (Québec-Canada) G1H 6R1, adresse du CCCOT: 580, rue Booth, Ottawa, ON K1A 0G1, Canada. Site web: http://www.rncan.gc.ca.

CHIBANE, B., BENTCHAKAL M., MEDJERAB A. et BENFARES B. (2015). Etude de la variabilité et de la structure de pluie annuelle dans une région semi-aride : cas du bassin versant de la Macta (Nord-Ouest algérien) ». *Larhyss Journal* PISSN 11123680/EISSN 26027828, no 23 (2015) : pp : 213-229.

CISSE, M., SAMBOU S., DIEME Y., DIATTA C. et BOP M. (2014). Analyse des écoulements dans le bassin du fleuve Sénégal de 1960 à 2008. *Rev. Sc. eau/Journal of Water Science* 27, n° 2 (2014): 167-187.

CISSÉ, A., EYMARD L., NDIONE J.A. et GAYE A.T. (2015). Analyse statistique des relations pluie-végétation au Ferlo (Sénégal) ». *XXVIIIème Col. Assoc. Internationale de Climatologie, Liège (France)* : 307-312, 2015.

CISSE, S. (2016). « Etude de la variabilité intra saisonnière des précipitations au Sahel : impacts sur la végétation (cas du Ferlo au Sénégal) », 2016.Thèse de doctorat, spécialité : Climatologie, Université Pierre et Marie Curie (Sorbonne-Universités), Université Cheikh Anta Diop de Dakar (Sénégal). 160 pages.

CONSTANTIN (OPREA) D.M., MARINICĂ A.F., MARINICĂ I., ZAHARIA L., GRIGORE E., ONȚEL I., IOANATOROIMAC G. and TIȘCOVSCHI A.A. (2019): The thermal and pluvial features of the agricultural year 2018 in the southwestern Romania and the effects on the crops. Communication au 32ème colloque Internationale de l'AIC sous thème : « Le changement climatique, la variabilité et les risques climatique » ; Thessaloniki-Grèce, 29 mai au 1 juin 2019 ; Département de Météorologie et de Climatologie-Collège de Géologie Université d'Aristote de Thessaloniki, Grèce, pp : 543-548

DAGNELIE P. (1970). Théorie et méthodes statistiques. Volume1 et 2. Editions Duculot. Gembloux (Belgique).

DAGNELIE P. (2013). Statistiques théorique et appliquée. T 1 et 2. Ed. Boeck Services. Bruxelles (Belgique).

DAKI, Y. ZAHOUR G., LACHGAR R. et EL HADI H., (2016). Caractérisation de la sécheresse climatique du bassin versant d'Oum Er Rbia (Maroc) par indice standardisé des précipitation (SPI) ». *Eur. Sc. J.* 12, n° 14 (2016).

DALLEL J., SAKKA M. et HENIA L (2011). Evolution de la biomasse végétale en Tunisie entre 1981 et 2010. Communication au 24ème colloque de l'Association Internationale de Climatologie, du 6 - 9 Septembre 2011 sous le thème : Climat montagnard et risques. Université de Ferrara (Rovereto-Italie) pp: 159-164
DAYON G. (2015). Evolution du cycle hydrologique continental en France au cours des prochaines décennies. Thèse de doctorat, Ecole doctorale et spécialité : SDU2E : "Océan, Atmosphère et Surfaces Continentales, Université de Toulouse 3 Paul Sabatier, Toulouse (France). 223 p.

DELANY, F. (1965). Mission hydrogéologique des régions à substratum schisteux (Côte d'Ivoire), Rapport de mission SODEMI Côte d'Ivoire.

DE PARCEVAUX S., PAYEN D., BROCHET P., SAMIE CH., HALLAIRE M. et MERIAUX S., 1990. Dictionnaire encyclopédique d'Agro-météorologie. Edition Météo-France, INRA, 323 p.

DIONNET M., KUPER M., GARIN P., HAMMANI A., ELIAMANI A., et SAAF M. (2006). Accompagner les acteurs dans le changement de leur système : Un jeu de rôles pour des projets collectifs d'irrigation au Tadla, Maroc. In *L'avenir de l'agriculture irriguée en Méditerranée : Nouveaux arrangements institutionnels pour une gestion de la demande en eau*, 2006. Actes du séminaire Wademed, Cahors, France, 67 novembre 2006. CIRAD, Montpellier, France.

DOORENBOS, J. and PRUITT, W.O. (1977). Crop water requirements. Ir.& dr. No. 24, (rev.) FAO, Rome.

DOORENBOS, J. and KASSAM, A. H. (1979): Yield response to water. FAO Ir& dr. No. 33, FAO, Rome (Italy).

DOUKPOLO B. 2014. Changements climatiques et productions agricoles dans l'Ouest de la République Centrafricaine. Thèse de Doctorat en Géographie et Géosciences de l'Environnement (Agro climatologie et Développement), Université de Abomey-Calavi. https://tel.archivesouvertes.fr/tel01053261.

DRIOUECH F. (2010). Distribution des précipitations hivernales sur le Maroc dans le cadre d'un changement climatique : descente d'échelle et incertitudes ». PhD Thesis, Institut National Polytechnique de Toulouse (INP Toulouse), spécialité : Sciences de l'Univers, de l'Environnement et de l'Espace, Université de Toulouse.

EDWARDS D.C. and Mc KEE T.B. (1997). Characteristics of 20th century drought in the United States at multiple time scales. Climatology Report Number 97-2, Colorado State University, Fort Collins, Colorado.

EL AJHAR L., El KHACHINE D., EL BAKOURI A., EL KHARRIM K. et BELGHYTI D. (2018). Evolution de la pluviométrie de 1960 à 2015 au Maroc. *Int. J. of Res. Sc. and Manag.* 5(10) : Octobre, 2018.pp : 47-56.

EL JIHAD M. D. et TAABNI M. (2017). L'eau au Maghreb : Quel « mix » hydrique face aux effets du changement climatique. Contribution au rapport : « Eau et climat en Afrique du Nord et au Moyen Orient », coordonné par Zeineddine NOUACEUR. Université de Rouen-Normandie et organisme Campus-France ; éditions-Transversal, Târgovişte (Roumanie), 2017, pages : 11-25.

EL HAFID D., ZERROUQI Z. et AKDIM B. (2017). Etude des séquences de sècheresse dans le bassin d'Isly (Maroc oriental). *Larhyss Journal*, 31 (2017), pp: 83-94.

EL HARI A., CHAIK M., LEKOUCH N., SEDKI A. and LAHROUNI A. « Water needs in citrus fruit in a dry region of Morocco ». *Journal of Agriculture and Environment for International Development* (JAEID) 104, no 3/4 (2010): 91-99.

EL KSABI R. et LAGHZALI M. (2011). Approches d'ingénierie et démarches de conception des projets collectifs de reconversion à l'irrigation localisée : Cas du périmètre du Tadla. Communication des ingénieurs de l'entreprise des études NOVEC, adresse : Immeuble NOVEC, Park Technopolis,11 100, Sala El Jadida à Salé (Maroc). Exposé au Salon international de l'eau et de l'irrigation de la région Tadla-Azilal (SIEITA), dates du 19 au 22 octobre 2011, Beni Mellal (Maroc).

EL KHOUDRI K. et DAMNATI B. (2015). Les changements climatiques dans le Tafilalet (Sud-Est marocain) : Analyse des tendances. Climate change on Tafilalet region. *Larhyss Journal*, 3 (2015), pp: 16-25.

FAO (Food and Agriculture Organization Nations) (1998). *World reference base for soil resources*. Vol. 3. Food & Agriculture Org., 1998.

FATY A., KANE A. et NDIAYE A.L. (2017). Influence de la manifestation climatique sur les régimes pluviométriques saisonniers dans le haut bassin versant du Sénégal. Revue des sciences de l'eau / Journal of Water Science 30, no 2 (2017) : 79-87.

FEKI M. (2015). Variation et variabilité de la dynamique de la végétation au Nord de la Tunisie. 28ème colloque de l'Association Internationale de Climatologie du 01 au 04/07/2015 sous le thème : modélisation & variabilités organisé par Laboratoire de Climatologie et Topo-Climatologie, Univ. Liège (Belgique). 2015. pp: 336-342

FAO (Food and Agriculture Organization of United Nations) (2002). World agriculture; towards 2015/2030, report abstract, Rome, Italy, ISBN 92-5-204761-1.

FAO et INRA du Maroc. (2015). MOSAIC : Système de modèles pour évaluation des impacts du changement climatique sur l'agriculture au territoire marocain, 26 pages. Website: www.changementclimatique.ma

FERCHICHI, A. (2017). Evaluation des services de l'eau agricole et renforcement des interactions des acteurs de la gestion de l'eau : Cas d'un système à multi-niveaux de gestion (périmètre public irrigué de Zaouiet Jedidi, Cap Bon-Tunisie. Thèse de doctorat en sciences agronomiques ; spécialité: Génie Rural Eau et Forêt ; Institut National Agronomique de Tunisie, Ecole Doctorale Sciences et Techniques de l'Agronomie et de l'Environnement (Tunisie). Page :81.

FOSSOU, R. M. N., LASM T., SORO N., SORO T., SORO G., DE LASME O.Z., BAKA D., ONETIE O.Z. et OROU R. (2015). Variabilité climatique et son impact sur les ressources en eaux souterraines : Cas des stations de Bocanda et de Dimbokro au Centre-Est de Côte d'Ivoire. *Larhyss J.* n°21 (2015): 97-120.
FRENKEN K. and GILLET V. (2012). Irrigation water requirement and water withdrawal by country, report of AQUASTAT, FAO, Rome, Italy.

GAALOUL N. et LAIGNEL B. (2012). Variabilité et effets de changement climatique sur les ressources en eau souterraines en Tunisie : Une synthèse globale des résultats et recommandations. Actes du colloque international de Rouen sous thème : « Regards croisés Nord-Sud », Univ. de Rouen (France), 25-26 Septembre, 2012 ; pp : 71-72.

GAMMAR A. M. et HLAOUI Z. (2012). Utilisation des indicateurs du déficit hydrique dans les synthèses bioclimatiques méditerranéennes. Communication au Colloque de l'Association Internationale de Climatologie, du 05 au 08 Septembre 2012, sous thème : Les climats régionaux : Observation et modélisation. Université Joseph Fourier (Grenoble 1), Grenoble (France), 2012, pp : 321-326

GAUSSEN H. ET BAGNOULS F. (1957) « Les climats biologiques et leur classification », *Annales de Géographie*, vol. 6, no 355, 1957 ; pp : 193-220.

GAUTIER F, LUBESNIEL H., SABATIER R., MASSON J. M., PATUREL J.E. et SERVAT E. Variabilité du régime pluviométrique d'Afrique Ouest non sahélienne 1950-1989 ». *Hydr. Sc. J.* 43, n° 6 (1998): 921-935.

GHENIM A. N. MEGNOUNIF A. SEDDINI A. et TERFOUS A. (2010). « Fluctuations hydro-pluviométriques du bassin-versant Oued Tafna-Beni Bahdel (Nord-Ouest algérien) ». *Rev. Sèch*. vol. 21, n°2, Avril-Juin 2010

GIEC (I.P.C.C.) (2001). Climate Change 2001 : The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change[Houghton, J.T., Y. Ding, D.G. Griggs, M. Noguer, P.J. vander Linden, X. Dai, K. Maskell, et C.A. Johnson (éds.)].Cambridge University Press, Cambridge, Royaume Uni et New York, NY, Etats-Unis, 881 pp.

GIEC (I.P.C.C.) (2007). Bilan 2007 des changements climatiques. Rapport de synthèse du groupe d'experts intergouvernemental sur l'évolution du climat (GIEC) n°4, 114 pages, 2007.

GIEC (I.P.C.C.) (2013). Les changements climatiques, les éléments scientifiques. Publication du groupe d'experts inter-gouvernemental sur l'évolution du climat (GIEC), 5ème rapport du GIEC, résumé pour décideurs. 2013. Website: https://www.ipcc.ch.

GIEC (I.P.C.C.). 2014 : 2ème partie du 5ème rapport d'évaluation intitulée : « Changements climatiques en 2014 : Incidences, adaptation et vulnérabilité ». Résumé pour décideurs. Website : https://www.ipcc.ch

GIEC (I.P.C.C.). (2019). Climate Change and Land. An IPCC Special Report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems. 2019. Intergovernmental Panel on Climate Change. Electronic copies of this report are available from the IPCC. Website: https://www.ipcc.ch.

GOOSSENS Ch. R. (1985). « Principal component analysis of Mediterranean rainfall ». *J. Clim*.5, no 4 : 379-388

GUERMAZI E., MILANO M. et REYNARD E (2018). Evolution des ressources en eau dans un bassin versant aride en Tunisie centrale. 31ème *Col. Assoc*. Intern. Climatologie, du 02 au 03/07/2018 sous le thème: Echelles spatiales et temporelles fines. Université Nice, Sophia Antipolis. Nice (France). 2018 ; pp: 169-174.

GUTTMAN N.B. (1994). On the sensitivity of sample L moments to sample size. Journal of Climatology 7: 1026-1029. Guttman NB. 1994. *Journal of Climatology* 7: 1026-1029.

HABIBA K. NASR Z., KHALDI A., WOO S.Y., NOURI M. et REJEB M. N. (2012) « Efficience d'utilisation de l'eau et de l'azote par les plants de chêne-liège ». IOBC/wprs Bulletin 76 (2012): 35-42.

HADJOUT N. et LARFI B. (2012). L'élaboration d'un bilan hydrogéologique de la nappe côtière de la plaine du Gharb, communication au colloque international de Rouen sous thème : « Regards croisés Nord-Sud », Université de Rouen, Septembre, 2012, page: 83.

HAFIDI A., SALGHI R., ABOUATALLAH A., AYOUB M., AYOUB F. and AFFI N. (2017). Calibration of GFS model by using real time weather data for water requirement forecasting of citrus crop at the Moroccan Souss basin, s. d. 2017 *Int. Journal of Engineering Research and Application* (Part 2) Nov. 2017, pp.0111

HAFIDI A. ABOUATALLAH A., AFFI N., CHAOUIK A. and SALGHI R. Yobeen forecast software reliability assessment using measured data of the Souss upstream region. *App. J. Env. Eng. Sc.* 5, n° 1 (s. d.): 5-1. 2019.

HAFIDI A., ABOUATALLAH A., AFFI N. and SALGHI R. Calibration of GFS model by using real time weather data for water requirement forecasting of vegetables in Chtouka region. *Mor.J.Chem.* 8, n° 1 (s. d.): 8-1. 2020.

HAIR J. F., BLACK W. C., BABIN B. J., ANDERSON R. E. and ATHAM R. L. (1998). *Multivariate data analysis*. Vol. 5. 3. Prentice hall Upper Saddle River, NJ, 1998.

HAJRI J. (2006). Bilans hydriques agro-climatiques au Cap Bon (Tunisie) : Le cas des agrumes. pp 220-225. Communication au 29ème Colloque international sous thème : « Les risques liés au temps et au climat ». Actes du colloque d'Épernay 6-9 septembre 2006, organisé par l'Université Denis Diderot Épernay (France).

HAKMAOUI, A.; ZAKI N., BAYE Y. et OUATMANE A. (2013). Techniques de production du piment rouge (Niora) au périmètre du Tadla. Bulletin mensuel d'information et de liaison : Transfert de technologie en agriculture, n°198, Mars 2013, IAV Hassan II (Rabat).

HAMMANI, A., KUPER, M., DEBBARH, A., BOUARFA, S., BADRAOUI, M., BELLOOUTI, A. (2004). Evolution de l'exploitation des eaux souterraines dans le périmètre irrigué du Tadla. *Revue HTE*, (536) : 10-43

HAMMANI, A., KUPER M., BEKKAR Y. et ZAZ H. (2006) « Exploitation des eaux souterraines dans le périmètre irrigué de Tadla (Maroc) : Etat des lieux et éléments de méthodologie pour contribuer à une réflexion sur une
gestion intégrée et durable des eaux souterraines et de surface ». *Deuxième atelier régional du projet Economies d'eau en Systèmes Irigués au Maghreb (SIRMA)*, 2006.

HAMMANI A. et KUPER M. (2007). Caractérisation des pompages dans le Tadla, Maroc. *Troisième atelier régional du projet Economies d'eau en Systèmes Irigués au Maghreb SIRMA (SIRMA)*, 10-p, CIRAR, 2007.

HANCHANE, M. (2013). Méthodologie de régionalisation spatio-temporelle pour une analyse des précipitations (1961-1992) : Application au Maroc atlantique, communication de la Faculté Polydisciplinaire de Taza de l'Université Sidi Mohamed Ben Abdellah. Actes du Deuxième Colloque international « Eau & Climat : Regards croisés Nord/Sud », Fès (Maroc), 27-28 novembre 2013

HARBOUZE, R., LE GRUSSE PH, BELABES K., RAKI M., BOUAZIZ A. et RUELLE P. (2008). Efficiences économiques comparées des systèmes de production dans différentes situations d'accès à la ressource en eau. Application dans le périmètre irrigué du Gharb, Maroc. *Econ. eau sys. Irrigués au Maghreb*, 16-p. Cirad, 2008.

HATFIELD J. L., ASRAR G. and KANEMASU T. E. (1984). Intercepted photo-synthetically active radiation estimated by spectral reflectance. Remote sensing of Environment 14, no 1-3 (1984): 65-75.

HEBAL A. et REMINI M. B. (2012). Study of spatiotemporal variability of rainfall in Northern Algeria. In: *Proceeding of the first symposium: « Water and climate in the Maghreb North / South Crossed views »*. University of Rouen, Rouen, France.

HIRSCH, R.M. and SLACK, J.R. (1984): A nonparametric trend test for seasonal data with serial dependence, *Water Resources Research* (20): 727-732.

HOLLING C. S. (1973). Resilience and stability of ecological systems. *An. Rev. of ecol. and syst.* 4, n° 1: 1-23.

HOUGHTON, J. 2009. Global warming: the complete briefing. Chapter 3: The greenhouse gases. pp: 34-68. Cambridge university press, 2009.

HSIAO, T. C., STEDUTO, P. and FERERES, E. (2007). A systematic and quantitative approach to improve water use efficiency in agriculture. *Irrigation Science* 25: 209-231 (2017).

HUBERT P, CARBONNEL JP et CHAOCHE A (1989). Segmentation des séries hydro-météorologiques : Application à des séries de précipitations et de débits en Afrique de l'ouest. *J. Hydrology* 110 : (3) 49-67.

ICRA (Centre International pour la Recherche Agricole orientée vers le développement), ORMVA du Tadla et Institut National de la Recherche Agronomique-C.R.R.A. du Tadla (2003). Reconvertir un système d'irrigation: le pari difficile du périmètre irrigué du Tadla (Maroc). Série de documents de travail N° 114-2003 (Maroc) : pp 5474. Edition d'ICRA, Agropolis International, Avenue Agropolis, 34394 Montpellier CX05 - France.

INIESTA F., TESTI L., ORGAZ F. and VILLALOBOS F. J. (2009). The effects of regulated and continuous deficit irrigation on the water use, growth and yield of olive trees. *Eur. J. Agron.* 30 n° 4 (2009): 258-265.

JARLAN L., KHABBA S., ERRAKI S., LE PAGE M., HANICH L., FAKIR Y. (2015). MERLIN O., MANGIAROTTI S., GASCOIN S. and EZZAHAR J. Remote sensing of water resources in semiarid Mediterranean areas: The joint international laboratory TREMA. *International Journal of Remote Sensing* 36, n° 1920 (2015): 4879-4917.

JOUILIL, I., BITAR K., SALAMA H., AMRAOUI A., MOKSSIT A. et TAHIRI M. (2013). Sécheresse météorologique au bassin hydraulique Oum Er Rbia durant les dernières décennies. *LARHYSS Journal*, ISSN 11123680/EISSN 26027828, no 12 (2013).

KAYA, Ü, F. Ö. G., GÖKHAN Ç., ERHAN A., ŞERAFETTİN A. and KÖSEOĞLU O. (2017). Effect of deficit irrigation regimes on yield and fruit quality of olive trees (cv. Memecik) on the Aegean Coast of Turkey. *Irrigation and drainage* 66, n° 5 (2017): 820-827.

KHARROU M. H., ERRAKI S., CHEHBOUNI A., DUCHEMIN B., SIMONNEAUX V., LEPAGE M., OUZINE L. et JARLAN L. (2011). Water use efficiency and yield of winter wheat under different irrigation regimes in a semi-arid region. *Agricultural Sciences in China* 2, n° 03 (2011): 273-282.

KHODJA H., LUBESNIEL H., SABATIER R., MASSON J.M., SERVAT E. et PATUREL J.E. (1998). Analyse spatio-temporelle de données pluviométriques en Afrique de l'Ouest. Recherche d'une rupture en moyenne. Une alternative intéressante : les tests de permutations. *Rev. stat. app.* 46, n° 1 (1998) : 95-110.

KOBRY A et ELIAMANI A. (2004). L'irrigation localisée dans les périmètres de grande hydraulique : Atouts et contraintes dans le périmètre du Tadla au Maroc. In : Actes du Séminaire Modernisation de l'Agriculture Irriguée Rabat, du 19 au 23 avril 2004, Projet Inco-Wademed, Maroc.

KENDALL, M. G. and STUART A. (1976). The Advanced Theory of Statistics. 3'd edition, Vol.3. Griffin and Company, London, 1976.

KOUDOU A., KOUAME K. A., NIAMKE K. H., KOUAME K. F., SALEY M.B. et ADJA M. G. (2015). Contribution de l'analyse en composantes principales à la régionalisation des pluies du bassin versant du N'Zi, Centre de la Côte D'Ivoire. *Rev. Ivoir. Sci. Technol* 26 (2015) : 156-172.

KRIMISA S., BOUCHAOUI L. et ABOURICHA S (2017). Apport des techniques statistiques à l'étude de la variabilité climatique : Exemple du bassin de Tadla-Azilal (Maroc) ». *Larhyss Journal*, no 30 (2017): 7-26.

LAURENT G. (1995). Variabilité des précipitations annuelles en Ardèche : Analyse factorielle en composantes principales. *La Météorologie*, 1995.

LEE A.F.S. and HEGHINIAN S.M. (1977). A shift of the mean level in a sequence of independent normal random variables A bayesiian approach. *Technometrics* J. (19): 50-36.
LENSSEN N. J. L., GAVIN A. S., HANSEN J. E., MENNE M.J., PERSIN A., RUEDY R. and ZYSS D. (2019). Improvements in the GISTEMP uncertainty model. *J. Geoph. Res.: Atmos.* 124, n° 12 (2019): 6307-6326.

LIONBOUI H. (2018). Modélisation agro-économique dynamique pour la gestion durable des ressources en eau dans les zones arides et semi-arides. Thèse de doctorat. FST-Beni Mellal. Université Moulay Slimane. Béni Mellal (Maroc).

LUBESNIEL H., MASSON J. M., PATUREL J. E. et SERVAT E. (1998). Variabilité climatique et statistiques : Etude par simulation de la puissance et de la robustesse de quelques tests utilisés pour vérifier l'homogénéité de chroniques. *Revue des sciences de l'eau*, (3) : 383-408.

LUBESNIEL H., MASSON J. M., SERVAT E., PATUREL J. E., KOUAMÉ B. et BOYER J. F. (1994). Caractérisation de fluctuations dans une série chronologique par application de tests statistiques. Etude bibliographique. ORSTOM, FRIEND AOC. Programme ICCARE, N° 3, 21p.

MANN H. B. (1945). Nonparametric Tests against Trend, J. Econometrical Soc. 1945, n°13, 245-259.

MARTIN-VERTEDOR A.I., RODRIGUEZ J.M.P., LOSADA M.H.P. and FERERES C.E.F. (2011). Interactive responses to water deficit and crop load in olive (*Olea europaea* L., cv. Morisca). II- Water use, fruit and oil yield. *Agricultural Water Management* 98(6): 950-958.

McKEE T.B., DOESKEN N.J. and KLEIST J. (1993). The relationship of drought frequency and duration times scales. American Meteorological Society. 8th *conf. App. Climatology*, 1722 Jan., Anaheim, CA, pp: 179-184.

McKEE T.B., DOESKEN N.J. and KLEIST J. (1995). Drought monitoring with multiple times scales. American Meteorological Society 9th conference on Applied Climatology, 1522 Janvier 1995, Dallas, TX, pp: 233-236.

MEDEJERAB, A. et HENIA L. (2011). Variations spatio-temporelles de la sécheresse climatique en Algérie Nord-Occidentale ». Courrier du savoir 11 (2011) : 71-79.

Ministère Marocain de l'Agriculture. (2020). Portail au site web du ministère marocaine de la pêche maritime, du développement rural et des eaux et forêts. Rubrique internet consultée : « Ressources Hydriques ». Website: http://www.agriculture.gov.ma/

MEDDE-France (2012). Hydrologie de surface. Projet Explore 2070. Rapport de synthèse A1. Ministère français de l'Ecologie, du Développement Durable et de l'Energie, 2012. 93p.

MJEJRA M., HENIA L. et DUBREUIL V. (2012). Apport de l'imagerie satellite à l'estimation de l'évapotranspiration réelle journalière dans le delta de la Mejerda en Tunisie ». Communication au 25ème Colloque internationale de l'Association Internationale de Climatologie du 05 au 08 Septembre 2012 sous thème : « Les climats régionaux : observation et modélisation ». Université Joseph Fourier (Grenoble 1), Grenoble (France), 2012, pp : 541-546.

MJEJRA M. et SAMAALI H. (2019). Essai de cartographie du risque de stress hydrique dans la plaine deltaïque de l'Oued Majerda (Tunisie). Communication au 32ème Colloque internationale de l'Association Internationale de Climatologie sous thème : « Le changement climatique, la variabilité et les risques climatique » ; Thessaloniki-Grèce. 29 mai au 1 juin 2019 ; Département de Météorologie et de Climatologie-Collège de Géologie Université d'Aristote de Thessaloniki ; Grèce ; pp : 229-233.

MONIOD F. (1973). Étude hydrologique de l'Oum Er Rbia. Cahiers. O.R.S.T.O.M. séries. Hwhl., 10(2) 197.

MOURADI M., BOUIZGAREN A., FARISSI M. and GHOULAM C. 2018. Assessment of deficit irrigation responses of moroccan alfalfa landraces grown under field conditions ». *Ir. & drainage* 67, n° 2 : 179-190.

MRAD D. (2012) Impact de la variabilité climatique sur les ressources en eau cas du bassin versant Seybousse. Université Mohammed Chérif Messadia, Actes du colloque International de Rouen ; sous le thème : « Regards croisés Nord-Sud », Université de Rouen, 25-26 Septembre, 2012 ; Rouen (France). Page: 122.

MURARESCU O., MURATOREANU G., NOUACEUR N. and FRINCULEASA M. (2017). Agro-Meteorological Drought in Boian Plain (Romania). Case Study: Years 2008-2013. *Annals* ” Valahia” University of Targoviste Agriculture 11, No 1 (2017): 62-65.

NAJIB S., FADILI A., MEHDI K., SALOUI A., ILMEN R. et SEBBAR A. (2018). Evaluation de l'effet des variations climatiques sur l'aquifère de la Chaouia côtière (Maroc). Actes du 31ème Colloque de l'Association Internationale de Climatologie, du 02 au 03/07/2018 sous le thème : Echelles spatiales et temporelles fines. Université Nice Sophia Antipolis. NICE. 2018. pp : 128-133.

NORRANT C. et DOUGUEDROIT A. (2004). Tendances des précipitations mensuelles et quotidiennes dans le Sud-Est méditerranéen français (1950-1951/1999-2000) ». *Annales* de l'AIC. Climatologie, no Volume 1 (1 janvier 2004) : 4564. Website : https://doi.org/10.4267/climatologie.987.

NICHOLSON, S. E. and PALAO I.M. (1993). A Re- evaluation of Rainfall Variability in the Sahel. Part I. Characteristics of Rainfall Fluctuations ». *International Journal of Climatology*, 01 Mai 1993.

OCDE (2012) Perspectives de l'environnement de l'OCDE à l'horizon 2020, résultats du Modèle IMAGE. Éditions de l'OCDE, 2012; pp: 12-13.

OMARI, S., BORD J. P., BOUAZIZ A., LE GRUSSE Ph, POUSSIN J. C. et RUELLE P. (2008). Le SIG, un outil de gestion de l'eau dans un périmètre irrigué: cas du secteur S7 au Gharb, Maroc. Actes du quatrième atelier régional du projet Sirma sous thème : Economies d'eau en systèmes irrigués au Maghreb, Mostaganem, Algérie, 2628 mai 2008. Cirad, Montpellier, France, colloques-cédérom. 11-p. Cirad, 2008.

OMM (2001). Bilan des changements climatiques 2001. Rapport de synthèse de l'Organisation Météorologique Mondiale, 204 pages.
OMM (2012). Guide d'utilisation de l'indice normalisé de précipitations. Rapport N°1090 de l'Organisation Météorologique Mondiale (OMM), page :1. Edition OMM-Genève (Suisse).

OMM (2020). Climat. Site web de l'Organisation Météorologique Mondiale : https://public.wmo.int/en/our-mandate/climate.

O.R.M.V.A du Tadla. (1996). Stratégie et plan d'action de l'amélioration de l'application de l'eau à la parcelle dans le périmètre de Tadla. Rapport de l'établissement ORMVA du Tadla; siège à Fkih Ben Salah (Maroc).

ORMVA du Tadla (2007). Bilan annuel de réalisations en équipement en irrigation localisée au périmètre du Tadla.

ORMVA du Tadla (2012). Analyses de faisabilité technicoéconomique et financière de la reconversion collective à l'irrigation localisée dans la zone des Béni Moussa ORMVA du Tadla; siège à Fkih Ben Salah (Maroc).

ORMVA du Tadla, 2015. Bilans annuels des essais culturaux en station des expérimentations hydroagricoles (SEHA) au périmètre du Tadla, pour la période de 2011-2014. Siège à Fkih Ben Salah (Maroc).

ORMVA du Tadla (2018). Monographie de zone d'action de l'ORMVA du Tadla; Fkih Ben Salah (Maroc).

ORMVA du Tadla (2019). Bilan des activités de la Station des Expérimentations Hydro-Agricoles d'Ouled Ganou (SHAOG) pour la période de 2003-2018. ORMVA du Tadla ; Fkih Ben Salah (Maroc).

ORMVA du Tadla (2020 a). Recueil de la base de données agro-météorologique du périmètre du Tadla pour période de 1958-2020. ORMVA du Tadla ; siège Fkih Ben Salah (Maroc).

ORMVA du Tadla (2020 b): Bilans des irrigations au périmètre du Tadla ; 1979-2020. ORMVA du Tadla (Maroc).

ORMVA du Tadla. (2020 c). Informations sur établissement ORMVA du Tadla (Fkih Ben Salah ; Maroc). Rubriques internet consultées : « Monographie », « Développement Agricole », « Gestion des Réseaux » et « Aménagement des Réseaux ». Site du web : https://ormvatadla.ma/

OUATTAB M., HAMMOUDY W., DAHMAN A., ILMEN R., SALOUI A., HSAINE M. et SEBBAR A. (2019). « Evolution des températures observées et projections futures en région Casablanca-Settat (Maroc) ». Communication au 32ème colloque internationale de l'AIC sous thème: « Le changement climatique, la variabilité et les risques climatique » ; Thessaloniki-Grèce, 29 mai au 1 juin 2019; Département de Météorologie et de Climatologie. Collège de Géologie Univ. d'Aristote de Thessaloniki; Grèce; pp: 33-38.

OUBICHA M. (2019). Cartographie des cultures à l'aide des images satellitaires dans le périmètre irrigué du Tadla. Rapport de projet pour examen d'aptitude professionnelle d'accès au grade Technicien 2ème grade, exercice 2019. ORMVA du Tadla. Fkih Ben Salah (Maroc).

OUSTOU B. (2015). Etude perspective de l'impact de la reconversion de l'irrigation gravitaire à l'irrigation localisé sur les ressources en eau souterraines au périmètre irrigué de Tadla. Projet de fin d'études, présenté pour l'obtention du diplôme d'Ingénieur en Génie Rural, option : Irrigation, Gestion de l'Eau et Environnement ; 2015. Institut Agronomique et Vétérinaire Hassan II-Rabat (Maroc); pp: 74-76.

OUZEMOU J.E., EL HARTI A., LHISSOU R., EL MOUJAHID A., BOUCH N., EL OUAZZANI R., BACHAOUI E. and EL GHMARI A. (2018). Crop type mapping from pan-sharpened Landsat 8 NDVI data: A case of a highly fragmented and intensive agricultural system. *Remote Sensing Applications: Society and Environment*, vol.11 (2018): 94-103.

OVIDIU M., MURĂTOREANU G., NOUACEUR Z. and FRÎNCULEASA M. (2017). Agro-Meteorological Drought in Boian Plain (Romania). Case Study: Years 2008–2013. *Annals "Valahia" University of Targoviste Agriculture* 11, n° 1 (2017): 62–65.

PATUREL J.E, SERVAT E., LUBES N.H. et DELATTRE M. O. (1997). Variabilité climatique et analyse de séries pluviométriques de longue durée en Afrique d'Ouest et Centrale non sahélienne. *Comptes-rendus d'Académie des SciencesSéries IIA Earth and Planetary Science* 325, n° 10 (01/11/1997) ; (7): 79-82.Website : https://doi.org/10.1016/S12518050(97)827565.

PETTITT A. N. (1979). A nonparametric approach to the change-point problem, *Appl. Stat.*, v 28, N°2, 126135.

PHOCAIDES A. Manuel des techniques d'irrigation sous pression. FAO, 2008. FAO, Rome (Italie).

PLANTON S (2013). Annex III: Glossary of IPCC. In: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Inter-governmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA.

PULINA M. A. (2012). Sécheresses et déficience hydrique des sols dans des zones agricoles de la Sardaigne durant la période 1951-2010. Colloque de l'AIC, 05 au 08 Septembre 2012 sous thème : Les climats régionaux : observation et modélisation. Université Joseph Fourier (Grenoble 1) Grenoble, 2012.

RAKOTOMALALA R. (2008). Tests de normalité-Techniques empiriques et tests statistiques. Université Lumière Lyon 2, Version 2 (2008) pp : 14-16.

RAMDI A. et EL KHALFI S. (2019). Diagnostic et évaluation de l'offre et la demande en eau d'irrigation dans le secteur pilote du périmètre irrigué du Tadla. Rapport du projet de Fin d'Etudes pour obtention du diplôme d'Ingénieur d'Etat en Génie Rural, 2019. IAV Hassan II-Rabat (Maroc) ; pp : 68-71.

REZIG M., SAHLI A. et LUC J. P. (2005). Impact de la distribution de l'eau sur la variabilité spatiale de la croissance et du rendement du Piment irriguée par système gravitaire. Colloque WATMED 2 en date du 14 au 17 Novembre 2005 à Marrakech (Maroc).

RODRIGUES G. C., PAREDES P., GONÇALVES J. M., ALVES I. and PEREIRA L.S. (2013). Comparing sprinkler and drip irrigation systems for full and deficit irrigated maize using multicriteria analysis and simulation modelling: Ranking for water saving vs. farm economic returns. *Agr.water manag.* (126): 85–96.
ROERINK, G. J., JACOBS C. and HAMMANI A. (2009). Assessment of groundwater extraction in the Tadla irrigated perimeter (Morocco) using the SSEBI remote sensing algorithm. 2009. Project Aquastress relating to the European Union. Web site: www.aquastress.net/index.php

RONCHAIL J., MORCEL C., COHEN M. et LABAT D. (2012). Précipitations et rendements des oliviers dans la province de Jaén (Andalousie), période de 1955-2009. Colloque de l'AIC05 au 08/09/2012 sous thème : Les climats régionaux : observation et modélisation. Université Joseph Fourier (Grenoble 1) Grenoble, 2012.

Royaume du Maroc (2019). Programme prioritaire pour l'approvisionnement en eau potable et l'irrigation. Rapport de la Commission Technique Interministérielle, représentante du Ministère de l'Intérieur, du Ministère de l'Agriculture, de la Pêche Maritime, du développement Rural, des Eaux et des Forêts, du Ministère de l'Economie et des Finances, du Ministère de l'Equipement, du Transport, de la Logistique et de l'Eau, du Ministère de l'Energie, des Mines et du Développement Durable et de l'Office National de l'Electricité et de l'Eau Potable ; pp : 23-35. Rabat (Maroc).

SALAMA, H. 2011. « Régionalisation pluviométrique du bassin Tensift (Maroc) ». *Larhyss Journal* PISSN 11123680/EISSN 26027828, n° 9 (2011).

SARR B. (2006). Les changements climatiques et la vulnérabilité des ressources en eau et de l'agriculture, Centre Régional AGRHYMET, DFR, Niamey. Communication scientifique, 22 pages.

SEBBAR A., BADRI W., FOUGRACH H., HSAINE M. et SALOUI A. (2011). Étude de la variabilité du régime pluviométrique au Maroc septentrional (1935-2004). Sciences et changements planétaires / Séch. 22, N° 3: 139-148.

SEBBAR, A., HSAINE M., FOUGRACH H. et BADRI W. (2012). Étude des variations climatiques de la région centre du Maroc ». 25ème Colloque de l'Association Internationale de Climatologie, Grenoble sous thème : « Climats régionaux : observation et modélisation », 2012, pp : 709-714.

SEBBAR, A., HSAIN M., FOUGRACH H. et BADRI W. (2013). Carte des précipitations annuelles au Maroc (1935/2006). 26ème Colloque de l'Association Internationale de Climatologie. Cotonou, Bénin sous thème : « Climat, agriculture, ressources en eau : D'hier à demain », 2013, pp : 37-43.

SEN, P.K. (1968). Estimates of the regression coefficient based on Kendall's Tau. *Journal of the American Statistical Association*. (63):1379-1389.

SERVAT E., PATUREL J.E., KOUAME B., TRAVAGLIO M., OUEDRAOGO M., BOYER J. F., LUBESNIEL H., FRITSCH J.M., MASSON J.M. et MARIEU B. (1998). « Identification, caractérisation et conséquences d'une variabilité hydrologique en Afrique de l'Ouest et centrale », *IAHS Publication*, vol. (252) : 323-337.

SHEPARD D. (1968). A two dimensional interpolation function for regularly spaced data. In, 23d National Conference of the Association for Computing Machinery, Princeton, NJ, ACM, pp: 517-524.

SHONGWE V. D., MAGONGO B.N., MASARIRAMBI M. T., ABSALOM M. and MANYATSI A. M. (2010). Effects of irrigation moisture regimes on yield and quality of paprika (*Capsicum annuum* L.). *Physics and Chemistry of the Earth, Parts A/B/C* 35, n° 1314 (2010): 717-722.

SIKAOUI L., BOUIZGAREN A., KARROU K., BOULAL H., BOULOUHA B., KARAMA K., NANGIA V. and OWEIS T. (2014). Effect of deficit irrigation on vegetative growth and fruit yield parameters of young olive trees (Olea europaea L.) in semi-arid area of Morocco. Paper presented at the 7th International Conference on Water Resources in the Mediterranean Basin, 10-12 October, Marrakech, Morocco.

SKIREDJ A., BEKRAOUI A. et EL ALAOUI M. (2001). Avertissement à l'irrigation pour la culture de Pomme de Terre dans le périmètre du Loukkos (Maroc). *Rev. Marocaine Sc. Agr. et Vété*. 21, n° 3 (2001): 169-176.

SPECTY, R. et ISBERIE C. (1996). L'avertissement à l'irrigation à l'épreuve de la pratique agricole. *La Météorologie*, 8ème série N° 16 Décembre 1996.

SPF (2020). Le Service Publique Fédéral de Santé publique, Sécurité de la Chaîne alimentaire et Environnement SPF SPSCAE, Gouv. belge. 2020. Site web: https://www.climat.be/frbe/servicefederal/

SOUDI B. (1994). Qualité des sols et des eaux : situation actuelle et méthodes de suivi de paramètres de qualité. Rapport n° 15 du Projet MRT, ORMVAT. Maroc.

STOUR L. et AGOUMI A. (2008). Sécheresse climatique au Maroc durant les dernières décennies. *Hydroécologie appliquée* 16 (2008): 215-232.

TAIBI S., MEDDI M., SOUAG D. et MAHE G. (2013). Évolution et régionalisation des précipitations au Nord de l'Algérie (1936-2009). Climate and land surface changes in hydrology, *IAHS Publ*. 359 (2013):191-197.

TANOUTI O. (2017). La gestion intégrée des ressources en eau à l'épreuve du bassin versant: cas du bassin du Tensift au Maroc. PhD Thesis, Université Nanterre de Paris, Paris (France) 10, 2017. Page: 247.

THIEME M.L., LEHNER B., ABELL R. and MATTHEWS J. (2010): Exposure of Africa's freshwater biodiversity to a changing climate. *Conservation Letters*, 3(5), 324331.

TOPAK R., SÜHERI S. and ACAR B. (2011). Effect of different drip irrigation regimes on sugar beet (Beta vulgaris L.) yield, quality and water use efficiency in Middle Anatolian, Turkey. *Irr. Science* 29, n° 1: 79-89.

TOUHAMI M. (2017). Régionalisation et variabilité pluviométrique dans le Nord Centre-Ouest algérien (Approche statistique). 2017. Mémoire pour obtention du diplôme de Master en Sciences de la Nature et de la Vie, Filière : Biologie, Spécialité : Eau et bioclimatologie Université Djilali Bounaama-Khemis Miliana, Faculté des Sciences de la Nature et de la Vie et des Sciences de la terre, Département de Biologie.
TRAMBLAY Y., GACHON P., STHILAIRE A. et CHAUMONT D. (2005). Variabilité et extrêmes de température et de précipitation entre 1941 et 2000 : le cas du golfe du fleuve Saint-Laurent. Rapport de recherche R805. INRS, Centre Eau, Terre et Environnement, 2005.

TRAORE V., NDIAYE M. L., MBOW C., GIOVANI M., SARR J., BEYE A. and DIAW A. (2017). Khronostat Model as Statistical Analysis Tools in Low Casamance River Basin, Senegal. 2017 (01/01/2017): 10-22.

TUCKER, C.J. (1979). Red and Photographic Infrared Linear Combinations for Monitoring Vegetation', Remote Sensing of Environment, 8(2): 127-150.

VAN HALSEMA, G. E. and VINCENT, L. (2012). Efficiency and productivity terms for water management: A matter of contextual relativism versus general absolutism. Agricultural Water Management 108: 9 - 15 (2012).

VERMEIREN, L. et JOBLING G. A. (1983). L'irrigation localisée. *Bulletin FAO d'irrigation et de drainage*, n° 36 (1983). FAO, Rome (Italie).

WIGLEY T. M. L., LOUGH J. M. and JONES P. D. (1984). Spatial patterns of precipitation in England and Wales and a revised, homogeneous England and Wales precipitation series. *J. of Clim*. 4, n° 1 (1984): 1-25.

WITTLING, C. S. et MOLLE B. (2017). Evaluation des économies d'eau à la parcelle réalisables par la modernisation des systèmes d'irrigation. PhD Thesis, IRSTEA, Université de Montpellier en collaboration avec l'institution UMR GEAU, adresse : 361, rue JF Breton BP 5095 – 34196, Montpellier Cedex 5 2017.

WU, H., HAYES, M.J., WILHITE, D.A and SVOBODA, M.D. (2005). The effect of the length of record on the standardized precipitation index calculation. *International Journal of Climatology*, (25): 505-520.

YAMAUCHI H. (1984). Impact on groundwater resources of conversion from furrow to drip irrigation. *JAWRA J. American Water Resources Association* 20, n° 4 (1984): 557-563.

YAO F. Z., REYNARD E., FALLOT J.M., OUATTARA I., KOUAKOU E. et SAVANE I. (2018). Analyse de la variabilité climatique dans le bassin versant agricole du Bandama-Blanc, Nord de la Côte d'Ivoire. Actes du 31ème Colloque de l'Association Internationale de Climatologie, du 02 au 03/07/2018 sous le thème : « Echelles spatiales et temporelles fines ». Université Nice-Sophia Antipolis ; Nice ; 2018 ; pp : 399-407.

YAPIA.F., OGAY.M.S., N'GNESSANBIV. H., KPAN O.J.G. et BIEMI. J. (2018). Caractérisation des séquences de sécheresse dans la région de la Marahoué partir des indices SPI et SPEI. Actes du 31ème Colloque de l'Association Internationale de Climatologie, du 02 au 03/07/2018, sous le thème : « Echelles spatiales et temporelles fines », Université Nice-Sophia Antipolis, NICE. 2018, pp : 158-162.

YERGEAU E. et POIRIER M. (2013). SPSS à l'Université de Sherbrooke : Régression logistique. Site francophone d'information sur SPSS. Accès au site web : http://spss.espaceweb.usherbrooke.ca/pages/stat

ZAIRI A., BEN NOUNA B., RUELLE P., NASR Z., OUESLATI T. et AJMI T. (2000). Irrigation et effet du déficit hydrique sur les productions: cas du Blé et de la Pomme de Terre ». *Actes du séminaire : « Economie de l'eau en irrigation », Hammamet-Tunisie : 1416 Novembre 2000*, 2000, 93-113.

ZITOUNI (2016). Apport de la télédétection dans l'évaluation de la variation des surfaces d'eau, du sol et du couvert végétal dans la Plaine du Tadla entre 1984 et 2015. Mémoire de stage de fin d'études pour l'obtention du diplôme du Master. Faculté des Sc.et Tech.de Marrakech (FSTM) Université Cadi Ayad (Maroc).

ZAMRANE Z., LAFTOUHI N.E., JARLAN L., LAIGNEL B. et MASSEI N. (2012). Caractérisation des relations précipitations/NDVI sur le Tensift : Variabilité inter et intra-annuelle et lien avec la dynamique atmosphérique et océanique à grande échelle. Communication au colloque international de Rouen sous thème : « Regards croisés Nord-Sud », 25-26 Septembre ; 2012. Université de Rouen (France) ; p 146.

ZHENG J., HUANG G., WANG J., HUANG Q., PEREIRA L.S., XU X. and LIU H. (2013). Effects of water deficits on growth, yield and water productivity of drip irrigated onion (*Allium cepa* L.) in an arid region of Northwest China. *Irrigation Science* 31, n° 5 (2013): 995-1008.
ANNEXES
# Annexe 1

## METHODES D'ETUDE DE VARIABILITES PLUVIOMETRIQUES ET DU PHENOMENE DE SECHERESSE

### A-Méthodes d'étude des variabilités des précipitations et des sécheresses

1-Méthode de l'intervalle de confiance

2-Méthode par indice pluviométrique de l'écart à la normale

3-Méthode par indice des précipitations standardisé (indice de Nicholson)

4-Méthode du filtre passe-bas de Hanning d'ordre 2

5-Méthodes d'étude de tendances sur séries temporelles agro-météorologiques

5-1-Test de Mann-Kendall

5-2-Méthode de Sen's

5-3-Tests d'homogénéité sur séries temporelles (Tests: Pettitt; Lee et Heghinian; Buishand; segmentation de Hubert)

6-Méthodes d'étude de tendances par les analyses multivariées des variances

### B-Méthodologie d'étude de variabilité thermique et évapo-transpiration de référence

1-Paramètres étudiés thermique et de la demande climatique en évaporation

2-Tests de variabilité thermique et de la demande climatique en évaporation

3-Etude de relation entre la température et l'évapo-transpiration de référence

4-Evaluation du bilan hydrique au sol dans un espace agro-écologique

### C-Méthode d'étude de l'évolution de l'aridité en zone semi-aride (Diagrammes de Gaussen ; Indice de De Maritonne ; Quotient d'Emberger)

### D-Techniques de la cartographie spatio-temporelle des tendances agro-climatiques

### E-Procédés d'étude d'impact des changements climatiques sur les besoins en eau d'irrigation des cultures en fonction du système d'irrigation

1-Expérimentations agricoles

2-Enquêtes agricoles chez des agriculteurs

3-Evaluation des besoins nets et bruts en eau irrigation par culture et par système d'irrigation

3-1-Détermination des efficiences des systèmes usuels d'irrigation

3-2-Evaluation de performance agronomique des cultures en fonction de l'eau

3-3-Modalités des projections à long-terme des besoins en eau des cultures

### F-Méthodes adoptées pour les traitements statistiques et numériques des données