Influence de La Lithologie Et Des Condit
Influence de La Lithologie Et Des Condit
Influence de La Lithologie Et Des Condit
وزارة ا ﺘ ﻴ ا ﺎ ﻲ و ا ﺤﺚ ا
BADJI MOKHTAR-ANNABA UNIVERSITY ﻨﺎ ﺔ-ﺟﺎ ﺔ ﺎﺟﻲ ﺨﺘﺎر
UNIVERSITE BADJI MOKHTAR-ANNABA
Thèse
Présentée en vue de l’obtention du diplôme de Doctorat en sciences
Option
Hydrogéologie
Par
GOUAIDIA Layachi
(Ingénieur d’Etat en Hydrogéologie. Magister en Chimie et environnement)
DEVANT LE JURY
PRESIDENT : Pr DJABRI L. Université de Annaba
RAPPORTEUR : Pr BOUDOUKHA A. Université de Batna
Meskiana plain, whose current hydrous resources are relatively limited consist of underground
water of the Mio-Plio-Quaternary aquifer. In recent years a high level of water demand was felt due to
remarquable increase water for agriculture use. On one hand, this has induced a shortage of water and
a deterioration of its quality on the other.
This research project comes to shed more light on the way the aquifer system works, i.e its
boundary conditions, influence of the lithological parameters and the hydrodynamic and physico-
chemical parameters. This has been carried out through several field excursions between Mars 2002
and November 2005.
Data analysis and graphical representation were done using statistical methods such as
principal component analysis by way of modern statistical softwares (WATEC, Diagrams and Surfur
version 8).
Data interpretation and analysis has revealed the following:
-Meskiana alluvial aquifer is mainly composed of sandy, sandstone and conglomeratic materials with
some passages of relatively thin clay layers.
It is limited to the north east by diapiric formations. The most important part of it’s recharge comes
from the carbonate formations of the borders.
-The chemical quality of water is poor resulting from the strong salinity influenced by the lithology of
the aquifer and the climatic factors (precipitations and the temperatures). Pollutants originating from
agricultural activities and urban wastes are also present.
The results obtained from this study agree perfectly with the questions raised in the proposal.
ﺨﺺ:
ﺎدر ا ﺎ ا ﺎ ﺔ و ﺔ ,إ أن ﺎع ا ﻮﻇﺎ ,ﺎ ﺔ ﻜ ﺎ ﺔ ل ا ﻮات اﻷ ﺮة ﻮا ا ﺎد ﺎ ﺮف ﻬ
ﺔ ﺎه ﻰ هﺬ ا ﺎدة ا ﻮ ﺔ.آ ﺎ أن هﺬ ا ﺎ ﺎت ا ﺰا ﺪ ا ﺔ ا ﺎ ا ﻮ ﺔ ﺮ ﺪ ﺪودة و ذ ﻚ را ا
ﺪ ﻮ ﺔ هﺬ ا ﺎ . ﺪر آ ﺮ
ﻜ ﺎ ﺔ ﺬ ﺪودهﺎ إ ﻰ ﺁ ﺎت ا ﺎل ﺎ ﺪة ا ﺎ ا ﻮ ﺔ ﻬ ﻮ أ هﺬا آﺎن ﺎر ﻮ ﻮع ا إن ا
ا درا ﺔ ا ﺎ ا ا ﻮا ا ﺎ ﺔ و ﻮ ﺔ ا ﺮ ﺔ.آ ﺎ ﺎ ﻬﺎ ا ﻬ ﺪرود ﺎ ﻜ ﺔ و ا ﺰ ﻮآ ﺎ ﺔ ا ﺎ ﺔ
ﺎدر هﺬا ا ﻮث. در ﺔ ﻮث ﻚ ا ﺎ و ا ﻜ
ﺎ ﺔ و ﺮا ﺔ ﺮات أ ﺪﺪ ﺮﺎ ا ﺎ ﺔ ا ﺎ ﺔ ﺎ ا ﻮ ﺔ وآﺬا در ﺔ ﻮ ﻬﺎ ﺪ ﺎ ﺔو
ﻬﺮ ﺎرس 2002إ ﻰ ﺎ ﺔ ﻬﺮ ﻮ ﺮ .2005 ا ا ﻬ ﺪرود ﺎ ﻜ ﺔ و ا ﺰ ﻮآ ﺎ ﺔ .ﻚ ا ﺎ ﺔ ا ﺮت ﺎ ا
اﻹ ﺎ و ا ﺮ ا ا ﻮﺎ ﺔ ﺪ ﺔ ﺎل و ﺎ و ﺮا ا ﺎت ا إن ﺎ ﺔ و
ﺎت ﺎ : ﻚا أﻇﻬﺮت ﺎ ﺮا و اﻷ ﻜﺎل
ﺬ ﺔ ا ﺎ ا ﻮ ﺔ. آ ﺮ ﻜ ﺎ ﺔ ﺎه ﻮر ا ﻜ ـــــﺔ ا ﺔ ﺎ ﺎ ﻰ ﺎ ﺔ ﻬ -إن و ﻮد ا
ﺎ ﺔ ,ﺪودهﺎ ا ﺎ ﺔ ا ﺮ ﺔ ا ﺎت ا ﻮي ا اﺮ ا ﻜﺎﺔ -ﻜﻮن ﺎ ﺪة ا ﺎ ا ﻮ ﺔ ﻬ
ﺎ ﺆﺮ ﻰ ﻮ ﺔا ﺎ. ﻮر ا ﺔ ﺎ ا ﺎ ﻮا
ﻮر وآـــﺬا ﻮا ا ﺎ ــــﺔ ) ا ﺮ و ﻮ ﺔا ﻬﺔ -إن ﻮ ـــﺔ ا ﺎ رد ﺔ ﺪا وذ ﻚ ﻜ ﺮة ﻮ ﻬﺎ ا ﺎ ــﺔ
ﻮق ﺮاآ ﺰهﺎ / 100ل ﺔ ﺎ ﺔ ا ﺮات ا ا ﺎ ﺎت ا ا ﻮ ﺎت ا ا ﺎ ( .آ ﺎ ﺬآﺮ ﻮ ﻮد
ﺔ. ﺎت ا ﺎ ا اﻷ ﺎن و آﺬا
. ﺎءت ﻮا ﺔ إ ﻜﺎ ﺔ ﻮ ﻮع ا آ ﺎ ﺬآﺮ أن هﺬ ا ﺎ
REMERCIEMENTS
L’élaboration de cette thèse de doctorat est le fruit de six ans de recherches, alternés
avec mes activités pédagogiques d’enseignement au Centre Universitaire de Tébessa. Son
achèvement n'aurait pu voir le jour sans la collaboration de nombreuses personnes qu'il m'est
agréable de remercier.
Je tiens dans un premier temps à remercier le Pr BOUDOUKHA Abderahmane, qui a
accepté de diriger ce travail et a veillé au bon déroulement de ces années, en m’apportant des
critiques constructives et des conseils pertinents, et ce suite à l’intérêt qu’il a accordé à mes
travaux de recherches.
Toute ma reconnaissance va également à DJABRI Larbi, Professeur au département de
Géologie de l’Université de Annaba, qui m’a prodigué encouragements et conseils, tout en me
permettant de profiter de ses connaissances. Il a toujours manifesté intérêt et appui à mes
travaux.
Je remercie le Dr Salim Hamid, enseignant au département de géologie de l’Université
de Damas (Syrie), qui s’est joint au Pr BOUDOUKHA pour encadrer une partie de ce
travail. Je lui adresse de sincères remerciements pour son implication dans mes recherches et
son aide très précieuse.
Mes vifs remerciements sont adressés au Dr François M., directeur du laboratoire
central de chimie à l’Université de Damas pour son bon accueil et les multiples analyses qu’il
m’a faites.
Mes remerciements s’adressent au Pr BEN MAMMOU A., Directeur du Laboratoire
des ressources minérales et de l’environnement au Département de Géologie de la Faculté des
Sciences de Tunis, qui m’a bien accueilli dans son laboratoire et m’a permis de contacter des
professeurs de géologie, tels que Mrs BOUHELEL S, TRIKI M. et bien d’autres, avec qui j’ai
eu des échanges scientifiques fructueux. Qu’ils soient assurés de ma sincère reconnaissance.
J’ai l’agréable tâche de témoigner ma grande reconnaissance aux différents services
techniques auprès desquels j’ai bénéficié d’un aimable appui, notamment ceux qui m’ont
facilité l’accès aux données de mesures et à la documentation, en particulier :
- l’antenne de l’A.N.R.H. de Tébessa (particulièrement Mr AGGOUNE Derradji) ;
- la Direction de l’Hydraulique de la Wilaya d’Oum El Bouaghi ;
- la Direction des Services Agricoles de la Wilaya d’Oum El Bouaghi ;
- l’Office National de la Météorologie (station de Tébessa) ;
- la Direction de CERAD (filiale de l’entreprise Ferphos) ;
Mes sincères remerciements s’adressent aux membres du jury qui ont accepté de lire et
juger ce travail. Le Pr DJABRI Larbi, qui a bien voulu m’honorer en présidant ce jury, malgré
ses énormes préoccupations, de même que le Pr HANI Azzeddine, les Dr HEMILA Mohamed
Laid, CHAFFAAI Hichem et CHABOUR Nabil, qui ont bien voulu juger ce travail et enrichir
le débat aux cotés du professeur BOUDOUKHA, mon directeur de thèse.
Un grand merci à Mme DJERRAB pour la correction orthographique, qui a permis de
faciliter la lecture et l’exploitation de cette thèse.
Je témoigne mon amitié à tous mes collègues de l’Institut des Sciences de la Terre de
Tébessa pour les bons moments passés ensemble et les échanges scientifiques.
Avant de terminer, je dois remercier ma famille pour son encouragement, son soutien et
sa patience.
Liste des Tableaux
INTRODUCTION GENERALE 01
Cadre physique du bassin
1 Situation géographique……………………..………………………………. 03
2 Réseau hydrographique ……………………………………………………. 03
3 Paramètres morphométriques……………………………………….……… 03
4 Activités économiques ………………………………………..……………. 03
Chapitre I : GEOLOGIE
I.1 Introduction ………………………………………………………...……… 05
I.2 Description lithostratigraphiques des 05
formations ……………………..……….
I.2.1 TRIAS ……………………………………………………………..………. 05
I.2.2 CRETACE …………………………………………………………………… 08
I.2.3 PALEOGENE…………………………………………………………………. 09
I.2.4 NEOGENE……………………………………………………………………. 10
I.2.5 QUATERNAIRE ……………………………………………………………… 10
I.3 Hydrogéologie et 10
stratigraphie………………………………………...………
I.3.1 Formations perméables ………………………………………………...……. 10
I.3.2 Formations imperméables …………………………………………………… 10
I.4 Géomorphologie……………………………………………………...……… 11
I.5 Aperçu tectonique…………………………………………………………….. 11
I.6 Conclusion…………………………………………………………………… 13
Chapitre II : GEOPHYSIQUE
II.1 Introduction ………………………………………………………………….. 14
II.2 Prospection électrique………………………………………………...…….. 14
II.2.1 Interprétations des résultats de la prospection électrique…………………….. 15
II.2.2 Apport de la prospection électrique……………………………………...….. 15
II.3 Apport des forages 15
mécaniques…………………………………………….…….
II.4 Conclusion……………………………………………………..…………….. 19
Chapitre III : HYDROCLIMATOLOGIE
III.1 Introduction…………………………………………………………………. 20
III.2 Acquisition et correction des données 20
……………………...……………………
III.2.1 Méthode des rapports …………………………………………….………….. 20
III.3 Caractérisation du climat………………………………………….…………. 22
III.3.1 Diagramme ombrothermique …………………………………...…………… 22
III.3.2 Indice d’Emmanuel De Martonne…………………………………………… 23
III.3.3 Indice d’aridité mensuel…………………………………..…………………. 24
III.3.4 Climagramme de Louis Emberger …………………….…………………….. 25
III.3.5 Indice de Paul Moral…………………………………………………………. 26
III.3.6 Comparaison des indices climatiques ……………………..………………… 26
III.4 La température………………………………………….………………. 27
III.4.1 Variations des températures moyennes annuelles ………………..…………. 27
III.4.2 Variations de la température moyenne mensuelle………………..………….. 27
III.5 La pluviométrie………………………………………………….………….. 28
III.5.1 Variations interannuelles des précipitations…………………………………. 28
III.5.2 Variations moyennes mensuelles des précipitations…………….……….. … 28
III.6 Estimation de l’évapotranspiration…………………………..…………… 29
III.6.1 Introduction ……………………………………………………..…………… 29
III.6.2 Estimation de l’ETP……………………………………………..…………… 29
III.6.2.1 Formule de L.Turc ……………………………………………...…………… 29
III.6.2.2 Formule de Thornthwaite ……………………………………….…………… 30
III.6.2.3 Discussion des résultats………………………………………….………..…. 31
III.6.3 Estimation de l’ETR……………………………………………..…………… 31
III.6.3.1 Formule de L. Turc……………………………………………..…………… 31
III.6.3.2 Formule de M. Coutagne…………………………………………………… 32
III.7 Le bilan hydrologique………………………………………...…………….. 32
III.7.1 Calcul du bilan hydrique par la méthode de Thornthwaite ……………..…… 32
III.7.2 Interprétation du bilan hydrologique………………………….……….……. 33
III.7.3 Répartition de l’excédent …………………. ………………………...……… 34
III.7.4 Commentaire sur le bilan hydrologique …………………………….……… 35
III.8 Les eaux de surface…………………………………………………………. 35
III.8.1 Estimation des eaux de surface …………………………………...…………. 35
III.8.2 Ecoulement de surface et régime hydrologique …………………….……….. 35
III.9 Conclusion ………………………………………………………….………. 37
Chapitre IV : HYDROGEOLOGIE
IV.1 Introduction……………………………………………………….…………. 38
IV.2 Géométrie de l’aquifère………………..………………………….………….. 38
IV.3 Piézométrie…………………………………………………………..………. 38
IV.3.1 Inventaire des points d’eau ……………………………………...…………… 38
IV.3.2 Evolution de la piézométrie ……………………………………….………… 38
IV.3.3 Carte piézométrique mars 2002 ……………………………………………… 39
IV.3.4 Carte piézométrique juillet 2002………………………………...…………… 40
IV.3.5 Carte piézométrique mars 2003…………………………………….………… 40
IV.3.6 Carte piézométrique juillet 2003………………………………...…………… 40
IV.3.7 Carte piézométrique mars 2004…………………………………….………… 42
IV.3.8 Carte piézométrique juillet 2004………………………………...…………… 42
IV.4 Les conditions aux limites …………………………………………………… 44
IV.5 Battement de la nappe ………………………………………………..……… 45
IV.6 Etude des paramètres hydrodynamiques…….............................................. 45
IV.6.1 Détermination des paramètres hydrodynamiques de la nappe ………………. 45
IV.6.2 Répartition spatiale de la perméabilité ……………….……………………… 47
IV.6.3 Répartition spatiale de la transmissivité ……..……………………………… 47
IV.7 Conclusion………………………………………..………………………….. 49
Chapitre V : HYDROCHIMIE
V.1 Introduction………………………………………………….………………. 50
V.2 Echantillonnage et analyse hydrochimique…………...……………………… 50
V.2.1 Moyens et méthodes d’analyse………………………….……………………. 51
V.2.1.1 Mesures in situ………………………………………………….…………….. 51
V.2.1.2 Analyses aux laboratoires…………………………………………………… 51
V.3 Description sommaire des paramètres physico-chimiques……………...... 52
V.3.1 Paramètres physiques………………………………………………………… 52
V.3.1.1 Potentiel d’Hydrogène (pH) ……………………………….………………… 52
V.3.1.2 La conductivité électrique (CE) ……………………………………………… 52
V.3.2 Les éléments majeurs……………………………………………..…………. 53
V.3.2.1 Le calcium……………………………………………………………………. 53
V.3.2.2 Le magnésium………………………………………………...………………. 53
V.3.2.3 Le sodium……………………………………………………….……………. 53
V.3.2.4 Potassium……………………………….…………………….………………. 53
V.3.2.5 Les chlorures……………………………….…………………………………. 54
V.3.2.6 Les sulfates…………………………………………………………………… 54
V.3.2.7 Les bicarbonates……………..………………………………………………. 54
V.4 Détermination des principaux faciès chimiques……………….………….. 54
V.4.1 Représentation graphique des faciès hydrochimiques ……………..………… 55
V.4.1.1 Diagramme de PIPER………………………………………………...………. 55
V.4.1.1.1 Période de mars 2002…………………………………………………………. 56
V.4.1.1.2 Période de juillet 2002………………………………………………..………. 56
V.4.1.1.3 Période de mars 2003…………………………………………………………. 56
V.4.1.1.4 Période de juillet 2003………………………………………………..………. 56
V.4.1.1.5 Période de mars 2004…………………………………………………………. 60
V.4.1.1.6 Période de juillet 2004…………………………………………..……………. 60
V.4.1.2 Diagramme de SCHOELLER-BERKALOFF………………...……………… 64
V.4.1.2.1 Période de mars 2002…………………………………………..……………. 64
V.4.1.2.2 Période de juillet 2002……………………………………….………………. 64
V.4.1.2.3 Période de mars 2003…………………………………………………………. 64
V.4.1.2.4 Période de juillet 2003………………………………………..………………. 68
V.4.1.2.5 Période de mars 2004…………………………………………………………. 69
V.4.1.2.6 Période de juillet 2004……………………………………………..…………. 68
V.4.1.3 Répartition des faciès chimiques…………………………………………….. 68
V.5 Cartographie 72
hydrochimique……………………………….…………………
V.5.1 Cartes des Bicarbonates………………………………………………………. 72
V.5.2 Cartes des chlorures…………………………………………….……………. 73
V.5.3 Cartes des sulfates……………………………...…………..……….…..……. 73
V.5.4 Cartes du calcium………………………………...………………….………. 77
V.5.5 Cartes du sodium………………………………………...……………...……. 77
V.6 Etude de l’origine du chimisme des eaux souterraines 78
………..……………….
V.6.1 Faciès caractéristiques ……………………………………….………………. 80
V.6.2 Origine des éléments chimiques …………………………...………………… 81
V.6.2.1 Les éléments Ca2+ – HCO3- – SO42-……………………………..…………… 81
V.6.2.2 Les éléments Na+ - Cl- ………………………………………………………. 82
V.6.2.3 Les éléments HCO3-, Cl- et SO42- ………………………..…………………… 83
V.6.2.3 Les éléments Cl--SO42- ………………………………………………………. 84
V.6.3 Etude de l’acquisition de la minéralisation ……………...…………………… 84
V.7.1 Introduction …………………………………………………..………………. 87
V.7.2 Principe……………………………………………………….………………. 87
V.7.3 ACP globale………………………………………………...…………………. 87
V.7.3.1 Matrice de corrélation ……………………………………...………………… 87
V.7.3.2 Analyse de l’espace des variables…………………………………………….. 88
V.7.3.3 Analyse de l’espace des individus……………………………………………. 88
V.8 Aptitude des eaux à 89
l’irrigation…………………………………..……………
V.8.1 Introduction…………………………………………..………….……………. 89
V.8.2 Caractérisation des eaux pour l’irrigation ……………………...…………….. 90
V.8.2.1 Classification des eaux par la méthode de Richards …………………………. 90
V.8.2.1.1 Cartographie d’aptitude des eaux à l’irrigation selon la classification de
Richards………..……………………………………………………..……….. 93
V.8.2.2 Classification des eaux par la méthode de Wilcox ……………………...……. 94
V.9 Conclusion ……………………………………………………………….…… 96
Chapitre VI : POLLUTION DES EAUX SOUTERRAINES
VI.1 Introduction………………………………….…….…………………………. 98
VI.2 Les principales sources de 98
pollution………………………………………………
VI.3 Impact des activités agricoles 100
……………………….…………………………
VI.3.1 Production des nutriments…………………………...……………………….. 100
VI.3.1.1 L’élevage ………………………………………………………………………. 100
VI.3.1.2 Les engrais azotés ………………………………………..……………………. 101
VI.3.1.3 Le fumier ………………………………………………………………………. 101
VI.3.2.1 Répartition spatiale des teneurs en nitrates…………...……………………… 101
VI.3.2.2 Evolution des nitrates dans le temps …………………….…………….……… 102
VI.3.3 Origine agricole des nitrates…………………………………………………… 103
VI.3.4 Evolution des nitrates en fonction des paramètres physico-chimiques .. 104
VI.3.4.1 Relation nitrate-nitrite …………………………………………………………. 104
VI.3.4.2 Relation nitrate-ammonium …………………………………………………… 104
VI.3.4.3 Relation nitrate-Oxygène………………………………………..…………….. 104
VI.3.4.4 Relation nitrate-profondeur …………………………………………………… 105
VI.3.5 Les phosphates ………………………………………………………………… 105
VI.4 L’influence des rejets urbains ……………………………………………….. 108
VI.4.1 Apport des nitrates………………………..……………………………………. 108
VI.4.2 Le bore………………………………………………………..……………….. 111
VI.5 Les métaux 111
lourds……………………………..……………..……….…………
VI.5.1 Le fer …………………………………………………….…………….………. 111
VI.5.2 Le strontium………………………………………………..………………….. 113
VI.5.2.1 Répartition du strontium ………………………………….…………………… 114
VI.5.2.2 Le rapport caractéristique Sr2+/Ca2+ …………………………..……..………… 115
VI.6 Etude de la vulnérabilité de la nappe de Meskiana ………….…………….. 116
VI.6.1 Introduction ……………………………………………………………………. 116
VI.6.2 Vulnérabilité à la pollution chimique des eaux souterraines ……….……..…... 116
VI.6.3 Présentation de la méthode DRASTIC………………………………….……. 116
VI.6.4 Validité de la méthode …………………………………………………….….. 118
VI.6.5 Etablissement de la carte de vulnérabilité de la nappe de Meskiana ……….… 118
VI.6.5.1 Profondeur de la nappe ……………………………………………………….. 119
VI.6.5.2 La recharge efficace ……………………………………………………..……. 120
VI.6.5.3 Nature de la zone saturée ……………………………………………….…….. 120
VI.6.5.4 La nature du sol ……………………………………………………….………. 121
VI.6.5.5 La topographie …………………………………………………….…….…….. 122
VI.6.5.6 Nature de la zone non saturée ………………………………………….……… 123
VI.6.5.7 La perméabilité ………………………………………………………….…….. 124
VI.6.5.8 Carte de synthèse …………………………………………………..………….. 125
VI.6.6 Conclusion………………………………………………………………….….. 127
CONCLUSION GENERALE et RECOMMANDATIONS 128
REFERENCES BIBLIOGRAPHIQUES 131
ANNEXES
Liste des abréviations et des paramètres
************ *********
-1-
Les explications des phénomènes et le détail des résultats seront présentés dans les
chapitres qui constituent cette thèse.
La première partie concerne l’étude des caractéristiques physiques de la plaine de
Meskiana. Elle comprend trois chapitres :
Le premier expose les caractéristiques physico-géologiques de la plaine.
Le deuxième détaille la géométrie de l’aquifère basée sur la prospection électrique.
L’élaboration de cette étude confirmera la disposition des formations géologiques et leur
susceptibilité en tant qu’aquifère. A la fin de ce chapitre, on a pu acquérir une vision plus
précise sur la géométrie et l’extension des couches qui pourraient être aquifères ainsi que sur
leur composition lithologique (qui influe sur le chimisme de la nappe).
Le troisième chapitre regroupe enfin une étude des différents paramètres climatiques,
tels que le type de climat, l’évolution des précipitations, l’évolution des températures,
l’estimation de l’évaporation et le calcul du bilan hydrologique par la méthode de
G.W.Thornthwaite. L’estimation des écoulements de surface de l’Oued Meskiana a également
été abordée au cours de ce chapitre.
-2-
I.CADRE GENERAL
I.1-Situation géographique :
Le bassin de Meskiana s’étend sur une superficie de 1680 km² et se localise au Nord-
Est du territoire Algérien, à 80 km de la frontière Algéro-Tunisienne et à 250 km au Sud de la
ville de Annaba. Autrement dit, il se situe à mi-chemin entre la mer Méditerranée et le Sahara
algérien (figure 01) et fait partie du bassin de l’Oued Mellègue. Plusieurs agglomérations y
siègent, à savoir Bellala, Dalâa, Meskiana, Rehia et Zebar, comptant plus de 37000 habitants
selon un recensement fait en 1998. Cette population serait estimée à 48000 en 2004.
Les limites géographiques du bassin de Meskiana sont présentées comme suit :
-Limite Nord : Djebel Mesloula, Argoub El Mnachir et Djebel Khannaga ;
-Limite Nord-Ouest : Djebel Chettaya, Djebel Bou Thokhma, et Djebel Ahmar ;
-Limite Nord-Est : Djebel Belkfif ;
-Limite Sud-Est : Djebel Gouriguer, Djebel Es Stih, et Djebel Khemalal ;
-Limite Sud-Ouest : Djebel Boutoukhma et Chott Esbikha.
I.2-Réseau hydrographique :
Le bassin versant de Meskiana est traversé par un seul cours d’eau principal, qui est
l’Oued Meskiana. Il est considéré comme un collecteur des eaux de ruissellement qui affluent
sur ses deux rives.
Sur la rive droite, les principaux affluents (temporaires) prenant naissance à partir des
montagnes et se déversant dans l’Oued Meskiana sont : l’Oued Rbâa, l’Oued Hmadjra, l’Oued
Smar El Houd et l’Oued El Mallah.
Sur la rive gauche, on compte l’Oued Djdida et l’Oued Ain Sedjra.
I.3-Paramètres morphométriques :
En se basant sur l’assemblage des cartes topographiques au 1/50000 (Dalâa feuille
n°204, Meskiana feuille n°177, El Aouinet feuille n°150, Youkous Les Bains feuille n°205,
F’kirina feuille n°176 et Ain Beida feuille n°149), on a pu délimiter le bassin versant de
l’Oued Meskiana, dont les caractéristiques morphométriques se résument dans le tableau 1.
I.4-Activités économiques :
La plaine de Meskiana appartient à un domaine sub-saharien, qui bénéficie d’un climat
semi-aride à faibles précipitations et par conséquent le développement du couvert végétal est
limité. Les agglomérations n’ont été dotées d’aucune infrastructure industrielle, à l’exception
d’une usine de laines et textiles actuellement à l’arrêt, alors que l’agriculture et l’élevage se
sont bien développés.
-3-
En effet, cette région appartient, agro-écologiquement, aux hautes plaines intérieures
orientales céréalières, qui ont des microclimats généralement contraignants (hiver froid à très
froid, été chaud à très chaud et sec). La pluviosité est généralement limitée et l’eau constitue
l’élément clé du système de culture mis en place. Durant les années pluvieuses où l’eau
d’irrigation est disponible, on rencontre les cultures maraîchères, les cultures fourragères et
les céréales avec un assolement biennal et parfois triennal, alors que durant les années de
sècheresse, les cultures maraîchères et fourragères seront réduites et on se contente des
céréales pour seulement équilibrer les revenus des agriculteurs.
Ces dernières années, le secteur de l’agriculture et de l’élevage a connu un
épanouissement important suite au lancement des programmes de PNDA en 2000 et celui du
FNRDA en 2001. Ces programmes ont provoqué l’extension des terres irriguées et multiplié
les cycles culturaux.
Rehia
Bellala
Meskiana
Zebar
Dalâa
0 10 20 Km
-4-
Chapitre I :
GEOLOGIE
I.1. INTRODUCTION :
I.2.1 TRIAS :
Le Trias est représenté par des masses chaotiques sans stratification (J.M. VILA,
1977), dans lesquelles divers éléments sont reconnaissables tels que des cargneules jaunâtres,
des dolomies noires, des calcaires, des gypses marneux et fibreux. Tous ces éléments sont
couverts par des argilites vertes, marquant ainsi le Trias du constantinois. Ces roches
triasiques forment des diapirs, qui sont partout en contact anormal avec les couches plus
récentes. :
-Au Nord, à proximité d’El Aouinet, au niveau de Mesloula, Argoub Ez Zitoun et Argoub El
Houra, le Trias se trouve en contact anormal avec les récifs albo-aptiens.
-Au Nord-Est (Djebel Belkefif), le Trias est en contact anormal avec les mêmes récifs et
également avec le Cénomanien et le Turonien (figure 03, coupe AA).
-5-
925 930 935 940 945 950 955 960 965
290
285
Rehia
280
Bellala
if
kef
Bel
275
Dj
ar
hm
270
Meskiana
eA
ern
gu
265
Dj
260
Zebar
ma
er
igu
kh
255
er
To
Gu
u
Bo
Dj
Dj
250
Dalâa
245
240
235 Echelle
0 10 20 Km
Légende
-6-
In Hemila 1988 modifié
Forage hydraulique
Légende:
Faille
Quaternaire : Alluvions actuelles et récentes. Contact stratigraphique
Santonien : Marne grise à rare intercalation de marno-calcaires. Trias : Cargneule, gypse marneux, dolomie et grès.
-7-
I.2.2 CRETACE :
I.2.2.1 Aptien :
Il s’agit de formations de dimensions assez réduites affleurant au Nord, qui sont en
contact anormal avec des pointements triasiques, au niveau de Mesloula, Djebel Belkefif,
Argoub Ez Zitoun et Kodiat El Hadjel. Ces formations sont représentées par des assises
carbonatées qui sont recouvertes par le Miocène (figure 3, coupe AA’).
Les formations aptiennes comprennent deux faciès, le premier marno-calcaro-gréseux
et le second à calcaires récifaux (DOZET, 1985). Dans le premier faciès, on rencontre des
marnes jaunâtres et verdâtres, des calcaires gris à orbitolines et des grès quartzeux jaunes. On
distingue la partie sous récifale qui se dispose au dessous des calcaires récifaux et la partie
post-récifale qui se trouve au dessus des calcaires récifaux. Les couches de l’Aptien sous
récifal ne sont nulle part entièrement exposées. Les couches de l’Aptien post récifal de
Mzouzia ont une puissance de 55 m.
Le faciès des calcaires récifaux comprend des calcaires gris à grains fins, des calcaires
à Huîtres, des calcaires à orbitolines, des calcaires coralliaires, des calcaires bréchiques, des
calcaires oolithiques et des dolomies. La puissance des calcaires récifaux est de 10 m à 180 m.
Les fossiles les plus rencontrés dans l’Aptien sont les huîtres, divers lamellibranches,
des gastéropodes, des ammonites et des coraux
I.2.2.4 Turonien :
Il est observé au Nord-Est à Douar El Machtab et il est représenté par :
-Turonien basal : est formé d’une série épaisse de marnes grises, des marnes calcaires en
plaquettes, des marno-calcaires, des calcaires marneux brunâtres avec un épaisseur pouvant
atteindre 270 m (figure3, coupe BB’). Dans les marnes calcaires et les calcaires marneux, de
nombreuses empreintes d’inocérames ont été déterminées. Les marnes argileuses renferment
de grosses globigérines.
-Turonien moyen et supérieur: sont représentés par des marnes beiges à grises, qui alternent
avec des calcaires zoogènes noduleux.
-8-
partie basale du Campanien représentent du point de vue lithologique un complexe homogène
sans repères.
- Le Coniacien est formé par des alternances de marnes beiges ou grises et des calcaires gris
biodétritiques et grumeleux. L’épaisseur du Coniacien est de 90 m (figure3, coupe CC’).
- Le Santonien se compose de marnes grises et de calcaires marneux à inocérames, avec une
épaisseur d’environ 100 m.
I.2.2.6 Campanien :
Il affleure sur les deux bordures occidentale, à Djebel Bou Thokma, et orientale, à
Djebel Gourigueur.
- Le Campanien inférieur : à sa base, on rencontre une unité lithologique de couches
marneuses composées essentiellement de marnes argileuses. Cette dernière est d’une
puissance de 100 m. Le sommet du Campanien inférieur est représenté par une unité
lithologique la plus caractéristique et la plus marquante du Campanien, qui se compose de
calcaires crayeux gris et biomicrites, bien stratifiée avec de nombreuses intercalations de
marnes. Des empreintes d’inocérames sont renfermées dans ces calcaires et ces marnes. Les
couches constituant le sommet du Campanien inférieur ont une épaisseur qui varie entre 45 m
et 100 m.
- Le Campanien supérieur : dans la partie inférieure, sur une puissance de 25 m,
alternent des couches de marnes calcaires et de marnes. Ces dernières sont surmontées par une
série d’environ 130 m d’épaisseur constituée de marnes argileuses gris verdâtres
I.2.2.7 Maestrichtien :
Les calcaires maestrichtiens forment les monts du Djebel Gourigueur, Djebel Bou
Tokhma et Djebel Garne Lahmar. Ces calcaires sont gris clairs, crayeux, stratifiés et
contiennent assez souvent des rognons et des intercalations de silex.
- Le Maestrichtien inférieur est marqué par une dominance marneuse : Ce sont des marnes
gris-beiges, suivies de marnes beiges et de biomicrites argileuses, en bancs décimétriques. La
série ne dépasse pas 40 m d’épaisseur.
- Le Maestrichtien supérieur est représenté par une formation calcaire de 40 m de
biomicrosparites grises surmontée par 40 m de biomicrites. Signalons que l’épaisseur du
Maestrichtien à Djebel Gourigueur est de l’ordre de 460 m.
I.2.3 PALEOGENE
I.2.3.1 Paléocène :
Il est observé seulement à la limite orientale, entre Djebel Guerigueur et Djebel
Serdies. Il est représenté par une série de 50 m de puissance, constituée de marnes gris
foncées qui sont parfois argileuses.
I.2.3.2 Eocène :
Il a été reconnu sur les flancs Est et Ouest du vaste anticlinal de Meskiana, notamment
à Djebel Gourigueur et au synclinal de Ain Dalâa. Il est caractérisé par :
-des formations montiennes, constituant la base de l’Eocène, se composant de marnes gris
verdâtres ;
- des formations thanétio-yprésiennes qui se disposent au dessus des formations précédentes.
Elles se composent d’alternances de marnes gris beiges et de biomicrites finement gréseuses
et noduleuses et en bancs décimétriques. Ces dernières sont surmontées par des alternances de
micrites blanchâtres à nodules phosphatés et de biomicrites interclastiques en bancs métriques
ou décimétriques à grains phosphatés et renferment des gastéropodes et des lumachelles.
L’épaisseur des couches thanétio-yprésiennes est d’environ 200 m.
-9-
I.2.4 NEOGENE :
I.2.4.1 Miocène :
Les affleurements du Miocène sont de dimensions réduites et se répartissent dans la
partie Nord de la plaine surtout à Mesloula, Argoub El Manchar, Argoub Ez Zitoun, Kodiat El
Agab, Kodiat Rehia et Kodiat Naga.
-Le Miocène inférieur est représenté par des conglomérats avec des galets de calcaires
maestrichtiens à ciment sableux. Dans ce ciment, il y a de nombreuses coquilles d’huîtres.
-Le Miocène supérieur est représenté par le Tortonien qui affleure particulièrement à Kodiat
Mami. Il est marqué par une série de couches épaisses, constituées de grès quartzeux mal
cimenté à grains grossiers. Ces dernières se trouvent en discordance avec les couches
santoniennes, campaniennes et maestrichtiennes.
I.2.5 QUATERNAIRE :
I.2.5.1 Villafranchien :
L’aspect le plus fréquemment rencontré est celui des croûtes calcaires massives et
blanches. Ces formations perchées par rapport aux surfaces précédentes sont faciles à repérer,
probablement karstifiées ou dissoutes. Ces faciès représentent habituellement le
Villafranchien.
I.2.5.2 Glacis polygéniques nappant les reliefs:
Ces glacis s’étendent sur de grandes surfaces de la plaine en pente douce, ils
recouvrent le cycle antérieur de glacis, qui est facile à reconnaître car il est constitué de
croûtes calcaires massives.
I.2.5.3 Alluvions actuelles et récentes:
Ces alluvions ont une extension limitée à la vallée de Meskiana. Elles sont constituées
de limons fins, de sables fins, de graviers et d’alluvions argileuses qui se déposent des deux
cotés de l’Oued Meskiana.
Les différentes formations géologiques du bassin de Meskiana se résument dans la figure 04.
- 10 -
Etage Epaisseur (m) Lithostratigraphie Lithologie
QUATERNAIRE
Coniacien 150
Marne grises et alternance de marne beige et calcaire
gris .
T uronien moy et sup 150
Alternance de marne gris e à beige et calcaire
zoogène noduleux.
Marne gris e ou noire avec une intercalation de
CRETACE
I.4 GEOMORPHOLOGIE :
La géomorphologie est caractérisée par l’existence :
- des massifs de bordures limitant la plaine ;
- d’une vallée se situant en bordure Nord d’un anticlinal crétacé, dont le cœur est constitué de
calcaires aptiens. Au Nord-Est de la dite vallée, le Trias apparaît en diapirs, en recoupant les
marnes cénomaniennes et les marno-calcaires turoniens. Un recouvrement alluvial s’étale au
niveau de l’axe de l’anticlinal, alors que les éboulis recouvrent les superficies aux pieds des
reliefs.
- 11 -
I.5 APERCU TECTONIQUE
La plaine de Meskiana appartient à l’Atlas saharien (J.M. VILA, 1980). Elle est
caractérisée par des éléments tectoniques simples, tels que les plis de direction NE-SW et les
dislocations qui ont souvent un cours transversal à la direction de ces plis. On rencontre
également des structures particulières, qui se sont survenues sous l’influence de la tectonique
diapirique. La période crétacée a été, du point de vue tectonique, assez calme (DOZET,
1985). Des grandes activités qui se sont déroulées durant cette période, on cite une subsidence
continue au Sud du bassin et les soulèvements relatifs de quelques parties septentrionales de
ce dernier (figure 05).
Légende
Axe anticlinal
Allochtone indifférencié
Axe synclinal
Principales failles Mio-Pliocène des fossés
Charriage frontal de Autochtone post-triasique
l’allochtone
Trias
Ecaillage
- 12 -
Entre le Crétacé et le Miocène, deux principales phases tectoniques ont eu lieu pour
former de nombreuses structures. Durant la première phase, de vastes plissements se sont
produits en donnant naissance aux anticlinaux et synclinaux suivants : les plis de l’Oued
Smara, El Houd, le synclinal de Djebel Belkefif, le synclinal de Kodiat Zennbaî, le synclinal
de Reguibet El Assel, le synclinal de Djebel Gourigueur et l’anticlinal de Kessar El Tir. La
phase des plissements fut suivie de la phase de la tectonique cassante où de nombreuses fortes
failles et dislocations longitudinales et transversales sont survenues, telles que celle de Djebel
Gourigueur Ain Dalâa, dont le rôle du point de vue hydrogéologique semble très important.
I.6 CONCLUSION :
* A partir de la lithologie des formations géologiques, on peut avoir une idée sur la
minéralisation des eaux souterraines. Ainsi,
-Les calcaires donnent un faciès généralement bicarbonaté calcique ou magnésique ;
-Les marnes sont à l’origine des faciès sulfatés ;
-Les alluvions donnent plusieurs faciès non salés.
On note que la présence du Trias peut provoquer la salure des eaux souterraines.
- 13 -
Chapitre II :
GEOPHYSIQUE
II.1 INTRODUCTION :
L’exploration et l’exploitation des couches aquifères d’une formation géologique donnée
sont étroitement liées à la connaissance de la nature lithologique de cette formation.
L’identification de la géométrie de ces couches aquifères nécessite des moyens de reconnaissance,
tels que la prospection électrique et les forages mécaniques à titre d’exemple. L’utilisation de ces
moyens de reconnaissance va permettre la connaissance des dimensions de ces aquifères, et par
conséquent aider à la bonne gestion de ces ressources hydriques.
290
285
280
Rehia
Bellala
275
BV
A
270 Meskiana
B
C
265 CD
260 E
Zebar
F
G
255 Légende
H
K
250 X A Profil de sondages électriques.
L
Figure 06: Carte de position des sondages électriques verticaux et des profils.
- 14 -
II.2.1 Interprétations des résultats de la prospection électrique
Ces deux campagnes géoélectriques, basées sur 25 profils de différentes longueurs AB ( 4000
m, 3000 m, 2000 m, 600 m et 100 m ), visaient à atteindre les objectifs suivants :
- Reconnaître la profondeur d’un éventuel substratum calcaire résistant du Crétacé inférieur sous
le recouvrement marneux et marno-calcaire du Cénomanien ;
- Identifier la géométrie du remplissage alluvial et déterminer les zones de meilleures
transmissivités dans ce dernier ;
- Reconnaître l’allure de l’anticlinal crétacé de Meskiana et celle du synclinal éocène perché de
Dalâa.
L’examen des coupes géoélectriques (figure 07a et 07b) confirme la structure anticlinale
dans la vallée de Meskiana. En effet, cette structure a été mise en évidence à partir des profils G,
H et K (leurs coupes géoélectriques sont en annexeI), où on observe la forme de plis pincés d’axe
orienté vers le Sud-Est de la plaine. A partir du profil L, cette structure est perturbée par la
présence d’une couche très conductrice (résistivité inférieure à 2 ﺟ.m).Sous le recouvrement Mio-
Plio-Quaternaire, la structure anticlinale se caractérise par la présence de faciès conducteurs très
épais, qui s’intercalent avec des niveaux résistants de fortes résistivités (70 à 300 ﺟ.m).
Les faciès conducteurs représentent des formations marneuses et argileuses dont l’âge
s’étale du Cénomanien jusqu’au Campanien. Les niveaux résistants représentent probablement
des calcaires et calcaires marneux du Turonien inférieur et supérieur, et des calcaires du
Coniacien qui affleurent en bordures et se prolongent sous le recouvrement alluvial (Rapport
C.G.G, 1974).
Le synclinal perché de Dalâa est présenté par deux niveaux très résistants, dont la
résistivité dépasse les 300 ﺟ.m, et qui sont séparés par un niveau conducteur, de 4 à 20 ﺟ.m. Le
premier niveau conducteur est très épais (environ 350 m) et se prolonge sous le recouvrement
alluvial, on pense qu’il est constitué par des calcaires thanétiens-yprésiens. Le second niveau
résistant représente les calcaires maestrichtiens. Le niveau conducteur marneux, qui sépare ces
derniers, serait d’âge Dano-Montien, son épaisseur est de l’ordre de 100 m (Rapport C.G.G,
1974).
- 15 -
Echelle
Echelle
Légende
11 12 13 1 8 2 100 3 4
- 17 -
930 935 940 945 950 955 960
290
285
Rehia
280 Bellala
275
270 Meskiana
265
260
Légende
Zebar
255
85m
250 Courbe isopaque des alluvions.
Dalâa Agglomération.
245
240
Echelle
0 10 20 Km
290
285
Rehia
280
Bellala
275
Meskiana
270
265
260 Légende
Zebar 850m
255 Courbe isosubstratum.
Agglomération.
250
Dalâa
245
Echelle
240 0 10 20Km
représentés d’une part par des couches conductrices constituées par des marnes et des
argiles, et d’autre part par des couches résistantes formées de sables, de graviers, de
galets et de cailloutis. Ces dernières couches constituent la nappe phréatique du
• une structure synclinale, qui se manifeste au niveau de la zone de Dalâa, est constituée
jusqu’au Campanien. Le cœur de cet anticlinal est d’âge Turonien.
- 19 -
Chapitre III :
HYDROCLIMATOLOGIE
III.1 INTRODUCTION :
- 20 -
U : Moyenne des précipitations mensuelles de la station à ajuster (Meskiana) pour la période
(1972/1973 – 1982/1983), où elle possède une série complète et un écart type relativement
faible.
La procédure consiste à estimer le rapport b à partir des données homogènes d’une même
période, puis à appliquer l’équation (1) à la période à vérifier (1983/1984- 2004/2005).
L’estimation du coefficient de correction est consignée dans le tableau 04 :
14000
Cumul des précipitationsà la station
12000
de Meskiana (mm)
10000
8000
6000
4000
2000
0
0 2000 4000 6000 8000 10000 12000 14000
Cumul des précipitations à la station de Tébessa (mm)
Les calculs statistiques montrent que la précipitation moyenne annuelle est de 341 mm.
Ceci est vérifié en se référant à la carte des précipitations moyennes annuelles, établie par
- 21 -
l’ANRH sur une période de 60 ans (figure 11). La précipitation maximale est enregistrée aux
mois de décembre et mars. La précipitation minimale est enregistrée aux mois de juin et
juillet. La saison printanière est la plus pluvieuse avec 103 mm suivie par la saison automnale
avec une moyenne de 99 mm.
MESKIANA
Zone
d’étude
- 22 -
P (mm) T (°C)
50 25
40 20
Saison sèche
T en (°C)
30 15
P (mm)
Saison humide
20 10
10 5
0 0
v
ai
nv
ut
v
t
ril
il
pt
in
rs
oc
no
dé
fé
ju
m
ao
ju
m
se
av
ja
Mois
- 23 -
260 8 70 60
0 0 50
240
Précipitations (mm)
0
I=
220 P Ecoulement abondant
T + 10
0
200 40
0
180
0 Ecoulement Exoréique
160
0 30
140
0 Zone tempérée, drainage extérieur,
120 irrigation indispensable
0
100
20
0
80 Régime semi-aride
0 Ecoulement temporaire
60 Formations herbacées
0 10
40
0 Meskiana Régime désertique, écoulement
20 temporaire, drainage intérieur
0 Endoréisme, hyperaridité, aréisme
Températures (C°)
0 2 4 6 8 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4
0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0
Figure 13 : Abaque de l’indice d’aridité annuel de De Martonne.
I=
12.P
T ° + 10
(3) où
- 24 -
Selon la classification précédente par De Martonne, les mois de septembre, octobre,
novembre, décembre, janvier, février et avril présentent un climat semi-aride. Le mois de mars a
un climat tempéré. Les mois de juin, et août sont caractérisés par un climat steppique, alors que le
mois de juillet devient hyperaride.
260 8 70 60
0 0 50
240
Précipitations (mm)
0
220
i=
12 P Ecoulement abondant
0
t + 10
200 40
0
180
0 Ecoulement Exoréique
160
0
140 30
0 Zone tempérée, drainage extérieur,
120 irrigation indispensable
0
100
20
0 Régime semi-aride
80
0 Ecoulement temporaire
60 Formations herbacées
0 Mar Nov Oct Mai Sept 10
40
Jan Fev λ
Avr λ λ
Juin Régime désertique, écoulement
0 λ λ
20 λ
Août temporaire, drainage intérieur
λ λ λ λ λ
0 Déc. Juil Endoréisme, hyperaridité, aréisme
Températures (C°)
λ
0 2 4 6 8 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4
0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0
Figure 14: Abaque de l’Indice d’aridité mensuel de De Martonne (1923).
2
P : précipitation moyenne annuelle (mm)
M : moyenne des maxima du mois le plus chaud (degrés/Kelvin)
m : moyenne des minima du mois le plus froid (degrés/Kelvin)
L’application numérique donne les résultats suivants :
P = 340.99 mm
M= 306.51 °K
m = 273,94 °K
Ce qui donne :
Q2 = 36,07
Le repport de notre station sur la figure 15 permet de dire que le climat de la plaine de
Meskiana est de type semi-aride. Figure 15.
- 25 -
Humide
Subhumide
Semi-aride
Meskiana
Aride
Saharien
Pour
P = 341 mm
T = 15.58 °C
Ce qui donne
IM = 1.16
On peut dire que le climat de la plaine de Meskiana est à la limite sec-humide, c’est-à-dire qu’il peut
être qualifié de semi-aride.
III.3.6 Comparaison des indices climatiques :
Sachant que les éléments climatiques les plus importants pour caractériser un climat sont les
précipitations, les températures et l’évaporation ; les méthodes combinant ces éléments sont plus
intéressantes que celles qui incluent un seul élément. Toutes ces méthodes ont qualifié le climat de la
plaine de Meskiana comme étant semi-aride, ce qui est en conformité avec la carte des étages
bioclimatiques, établie par M. Côte en 1998 (figure16).
- 26 -
Mer Méditer ranée
Skikda Annaba
Jijel
Bejaia
Constantine Guelma
U N I S I E
S.Ahras
Setif
Bordj Bou Arreridj
Ain Beida
T
M'sila Batna Meskiana
Khenchela Tébessa
Légende Chott el Hodna
Domaines bioclimatiques.
Domaine humide
Domaine sub-humide
Domaine semi-aride Biskra Zone d’étude
Domaine semi aride inférieur
Domaine sub-aride 0 50 Km
Domaine aride
Figure 16: Carte simplifiée des zones bioclimatiques de l’Est algérien Côte M., 1998).
III.4 LA TEMPERATURE :
La température est un élément fondamental du climat. Elle est liée à la radiation solaire. Sa
variation influe sur la transformation des eaux en vapeur, que ce soit à la surface ou dans le sous-
sol. De ce fait, elle influe sur le degré d’évapotranspiration et par conséquent elle agit sur le taux
de salinité des eaux. Toutefois, la température a un rôle important dans la variation des
composantes du bilan hydrologique.
- 27 -
20
15
T en °C
10
0
1972/1973
1973/1974
1974/1975
1975/1976
1976/1977
1977/1978
1978/1979
1979/1980
1980/1981
1981/1982
1982/1983
1983/1984
1984/1985
1985/1986
1986/1987
1987/1988
1988/1989
1989/1990
1990/1991
1991/1992
1992/1993
1993/1994
1994/1995
1995/1996
1996/1997
1997/1998
1998/1999
1999/2000
2000/2001
2001/2002
2002/2003
2003/2004
2004/2005
Années
30
25
20
T en C°
15
10
5
0
octobre
novembre
aout
mars
décembre
janvier
septembre
mai
avril
juin
février
juillet
Mois
III.5 LA PLUVIOMETRIE :
L’étude pluviométrique présente un intérêt considérable en hydroclimatologie, car elle permet
de décrire le régime d’écoulement.
- 28 -
600
500
400
P (mm)
300
200
100
0
1970/1971
1971/1972
1972/1973
1973/1974
1974/1975
1975/1976
1976/1977
1977/1978
1978/1979
1979/1980
1980/1981
1981/1982
1982/1983
1983/1984
1984/1985
1985/1986
1986/1987
1987/1988
1988/1989
1989/1990
1990/1991
1991/1992
1992/1993
1993/1994
1994/1995
1995/1996
1996/1997
1997/1998
1998/1999
1999/2000
2000/2001
2001/2002
2002/2003
2003/2004
2004/2005
Figure 19 : Evolution des précipitations annuelles à la station de Meskiana pour la
période (1972/1973-2004/2005).
50
40
P en (mm)
30
20
10
0
octobre
novembre
aout
mars
septembre
décembre
janvier
avril
mai
juin
février
juillet
Mois
- 29 -
ETP mm/mois = 0.40[T/T +15](Ig + 50) (6)
Dans laquelle on a
T : température moyenne mensuelle en °C, mesurée sous abri,
Ig : radiation solaire globale exprimée en cal/cm2/jour.
Ig = Iga (0.18 + 0.62 h/H).
Iga : radiation maximale théorique tirée des abaques (voir annexes).
La constante 0.4 est remplacé par 0.39 pour les mois de 30 jours et par 0.37 pour le mois de
février.
H : durée astronomique du jour, définie à partir des abaques (voir annexes).
h :durée d’insolation effective mesurée par l’héliographe.
Quand l’humidité relative mensuelle de l’atmosphère est inférieure à 50%, on multiplie la valeur
de l’ETP par le coefficient de correction suivant :
k = 1 + [(50 – hr)/70] avec
hr :humidité relative de l’air durant le mois en question.
Les résultats de l’application de cette formule sur la plaine de Meskiana sont présentés dans le
tableau 07.
Tableau 07: Valeurs de l’ETP estimées par la méthode de Turc.
Sept Oct Nov déc Jan Fév Mar Avr Mai Juin Juil Août Année
T (°C) 21.47 16.65 10.92 7.55 6.27 7.54 10.06 13.04 18.19 23.16 26.45 25.7
Iga 630 593 456 393 430 556.2 714.5 852.6 943.7 973.5 955.5 864.6
H en heures 372 352.3 334.4 304.2 315.3 308.4 373.3 390.4 432 432.2 442.3 416
h en heures 256 215 172 138 136.1 136.7 169.3 204 285.7 277.3 333.2 305.
5
Ig 342.8 329 196.1 181.3 192.4 253 329.1 429.6 556.6 562.5 618.3 549.
3
Hr 50 56.4 65.7 67.5 65.3 67.8 56.7 45.8 49.5 41.2 35.6 39.5
ETP en mm 101.68 79.75 41.46 30.9 28.5 37.5 60.8 89.2 132.7 148.2 170.6 151. 1064.6
3 8
I = Σi et i = (t/5)1.514
a = (1.6.I/100) + 0.5 (D’après Serra L. 1954)
Les valeurs obtenues sont corrigées par un coefficient de correction k qui dépend de la latitude du
bassin. Tableau 08.
- 30 -
Tableau 08 : Valeurs de l’ETP estimées par la formule de Thornthwaite.
Sept Oct Nov Déc Jan Fév Mar Avr Mai Juin Juil Août Total
T (°C) 21.47 16.65 10.92 7.55 6.27 7.54 10.06 13.04 18.19 23.16 26.45 25.69
ETPnc (mm) 96.21 63.2 31.52 17.95 12.67 17.26 27.47 41.81 73.30 108.2 134.7 129.3
K 1.03 0.97 0.85 0.83 0.86 0.84 1.03 1.1 1.21 1.22 1.24 1.16
ETPc (mm) 99.1 61.3 26.8 14.9 10.9 14.5 28.3 46 88.7 132 167 150 839.5
Il est utile de comparer les résultats obtenus par les des deux méthodes utilisées afin de
connaître celle qui s’adapte le mieux à la région. Le graphe relatif (figure 21) à la station de
Meskiana a permis de conclure que les résultats auxquels nous avons aboutis sont relativement
égaux. Pour notre région, nous allons appliquer la méthode de Thornthwaite, qui est
universellement appliquée pour le calcul des bilans hydrologiques.
180
160
140
ETP en (mm)
120
100
80
60
40
20
0
S ept Oct Nov Déc Jan fév Mar Avrl Mai Juin Juil Août
Mois
- 31 -
λ = 1/ (0.8 + 0.14T)
P :précipitation moyenne annuelle en (m).
L’application de cette formule est conditionnée par la relation suivante :
1/8λ ≤ 1/2λ
Soit ETR = 303 mm
Sept Oct Nov Dec Jan Fév Mar Avr Mai Juin Juil Août Total
P(mm) 38.86 28.21 31.88 25.86 26.02 24.9 39.3 28.6 35.48 25.15 9.81 26.81 340.99
6 5
T(C°) 21.47 16.65 10.92 7.55 6.27 7.54 10.0 13.0 18.19 23.16 26.45 25.69
6 4
ETP(mm) 99.1 61.3 26.8 14.9 10.9 14.5 28.3 46 88.7 132 167 150 839.5
ETR(mm) 38.86 28.21 26.88 14.9 10.9 14.5 28.3 46 44.08 25.15 9.81 26.81 314.32
EX(mm) 00 00 00 00 5.16 10.4 11.0 00 00 00 00 00 26.67
6 5
DA(mm) 60.24 33.09 00 00 00 00 00 00 44.62 106.85 157.19 123.19 525.18
- 32 -
Tableau 10: Bilan hydrologique pour l’année la plus humide à la station de Meskiana, par la méthode de
Thornthwaite (2002/2003).
Sept Oct Nov Dec Jan Fév Mar Avr Mai Juin Juil Août Total
P(mm) 32.6 77.9 59.3 36.5 53.5 25.5 31 64.1 83.7 30.9 16.2 36.8 548
T(C°) 21.6 17.2 11.2 8.4 6.4 6.2 10.6 14.1 18.5 24.4 27.6
ETP(mm) 100.81 63.24 31.64 18.0 11.3 10.9 28.1 47.7 74.72 124.8 177.06 150.48 839.01
4 2 2 9 9
P-ETP - + + + + + + + + - - -
68.21 14.66 27.66 18.4 42.1 14.5 2.81 16.3 08.98 93.90 160.86 113.68
6 8 8 1
RFU26(mm) 00 14.66 26 26 26 26 26 26 26 00 00 00
ETR(mm) 32.6 63.24 31.64 18.0 11.3 10.9 28.1 47.7 74.72 56.9 16.2 36.8 428.36
4 2 2 9 9
EX(mm) 00 00 16.32 18.4 42.1 14.5 2.81 16.3 8.98 00 00 00 119.64
6 8 8 1
DA(mm) 68.21 00 00 00 00 00 00 00 00 67.9 160.86 113.68 410.65
Tableau 11: Bilan hydrologique pour l’année la plus sèche à la station de Meskiana, par la méthode de
Thornthwaite (1993/1994).
Sept Oct Nov Dec Jan Fév Mar Avr Mai Juin Juil Août Total
P(mm) 16.05 26.2 9.9 19.5 44.2 20.9 9.4 13.3 31.7 50.6 14.5 6.7 262.95
T(C°) 22.3 19.1 11.1 7.7 7.5 8.9 11.9 11.8 21.9 24.2 27 28.6
ETP(mm) 105.61 74.69 29.80 14.5 14.0 19.0 33.1 33.5 99.36 122.36 139.26 160.65 846.08
5 2 8 2 8
P-ETP - - - + + + - - - - - -
89.56 48.49 19.90 4.95 30.1 1.82 23.7 20.2 67.66 71.76 124.76 153.95
8 2 8
RFU26(mm) 00 00 00 4.95 26 26 2.28 00 00 00 00 00
ETR(mm) 16.05 26.2 9.9 14.5 14.0 19.0 33.1 15.5 31.7 50.6 14.5 6.7 252
5 2 8 2 8
EX(mm) 00 00 00 00 9.13 1.82 00 00 00 00 00 00 10.95
DA(mm) 89.56 48.49 19.9 00 00 00 00 18 67.66 71.76 124.76 153.95 594.08
- 33 -
180
160
140 P (mm)
Lame d'eau (mm 120 ETP (mm)
100 ETR (mm)
80 RFU (mm)
60 EX (mm)
40 DA (mm)
20
0
Nov
Mar
Avrl
fév
Août
Jan
Sept
Oct
Juin
Déc
Juil
Mai
(Mois)
R = P3/3.ETP2 (11) où
R : Ruissellement (m).
P : Précipitations moyennes mensuelles (m).
ETP : Evapotranspiration potentielle (m).
L’application numérique donne
R = 22 mm
l’infiltration (I) :
I = EX – R (12) où
I = 5 mm
Il est à noter que cette valeur ne représente pas réellement l’infiltration dans tout le bassin, car son
estimation est basée sur les données d’une seule station pluviométrique, qui est celle de Meskiana.
Les mesures d’infiltration in situ, réparties sur toute superficie de la plaine, pourrait nous donner
une approche fiable de ce paramètre.
- 34 -
III.7.4 Commentaire sur le bilan hydrologique :
Le calcul du bilan hydrologique permet de faire les constatations suivantes :
-La pluie ne permet pas un stockage permanent de l’eau dans le sol, ce qui induit un déficit
agricole plus ou moins long. Ainsi, la saturation de la RFU, pour une année moyenne, ne dure que
03 mois de janvier à mars. Cette saturation ne dure que 02 mois pour une année sèche et peut
atteindre 07 mois pour une année humide. Cette situation induit un déficit agricole pendant 06
mois (de mai à octobre). Pour l’année la plus sèche, le déficit agricole est enregistré pendant 07
mois (de mai à novembre). Les sols agricoles demandent donc une irrigation par les eaux de
surface, qui vont produire une recharge supplémentaire de la nappe et une modification de la
composition chimique des eaux.
L’enregistrement de ce déficit agricole coïncide avec la période sèche.
-L’ETR ne couvre que 37% de l’ETP. Elle consomme 92% des précipitations. Ceci induit un
faible excédent et par conséquent la recharge de la nappe est très limitée.
III.8 LES EAUX DE SURFACE
III.8.1 Estimation des eaux de surface :
Le bassin de Meskiana, d’une superficie de 1680 km2 , est drainé par l’Oued Meskiana et
ses affluents. Ce dernier s’écoule du Sud vers le Nord et traverse ainsi la totalité de la plaine. La
confluence, au Sud d’El Aouinet, de l’Oued Meskiana et de l’Oued Ksoub constitue l’Oued
Mellègue. Il faut signaler que l’Oued Meskiana prend naissance au Sud d’Ain Dalâa et est
alimenté, principalement, par plusieurs sources émergeant des calcaires de bordures telles que Ain
Sedjera et Ain Bedjene. Le régime d’écoulement est semi-permanent, en relation avec le climat de
la région.
Malheureusement, ce cours d’eau n’est doté d’aucune station hydrométrique qui permet de
suivre son régime d’écoulement. L’estimation de l’écoulement s’est basée sur les deux stations
hydrométriques d’El Aouinet et celle de Morsott par différence de débit (figure 23).
III.8.2 Ecoulement de surface et régime hydrologique :
En se basant sur les données illustrées dans le tableau 13, nous avons réalisé la courbe d’évolution
des débits instantanés (figure 24). L’interprétation de la dite courbe montre que les fortes valeurs
sont enregistrées pendant la période des hautes eaux (0.529 m3/s au mois de novembre, 0479
m3/s au mois d’avril) et les plus faibles valeurs marquent la période des basses eaux avec 0.208
m3/s au mois de juin et 0.186 m3/s au mois de juillet.
Tableau 13: Estimation des débits instantanés en (m3/s) de l’Oued Meskiana pour la période
(2000/2001-2004/2005).
Sept Oct Nov Déc Jan Fév Mar Avr Mai Juin Juil Août Total
00/01 0.03
0.462 0.454 0.296 0.315 0.383 0.335 0.325 0.234 0.423 0.123 6 0.041 3.427
01/02
0.474 0.045 0.520 0.158 0.305 0.341 0.463 0.415 0.224 0.048 0.32 0.245 3.558
02/03
0.441 0.910 0.875 0.467 0.653 0.403 0.440 0.880 0.934 0.435 0.30 0.381 7.127
03/04
0.545 0.471 0.922 1.126 0.302 0.290 0.495 0.335 0.378 0.037 0.04 0.292 5.234
04/05
0.425 0.450 0.032 0.430 0.446 0.333 0.453 0.531 0.393 0.372 0.22 0.731 4.821
Moy
0.469 0.466 0.529 0.499 0.418 0.340 0.435 0.479 0.470 0.203 0.18 0.338
- 35 -
S.Ahras
S.Ain Dalia
S.Ouenza
S.El Aouinet
Légende
S.Ain Erkel
S.Morsott
Echelle
0 10 20 Km
On note que l’évolution des débits montre un certain parallélisme avec celle des
précipitations (figure 25), ceci montre l’influence de ces dernières sur les écoulements de surface.
En se basant sur l’évolution du coefficient mensuel de débit, on remarque que la période des
hautes eaux s’étale du mois de septembre jusqu’au mois de mai et celle des basses eaux
commence au mois de juin et se termine au mois d’août (figure 26).
600
500 Légende
400
Q (l/s)
300 Q(l/s)
200
100
0
sep oct nov dec jan fev mar avr mai juin juil aout
Mois
50 600
40 500 Légende
400
P (mm )
Q (l/s)
30 P ( mm )
300
20 Q(l/s)
200
10 100
0 0
sep oct nov dec jan fev mar avr mai juin juil aout
Mois
1,4
1,2
Légende
1
1
CMD
0,8
0,6 CMD
0,4
0,2
0
sep oct nov dec jan fev mar avr mai juin juil aout
Mois
III.9 CONCLUSION
La plaine de Meskiana possède un climat méditerranéen de type semi-aride, avec un hiver
froid et pluvieux, et un été chaud et sec, influencé par des facteurs géographiques (éloignement
de la mer, altitude et orientation des reliefs). La répartition des précipitations et des températures
influe sur la chimie des eaux par dilution pendant la période pluvieuse, ou concentration lors de
l’évaporation durant la période sèche.
La température moyenne annuelle qui y règne est de 15,58°C.
Les précipitations moyennes annuelles, sur une période de 33 ans (72/73 – 2004/2005), sont
estimées à 341 mm. L’évapotranspiration réelle est estimée à 314 mm, mobilisant 92 % des
précipitations.
Le ruissellement ne représente que 6.45 % des précipitations, ce qui représente 22 mm.
L’infiltration ne laisse percoler que 5 mm, ce qui ne représente que 1.55% des précipitations.
Ceci laisse penser que l’alimentation de la nappe ne se fait pas exclusivement à l’aide de
l’infiltration efficace, mais aussi sur les calcaires fissurés qui bordent la nappe.
Le régime des écoulements de l’Oued Meskiana dépend de celui des précipitations.
- 37 -
Chapitre IV :
HYDROGEOLOGIE
IV.1 INTRODUCTION:
Les formations aquifères, tant par leur nature lithologique que par leur extension dans
l’espace, occupent une place importante dans notre présente étude, car elles influent sur le
fonctionnement hydrogéologique de celles-ci. Ainsi, ce travail a pour but de comprendre
l’impact des variations climatiques, dans le temps, et la variation de la lithologie dans
l’espace, sur le comportement hydrodynamique de la nappe alluviale de Meskiana.
Les relevés piézométriques, effectués pendant 06 campagnes (03 en basses eaux et 03
en hautes eaux) et la détermination des paramètres hydrodynamiques de l’aquifère (la
transmissivité et la perméabilité) renseignent sur les caractéristiques de l’écoulement et le
comportement de la nappe, permettent de cerner le risque de pollution et de protéger la nappe
aquifère.
IV.2GEOMETRIE DE L’AQUIFERE :
L’interprétation des résultats de la prospection électrique (CGG, 1974) et des données
acquises suite aux sondages mécaniques a permis de visualiser la répartition du matériel du
remplissage dans la nappe alluviale de Meskiana :
La zone de Dalâa est en réalité une dépression alluvionnaire, où l’épaisseur des alluvions peut
atteindre 155 m. Dans la zone de Zebar, qui se situe entre Dalâa et Meskiana, les dépôts du
remplissage ne dépassent pas 100 m d’épaisseur.
Il faut signaler que ces alluvions reposent sur un substratum marneux qui remonte en surface
au Nord de la plaine.
IV.3 PIEZOMETRIE :
- 38 -
925 930 935 940 945 950 955 960 965
290
285
Rehia
280
68
70
67 69
Bellala
65
64
62 63
275
61
MK1 Meskiana
60 59 58
A3BV5
57 56
54 F1 55
270 49 52 51 53
48 50
47 46
45 B8
43 44
41 B1042
265 39Zebar
CD7
40
38
36 37
MKN-1
35 34 E5 33
31 32
27 29 28 30
260
23 MK3
26 24
25 Légende
22
21 MK4
18 19 20
17
MK2
255 16 K5 13 14 15 MKN-1
Forage de reconnaissance pétrolière n° MKN-1
AL5 12
11
9 10 AL5 Forage hydraulique n° AL5
250 7 8
6 Dalâa 43 Puits domestique et agricole n° 43
5
M13 X14
4 3
Agglomération
245 D1
2
1
240 0 10 20 Km
235
Figure 27: Carte d’inventaire des points d’eau de la nappe de la plaine de Meskiana.
- 39 -
constitue également une limite imperméable, car les courbes piézométriques sont perpendiculaires
aux bordures.
930 940 950 960
290
280 Rehia
Bellala
270 Meskiana
260 Légende
Zebar
980m Courbe piézométrique
250 Sens d’écoulement
Agglomération
Dalâa
240
Echelle
0 10 20Km
- 40 -
930 940 950 960
290
Rehia
280
Bellala
270 Meskiana
260
Zebar Légende
0 10 20Km
290
280 Rehia
Bellala
270 Meskiana
260 Légende
Zebar
960m Courbe piézométrique
250 Sens d’écoulement
Agglomération
Dalâa
240
Echelle
0 10 20Km
- 41 -
930 940 950 960
290
Rehia
280
Bellala
270 Meskiana
260
Légende
Zebar
910m Courbe piézométrique
250 Sens d’écoulement
Dalâa
Agglomération
240
Echelle
0 10 20Km
Croûte calcaire Grès grossier Marne verdâtre Calcaire
Marne grise Marne et calcaire Gypse
Marne
- 42 -
930 940 950 960
290
Rehia
280
Bellala
270 Meskiana
260
Légende
Zebar
980m Courbe piézométrique
250 Sens d’écoulement
Dalâa
Agglomération
240
Echelle
0 10 20Km
Croûte calcaire Grès grossier Marne verdâtre Calcaire
Marne grise Marne et calcaire Gypse
Marne
Figure 32: Carte piézométrique de la nappe de Meskiana, mars 2004.
290
Rehia
280
Bellala
270 Meskiana
260
Légende
Zebar
910m Courbe piézométrique
240
Echelle
0 10 20Km
Croûte calcaire Grès grossier Marne verdâtre Calcaire
Marne grise Marne et calcaire Gypse
Marne
Figure 33: Carte piézométrique de la nappe de Meskiana, juillet 2004.
- 43 -
IV.4 CONDITIONS AUX LIMITES :
Dans le but de mieux comprendre le fonctionnement hydraulique de la nappe, on a établi
la carte des conditions aux limites de la nappe aquifère de Meskiana (figure 34). Cette carte a été
établie en se basant sur les caractéristiques lithologiques des formations géologiques, les
propriétés géophysiques de celles-ci et les relevés piézométriques.
La bordure orientale de Dj Guerigueur, formée essentiellement par des calcaires fracturés
d’âge Maestrichtien, constitue une limite perméable, car elle alimente latéralement la nappe
alluviale. La bordure occidentale de Dj Bou Tokhma et Dj Guern Ahmar, formée de calcaires
perméables d’âge Maestrichtien, constitue une limite à flux entrant. Vu leur nature lithologique,
les marnes cénomaniennes, affleurant au Nord-Est de la plaine, forment une limite étanche.
Les marnes et les argiles emscheriennes constituent le substratum imperméable,
empêchant ainsi les eaux d’infiltration de pénétrer en profondeur. Ces dernières s’écoulent
gravitairement vers l’aval. Au Nord, le substratum imperméable surgit en surface, ce qui rend la
quasi-totalité de cette limite imperméable, avec une zone centrale à flux imposé sortant.
L’Oued Meskiana draine la nappe alluviale en période des basses eaux, comme il peut l’alimenter
en période de hautes eaux, et correspond par conséquent à une limite à potentiel imposé. Les
pluies efficaces sont pratiquement nulles, comme on l’a vu lors de l’étude du bilan hydrologique.
290
285 Rehia
280
Bellala
275
Meskiana
270
265
Légende
260 Zebar
Limite à potentiel imposé
255
Limite à flux imposé sortant
Figure 34: Carte des conditions aux limites de la nappe de la plaine de Meskiana.
- 44 -
IV.5 Battement de la nappe:
La comparaison de la piézométrie entre la période de basses eaux et celle de hautes eaux a
permis de dresser une carte de battement de la nappe pour les années 2002, 2003 et 2004. (Figures
35 à 37).
Cette opération a montré que les fluctuations les plus importantes se manifestent en bordures de
nappe et peuvent atteindre 0.7 m. Au centre de la nappe, celles-ci varient entre 0 et 0.4 m
930 935 940 945 950 955 960
290
285
Rehia Bellala
280
275
265
Courbe de battement de la nappe (m)
0.4
260 Zebar
Limite du bassin
255
Agglomération
250
Dalâa
245
Echelle
240 0 10 20 Km
s=
0.183.Q 2.25T.t
log (13) où :
T R 2 .S
s: rabattements en (m) ;
T : transmissivité en (m2/s) ;
Q : débit de pompage en (m3/s) ;
S : coefficient d’emmagasinement ;
t : temps écoulé depuis le début du pompage en (secondes) ;
R : distance entre le piézomètre et le forage en (m).
Les graphes de cette interprétation, sont représentés en annexe III.
- 45 -
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
Meskiana
270
Légende
265
Courbe de battement de la nappe (m)
0.8
260 Zebar Limite du bassin
255 Agglomération
250
Dalâa
245
Echelle
240 0 10 20 Km
Figure 36: Carte de battement de la nappe de Meskiana en (m) durant l’année 2003.
930 935 940 945 950 955 960
290
285
275
270 Meskiana
265
260
Légende
Zebar
- 46 -
La détermination de la perméabilité a été faite en se basant sur la transmissivité
et l’épaisseur de la couche mouillée de l’aquifère. Cette dernière a été estimée en se
basant sur les profils géoélectriques et les colonnes lithostratigraphiques des forages.
L’application de la relation suivante a permis la détermination de la perméabilité, dont
les résultats sont consignés dans le tableau 14.
K=
T
(14) où
e
K : perméabilité (m/s) ;
T : transmissivité (m2/s) ;
e : épaisseur de la couche mouillée aquifère (m).
- 47 -
930 940 950 960
290
Rehia
Bellala
280
Meskiana
270
Légende
Zebar
260
290
Rehia
Bellala
280
Meskiana
270
Zebar Légende
260
- 48 -
IV.7 CONCLUSION :
- 49 -
Chapitre V :
HYDROCHIMIE
V.1 INTRODUCTION :
Il est évident que le chimisme des eaux souterraines dépend, principalement, de la
composition lithologique des couches traversées et du temps de séjour des eaux. Cette interaction
influe sur la teneur des éléments majeurs (Ca2+, Mg2+, Na+, K+, Cl-, SO42-, HCO3- …..). Les
concentrations de ces éléments naturels sont conditionnées par divers facteurs tels que les
paramètres climatiques, l’activité anthropique, les échanges entre aquifères et les eaux de surface
Les activités anthropiques ont une influence non négligeable sur la qualité physico-
chimique des eaux souterraines. On observe principalement deux types d'influence:
- En agriculture : apports de K+, Ca2+, NO3-, SO42-, NH4+ et PO43- .
- Rejets des eaux usées : apports de Na+, K+, NO3-, Cl- et quelques métaux lourds.
Le présent chapitre traite du chimisme des eaux souterraines, et nous permettra d’essayer
tentons d’expliquer le comportement, l'origine et l’évolution des éléments chimiques décrits pour
expliquer le fonctionnement hydraulique du système aquifère.
L'étude se base principalement sur les informations de six campagnes de prélèvement à
partir de 60 points d'eau, répartis sur le terrain. Comme la plaine de Meskiana est à vocation
agricole, un intérêt particulier sera accordé à cet effet.
V.2 ECHANTILLONNAGE ET ANALYSE HYDROCHIMIQUE:
Un réseau d’échantillonnage a été choisi pour permettre d'acquérir des données représentatives sur la variabilité
spatiale et temporelle de la qualité des eaux souterraines de la nappe étudiée. Ce réseau se compose de quatre
ensembles de points d’eau (figure 40) et couvre la totalité de la plaine de l’amont vers l’aval:
930 935 940 945 950 955 960
290
285 Rehia
Bellala
280
70
68 67 69
65
64
6263
275 61
Meskiana
60 59 58
5756
54 52 5553
270 49 51
48 50
47
45 46
43 4144 42
265 39 38 40
36 37
35 34 33
31 32 30 Légende
260 272928
26
222325 24Zebar
21 20 11
17 18 19 Puits échantillonné n°11
255 16 13 1415
11 12 Cours d’eau
9 10
250 7 8 Agglomération
5 6
Dalâa
4 3
245 Echelle
2
1
240 0 10 20 Km
- 50 -
-Puits proches des bordures calcaires ;
-Puits proches des cours d’eau (pour révéler les échanges entre les eaux de surface et les eaux
souterraines) ;
-Puits se répartissant dans les sédiments mio-plio-quaternaires ;
-Puits se répartissant dans les terrains cultivés ;
-Points de prélèvement à proximité de l’Oued Meskiana.
Les prélèvements ont été effectués selon une fréquence saisonnière, pendant six
campagnes où des analyses complètes ont été effectuées. Deux campagnes ont été réalisées
respectivement en décembre 2004 et novembre 2005 pour analyser les métaux lourds. Tableau
15.
Tableau 15: Les campagnes d’échantillonnage effectuées durant l’étude.
Campagne Période Nombre Eléments Laboratoire d’analyse
d’échantillons dosés
1 Mars 2002 43 Éléments CERAD,S.N Ferphos de Tébessa
majeurs, NO3-
2 Juillet 2002 39 Éléments CERAD,S.N Ferphos de Tébessa
majeurs,NO3-
3 Mars 2003 39 Éléments CERAD,S.N Ferphos de Tébessa.
majeurs,NO3-
4 Juillet 2003 39 Éléments CERAD,S.N Ferphos de Tébessa
majeurs,NO3-
5 Mars 2004 31 Éléments CERAD,S.N Ferphos de Tébessa
majeurs,NO3-
6 Juillet 2004 33 Éléments CERAD,S.N Ferphos de Tébessa
majeurs,NO3-,B,
Sr2+,NH4+
7 Décembre 2004 30 NO3-, NO2-, O2-, Laboratoire central,Univ de Damas (Syrie)
PO43-, fer, Sr2+
8 Novembre 2005 30 NO3-, PO43-, fer Laboratoire central,Univ de Damas (Syrie)
Cu2+) ont fait l’objet d’analyse pour quelques campagnes seulement. Les polluants résultants des
activités agricoles ont été l’objet d’un suivi particulier (NO3-, NO2-, NH4+ et PO43-).
Les méthodes analytiques utilisées sont :
-Méthode colorimétrique pour les éléments (Cl-, SO42-, HCO3-, NO3-, NO2-, NH4+ ,PO43- ,F et Fe
total ).
-Méthode de spectrophotométrie d’absorption atomique : (Ca2+, Mg2+, Na+, K+, Sr2+, Mn2+, Al3+,
Zn2+, Cd2+ et Cu2+). On signale que l’analyse spectrophotométrique a été faite en utilisant un
spectrophotomètre DR 2000 (HACH), aux laboratoires du Centre des Etudes et de Recherches
Appliquées et du Développement (CERAD), filiale de l’Entreprise Nationale du Fer et du
Phosphate.
- 51 -
Pour la fiabilité des résultats d’analyse, on a procédé à l’application de la méthode de la balance
ionique et une erreur de 10 % a été acceptée.
Eléments Norme Max Min Moy Ecart-type Max Min Moy Ecart-
OMS type
HCO3- (mg/l) 250 263 102 185.6 39.96 260 88 168.85 39.49
2-
SO4 (mg/l) 250 1098 75 463.6 278.12 1151 95 472.1 267.1
-
Cl (mg/l) 250 933 38 448 263.75 949 47 456.8 265.94
Ca2+ (mg/l) 100 805 76 267.2 157.6 845 70 280.42 154.4
Mg2+ (mg/l) 50 73 10 34.15 18.45 61 8 26.93 14.28
Na+ (mg/l) 150 425 23 235.1 121.1 506 32 253.24 129.73
K+ (mg/l) 12 18 06 12.35 3.47 21 7 13.44 2.95
NO3- (mg/l) 50 91 20 56.9 20.8 125 24 60.97 21.42
pH 6.5<pH<9.5 8.7 7.4 8.3 0.38 8.7 7.3 8.34 0.34
-1
Cond (µS.cm ) 2500 5050 380 2571.4 1171.1 5230 460 2616.5 1173.9
Les valeurs de la conductivité électrique montrent que la plupart des puits dépassent la
norme recommandée relative à la potabilité des eaux, surtout pendant la période des basses eaux.
Les valeurs minimales sont enregistrées à l’amont de la nappe.
Les valeurs maximales sont de l’ordre de 5230 (µS.cm-1), on les rencontre dans la zone
d’accumulation des eaux souterraines, au Nord-Est de la nappe, elles sont dues aux dépôts
évaporitiques issus de Djebel Belkefif.
- 52 -
V.3.2 Les éléments majeurs :
- 53 -
V.3.2.5 Les chlorures (Cl-) :
L’ion chlorure possède des caractéristiques différentes de celles des autres éléments, il
n’est pas adsorbé par les formations géologiques, ne se combine pas facilement avec les éléments
chimiques et reste très mobile. Il constitue un bon indicateur de la pollution. Les teneurs
enregistrées, pendant toute la période d’observation, montrent des valeurs importantes qui
dépassent la norme de potabilité de l’OMS fixée à 250 mg/l. On signale que 74 % des puits
dépassent cette norme. Les chlorures peuvent avoir plusieurs origines (écart type 265.94: mg/l) :
-Les argiles sableuses gypsifères du Mio-Plio-Quaternaire qui couvrent la totalité du bassin ;
-Les dépôts évaporitiques provenant de Dj Belkefif ;
-Le déversement des eaux usées dans la plaine ;
-Les marnes emschériennes formant le substratum.
- 54 -
Tableau 17:Faciès chimiques des eaux souterraines de l’aquifère de Meskiana.
Faciès chimiques
- 55 -
V.4.1.1.1 Période de mars 2002 (fig. 41a et fig. 41b) :
A l’amont de la nappe aquifère, l’ensemble des nuages de points ne présente pas de
dominance anionique ou cationique. On signale cependant la présence d’un faciès chloruré
calcique, relativement dominant, et un faible pourcentage (15%) de puits ayant acquis un faciès
bicarbonaté calcique dû à l’alimentation par les bordures. La présence de quelques cations
sodiques s’explique par les marnes emschériennes gypsifères qui constituent le substratum de la
nappe.
On note que les ions ne connaissent pas encore une certaine stabilité à cet endroit, vu les
conditions hydrodynamiques de la nappe, et ont une tendance à se concentrer.
Les puits qui se situent à l’aval présentent un état différent, avec des anions chlorurés et
sulfatés dominants, alors que les cations calciques marquent la totalité des échantillons et sont
secondés par les cations sodiques. Ceci révèle la dominance de trois faciès hydrochimiques :
chloruré-calcique, chloruré-sodique et sulfaté-calcique, avec l’absence du faciès bicarbonaté
calcique.
Le diagramme global confirme la présence des faciès signalés. Les conditions climatiques
(les précipitations), la lithologie et les conditions hydrodynamiques (vitesse d’écoulement) ont
favorisé le lessivage des couches traversées, ce qui a engendré ces faciès.
- 56 -
Figure 41a: Diagramme de Piper, campagne mars 2002 (amont de la nappe).
- 57 -
Figure 42a : Diagramme de Piper, campagne juillet 2002 (amont de la nappe).
- 58 -
Figure 43a: Diagramme de Piper, campagne mars 2003 (amont de la nappe).
- 59 -
A l’aval de la nappe, les puits qui se trouvent à proximité des bordures calcaires (P3, P8
et P15) ont perdu leur faciès bicarbonaté calcique, les nuages de points se concentrent autour de
deux pôles : le premier calcique et le second chloruré à sulfaté.
La combinaison de ces ions a donné naissance aux faciès chloruré calcique et sulfaté
calcique. Comme on l’a signalé auparavant, la lithologie, les conditions hydrodynamiques de la
nappe et l’évaporation sont à l’origine de ces faciès.
- 60 -
Figure 44a:Diagramme de Piper, campagne juillet 2003 (amont de la nappe).
- 61 -
Figure 45a:Diagramme de Piper, campagne mars 2004(amont de la nappe).
- 62 -
Figure 46a:Diagramme de Piper, campagne de juillet 2004 (amont de la nappe).
- 63 -
V.4.1.2 Diagramme de SCHOELLER-BERKALOFF:
Le diagramme de Schoëller- Berkaloff permet de représenter le faciès chimique de
plusieurs échantillons d’eaux. Chaque échantillon est représenté par une ligne brisée. La
concentration de chaque élément chimique est figurée par une ligne verticale en échelle
logarithmique. La ligne brisée est formée en reliant tous les points qui représentent les différents
éléments chimiques.
Un groupe d'eau de minéralisation variable mais dont les proportions sont les mêmes pour
les éléments dissous, donnera une famille de lignes brisées parallèles entre elles. Lorsque les
lignes se croisent, un changement de faciès chimique est mis en évidence.
- 64 -
Lieu
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p4
10
p6
10 p8
10
10 p10
10 p12
10 10
p13
0,1 0,1
1 p15
1 p16
1
Amont de la nappe p17
1
1 p18
1 1
p20 Lieu
0,01 0,01
0,1 p21
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L 0,1
mg/L mg/L mg/L
0,1 mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000
0,1
300 Berkaloff
10000
0,1
0,02 10000
0,1 10000
0,1
0,03
0,03
100 0,04 0,05 0,07 0,07 100
0,001 1000 0,001
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p34
100
100 100 p35
1 1 p38
10
p41
10 p42
10
10 p43
10 p44
10 10
p45
0,1 0,1
1 p46
1 p47
1
p48
1
1 p49
1 1
p50
0,01 0,01
0,1 p51
Aval
0,1 de la nappe0,1
0,1
0,1
0,03 0,02 0,1 0,1
0,03
0,04 0,05 0,07 0,07
0,001 0,001
Figure 47 : Diagramme de Schoëller- Berkaloff, mars 2002.
- 65 -
Lieu
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p4
10
p6
10 p8
10
10 p10
10 p12
10 10
p13
0,1 0,1
1 p15
1 Amont de la nappe1 p16
p17
1
1 p18
1 1
p20
0,01 0,01
0,1 p21
0,1
0,1 Lieu
Ca Mg Na+K Cl
0,1 SO4 HCO3+CO3 NO3
meq/L mg/L mg/L
0,02 mg/L mg/L 0,1
mg/L mg/L
0,1 mg/L
0,1 meq/L Schöeller
0,03
300
0,001
6000
3000
0,03
6000 10000
0,04 0,05 0,07 0,07 300
0,001
Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p34
100
100 100 p35
1 1 p38
10
p41
10 p42
10
10 p43
10 p44
10 10
p45
0,1 0,1
1 p46
1 p47
1
p48
1
1 p49
1 1
p50
0,01 Aval de la 0,01
0,1nappe p52
0,1
0,1
0,1
0,1
0,03 0,02 0,1 0,1
0,03
0,04 0,05 0,07 0,07
0,001 0,001
Figure 48: Diagramme de Schoëller-Berkaloff, juillet 2002.
- 66 -
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p4
10
p6
10 p8
10
10 p10
10 p12
10 10
p13
0,1 0,1
1 p15
Amont de la nappe
1 p16
1
p17
1
1 p18
1 1
p22
0,01 0,01
0,1 p23
0,1
Lieu
0,1
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
0,1 Schöeller
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L
0,1
300 6000
0,03 0,02
3000 6000 10000 0,1 0,1 300 Berkaloff
0,03 10000
0,04 0,05 0,07 0,07
0,001 10000 10000 0,001
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p38
100
100 100 p41
1 1 p42
10
p43
10 p45
10
10 p46
10 p47
10 10
p48
0,1 0,1
1 p49
1 p50
1
p51
1
1 p52
1 1
p54
0,01 0,01
0,1 p58
0,1
0,1
0,1
Aval de la nappe 0,1
0,03 0,02 0,1 0,1
0,03
0,04 0,05 0,07 0,07
0,001 0,001
Figure 49: Diagramme de Schoëller- Berkaloff, mars 2003.
- 67 -
V.4.1.2.4 Période de juillet 2003 (fig.50) :
La plupart des puits ont acquis le faciès chloruré calcique, tandis que les quatre puits (P1,
P3, P8 et P15) gardent toujours le faciès bicarbonaté calcique. Les concentrations sont
relativement fortes par rapport aux campagnes précédentes.
Au centre et au Nord de la nappe, les teneurs deviennent plus importantes. Le puits P52,
qui se situe au Nord de la ville de Meskiana, à proximité de l’oued, est fortement minéralisé et
atteint les concentrations maximales. On pense que l’évaporation a causé l’augmentation de ces
teneurs.
- 68 -
Lieu
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p4
10
p6
10 p8
10
10 p10
10 p12
10 10
p13
0,1 0,1
1 p15
1 p16
1
Amont de la nappe 1
p17
1 p18
1 1
p22
0,01 0,01
0,1 p23
0,1
0,1 Lieu
Ca Mg Na+K Cl0,1 SO4 HCO3+CO3 NO3
meq/L mg/L mg/L
0,02 mg/L mg/L 0,1
mg/L mg/L
0,1 mg/L
0,1 meq/L Schöeller
0,03
300
0,001
6000
3000
0,03
6000 10000
0,04 0,05 0,07 0,07 300
0,001
Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p41
100
100 100 p42
1 1 p43
10
p44
10 p46
10
10 p47
10 p48
10 10
p49
0,1 0,1
1 p50
1 p51
1
p52
1
1 p54
1 1
p58
0,01 0,01
0,1 p60
Aval de la nappe
0,1
0,1
0,1
0,1
0,03 0,02 0,1 0,1
0,03
0,04 0,07 0,07
Figure
0,001 50:Diagramme de Schoëller-Berkaloff
0,05 , juillet 2003.0,001
- 69 -
Lieu
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p6
10
p8
10 p10
10
10 p13
10 p15
10 10
p16
0,1 0,1
1 p17
1 p18
1
p22
1
Amont de la nappe 1 p23
1 1
p29
0,01 0,01
0,1 p31
0,1
0,1
0,1
0,1
0,03 0,02 0,1 0,1 Lieu
0,03
Ca Mg Na+K Cl 0,04 SO4 0,05 0,07
HCO3+CO3 NO3
0,07
0,001 0,001
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p41
100
100 100 p42
1 1 p46
10
p48
10 p49
10
10 p50
10 p51
10 10
p52
0,1 0,1
1 p54
1 p58
1
p60
1
Aval de la nappe 1 p64
1 1
p65
0,01 0,01
0,1 p68
0,1
0,1
0,1
Figure 51: Diagramme de Schoëller-Berkaloff,
0,1
mars 2004.
0,03 0,02 0,1 0,1
0,03
0,04
-0,07
70 - 0,07
0,001 0,05 0,001
Lieu
Ca Mg Na+K Cl SO4 HCO3+CO3 NO3
meq/L mg/L mg/L mg/L mg/L mg/L mg/L mg/L meq/L Schöeller
300 6000
3000 6000 10000 300 Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p1
100
100 100 p3
1 1 p6
10
p8
10 p10
10
10 p12
10 p13
10 10
Amont de la nappe p15
0,1 0,1
1 p16
1 p17
1
p22
1
1 p23
1 1
p27
0,01 0,01
0,1 p31
0,1
0,1 Lieu
Ca Mg Na+K Cl0,1 SO4 HCO3+CO3 NO3
meq/L mg/L mg/L
0,02 mg/L mg/L 0,1
mg/L mg/L
0,1 mg/L
0,1 meq/L Schöeller
0,03
300
0,001
6000
3000
0,03
6000 10000
0,04 0,05 0,07 0,07 300
0,001
Berkaloff
10000
10000 10000
100 100
1000
1000
1000
1000
1000
1000 1000
10 10
100
100
100
100 p41
100
100 100 p42
1 1 p44
10
p46
10 p47
10
10 p48
10 p49
10 10
p50
0,1 0,1
1 p52
1 p54
1
p58
1
1 p60
1 1
Aval de la nappe p62
0,01 0,01
0,1 p64
0,1
0,1
0,1
0,1
0,03 0,02 0,1 0,1
0,03
0,04 0,05 0,07 0,07
0,001 0,001
Figure 52:Diagramme de Schoëller-Berkaloff, juillet 2004.
- 71 -
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
P68 P69
P65
275 P64
P62
P60 P58
270
P49
P54P52 Meskiana
P51
P48 P50
P43
P47
P45 P46
P44 P42
Légende
265 P41
P38
P35P34 Faciès bicarbonaté calcique
260 P31 P32 P29
P27
P23
P22
P21 Zebar
P18 P20 Faciès chloruré calcique
255 P17
P16 P13 P15
P12 Faciès chloruré sodique
P10
250 P8
P6 Dalâa Faciès sulfaté calcique
245 P4 P3
P51
Puits échantillonné n° 51
240 P1 Agglomération
0 5 10 Km
- 72 -
La cartographie hydrochimique a été illustrée par interpolation triangulaire, en utilisant le
logiciel Surfer 8 (Golden software, 2002).
- 73 -
930 935 940 945 950 955 960
290
285
Rehia
280 Bellala
275
270 Meskiana
265
255
350 Courbe d’isobicarbonate en (mg/l)
250 Agglomération
Dalâa
245
Echelle
240 0 10 20 Km
290
285
Rehia
280
Bellala
275
270 Meskiana
265
260 Légende
255 Zebar
350 Courbe d’isobicarbonate en (mg/l)
250 Agglomération
Dalâa
245
Echelle
240
0 10 20 Km
290
285
Rehia
280
Bellala
275
270 Meskiana
265
260 Légende
Zebar
380 Courbe isochlorure en (mg/l)
255
Agglomération
250
245 Dalâa
240 Echelle
0 10 20 Km
290
285
Rehia
280
Bellala
275
270 Meskiana
265
260 Légende
Zebar
255 390
Courbe isochlorure en (mg/l)
250 Agglomération
Dalâa
245 Echelle
240 0 10 20 Km
290
285
Rehia
280 Bellala
275
Meskiana
270
265
260
Zebar
255
Légende
300 Courbe isosulfate en (mg/l)
250
Dalâa Agglomération
245
240 Echelle
0 10 20 Km
290
285
Rehia
280
Bellala
275
Meskiana
270
265
260 Zebar
Légende
255
700 Courbe isosulfate en (mg/l)
250
Agglomération
Dalâa
245
240 Echelle
0 10 20 Km
- 76 -
V.5.4 Cartes du calcium (fig.60 et fig.61) :
Les concentrations oscillent en général entre 75 et 850 mg/l. Ces dernières sont
relativement faibles à l’amont de la nappe par rapport aux autres régions, car les eaux souterraines
ne sont pas encore minéralisées et proviennent de la dissolution des formations carbonatées de
bordures. On remarque que les teneurs croissent dans le sens de l’écoulement. Les fortes teneurs
apparaissent au Nord-Est, on pense qu’elles résultent de la dissolution des gypses contenus dans
les marnes emschériennes et les dépôts évaporitiques provenant de Djebel Belkefif.
La variation des teneurs en calcium est semblable à celle des chlorures, ce qui confirme la
prédominance du faciès chloruré-calcique.
Les ions Ca2+ proviennent de la dissolution des formations carbonatées (CaCO3) et des gypses
contenus dans les marnes emschériennes (CaSO4) selon les équations d’équilibre suivantes :
CaCO3 = Ca2+ + CO32-
CaSO4, 2H2O = Ca2+ + SO42- + 2H2O
On remarque que les fortes concentrations coïncident avec la période des basses eaux, ce
qui nous laisse penser que l’élévation des concentrations est due au phénomène d’évaporation. La
dissolution des calcaires ne contribue qu’avec un faible degré.
L’analyse des données chimiques à l’aide de la cartographie a pu donner une idée sur le
mécanisme du chimisme des eaux souterraines de la nappe aquifère.
L’influence des formations géologiques sur la nature des faciès est très remarquable : les calcaires
ont donné le faciès bicarbonaté calcique, confirmant ainsi l’alimentation latérale de la nappe, à
partir des bordures calcaires.
Les marnes emschériennes, qui sont gypsifères, ainsi que les dépôts triasiques, sont à
l’origine des faciès chloruré calcique, chloruré-sodique et sulfaté-calcique, par les phénomènes de
lessivage et d’évaporation.
Les conditions hydrodynamiques de la nappe ont favorisé l’acquisition de ces faciès :
-La surface piézométrique n’est pas assez profonde ce qui augmente l’effet de l’évaporation dans
la zone septentrionale, où le gradient hydraulique diminue, donnant ainsi suffisamment de temps à
l’interaction entre les eaux souterraines et la lithologie ;
-Pendant la période des hautes eaux, les précipitations provoquent la dissolution des
calcaires de bordures et le lessivage des dépôts de remplissage.
- 77 -
930 935 940 945 950 955 960
290
285
275
270 Meskiana
265
Légende
260 Zebar
550
Courbe d’isocalcium en (mg/l)
255
Agglomération
250
245 Dalâa
Echelle
240 0 10 20 Km
290
285
280 Rehia
Bellala
275
270 Meskiana
265
260 Légende
Zebar
255
550 Courbe d’isocalcium en (mg/l)
250
Agglomération
245 Dalâa
Echelle
240
0 10 20 Km
- 78 -
930 935 940 945 950 955 960
290
285
Rehia
280 Bellala
275
Meskiana
270
265
Légende
260 Zebar
245
Echelle
240 0 10 20 Km
290
285
280 Rehia
Bellala
275
Meskiana
270
265
Zebar
260
Légende
255
550 Courbe d’ isosodium en mg/l
250
Dalâa Agglomération
245
Echelle
240 0 10 20 Km
- 79 -
- Les eaux souterraines s’écoulent suivant la direction SW-NE et se chargent progressivement au
fur et à mesure que l’on se déplace vers la zone d’accumulation, ce qui explique les fortes teneurs
observées dans cette région.
Quelques éléments chimiques proviennent de plusieurs origines différentes tels que le
calcium et le sodium, celle qui influe sur la prédominance des faciès doit être mise en évidence.
Ces faciès se sont-ils stabilisés ou évoluent-ils vers d’autres tendances ?
Ces éléments chimiques ont-ils la même évolution dans le temps et dans l’espace ?
Pour répondre à ces questions, une étude détaillée s’impose pour mettre en relief l’origine,
l’évolution et l’influence de chaque élément sur l’acquisition des faciès. La réalisation d’une telle
étude se basera sur les propriétés de chaque élément chimique et les outils statistiques.
V.6.1 Faciès caractéristiques :
A partir de la représentation graphique des faciès chimiques par les diagrammes de Piper
et celui de Schoelleur-Berkaloff, on a constaté la dominance de 04 faciès qui sont : le chloruré-
calcique, le chloruré-sodique, le sulfaté-calcique et le bicarbonaté-calcique. Ceci est vérifié par les
teneurs dominantes des ions suivants : Cl- (26.73 méq/l), SO42- (23.97 méq/l), HCO3- (4.31 méq/l),
Ca2+ (42.25 méq/l) et Na+ (22 méq/l).
10
Chloruré sodique
2+
Na /Ca
1
+
10
Chloruré sodique
2+
Na /Ca
1
+
La réalisation d’une représentation graphique des faciès, qui se base sur les quantités en réaction, est
Sulfaté calcique Chloruré calcique
0,1
0,1 1 10 100
2-
Cl-/SO4
- 80 -
très utile pour la confirmation de la présence de ces faciès caractérisant les eaux
souterraines de la nappe. Etant donné la dominance des ions précédents, on a réalisé deux
diagrammes (figure 64, A et B) qui utilisent les rapports Na+/ Ca2+– Cl-/ HCO3- et Na+/ Ca2+ -
Cl-/ SO42-. L’interprétation de ces deux diagrammes a montré la concentration des rapports
d’ions en 04 pôles indiquant les 04 faciès, déterminés précédemment, qui caractérisent la
nappe. On note l’absence des faciès : bicarbonaté sodique et sulfaté sodique.
V.6.2 Origine des éléments chimiques :
Les eaux souterraines qui circulent dans la nappe alluviale de Meskiana proviennent
des formations carbonatées de l’Eocène et du Maestrichtien qui bordent la plaine. Ensuite,
elles s’écoulent dans un matériel très hétérogène (des alluvions, des sables, des graviers, des
dépôts évaporitiques,,,etc.). Le parcours d’écoulement et les couches traversées influent sur
l’acquisition des éléments chimiques. En se basant sur les mesures et les analyses chimiques
de toutes les campagnes, on va essayer de déterminer l’origine probable de tous les ions.
100
Excès en Ca2+
A
Ca 2+ (méq/l )
10
1
1 10
HC O 3 - (mé q/l)
50
B
40 Origine gypsifère
Ca en méq/l
30
20
2+
Afin10
de différencier le calcium d’origine carbonatée duOrigine
calciumcarbonatée
d’origine évaporitique, nous avons
0
0 5 10 15 20 25 30
- 2-
(HCO3 + SO4 ) en méq/l
- 81 -
établi un graphe où nous avons représenté le calcium en fonction des bicarbonates et des
sulfates (figure 65B). L’analyse de ce graphe montre que 35 échantillons sur 225 ont une origine
carbonatée soit 15% de l’ensemble des échantillons.
Na+/cl
30
A
Excès en Na+
Na en méq/l
20
+
10
0
0 10 20 30
-
Cl en méq/l
2
Ca / (HCO 3 + SO 4 )
Excès en Ca2+
2-
Etat naturel
1,5
B
-
0
0 0,5 1 1,5 2
+ -
Na /Cl
- 82 -
V.6.2.3 Les éléments HCO3-, Cl- et SO42- :
Après avoir fait le point sur les origines des cations, on a essayé de représenter l’évolution
des anions en fonction de la conductivité. L’interprétation de la figure 67A montre que la
conductivité électrique n’évolue pas avec les bicarbonates, alors que les chlorures et les sulfates
croissent avec la salinité. Cette dernière est donc due aux éléments salifères.
30
Anions en (méq/l)
25
20
15
10
5
0
0 1000 2000 3000 4000 5000 6000
Conductivité électrique en (uS/cm)
B
Figure 67a : Evolution des anions en fonction de la conductivité électrique.
2,5
Pôle carbonaté
HCO 3 -/(Cl - + SO 4 2-)
1,5
0
0 1000 2000 3000 4000 5000
Conductivité éle ctrique e n (uS/cm)
- 83 -
V.6.2.3 Les éléments Cl--SO42- :
Ces deux ions sont issus de deux formations différentes, le SO42- provient du gypse et le Cl-
provient de la dissolution de la halite et probablement d’une pollution anthropique. On a associé
l’évolution de ces ions à la conductivité électrique pour mettre en relief la façon selon laquelle ils
évoluent et déterminer lequel entre eux est à l’origine de la salinité des eaux de la nappe.
L’interprétation du graphe (figure.68A) a révélé que la majorité des points d’eau présente un
rapport (Cl-/SO42-) supérieur à 1, indiquant ainsi la dominance des ions Cl- par rapport à SO42-, ce
qui confirme la dominance du faciès chloruré calcique et sodique. Quelques points d’eau
présentent des valeurs du rapport (Cl-/SO42-) inférieures à 1, indiquant une origine gypsifère
(figure 68B).
3,5
2,5
Cl -/SO 4 2-
2
A
1,5
1
B
0,5
0
0 1000 2000 3000 4000 5000 6000
Conductivité é le ctrique e n (uS/cm)
- 84 -
-les minéraux carbonatés présentent des degrés de saturation différents. En effet, celui de la
calcite varie de - 0.36 à + 0.19, suivi par celui de l’aragonite qui varie de - 0.5 à + 0.12 et enfin
celui de la dolomite qui varie de -1.56 à - 0.29.
-Si l’on admet que l’état d’équilibre se fait dans l’intervalle de - 0.5 à + 0.5, on peut dire que les
trois minéraux ont atteint leur état d’équilibre et non pas celui de la sursaturation.
La répartition dans l’espace de ces indices montre une augmentation des concentrations
dans deux sens, pour la calcite et l’aragonite : la première à partir des bordures vers l’axe de la
nappe et la seconde dans le sens de l’écoulement du Sud-Ouest vers le Nord-Est (figure 69 et
figure 70). Par contre, la concentration de la dolomite montre une évolution inverse à celles des
deux minéraux précédents, ce qui confirme l’origine gypsifère et évaporitique du calcium et non
pas celle carbonatée.
- Les minéraux évoporitiques montrent des degrés de saturation plus faibles que les minéraux
carbonatés (figure 69 et figure 70). Le gypse enregistre des indices qui varient de -1.63 à -0.08
(où 70 % des points d’eau analysés n’ont pas atteint le degré de saturation), suivi de l’anhydrite
dont l’indice varie de - 1.85 jusqu’à - 0.3 (83% des puits sont sous-saturés) et enfin l’halite dont
les indices varient de -7.64 à - 4.10 (donc très sous saturés)
P10
P12
P20
P21
P22
P23
P27
P29
P32
P34
P35
P48
P49
P50
P52
P54
P60
P64
P65
P1
P3
P4
P8
1
0 calcite
-1 dolomite
-2 aragonite
anhydrite
-3
gypse
-4
halite
-5
-6
-7
-8
-9
Figure 69: Variation de l’indice de saturation des minéraux, mars 2002.
P10
P12
P20
P21
P22
P23
P27
P29
P32
P34
P35
P48
P49
P50
P52
P54
P60
P64
P65
P1
P3
P4
P8
1
0 calcite
-1 dolomite
-2 aragonite
anhydrite
-3
gypse
-4
halite
-5
-6
-7
-8
-9
Figure 70: Variation de l’indice de saturation des minéraux, juillet 2002.
- 85 -
Tableau 18 : Evolution des indices de saturation des minéraux des campagnes de l’année 2002.
P 29 CS - 0.03 - 0.58 - 0.18 - 0 1.0 -0.78 - 5.24 CS - 0.19 - 1.22 - 0.32 - 0.80 - 0.58 -5.15
P 31 CC - 0.35 - 1.13 - 0.49 - 1.65 - 1.87 - 6.88 CC - 0.33 - 1.27 - 0.47 - 1.71 - 1.49 - 6.66
P 32 CC - 0.04 - 0.15 - 0.10 - 0.90 - 0.70 - 5.64 CC + 0.01 - 0.70 - 0.13 - 0.84 - 0.62 - 5.51
P 34 CS - 0.08 - 0.37 - 0.23 - 0.93 - 0.71 - 5.09 CS - 0.21 - 1.29 - 0.35 - 0.76 - 0.54 - 5.02
P 35 CC + 0.08 - 0.14 - 0.10 - 0.90 - 0.68 - 5.61 CC + 0.01 - 0.55 - 0.13 - 0.81 - 0.59 - 5.48
P38 CS - 0.12 - 0.71 - 0.26 - 1.49 - 1.71 - 6.66 CS - 0.13 - 0.70 - 0.26 - 0.83 - 0.62 - 5.67
P 48 CC + 0.04 - 0.31 - 0.11 - 0.86 -0.64 - 5.57 CC - 0.03 - 0.62 - 0.17 - 0.80 - 0.58 - 5.44
P 49 CC + 0.06 - 0.30 - 0.08 - 0.97 - 0.75 - 5.74 CC + 0.05 - 1.04 - 0.09 - 0.88 - 0.66 - 5.59
P 50 SC + 0.07 - 0.29 - 0.07 - 0.64 - 0.42 - 5.45 SC + 0.01 - 0.62 - 0.16 - 0.61 - 0.39 - 5.39
P 51 SC - 0.07 - 0.41 - 0.07 - 0.51 - 0.29 - 5.31 SC - - - - - -
Meskiana
P 52 SC + 0.04 - 0.50 - 0.10 - 0.56 - 0.34 - 5.36 SC - 0.01 - 0.76 - 0.15 - 0.53 - 0.31 - 5.33
P 54 CC + 0.02 - 0.65 - 0.70 - 0.80 - 0.51 - 5.52 CC - 0.05 - 0.94 - 0.21 - 0.74 - 0.52 - 5.39
P 58 SC + 0.022 - 1.02 + 0.12 - 0.34 - 0.12 - 5.29 SC - - - - - -
P 60 CC + 0.026 - 0.93 - 0.23 - 0.74 - 0.52 - 5.41 CC - 0.06 - 1.05 - 0.22 - 0.78 - 0.48 - 5.59
P 64 CC + 0.03 -1.17 - 0.22 - 0.60 - 0.37 - 5.25 CC - 0.08 - 1.40 - 0.28 - 0.63 - 0.41 - 5.22
P 65 SC + 0.032 - 1.11 + 0.08 - 0.32 - 0.32 - 5.15 SC + 0.19 - 1.42 + 0.05 - 0.30 - 0.08 - 5.19
P 69 - - - - - - - CC - 0.15 - 1.56 + 0.10 - 0.60 - 0.37 - 5.17
86
Les concentrations les plus importantes ont été enregistrées dans les régions Nord-Est de
Meskiana où affleurent les dépôts évaporitiques et celle de Zebar où se trouvent les marnes
emscheriennes gypsifères.
Le calcul de l’indice de saturation des différents minéraux dans l’eau indique que seuls
les minéraux carbonatés ont tendance à atteindre l’équilibre, surtout sous la forme de calcite et en
particulier pendant la période des hautes eaux. Par contre, les minéraux évaporitiques sont
toujours à l’état de sous-saturation, malgré les fortes concentrations qu’ils acquièrent.
V.7.1 Introduction :
La représentation graphique des éléments chimiques et l’étude des rapports
caractéristiques ont montré qu’un grand nombre de paramètres chimiques et physico-chimiques
étaient significatifs, pris individuellement. Pour expliquer l’évolution du chimisme, la plupart du
temps une (ou plusieurs) liaison plus ou moins forte existe entre ces paramètres. Il est donc apparu
intéressant d’essayer de traiter les données statistiquement, en utilisant une méthode d’analyse
multidimensionnelle comme l’Analyse en Composantes Principales (ACP).
V.7.2 Principe:
Chaque échantillon d’eau d’une campagne de prélèvement constitue une unité statistique.
Les différents paramètres soumis à l’analyse constituent les variables qui caractérisent ces unités
statistiques. Chacune des six campagnes analysées au cours des trois années (2002,2003 et 2004)
a donc été traitée statistiquement grâce à l’utilisation du logiciel SPAD (voir annexe IV).
Notre présente étude porte sur les éléments majeurs, qui influent sur l’évolution de la
minéralisation. On a inclus la conductivité électrique, dont l’évolution dépend de cette dernière.
Cl- 1,00
SO42- 0,75 1,00
HCO3- -0,55 -0,52 1,00
NO3- 0,62 0,48 -0,63 1,00
Ca2+ 0,63 0,92 -0,57 0,51 1,00
Mg2+ 0,37 0,30 0,25 -0,10 -0,01 1,00
Na++K+ 0,95 0,77 -0,61 0,68 0,65 0,29 1,00
Cond. 0,88 0,90 -0,55 0,57 0,82 0,29 0,88 1,00
87
D’autres corrélations moins significatives entre la conductivité électrique d’une part et les
bicarbonates permettent de dire que la salinité de l’eau est beaucoup plus due aux
formations salifères telles que le gypse, l’anhydrite, le sel gemme, la sylvinite, …
La forte liaison qui existe entre les sulfates et le calcium confirme cette tendance.
Les nitrates sont bien corrélés avec le sodium suite à l’utilisation abusive des engrais type
N P K.
88
Légende
5P49 5 : n° de campagne
P49 : n° du puits
V.8.1 Introduction :
Dans les zones arides et semi-arides, l’approvisionnement en eau d’irrigation constitue
l’un des facteurs déterminants dans la production agricole, aussi bien dans l’intensification
des cultures, que dans l’extension des surfaces irriguées. Pour les régions tempérées, les eaux
superficielles constituent la principale source d’eau d’irrigation ; alors que dans les zones
semi-arides, où cette ressource est rare ou inexistante, on fait appel aux eaux souterraines. Le
développement de l'agriculture dans ces zones rencontre actuellement, en dehors de la rareté
des ressources hydriques, de nouveaux problèmes tels que le risque de salinisation des sols
qui peut être apprécié par la conductivité électrique et celui de l'alcalinisation des sols. Cette
dernière est due aux échanges ioniques (du sodium, du calcium et du magnésium) entre l'eau
et le sol irrigué.
L’accumulation des sels hydrosolubles dans le sol d’enracinement influe négativement
sur la croissance des plantes, on a recensé deux séries d'effets de la salinité (M. De FORGES,
1972) : les uns portent sur les végétaux, les autres sur le sol. En effet, ces sels causent des
changements de la perméabilité et de l’aération du sol d’une part, et d’autre part ils
provoquent des perturbations du métabolisme des plantes et du processus osmotiques de ces
dernières.
Comme la région de Meskiana est à vocation agricole et les risques de salinisation des
sols se font sentir, une étude d’aptitude des eaux à l’irrigation s’avère nécessaire car cette
région se caractérise par une faible pluviométrie, une forte évaporation et une eau
souterraine trop chargée en chlorures et en sulfates. Les risques qu’encourt cette plaine
89
imposent une bonne connaissance des régimes hydriques et salins afin de pouvoir assurer une
durabilité d’utilisation des ressources en eau et des sols.
Dans la présente approche, on a fait l’estimation de quelques paramètres de salinité en
rapport avec l’usage agricole pour les six campagnes saisonnières.
( Na + + K + )
SAR =
(Ca 2+ + Mg 2+ ) / 2
(16)
-La classe C2S1 : Elle caractérise une eau bonne pour l’irrigation et pouvant être utilisée sans
contrôle particulier pour l’irrigation des plantes moyennement tolérantes aux sels. Cette classe
représente les eaux de faible minéralisation que l’on rencontre en amont de la nappe à
proximité des bordures. Cette classe n’a paru que dans les deux campagnes de l’année 2002 et
avec un faible pourcentage (7.2 %) (tableau .20) ;
90
4
30
4
C1S4 C2S4 C3S4 C4S4 30 C1S4 C2S4 C3S4 C4S4
Très élevé
Très élevé
3
3
Sodium Absorption Ratio (SAR)
Risque de sodicité
Risque de Sodicité
20 C2S3 20 C2S3
Moyen
Moyen
2
2
C3S3 C3S3
C1S2 C2S2 C1S2
C4S3 C4S3
C3S2 C2S2
10 C3S2
10
C4S2 C4S2
8 C1S1 8 . P34
C2S1 P18 .P29 . P23. C1S1 C2S1 ..P23 P29
Faible
Faible
6 P17 . . . 6 .. .
1
1
C3S1 P41 P38. C3S1 . P21
4 P46. C4S1 P60
P47 . . 4 . P13 C4SP32
. . 1 P52
2 P10. .
P13 P35 P27 P65 . 2 .. . .
P10 P12 P27 P35
.P1..P8P3 P31. . . . .
P12 P4 .
P6 P16 P32 P45 . .P1 .P8.P3 . . . .
P15P20 P4 P6 P16. P65
P15
2 3 5 6 7 8 9 1000 2 3 4 2 3 4 5 6 7 8 9 1000 2 3 4
Conductivité Spécifique en µS/cm (CE*10) à 25°c Conductivité Spécifique en µS/cm (CE*10) à 25°c
1 2 3 4 1 2 3 4
Faible Moyen Elevé T. Elevé Faible Moyen Elevé T. élevé
Risque de Salinité Risque de Salinité
Mars 2002
Juillet 2002
4
30
4
Très élevé
3
3
Sodium Absorption Ratio (SAR)
C1S3 C1S3
20 C2S3 20 C2S3
Risque de Sodicité
Moyen
Moyen
2
C3S3 C3S3
C1S2 C1S2
C4S3 C4S3
C2S2 C3S2 C2S2
10 10 C3S2
8 C1S1
C4S2
. .P23
P22 8
.C S
4 2
C2S1 P18. . P29 ..
C1S1 P22 P23
Faible
Faible
91
-Les classes C3S1, C3S2: Généralement, ces deux classes renferment des eaux admissibles pour
l’irrigation des cultures tolérantes aux sels, sur des sols bien drainés ou de bonne perméabilité et
la salinité doit être contrôlée. La classe C3S1 représente les puits qui se situent au Nord de la ville
de Meskiana et ceux qui se trouvent dans la région de Dalâa. Elle représente 35.7 % des puits des
six campagnes, alors que la classe C3S2 ne représente qu’un pourcentage de 6.5 % ;
-Les classes C4S1, C4S2: elles représentent les eaux médiocres fortement minéralisées,
susceptibles de convenir à l’irrigation de certaines espèces bien tolérantes aux sels et sur des sols
bien drainés et lessivés. La majorité des puits appartenant à la classe C4S1 se situent entre Zebar
et Meskiana et représentent 18.08% des points d’eau analysés. La classe C4S2 est d’une sodicité
plus importante que celle de la classe précédente. On la rencontre dans la région de Zebar et au
Nord-Est de la plaine de Meskiana, avec un pourcentage de 30. 82 %.
4
4
30 C1S4 C2S4 C3S4 C4S4 30 C1S4 C2S4 C3S4 C4S4
Très élevé
Très élevé
3
3
Sodium Absorption Ratio (SAR)
C1S3 C1S3
20 C2S3 20 C2S3
Risque de Sodicité
Risque de Sodicité
Moyen
Moyen
2
C3S3 C3S3
C1S2 C1S2
C4S3 C4S3
C2S2 C2S2
10 C3S2 10 C3S2
8
.C S.
4 2 8
.C S
4 2
C1S1 P22..P29
P23
C1S1 ..
P23 P22
Faible
Faible
... . .
P3 P41 P32 P52
P8. P1 P15 P1
2 3 4 5 6 7 8 9 1000 2 3 4 2 3 4 5 6 7 8 9 1000 2 3 4
Conductivité Spécifique en µS/cm (CE*10) à 25°c Conductivité Spécifique en µS/cm (CE*10) à 25°c
1 2 3 4 1 2 3 4
Faible Moyen Elevé T. Elevé Faible Moyen Elevé T. Elevé
Risque de Salinité Risque de Salinité
Mars 2004 Juillet 2004
Figure 72b: Diagrammes de Richards pour l’année, 2004
Les eaux de la nappe alluviale de la plaine de Meskiana présentent, en général, les mêmes
qualités pour l’irrigation au cours de toute la période d’observation (tableau 20). La dégradation
observée traduit l’effet de la lithologie (présence de gypses et d’argiles) d’une part et l’effet de
l’évaporation d’autre part (augmentation des concentrations).
92
Tableau .20: Evolution des classes d’irrigation selon la méthode de Richards.
Mars2002 Juillet 2002 Mars2003 Juillet 2003 Mars2004 Juillet 2004
290
285
Rehia
280
Bellala
69
64 65
275
60
270 54 52Meskiana
49
48 50
47
45 46
265 43 4144 42
38
35 34
260 31 322729 Légende
2223 Zebar
21 20 Eaux médiocres à l’irrigation
17 18
255
16 13 15
12 Eaux admissibles à l’irrigation
250 10
8 Eaux bonnes à l’irrigation
6
245 4 3
Dalâa Agglomération
Echelle
240 1 0 10 20 Km
93
V.8.2.2 Classification des eaux par la méthode de Wilcox :
Généralement, les plantes supportent mal les sols saturés en sodium. La classification de
Wilcox fondée sur la conductivité électrique et la teneur en sodium dans l’eau exprimée en
pourcentage. La représentation des différents échantillons sur ce diagramme permet la
caractérisation des eaux pour leur aptitude à l'irrigation (figure 74a et 74b). Le %Na+ est définie
par la relation :
( Na + + K + )
%Na + = × 100
(Ca + 2 + Mg + 2 + Na + + K + )
(17) Tous les ions sont exprimés en méq/l.
100 100
Mauvaise Mauvaise
90 90
Médiocre Médiocre
80 80
70 Admissible 70 Admissible
Excellente
Excellente
60 ..P23 60 .P18 .P34
% Na
% Na
P18 . P22
.P17 P41..P29
P17 . .P29P41. .P34 P23.. P38
50 . 50 P64. . P45
.P38
P69
Bonne P60..P54
P21
Bonne P35. . P21
40 P16 . . . .P45 40
P46 . P49..P32 . P43
P13 . .P46 P47
P62 ..P32 .P35 .P43 .P16 P27. .P48
P10 . P49. . P27 .P50 P10 . .P13
.
P12 P52. .P50
30 P6 . P51. 30 P31 . P4. P6 . .P42
P31 . ..
P12 P42 .P44 P65. . P20 P65.
P3 . P3 .
20 . P20 20
P15 . P8 .
P1 . . P8 P15 .
10 10 P1 .
0 0
400 1000 2000 3500 400 1000 2000 3500
Conductivité µs/cm à 25 °C Conductivité µs/cm à 25 °C
Mars 2002 Juillet 2002
100 100
Mauvaise Mauvaise
90 90
Médiocre Médiocre
80 80
70 Admissible 70 Admissible
P22. . P23
Excellente
Excellente
60 . 60 . . P18
% Na
% Na
P17 P17
. P23
50 . . ..P34
P38 P22
50 P29. P38.
P29
.P34
Bonne . P13 P64. P69 Bonne P64. .P69
. P12 P54. .P60 P12..P13 P68. P60. .P54
40 P16. P49. . P62.P43 40 .P10 . P46 P49. P43. . P48
P10. .P41 P32. .P27. P50 P31. P6. P32. . P27 . P47
P31. .P6 P51. P4. P42. . P44
30 P4.
30
P44. P68. P65.
P3. .P42 P58. .P65 .P15 P58.
P1..P8
20 20
..
P1 P8 P15
10 10
0 0
400 1000 2000 3500 400 1000 2000 3500
Conductivité µs/cm à 25 °C Conductivité µs/cm à 25 °C
Juillet 2003
Mars 2003 Figure 74a : Diagrammes de Wicox pour les années 2002 et 2003.
94
100 100
Mauvaise Mauvaise
90 90
Médiocre Médiocre
80 80
70 Admissible 70 Admissible
Excellente
Excellente
60 . . P23 60 .
% Na
% Na
P18 P17
. P17 .P22.P38 ..
P23 P22
P29. .P34
50 50 P38
P34 ..
Bonne P64. Bonne .
P69
P13. . P16 P68 .P60. .. .
P12 P13 P16
P60
. .
P54
40
P46 P49. P48. P52. 40 P10 . P49P32. . ..
P48
P47
.P10 . . P41 P50.
. P6 . .
P41 P29 . .P50
P31. P51.
P6.
..
P31
30 30 P44 P46
. P42 P3 . . P42 P65 .
P3. P58..P65 P58 .
20 P15 20 P8 . .P15
. .P8
P1. P1.
10 10
0 0
400 1000 2000 3500 400 1000 2000 3500
Conductivité µs/cm à 25 °C Conductivité µs/cm à 25 °C
Juillet 2004
Mars 2004
Figure 74b: Diagrammes de Wilcox pour l’année 2004.
Les diagrammes couvrant toute la période d’observation montrent que les eaux
souterraines de la nappe de Meskiana, selon la classification de Wilcox, appartiennent aux quatre
classes : Excellente, Bonne, Médiocre et Mauvaise (tableau 21).
- La classe Excellente : Cette classe n’apparaît que dans la campagne de mars 2002, elle
caractérise les eaux provenant des formations carbonatées de bordures, surtout dans la région de
Dalaâ ;
- La classe Bonne : Elle regroupe les eaux qui sont faiblement minéralisées, on la rencontre à
proximité des bordures dans les régions de Dalaâ et Zebar ;
- La classe Médiocre : Elle regroupe les eaux plus minéralisées que celles des classes précédentes,
elle se manifeste au centre, au sud Ouest et au Nord-Est de la plaine.
- La classe Mauvaise : Elle regroupe les eaux qui se situent au centre et au Nord de la plaine. Elle
est la plus répandue dans la nappe (54.84 %).
On constate que les eaux souterraines subissent une dégradation selon le sens de l’écoulement,
montrant ainsi l’effet de la lithologie sur leur qualité (figure 75).
95
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
69
64 65
275
60
270 54 52 Meskiana
49
48 50
47
45 46
265 43 4144 42
38
3534
260 31 322729 Légende
Zebar
2223 Eaux de mauvaise qualité pour l’irrigation
21 20
17 18
255
16 13 15 Eaux de qualité médiocre pour l’irrigation
12
250 10
8 Eaux de bonne qualité pour l’irrigation
6
Dalâa Eaux d’excellente qualité pour l’ l’irrigation
245 4 3 Echelle
Agglomération
240 1 0 10 20 Km
Remarque : En examinant les deux cartes d’aptitude des eaux à l’irrigation établies à l’aide des
deux méthodes de Richards et Wilcox, on remarque qu’elles présentent la même distribution de la
qualité. En effet, la classe ‘bonne’, donnée par la première méthode, correspond à celles
‘excellente’ et ‘bonne’, donnée par la deuxième méthode. La classe ‘admissible’ définie par la
méthode de Richards correspond aux classes ‘bonne’ et ‘médiocre’ définies par la méthode de
Wilcox. La classe ‘mauvaise’ est déterminée par les deux classifications.
V.9 CONCLUSION:
L’étude hydrogéochimique des eaux de la nappe du Mio-Plio-Quaternaire de Meskiana
permet de dire que :
-Les eaux sont caractérisées par une forte salinité et par la présence de quatre faciès dominants,
qui sont le chloruré calcique, le chloruré sodique, le sulfaté calcique et le bicarbonaté calcique ;
-La répartition spatiale des éléments chimiques montre que le chimisme est fortement lié à la
lithologie de la nappe, suite à la dissolution des formations carbonatées et des formations
évaporitiques ;
-Le suivi de la qualité des eaux dans le temps a montré l’effet des précipitations et de
l’évaporation sur la minéralisation des eaux souterraines ;
-L’évolution des éléments chimiques (Na+, Cl-, Ca2+, SO42-, HCO3-), qui caractérisent les
formations géologiques principales du bassin, a montré la dominance des ions salifères ((Na+, Cl-
) et des ions gypsifères (Ca2+, SO42-) par rapport à ceux carbonatés (Ca2+, HCO3-) dans
l’acquisition de la salinité. Ce bilan a également mis en évidence le phénomène d’échange de base
entre les eaux de la nappe et les passages argileux, ce dernier a été révélé par l’augmentation des
teneurs en Na+ par rapport à Cl- .
96
-Le calcul de l’indice de saturation des minéraux carbonatés et gypsifères a montré que les
premiers ont tendance précipiter avant les seconds, ceci a permis aux éléments chimiques
provenant de la dissolution des gypses d’acquérir des teneurs importantes.
- L’analyse en composante principale a confirmé le mode de minéralisation des eaux souterraines,
car les puits qui se situent en amont de la nappe se trouvent faiblement minéralisés au contraire de
ceux qui se positionnent en aval, qui sont trop minéralisés.
-La classification de Richards et celle de Wilcox ont prouvé la dégradation de la qualité des eaux
pour l’irrigation.
97
Chapitre VI :
Durant nos sorties sur terrain, on a pu faire un inventaire des différentes sources de
pollution (figure 76), qui se résument comme suit:
-L’élevage des animaux (ovins, bovins et volailles), qui produit environ 450.000 kg d’azote par
an, cette estimation ayant été faite par les services agricoles de la Wilaya d’Oum El Bouaghi en
2003;
-Les cultures maraîchère et céréalière, qui utilisent d’importantes quantités de fertilisants sous
forme de fumier, à raison de 280 Kg /ha par an, ou d’engrais azotés (100 à 300 Kg/ha);
-L’usine de l’ELATEX (actuellement en arrêt) ;
-Les rejets urbains, qui déversent des eaux usées domestiques dans la plaine, sans traitement
préalable ;
-Les décharges ménagères sauvages.
Ce recensement nous a permis de choisir un réseau d’une quarantaine de puits pour le
suivi de l’évolution de ce type de pollution. Le choix de ce réseau de surveillance s’est fait à
l’aide de critères lithologiques et hydrodynamiques. Huit campagnes de prélèvements ont été
menées à raison de deux par an.
Ces différentes campagnes ont permis le dosage de quelques polluants, mentionnés dans le tableau
22.
98
925 930 935 940 945 950 955 960 965
290
285
Rehia
280
Bellala
275
270
Meskiana
265
Légende
Industrie (textiles)
260
Elevage animal
Zebara
255 Cultures maraîchère et céréalière
Décharge ménagère sauvage
250
Rejet des eaux domestiques et industrielles
Dalâa
245 Agglomération
Echelle
240 0 10 20 Km
235
99
VI.3 IMPACT DE L’ACTIVITE AGRICOLE :
100
Tableau 23: Estimation de l’apport annuel en azote à partir de l’élevage.
101
930 935 940 945 950 955 960
290
125
285 120
Rehia 115
280 110
BellalaR
105
275 100
95
270 Meskiana 90
85
265 80
75
260 70
Zebar
65
255 60
55
250 50
Dalâa 45
245 Echelle 40
35
240 0 10 20Km
Figure 78: Carte des nitrates (mg/l) dans la nappe de Meskiana pour la
période de novembre 2005.
102
mars-02 juil-02 mars-03 juil-03
mars-04 juil-04 déc-04 nov-05
140
120
NO3 (mg/l)
100
80
-
60
40
20
0 P1 0
P1 3
P1 5
P1 6
P1 7
P2 2
P2 3
P3 1
P3 4
P3 8
P4 1
P4 2
P4 6
P4 8
P4 9
P5 0
P5 2
P5 4
P6 0
P6 4
P6 5
P6 9
P1
P3
P6
P8
Puits
140
120
NO3 dans l'eau (mg/l)
Origine différente A
100
80
Origine agricole
60
40
20
0
0 0,02 0,04 0,06 0,08 0,1
N dans le sol (mg/g)
Figure 80: Corrélation de la concentration des nitrates avec l’azote du sol.
103
VI.3.4 Evolution des nitrates en fonction des paramètres physico-chimiques :
0,3
NO2 (mg/l)
0,25
0,2
-
0,15
0,1
0,05
0
20 30 40 50 60 70 80 90 100
-
NO3 (mg/l)
Figure 81:Evolution des nitrates en fonction des nitrites, pendant le mois
de juillet 2004.
0,45
0,4
0,35
0,3
NH4 (mg/l)
0,25
+
0,2
0,15
0,1
0,05
0
20 40 60 80 100 120
-
NO 3 (mg/l)
Figure 82: Evolution des nitrates en fonction de l’ammonium,
pendant le mois de juillet 2004.
104
Quelques points sont observés avec des teneurs considérables en oxygène accompagnées
de faibles teneurs en nitrates (figure 83), on pense qu’il s’agit de puits qui se trouvent en dehors
des zones agricoles.
VI.3.4.4 Relation nitrate-profondeur :
La teneur des nitrates évolue d’une manière inverse avec la profondeur du niveau
statique des puits. Ainsi, la teneur de l’oxygène dissous diminue avec la profondeur, ce qui
entraîne la réduction des nitrates en nitrite et en ammonium (figure 84).
4
3,5
3
2,5
O2 (mg/l)
2
-
1,5
1
0,5
0
20 40 60 80 100 120
-
NO3 (mg/l)
Figure 83: Evolution des nitrates en fonction de l’oxygène
dissous pendant le mois de juillet 2004.
120
100
80
NO3 (mg/l)
60
-
40
20
0
0 2 4 6 8 10 12 14 16 18 20 22
profondeurs (m)
Figure 84: Evolution des nitrates en fonction de la profondeur du niveau
statique pour le mois de juillet 2004.
Les agriculteurs utilisent les engrais phosphatés d’une manière aléatoire avec une dose
qui varie de 15 à 25 Kg/ha, sans respecter les recommandations de l’ Institut National de la
105
Vulgarisation Agricole (INVA) et de l’Institut Technique des Grandes Cultures (ITGC) qui fixent
cette dose à 46 kg/ha. Les engrais utilisés sont fabriqués à partir des roches phosphatées. Le
phosphore présent dans ces roches n’est pas assimilable directement par les plantes, surtout dans
les sols basiques. La production industrielle de ce dernier se fait selon le processus suivant :
L’acide phosphorique constitue la base de fabrication de tous les types d’engrais phosphatés.
En ce qui nous concerne, nous avons procédé à deux campagnes d’analyses de cet élément
(décembre 2004 et novembre 2005).
L’évolution des phosphates est semblable à celle des nitrates (figure 85).En effet,
l’examen des cartes des phosphatés a montré une certaine ressemblance avec celle des nitrates, ce
qui laisse penser que leur origine est anthropique résultant ainsi de l’activité agricole et non pas
d’une origine naturelle (figures 86a et b).
Les teneurs en phosphates enregistrées varient de 0.25 mg/l à 0.5 mg/l, elles sont
inférieures à la norme de potabilité.
100
Teneurs (mg/l)
80
60
40
20
0
P1
P6
P10
P13
P16
P22
P27
P32
P41
P44
P47
P49
P54
P60
P64
Puits
Figure 85: Corrélation de la teneur des nitrates avec celle des
phosphates pour la période de décembre 2004.
106
930 935 940 945 950 955 960
290
0.45
285 0.44
Rehia 0.43
280 Bellala 0.42
0.41
275 0.4
0.39
270 0.38
Meskiana
0.37
0.36
265
0.35
0.34
260 0.33
0.32
Zebar
255 0.31
0.3
250 0.29
0.28
Dalâa 0.27
245 Echelle
0.26
0.25
240 0 10 20 Km
290
0.47
0.46
285 0.45
Rehia 0.44
280 Bellala 0.43
0.42
275 0.41
0.4
270 Meskiana 0.39
0.38
265 0.37
0.36
260 0.35
Zebar 0.34
255 0.33
0.32
250 0.31
Dalâa 0.3
0.29
245
Echelle 0.28
240 0 10 20 Km
107
VI.4 L’INFLUENCE DES REJETS URBAINS :
Les rejets urbains, qui se déversent dans la plaine sans subir de traitement préalable, sont
chargés de matières organiques et produisent des nitrates. Par le processus d’organisation, cette
transformation se fait selon les processus décrits par NICOLARDOT.B et al (2000). Figure 87.
En effet, après l’élaboration de l’évolution des éléments majeurs, on a remarqué une
certaine corrélation entre les nitrates et les chlorures (voir la matrice de corrélation au chapitre
précédent), où le coefficient de corrélation est de l’ordre de 0.62 pour une population de 160 puits
(ACP globale). En d’autre terme, les chlorures peuvent avoir plusieurs origines naturelles et
anthropiques (écart type 265.94: mg/l , tableau 15). On pense qu’une partie des chlorures tire son
origine des eaux usées domestiques qui constituent, en même temps une source de nitrates
d’origine organique.
L’évolution des nitrates et des chlorures, à travers les différentes campagnes, montre cette
dépendance. Dans le but de confirmer cette dernière, on a fait une comparaison entre l’évolution
des chlorures et celle des nitrates.
108
EFFLUENT N2O
+
NH4 N2
N organique NH3
Volatisation
AIR
SOL
Dénitrification
Minéralisation
Matière organique NH4+ Nitrification
du sol NO3-
Organisation
Lixiviation
Figure 87: Principaux processus de transformation des formes d’azote
contenues dans les effluents et déchets, (NICOLARDOT.B et al 2000).
EAU
NO3 Cl
1050
900
Concentrations (mg/l)
750
600
450
300
150
0
P1
P4
P8
P12
P15
P17
P20
P22
P27
P31
P34
P38
P42
P44
P46
P48
P50
P52
P58
P62
P65
P69
Puits
Figure 88a: Evolution des nitrates en fonction des chlorures, pour la période
de mars 2002.
NO3 Cl
1050
900
Concentrations (mg/l)
750
600
450
300
150
0
P1
P4
P8
P12
P15
P17
P20
P22
P27
P31
P34
P38
P42
P44
P46
P48
P50
P54
P64
P69
Puits
Figure 88b: Evolution des nitrates en fonction des chlorures, pour la
période de juillet 2002.
109
NO3 Cl
1050
900
Concentrations (mg/l)
750
600
450
300
150
0
P1
P4
P8
P12
P15
P17
P22
P27
P31
P34
P41
P43
P46
P48
P50
P52
P58
P62
P65
P69
Puits
Figure 88c: Evolution des nitrates en fonction des chlorures, pour la période
de mars 2003.
NO3 Cl
1050
Concentrations (mg/l)
900
750
600
450
300
150
0
P1
P4
P8
P12
P15
P17
P22
P27
P31
P34
P41
P43
P46
P48
P50
P52
P58
P62
P65
P69
Puits
Figure 88d: Evolution des nitrates en fonction des chlorures, pour la
période de juillet 2003.
NO3 Cl
1050
Concentrations (mg/l)
900
750
600
450
300
150
0
P1
P6
P10
P15
P17
P22
P29
P34
P41
P46
P49
P51
P54
P60
P65
P69
Puits
Figure 88e: Evolution des nitrates en fonction des chlorures, pour la
période de mars 2004.
110
NO3 Cl
1050
900
Concentrations (mg/l)
750
600
450
300
150
0
P1
P6
P10
P13
P16
P22
P27
P32
P38
P42
P46
P48
P50
P54
P60
P64
P69
Puits
Figure 88f: Evolution des nitrates en fonction des chlorures, pour la période
de juillet 2004.
VI.4.2 Le bore :
Il est présent à l'état de trace dans les eaux souterraines, son origine naturelle est très rare et
son existence ne peut qu’être anthropique. Il provient essentiellement des rejets industriels, urbains ou
domestiques, par sa présence en forte concentration dans les détergents sous forme de perborates (20
% de la composition des lessives en poudre); il est parfois utilisé comme traceur de pollution dans les
études d'assainissement.
L’observation de la cartographie du bore dans la nappe Mio–Plio-Quaternaire de Meskiana
(figure 89), permet de dire que cet élément constitue un indicateur de la pollution urbaine, car les
fortes concentrations sont liées aux points de rejets des eaux usées domestiques. Les teneurs
maximales enregistrées sont de l’ordre de 3.5 mg/l et se localisent en aval des agglomérations dans
l’axe de la nappe.
Le bore est un élément indispensable à la croissance des plantes et est directement appliqué au
sol comme engrais. Il est fortement mobile dans ce milieu et facilement entraîné par lessivage.
L’industrie des engrais, par le complexe ASMIDAL, n’utilise pas cet élément comme additif dans les
engrais, c’est pourquoi on a attribué sa présence à la pollution urbaine et non pas à l’activité agricole.
Il est à signaler que la dissolution des sédiments triasiques peut donner aussi naissance au bore.
VI.5 LES METAUX LOURDS :
Lors des deux campagnes de décembre 2004 et novembre 2005, on a effectué l’analyse de
quelques métaux lourds tels que le fer, le plomb, le zinc, le manganèse le cuivre, le cadmium et le
soufre. La concentration de ces éléments est insignifiante à l’exception de celle du fer qui présente de
fortes teneurs.
VI.5.1 Le fer :
Le fer est le quatrième élément abondant dans l’écorce terrestre et le premier parmi les métaux
lourds. On le trouve surtout sous forme Fe++ ou Fe+++. Le fer dissous dans les eaux souterraines est
sous forme réduite (Fe++), qui est soluble. Ce dernier est oxydé en Fe+++ par le contact avec l'oxygène
de l'air, c’est pourquoi la concentration de fer dans les eaux bien aérées est rarement élevée. Le
Fe+++précipite sous forme d'hydroxydes insolubles dans l'eau et sont précipités ou adsorbés par les
sédiments.
Dans la nappe de Meskiana, lors des deux campagnes du mois de juillet 2004 et du mois de
décembre 2004, les concentrations de fer varient de 0.1 à 1.2 mg/l (figure 90a et b).
Quelques puits, qui se situent au Nord de la ville de Meskiana, atteignent des teneurs qui
dépassent le niveau recommandé par l’OMS (0.3 mg/l). Sa présence est probablement due aux rejets
de l’usine d’ELATEX et/ou au lixiviat de la décharge ménagère de la ville de Meskiana. On pense
que sa présence pourrait être due également au matériel de remplissage alluvionnaire.
111
930 935 940 945 950 955 960
290
3.4
285 3.2
3
Rehia
280 Bellala 2.8
2.6
275 2.4
2.2
270 Meskiana
2
1.8
265
1.6
1.4
260
1.2
Zebar
255 1
0.8
250 0.6
Dalâa 0.4
245 Echelle 0.2
0
240 0 10 20 Km
290
1.2
285 1.15
1.1
Rehia 1.05
280 Bellala 1
0.95
275 0.9
0.85
0.8
270 Meskianae 0.75
hi 0.7
265 0.65
0.6
0.55
260 0.5
Zebar 0.45
255 0.4
0.35
0.3
250 0.25
Dalâa 0.2
0.15
245 0.1
Echelle
0.05
240 0
0 10 20 Km
112
930 935 940 945 950 955 960
290 1.1
1.05
285 1
0.95
Rehia 0.9
280 Bellala
0.85
0.8
275
0.75
0.7
270 Meskiana 0.65
0.6
265 0.55
0.5
0.45
260
0.4
Zebar
0.35
255 0.3
0.25
250 0.2
Dalâa 0.15
245 0.1
Echelle 0.05
240 0 10 20Km
VI.5.2 Le strontium :
Le strontium est un métal alcalino-terreux, et il est présent naturellement dans
l'environnement (les roches, le sol, l'eau). Les composés du strontium peuvent se déplacer dans
l'environnement assez facilement car beaucoup d'entre eux sont solubles. Généralement, on le
trouve dans les minéraux de strontianite (SrCO3) et de célestite (SrSO4). Cette dernière
accompagne les formations évaporitiques, car c’est un minéral qui est souvent associé au gypse
(CaSO42H2O), (Carre, 1975, cité par Hsissou 1996). Il peut également s’intégrer en faible
quantité dans le réseau des carbonates et est peu adsorbé par les argiles.
Le strontium peut s’intégrer dans la structure cristalline des minéraux carbonatés, car ces
deux éléments possèdent des rayons ioniques voisins.
Les concentrations excessives en strontium peuvent être provoquées par des procédés
industriels.
113
930 935 940 945 950 955 960
290
6.2
285 5.8
Rehia 5.4
280 Bellala 5
4.6
275
4.2
270 Meskiana 3.8
3.4
265
3
260 2.6
Zebar 2.2
255 1.8
1.4
250
1
290 5.6
5.4
5.2
285 5
4.8
Rehia 4.6
280 4.4
Bellala
4.2
4
275 3.8
3.6
270 3.4
Meskiana 3.2
3
265 2.8
2.6
2.4
260 2.2
2
Zebar 1.8
255 1.6
1.4
1.2
250 1
0.8
0.6
245 Dalâa 0.4
Echelle
0.2
240 0 10 20 Km
114
0,2
0,18
0,16
Teneur du Sr en (méq/l)
0,14
0,12
0,1
0,08
0,06
0,04
r = 0.77
0,02
0
0 1000 2000 3000 4000 5000 6000
Conductivité électrique en us/cm
Figure 92: Relation de la teneur du Sr avec la conductivité électrique pendant le
mois de juillet 2004.
Pour Meybeck (1984, cité par Hsissou Youssef 1996), les rapports des concentrations
molaires Sr2+/Ca2+ supérieurs à 1‰ sont caractéristiques des eaux gypsifères. De plus, comme
chaque ensemble évaporitique possède sa propre valeur du rapport Sr2+/Ca2+, les teneurs en Sr2+
et en Ca2+ des eaux ayant dissous les formations d’une région devraient permettre d’identifier les
circulations souterraines et de distinguer leurs origines.
Le calcul du rapport Sr2+/Ca2+, pour la campagne de juillet 2004, a révélé que la majorité
des points d’eau possèdent des valeurs supérieures à 1‰ (figure 93), ce qui traduit l’influence des
formations gypsifères sur la présence du strontium. Cette affirmation est en pleine conformité
avec les types de faciès rencontrés, où les puits dont le faciès est carbonaté se répartissent en
amont de la nappe et en faible pourcentage, alors que les puits ayant un faciès chloruré ou sulfaté
occupent la majorité de la nappe, surtout au Nord de cette dernière.
930 935 940 945 950 955 960
290 0.012
0.0115
285 0.011
0.0105
Rehia Bellala 0.01
280
0.0095
0.009
275 0.0085
0.008
270 0.0075
Meskiana
0.007
0.0065
265 0.006
0.0055
260 Zebar 0.005
0.0045
0.004
255
0.0035
0.003
250 0.0025
Dalâa 0.002
245 Echelle 0.0015
0.001
240 0 10 20 Km
VI.6.1INTRODUCTION
Après avoir achevé la partie précédente, qui a montré l’état de contamination de la nappe
de Meskiana, il s’est avéré intéressant de cartographier les zones vulnérables à la pollution. Cette
cartographie pourrait aider les gestionnaires du secteur de l’hydraulique pour mieux gérer cette
ressource hydrique en la protégeant contre la pollution et également les aider à implanter les
ouvrages de captage en dehors des zones contaminées.
L’objectif visé par l’établissement de la cartographie de la vulnérabilité est de faire
apparaître les possibilités de pénétration et de circulation des polluants dans la nappe. On peut
aussi évaluer la sensibilité des ouvrages de captage au risque de pollution et par conséquent le
limiter
Di=Dn.dp+Rn.Rp+An.Ap+Sn.Sp+Tn.Tp+In.Ip+Cn.Cp où (18)
116
n : notation accordée à chaque paramètre.
117
Chaque paramètre est affecté d’un poids qui traduit son influence dans l’estimation de
l’Indice de Pollution Drastic. Le tableau 27 montre les poids et les propriétés des sept paramètres.
VI.6.4Validité de la méthode :
La précision avec laquelle la méthode Drastic permet de distinguer les régions
vulnérables de celles qui sont protégées a été vérifiée dans différentes régions climatiques des
États-Unis (Maine, Kansas, Texas, Caroline du Sud, Floride, Idaho, Wisconsin, etc.), ainsi
qu'au Québec. Dans ce dernier cas, ISABEL et al (1990) ont vérifié la validité de la méthode
sur des unités hydrogéologiques différentes. Dans tous les cas, les zones réellement
contaminées correspondent à celles où les indices de vulnérabilité sont les plus élevés.
118
VI.6.5.1Profondeur de la nappe :
La profondeur de la nappe est un paramètre important, qui constitue un facteur
déterminant de la vulnérabilité. Plus cette profondeur est élevée, plus le contaminant met de
temps à atteindre la surface piézométrique.
Pour la nappe de Meskiana, ce paramètre a été déterminé pendant les relevés
piézométriques. On a remarqué que dans certains puits la profondeur est parfois nulle, ce qui
met en contact les eaux de surface avec les eaux souterraines et facilite la contamination de
ces dernières.
Le calcul des indices pour ce paramètre a été fait selon le tableau 28.
Tableau 28 : Notation et pondération du paramètre de la profondeur
Intervalle Dn Dp Dn.Dp
0 – 1.5 10 5 50
1.5 – 4.5 9 5 45
4.5 – 9 7 5 35
9 - 15 4 5 25
15 - 22 3 5 15
22.5- 30 2 5 10
Ces indices ont permis d’établir la carte thématique du paramètre profondeur (figure
94). L’observation de la dite carte montre que les plus faibles indices sont répartis dans la
zone de Dalâa et à proximité des bordures. Leurs valeurs ne dépassent pas 25, alors que les
plus importants indices se situent dans les zones centrale et septentrionale, avec des valeurs de
35 à 50.
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
Meskiana
270
265
260
Légende
Zebar
255
Profondeur ( de 0 à 9 m) Indice ( de 35 à 50 )
240 0 10 20 Km
119
VI.6.5.2La recharge efficace :
Elle concerne la quantité d’eau qui franchit la surface du sol pour renouveler les
réserves d’eau souterraine et entretenir le débit de l’écoulement souterrain, après circulation
dans l’aquifère. Les résultats du bilan hydrologique montrent que la valeur de ce paramètre est
estimée à 5 mm sur l’ensemble du bassin.
290
285
Rehia
280
Bellala
4
275
Meskiana
270
4
265
260 4
Zebar
255
Légende
4
250
Dalâa Recharge:5 mm Indice drastique:4
245
Echelle
240 0 10 20 Km
120
Tableau 30 : Notation et pondération du matériau aquifère.
Calasses An Ap An.Ap
Gravier et conglomérat 10 3 30
Sables et graviers 9 3 27
Sables argileux 4 3 12
L’observation de la carte des matériaux de la zone saturée montre que les indices les
plus faibles se situent à proximité des bordures et les plus importants au centre de la plaine.
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
Meskiana
270
265
260
Zebar Légende
121
L’examen de la carte des indices relatifs à la nature du sol montre que le sol argilo-
graveleux se répand sur la majorité de la surface de la plaine. Le sol argileux se répartit au
Nord de la plaine, alors que les sols sableux et conglomératique occupent des superficies
limitées au centre et au Sud de la plaine.
290
285
Rehia
280
Bellala
275
Meskiana
270
265
Légende
260
Zebar galets et conglomérats (I = 20)
255 Sables et gravier (I = 16)
Dalâa Argile (I = 2)
245
Echelle
240 0 10 20 Km
122
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
270 Meskiana
265
260
Légende
Zebar
255 Pente (0 à 2 %) (I = 30)
Dalâa
245
Echelle
240 0 10 20 Km
La carte des indices correspondant à la nature de la couche vadose montre que deux
zones sont très vulnérables selon ce paramètre : la première se situe au Sud de la plaine à
proximité de Dalâa et la seconde se localise au centre de la nappe. Au Nord, les indices sont
moins importants mais peuvent dépasser la valeur de 40. On signale la présence d’une
couverture argileuse qui occupe une surface limitée à coté des bordures occidentale et
orientale
123
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
Meskiana
270
265
Légende
Dalâa
245 Echelle
240 0 10 20 Km
VI.6.5.7La perméabilité
La perméabilité de la couche aquifère nous renseigne sur la vitesse de propagation des
polluants dans la nappe. Plus ce paramètre est important, plus le transfert des polluants est
rapide. L’estimation des indices de ce paramètre a été faite selon les classes présentées dans le
tableau 34.
Perméabilité en Cn Cp Cn.Cp
(10-4 m/s)
0.65 à 0.75 8 3 24
0.35 à 0.65 5 3 15
0.05 à 0.35 3 3 9
> 0.05 1 3 3
La carte des perméabilités montre que les plus importantes valeurs (0.65 à 0.75.10-4
m/s) se situent au centre de la plaine entre Zebar et Meskiana. Des valeurs moins importantes
(0.35 à 0.65.10-4 m/s) se rencontrent au Nord et au Sud de la plaine. Des perméabilités
moyennes (0.05 à 0.35.10-4 m/s) se trouvent sur la majorité de la plaine, alors que les plus
faibles perméabilités se localisent sur une zone limitée à l’Est de la nappe .Figure 100.
124
930 935 940 945 950 955 960
290
285
Rehia
280
Bellala
275
Meskiana
270
265
260 Légende
Zebar
255 I = 24
I = 15
250
I=9
Dalâa
245 I=3
Echelle
240 0 10 20 Km
125
925 930 935 940 945 950 955 960 965
290
285
Rehia
280
Bellala
69
65
275 64
MK1
60
A3BV5
F1
270
49 52 Meskiana
48
47 B8 46
45
265 43 CD7 4144 B10
42
38
E5 Légende
34
260 31 29 MK3
Zebar
150 < I Vulnérabilité forte
22 MK4
21
18
MK2 20 150 < I < 100 Vulnérabilité moyenne
255 17 K5
16 AL5 13 15
12 I < 100 Vulnérabilité faible
10
250
8 AL5 Forage hydraulique
6 M13 X14
Dalâa 43 Puits domestique ou agricole
245 4 3 D1
Agglomération
Echelle
240 1 0 10 20 Km
235
Croûte calcaire Grès grossier Marne verdâtre Calcaire
Marne grise Marne et calcaire Gypse
Marne
Remarque:
Cette approche présente une application intégrale de la méthode Drastic. Appliquée à
la plaine de Meskiana, qui se situe en milieu semi-aride, elle a permis de fournir une carte de
vulnérabilité des eaux souterraines de l'aquifère Mio-Plio-Quaternaire. À l'issue de cette
première approche, l'adéquation des zones très vulnérables avec celles où sont observés des
cas de pollution a permis de vérifier l’efficacité de la méthode.
VI.6.6 CONCLUSION :
Les différents points abordés le long de ce chapitre ont fait état de la pollution de la
nappe souterraine de Meskiana, mais il était difficile d’évaluer l’importance de cette
pollution. Notre approche a permis de préciser l’importance de ce problème et de préciser son
origine. Ainsi, la nappe superficielle du Mio-Plio-Quaternaire semble être fortement
contaminée par des nutriments d’origine agricole, suite à l’utilisation de fertilisants plus que
126
nécessaire pour les plantes. D’autre part, les rejets urbains, déversés sans aucun traitement
préalable dans le réseau hydrographique, sont une source importante de pollution.
L’examen de l’évolution des formes de l’azote en fonction de la lithologie et de
l’hydrodynamique des eaux souterraines a permis d’aboutir aux conclusions suivantes :
- La contamination de la nappe est accentuée surtout là où le niveau piézométrique est
relativement faible. Le pompage intensif et régulier des eaux de la nappe dans les zones
profondes pourrait augmenter significativement la circulation des eaux souterraines, par appel
d’eau, et accélérer la contamination des eaux de la nappe à partir de la surface ;
- La nature lithologique de la couche vadose et celle de la couche aquifère ont favorisé la
contamination de cette dernière et influencé la répartition des nitrates dans l’espace;
- La répartition du bore dans la nappe a mis en évidence la pollution par les rejets urbains.
-Cette étude a montré la présence de fer dans les eaux souterraines, probablement due à la fois
à l’usine d’ELATEX et à la décharge ménagère de la ville de Meskiana, ainsi qu’aux rejets
des eaux usées domestiques.
On note l’effet de la contamination naturelle, marquée par le taux de salinité élevé
avec de fortes concentrations en Na+, Cl-, SO42- et Ca2+.
L’établissement de la carte de vulnérabilité à la pollution des eaux souterraines de la
plaine de Meskiana a montré le champ de vision à l'ensemble de la plaine et ne se limite pas
seulement aux endroits où la pollution a déjà contaminé les eaux souterraines. Elle permet
ainsi la prise de mesures de protection, notamment dans la partie couvrant le champ de
captage.
127
CONCLUSIONS GENERALES
ET
RECOMMANDATIONS
CONCLUSION GENERALE :
Le suivi de la qualité des eaux dans le temps et dans l’espace a montré l’impact des
facteurs climatiques (les précipitations et l’évaporation) sur la minéralisation des eaux de la
nappe.
128
La comparaison des concentrations des éléments chimiques (Na+, Cl-, Ca2+, SO42-,
-
HCO3 ) qui caractérisent les formations géologiques principales du bassin, a mis en relief la
dominance des ions salifères ((Na+, Cl-) et les ions gypsifères (Ca2+, SO42-) par rapport à ceux
carbonatés (Ca2+, HCO3-) dans l’acquisition de la salinité.
Le calcul de l’indice de saturation des minéraux carbonatés et gypsifères a montré que
les premiers ont tendance à s’équilibrer avant les seconds, ce qui explique les teneurs
importantes des éléments chimiques provenant de la dissolution des gypses.
L’analyse statistique a pu montrer que les eaux souterraines se trouvant en amont de la
nappe sont faiblement minéralisées et celles qui se trouvent en aval sont fortement
minéralisées.
Les classifications de Richards et celle de Wilcox ont montré la dégradation de la
qualité des eaux pour l’irrigation.
Quant à l’aspect environnemental, les investigations sur les sources de pollution de la
plaine de Meskiana ont fait apparaître une forte contamination de celle-ci par des nutriments
d’origine agricole et par les rejets urbains, qui sont déversés sans aucun traitement préalable,
ce qui contribue à la dégradation de la qualité des eaux.
L’interprétation de l’évolution des formes de l’azote en fonction de la lithologie et de
l’hydrodynamique des eaux souterraines a montré que :
- La contamination de la nappe est très accentuée à de faibles profondeurs et l’éventualité de
la contamination des eaux profondes est envisageable sérieusement, surtout avec
l’exploitation systématique actuelle de cette nappe.
- La nature lithologique de la couche vadose et celle de la couche aquifère a favorisé la
contamination de cette dernière et a joué un rôle sur la répartition des nitrates dans la plaine.
-Les métaux lourds n’ont pas été détectés par les analyses des différentes campagnes
d’observation à l’exception du fer.
RECOMMANDATIONS :
Les recommandations ci-dessous visent à élucider certains problèmes qui ont été observés, à
préciser les incertitudes résultant de l’étude actuelle et à proposer des interventions en
conséquence :
-Installer un réseau de piézomètres à travers la nappe, pour contrôler la fluctuation de la
surface piézométrique, pour prendre les mesures de gestion conservatoires et appropriées de la
nappe.
-Rentabiliser et équiper les infrastructures d’observation et de mesure des paramètres
climatiques, et ce en évitant les ruptures durant les observations
-Comme la nappe Mio-Plio-Quaternaire de Meskiana constitue la seule réserve hydrique
actuelle pour la région, il est recommandé de la recharger artificiellement avec les eaux de
surface. Ceci ne pourrait avoir lieu qu’après la programmation d’études approfondies de
géologie, de géophysique et d’hydrogéologie.
-Encourager les agriculteurs à améliorer leurs pratiques agricoles, puisqu’un lien significatif a
été établi entre la présence de nitrates dans l’eau des puits et la zone d’agriculture intensive.
Ce résultat est en soi un avertissement, puisque le seuil qui a été révélé atteste de l’influence
indéniable des activités anthropiques.
-Maîtriser le stockage et l’utilisation des matières organiques résultant de l’élevage animal.
-Arrêter le déversement direct des rejets urbains et industriels au niveau de l’Oued Meskiana,
en réalisant un traitement préalable de ces effluents avant qu'ils atteignent les oueds et
surveiller en permanence la qualité des eaux au niveau de l’oued.
129
-Réétudier l’emplacement des décharges ménagères des agglomérations en les installant loin
des zones vulnérables à la pollution, et ce pour éviter leurs effets sur les ressources
hydriques.
-Réaliser une étude de modélisation des écoulements souterrains pour connaître les
mécanismes de transfert des polluants, ce qui peut beaucoup aider à la protection de la nappe
contre la pollution.
130
REFERENCES
BIBLIOGRAPHIQUES
AGENCE DES BASSINS HYDROGRAPHIQUES (2002) : Les cahiers de l’agence 120
pages.
AHKOUK S. (2004) : Impact des fertilisants agricoles et du mode d’irrigation sur la qualité
des eaux souterraines en zones irriguées sous climat semi-aride : cas de la plaine des Chtouka.
Bassin du Souss-Massa, Maroc. Thèse doctorat. Université Ibnou Zohr Agadir 134 pages.
ANONYME, (1973-1975) : Fiche technique des forages : A3, AL5, B10, B8, CD7 D1, E2,
E5, K5, MK2, MK4, MK5.
ARNAUD W. et al (2004) : Développement d’un système harmonisé de surveillance de la
qualité des terres agricoles en Région wallonne anticipant la future directive européenne sur
les sols. Biotechnol. Agron. Soc. Pages 70-83.
BOUZIANI M.(2005) :L’eau dans tous ses états. Editions DAR EL GHARB.239 pages.
BRGM (2000) : Gestion des sites (potentiellement) pollués. Version 2. Edition BRGM
CASTANY G. (1963) : Traité pratique des eaux souterraines. Dunod Edit. Paris, 641 pages.
CASTANY G. (1968) : Prospection et exploitation des eaux souterraines. Ed. Dunod Paris
France 683 pages.
CASTANY G. (1982) : Principes et méthodes de l’Hydrogéologie. Dunod Edit. Paris,
334pages.
CITEAU LAËTITIA (2004) : Etude des colloïdes naturels présents dans les eaux gravitaires de sols
contaminés : relation entre nature des colloïdes et réactivité vis-à-vis des métaux (Zn, Cd, Pb, Cu). Thèse
doctorat. INRA, Unité de Science du sol – Versailles.147 pages.
COOTE D.R. GREGORICH (2000) : La santé de l’eau – Vers une agriculture durable au
Canada. Direction de l’agriculture et agroalimentaire Canada, Ottawa (Ontario).
DEVEZ AURELIE (2004) : Caractérisation des risques induits par les activités agricoles sur
les écosystèmes aquatiques. Thèse de doctorat. ENGREF, Centre de Montpellier.231 pages.
DUROZOY G. (1956) : Carte géologique au 1/50 000 de Tébessa feuille 206, et sa notice
explicative.
FAO (2005) : Utilisation des engrais par culture en Algérie. Division de la mise en valeur des
terres et des eaux. Rapport de Première édition, publié par la FAO, Rome. 44 pages
FELLRATH M. (1995) : Qualité des eaux. Monographies en Ecologie des Polluants (2).
Institut de Génie de l’Environnement, Laboratoire de Génie Sanitaire. Ecole Polytechnique
Fédérale de LAUSANNE 220 pages
.
FORGES M. (1972) : Irrigation et salinité. CIHEAM, options méditerranéennes. Pages 40-
45.
HSISSOU Youssef et al. (1995) : Caractérisation des eaux de l’aquifère turonien du bassin
du Tadla (Maroc) par le rapport des concentrations molaires Sr²+/Ca²+. Journal of Hydrology
183. Pages 445-451.
HSISSOU Y., MUDRY J., MANIA J., BOUCHAOU L., CHAUVE P. (1997) : Dynamique
et salinité de la nappe côtière d’Agadir (Maroc), influence du biseau salé et des faciès
évaporitiques. Fifth Scien. Ass. of IAHS, Rabat, Maroc, IAHS publi. n° 244.Pages 73-82.
LECOMTE P. (2000) : Les sites pollués, traitement des sols et des eaux souterraines.
CNRCCP. Edition Lavoisier TEC.&DOC .194 pages
MARSILY DE G. (1993) : les différentes approches pour l’élaboration d’un modèle, Actes
du séminaire Eaux environnement : Modélisation du comportement des polluants dans les
hydrosystèmes, Paris ministère de l’environnement. Pages 9-23.
NAFAA B.et al (2004) : La salinité des eaux souterraines de la région Est de la plaine d’El
Outaya (Biskra) Algérie. Colloque international Terre et Eau. Université de Annaba 04,05&
06 décembre 2004.
PAIN J.-M et al (1996) : Le bore marqueur de pollution urbaine, le bore naturel : état des
connaissances et illustrations régionales dans le bassin Artois-Picardie. Techniques sciences
méthodes, génie urbain génie rural. ISSN 0299-7258.
REMON E. (2005): Soil characteristics, heavy metal availability and vegetation recovery at a
former metallurgical landfill: Implications in risk assessment and site restoration.
Environmental Pollution 137 (2005) 316e323.Pages 316-323.
RODIER J. (1996) : L’analyse de l’eau, eaux naturelles, eaux résiduaires, eau de mer. 8e
édition Dunod, Paris.
U.S. Environmental Protection Agency (1987): Health effects assessment for boron and
compounds. Environmental Criteria and Assessment Office.8 pages.
VILA J.M. (1973) : Notice explicative de la carte géologique, au 1/50.000, d’Ain Beida
(feuille n° 149).Editée par Sonatrach en 1977.6 pages.
VILA J.M. (1973) : Notice explicative de la carte géologique, au 1/50.000, de Dalâa (204).
Editée par Sonatrach en 1977. 6pages
VILA J.M. (1973) : La chaîne alpine d’Algérie orientale et des confins algéro-tunisiens.
Thèse de doctorat Es Sciences Naturelles. Université Paris VI. 665 pages.
Année sept oct nov déc jan fev mars avri mai juin juil aout
70/71 117,8 49,8 8,9 10,8 48,7 10,4 31,3 45,7 20,3 95,1 7,5 8,7
71/72 18,8 58,5 5,2 31 51,1 51 102 51,6 25 1,6 0 21,3
72/73 36,9 91,3 9 20,1 30,8 32,7 122,1 15,2 37,4 44,9 0,3 30,4
73/74 10,3 10,9 6 85,4 12,5 48,3 36,2 50,2 30,4 14,8 8,5 22,1
74/75 18,5 28,9 13,5 19,2 13,4 54,8 31,4 27,6 55,6 0 25,4 20,7
75/76 23,1 11 37,3 6,2 14,7 30,2 41,1 22,4 31,8 43,5 27,3 33,3
76/77 23,1 20,1 119,5 8,3 10,8 5,6 40,1 30,4 33,2 8,7 15 14,8
77/78 11,2 3,3 36,7 3,9 3,9 38,7 77,5 23 19,9 23,9 0 32,7
78/79 25,4 23 16,4 3,6 17,5 34,6 61,3 59,4 46,7 17,7 0 9,9
79/80 98,1 10,5 31,3 31,7 16,7 29,8 49,8 21,1 33 27,3 1 10,8
80/81 49,8 14 16,8 26,5 14,2 22,2 15,3 41,7 33,6 49,4 6,7 4,5
81/82 42,7 44,2 43,8 15,6 18,4 35,6 22,4 36,2 59,21 14,8 4,8 26,7
82/83 7 33,2 46,3 12,8 42,2 48,6 22,7 17,6 13,6 17,9 7,8 28,3
83/84 11,4 32,7 4,3 32,9 13,7 41,4 41,3 19 21 6,8 0 15,4
84/85 23,7 21,2 16,8 43 25,7 15,9 64,5 15,4 57,4 22,2 15,9 6
85/86 66,55 12,6 16,75 36,4 10,5 19 60,9 2,5 26,8 13,2 36,5 23
86/87 15,83 25,4 36,6 16,7 10,2 17,4 42,6 19,2 33,3 34,2 13,7 5
87/88 14,44 22,6 30,2 27,5 12,7 2 23,8 27,6 56,1 41,4 11,5 36,5
88/89 18 13,5 23,9 25,4 43,3 10,4 34,7 26,3 19 45,3 8,7 101,3
89/90 38,6 7,8 31,6 18,7 43 19,5 40,6 50,8 48,35 17,1 13,2 115,6
90/91 37,5 29,4 60,1 34,5 32,4 31,3 44,3 25 54,9 14,4 8,5 35,6
91/92 58,9 52,4 39,6 14,2 20,3 24,9 22,3 25,6 71,7 23,2 13,4 24,5
92/93 44,08 28,4 61,05 48,4 5,3 27,9 19,4 14,5 22,36 12,8 7,1 2,3
93/94 16,05 26,2 9,9 19,5 44,2 20,9 9,4 13,3 31,7 50,6 14,5 6,7
94/95 82 36,8 7,4 5,7 88 23,5 52,4 30,2 12,7 15,8 7,6 3,5
95/96 94,9 33,7 28 8,2 10,5 60,9 39,3 30,5 13,83 27,6 7,5 15,5
96/97 56 5 42,6 22,9 27,6 4 27,8 50,7 12,06 12 0 19,4
97/98 32,1 40,5 26,7 13,5 51,5 2 15,7 9,9 18,51 36,3 13 13,3
98/99 42,1 51,7 36,3 27,9 23,7 11,7 31,8 15,4 80 22 5 13,6
99/00 20,3 48,7 39,5 23,2 20 14,8 6,1 12 50,96 36,5 10 8,5
00/01 41,2 32,7 13,5 13,9 22,8 16 5,5 11,1 31,6 28 6 46,7
2001/2002 35,84 6 56 6,8 12,3 16 33,5 26 11,44 5 15 85
2002/2003 32,6 77,9 59,3 36,5 53,5 25,5 31 64,5 83,7 30,9 16,2 36,8
2003/2004 48,8 35,4 80,77 122,5 12,5 13,4 41,8 13,7 20,8 4,5 6,3 13,5
2004/2005 46,54 18,3 4,24 31,8 32,1 15 35,5 55,5 23,9 21 9,5 46,6
Les précipitations mensuelles enregistrées à la station de Tébessa :
Année sept oct nov déc jan fev mars avri mai jui juil aout
72/73 36,9 99,3 11,2 34,9 46 42,7 171,1 31,3 44,7 65,5 5,3 36,4
73/74 16,3 12,9 6 94,4 14 28,3 29,2 50,2 10,4 24,8 4,5 12,1
74/75 27,1 37,9 18,5 14,2 23,4 67,8 33,6 21,6 66,6 0 25,4 23,7
75/76 26,1 11 47,3 6,2 22,4 38,2 49,1 32,4 35,8 59 27,3 39,3
76/77 26,1 23,1 134,5 10,3 14,7 6,6 45,1 40,4 38,2 9,1 15 19,4
77/78 11,2 3,3 46,7 3,9 3,9 54,7 102,5 23 23,9 3,9 0 50,1
78/79 5,4 26 20,4 3,6 10,3 44,6 40,3 89,4 22,7 27,7 0 11,7
79/80 116 18,5 21,3 1,7 33,7 29,8 76,8 28,1 41 4,3 0,2 3,4
80/81 65,8 3,7 24,1 47,5 13,4 18,8 24,1 11,7 35,8 72,4 3,6 4,1
81/82 37,3 23 1,9 15,3 21,8 45,6 12,4 56,2 80,1 8,5 3,7 15,5
82/83 12 58,5 50,3 24,7 2,8 7,3 18,1 5,7 30,4 42,7 0,7 31,5
83/84 3,9 31,7 17,9 12,2 18,9 92,4 24 24,1 4,3 6,8 0,2 15,4
84/85 27,2 26,2 19,1 51 25,7 11,3 54,5 26,4 65,2 27,2 2,4 6
85/86 50,8 23,1 3,5 13,5 31,1 14,3 83,1 2,5 35,8 15,2 51 13,1
86/87 24,4 28,7 44,7 20,7 10,2 27,4 62,6 13,2 25,1 4,2 33,7 5
87/88 15,5 18,7 33,8 9,2 23,7 4,2 35,8 31,6 55,6 62,1 8,3 6,5
88/89 21,1 20,6 35,1 35,4 18,3 17,4 14 16,3 8,4 57,3 8,7 99,3
89/90 44,6 12 10,8 8,7 83 0,2 34,8 43,1 66,9 17,1 15,2 137
90/91 53,3 22,4 99,8 64,9 30,3 12,8 54 43 67,8 14,4 6,4 65,6
91/92 74,7 34,4 44,3 14,2 34 29,9 24,3 43,6 82 23,2 13,4 4,5
92/93 51,2 28,4 61,6 48,4 9,3 27,9 21,4 2,6 31,1 12,8 20,1 1,8
93/94 22,7 3,8 16,8 28,7 31 23,9 19,4 23,3 41 2,4 4,5 11
94/95 7,2 66,8 0,6 6,8 24,7 3 32,2 22,1 7,4 37,9 1,7 44,1
95/96 150 39,7 26,6 18,2 24,9 72,9 56,3 49,8 30,2 38,9 13,2 30
96/97 12,4 4,1 1,2 15 31,6 7,1 18,9 46,8 16,1 10,3 20,2 23,7
97/98 64 72,5 45,2 21,5 22,3 10,2 28,7 29,2 16,7 31 0 15,1
98/99 78,6 36,2 55,1 14,5 56,4 11,7 45,6 15,4 30,9 16,9 18,9 33,7
99/00 22,1 81,5 64,5 34,5 3,2 4,4 10,4 15 87 76 21 18,9
00/01 50,5 18 16,5 13,9 26,8 16 14,8 3 48,8 2 8 1
2001/2002 54,5 10,9 23,3 7,1 16,5 11,5 4,5 29,4 41 13 58 85
2002/2003 36,9 37,5 76,9 30 100 38,9 17,5 97,5 28,6 9,4 3 12
2003/2004 70,4 45,2 17,5 168,4 20 3,4 72,8 29 38,6 16 43,7 19,4
2004/2005 19,4 25,5 11 66,5 28,6 33,6 23,5 20 1 32 1 46,6
Les températures mensuelles enregistrées à la station de Tébessa :
Année sept oct nov déc jan fev mars avri mai jui juil aout
72/73 19,8 13,7 10,7 5,7 4,4 4,8 5,6 10,5 19,9 23,1 26,9 23,6
73/74 21,7 17,1 9,3 6,3 6,2 6,4 9,8 10,8 15,6 23,5 23,9 24,1
74/75 20,8 13,3 9,1 5,9 5,9 6 8,5 11,9 16,4 21,3 25,4 23,2
75/76 22,7 14,2 9,2 7,4 5,2 6,6 7,4 11,8 16 20 23,6 23,8
76/77 19,9 15,4 7,8 8,2 8 9,8 11,8 12,9 17,5 21,8 27,8 24,6
77/78 19,8 16,6 11 7,9 5,1 9,8 9,2 12,7 17 23,1 25,5 24,6
78/79 20,4 12,6 7,9 9,6 9,8 8,6 10,5 10,2 16,7 22 26,1 25,2
79/80 18,8 17,6 8,4 7,4 5,8 7,3 9 10,3 14,9 22,4 24,9 25,7
80/81 21,7 14,2 11 4,5 3,9 6,3 12,4 15,2 19 23,3 24,5 24,1
81/82 20,7 17,3 9,2 9,7 7,7 7,4 9,2 11,8 17 24,2 28 25,9
82/83 21,5 15,9 10,8 5,7 4,7 6,7 9,3 15,3 18,9 22,5 27,1 26
83/84 21,6 15 12,2 7,2 6,1 6 8,5 13 17 23,1 26,1 25,2
84/85 20,6 14,5 12 6,4 5,4 10,4 8,3 14,2 16,9 25,1 27,7 25,6
85/86 20,4 15,7 13 8 6,3 8,1 9,4 13,7 19,9 22 25 26,8
86/87 21 16,6 10,2 6,4 6,3 7,7 8,8 14,3 16,7 24,3 26,7 27,9
87/88 23,6 19,3 10,7 10,5 8,2 7,4 9,9 14,8 20,5 22,5 28 26,7
88/89 20,6 18,2 11,6 5,8 5,4 7,2 11,5 13,7 18,2 20,9 25,6 26,1
89/90 22,2 15,5 13,4 10,7 6,1 10,4 10,4 12,7 17,2 25,1 24,9 22,5
90/91 24,1 20 11,5 5,5 5,5 6,8 11,6 10,4 14,2 21,9 26,3 25,6
91/92 21,8 16,7 10,6 5,5 4,7 6,7 9,4 11,8 16,3 20,9 23,9 25,7
92/93 21,9 18 12,1 7,8 5,2 5,6 8,7 13,9 19,2 24,8 26,8 27
93/94 22,3 19,1 11,1 7,7 7,5 8,9 11,9 11,8 21,9 24,2 27 28,6
94/95 23,6 16,7 13,2 8,1 5,7 10,3 9,2 12,7 20,1 22,9 27,1 24,6
95/96 21,1 16,3 11,3 9,8 9 6,3 10,1 12,4 18,2 20,8 25,9 26,6
96/97 20,3 15 12,4 10,2 8,7 9,3 9,3 11,2 20,4 26,6 27,5 25,2
97/98 20,5 17 11,8 8,5 7,2 8,2 9,8 15,1 17,7 24,6 27,8 25,7
98/99 23,2 15 10,2 6,3 7,1 5,8 10,2 14,9 22,1 25,8 26,2 28,9
99/00 23,6 19,2 11,1 7,1 3,8 7,2 11,2 15,8 20,8 20 27 26,4
00/01 21,9 15,6 12,2 9,2 8 7,8 15,2 14,5 19,2 25,5 28,4 27,5
2001/2002 22,9 21,5 11,2 6,2 6,2 8,5 12,6 14,6 19,6 25,6 26,8 24,2
2002/2003 21,6 17,2 11,2 8,4 6,4 6,2 10,6 14,1 18,5 24,4 29,8 27,6
2003/2004 21,9 19,4 12,5 7,2 6,8 9,6 11,2 13 15,5 22,4 26,2 27
2004/2005 20,2 20,2 10,6 8,5 4,7 4,9 11,5 14,6 21,5 23,9 28,5 25,8
Les précipitations interannuelles cumulées aux stations de
Tébessa et Meskiana.
s=
0.183.Q 2.25T.t
log
T R 2 .S
où :
s: rabattements en (m) ;
T : transmissivité en (m2/s) ;
Q : débit de pompage en (m3/s) ;
S : coefficient d’emmagasinement ;
t : temps écoulé depuis le début de pompage en (secondes) ;
R : distance entre le piézomètre et le forage en (m).
T=
0.183.Q
C
où
T : transmissivité en (m2/s)
Q : débit de pompage (m3/s)
C : pente de la droite.
0,5
0,4
Rabattement (m)
0,3
0,2
0,1
T = 0,63.10-3 m2/s
0
1 10 100 1000
Temps (s)
Courbe de remontée du puits n° 12
0,6
0
1 10 100 1000
Temps (s)
0,6
0,5
Rabattement (m)
0,4
0,3
0,2
0,1 T = 1,66.10-3 m2/s
0
1 10 100 1000
Temps (s)
0,6
0,5
Rabattement (m)
0,4
0,3
0,2
0,1 T = 1,45.10-3 m2/s
0
1 10 100 1000
Temps (s)
Courbe de remontée du puits n° 38.
0,5
Rabbatement (m)
0,4
0,3
0,2
0,8
Rabattement (m)
0,6
0,4
0,2
T = 5,4.10-3 m2/s
0
1 10 100
Temps (s)
0,5
0,4
Rabattement (m)
0,3
0,2
0,1
T = 3,33.10-3 m2/s
0
1 10 100
Temps (s)
Courbe de remontée du puits n° 50
0,6
0,5
Rabattement (m)
0,4
0,3
0,2
0,1 T = 1,66.10-3 m2/s
0
1 10 100 1000
Temps (s)
0,6
0,5
Rabattement (m)
0,4
0,3
0,2
T = 1,2.10-3 m2/s
0,1
0
1 10 100 1000
Temps (s)
0,5
0,4
Rabattement (m)
0,3
0,2
0,1
T = 0,92.10-3 m2/s
0
1 10 100 1000
Temps (s)
ANNEXE IV
(HYDROCHIMIE)
Stratégie d’échantillonnage :
Les échantillons ont été prélevés en utilisant un échantillonneur dont le câble dépasse
200 m, puis placés dans des flacons de plastique de 1litre, préalablement rincés à l’eau de
l’échantillon. Ces derniers n’ont pas été acidifiés car leur analyse se fait dans l’immédiat au
laboratoire de CERAD. Les échantillons sont conservés dans une glacière portative pour
garder une température inférieure à 4°C, jusqu’à l’arrivée au laboratoire. On note que les
échantillons analysés au laboratoire central de chimie, université de Damas ont été acidifiés à
un PH inférieur à 2 avec de l’acide nitrique.
Mesures in situ :
Des paramètres physico-chimiques ont fait l’objet de mesure sur le terrain, en utilisant
une valise multi paramètre de marque Multi340i (WTW), tels que le PH, la conductivité
électrique et l’oxygène dissous.
Analyses aux laboratoires :
On a procédé à une analyse complète de tous les éléments majeurs (Ca2+, Mg2+, Na+,
K+, Cl-, SO42- et HCO3-).Les métaux lourds :Fe3+, Mn2+, Al3+, Zn2+, B ,Sr2+ ,Zn2+, Cd2+ et Cu2+
ont fait l’objet d’analyse pour quelques campagnes seulement. Les polluants résultants des
activités agricoles ont été le sujet d’un suivi particulier (NO3-, NO2-, NH4+ et PO43-).
Les méthodes analytiques utilisées sont :
-Méthode colorimétrique pour les éléments (Cl-, SO42-, HCO3-, NO3-, NO2-, NH4+ ,PO43- ,F et
Fe total ).
-Méthode de spectrophotométrie d’absorption atomique : (Ca2+, Mg2+, Na+, K+, Sr2+, Mn2+,
Al3+, Zn2+, Cd2+ et Cu2+).On signale que l’analyse spectrophotométrique a été effectuée, au
laboratoire du CERAD en utilisant un spectrophotomètre DR 2000 (HACH) et au laboratoire
centrale de chimie à l’université de Damas (Syrie) en utilisant l’appareil Varian Spectra 220.
Quand on a une réaction chimique quelconque l’équilibre de cette dernière se réalise selon la
constante K qui est calculée par la relation :
⎛ ΔG r0 ⎞
K = exp ⎜⎜ − ⎟⎟
⎝ RT⎠
La variation d’énergie libre ΔG° est la différence entre la somme des énergies des éléments
qui sont produits par la réaction moins la somme des énergies des éléments entrant dans la
réaction.
ΔG° = ΣG (produits)−ΣG(réaction)
Log K =
a CaSO4
Ca2+
a .a
SO42-
A zi μ
Log γ i = − + bi μ
2
1 + B ai μ
0
μ= ∑
1
zi 2 m i est la force ionique de la solution (<1)
2 i
zi est la charge de l’ion i ;
A et B sont constants et ne dépendent que de la température ;
ai0 est un paramètre représentant la taille de l’ion en solution ;
bi est une constante qui représente l’écart à la loi de debye-Hückel.
Etat de saturation :
Ensuite, armé des valeurs des coefficients d’activités, on peut calculer les activités des
espèces i puis l’état (ou l’indice) de saturation (IS) de la solution vis-à-vis d’un minéral. C’est
le rapport :
IS = Log Q/K
Carte des nitrates (mg/l), nappe de Meskiana mars 2003. Carte des nitrates (mg/l), nappe de Meskiana juillet 2003.
930 935 940 945 950 955 960 930 935 940 945 950 955 960
290 290
90
100
285 85 285
95
80
Rehia Rehia 90
280 Bellala 280 Bellala
75
85
275 70 275 80
65
75
270 Meskiana 270 Meskiana
60
70
55 65
265 265
50 60
260 45 260 55
Zebar Zebar
40 50
255 255
35 45
250 30 250 40
25 35
Dalâa 245 Dalâa
245 30
20
25
240 0 10 20Km 240 0 10 20Km
Carte des nitrates (mg/l), nappe de Meskiana mars 2004. Carte des nitrates (mg/l) , nappe de Meskiana juillet 2004.
930 935 940 945 950 955 960 930 935 940 945 950 955 960
290 290
95
285 90 285 90
85 85
Rehia Rehia
280 Bellala 280 Bellala
80 80
75 75
275 275
70 70
270 Meskiana 65 270 Meskiana 65
60 60
265 265
55 55
50 50
260 260
Zebar 45 Zebar 45
255 40 255 40
35 35
250 250
30 30
Dalâa
Dalâa 25 245 25
245
20 20
240 0 10 20Km 15 240
0 10 20 Km
Carte des nitrates (mg/l), nappe de Meskiana mars 2002. Carte des nitrates (mg/l), nappe de Meskiana juillet 2002.
Carte des nitrates (mg/l), nappe de Meskiana décembre Carte des nitrates (mg/l), nappe de Meskiana novembre 2005.
2004.