Cartographie des sols de Tafrata, Maroc
Cartographie des sols de Tafrata, Maroc
Résumé
L’objectif principal de cette étude était d’obtenir d’une carte numérique des sols en exploitant l’information de systèmes de
classification des images satellitaires et des bases de données géologiques et des variables environnementales qui expriment les
relations sol-paysage. Cette étude a été réalisée dans la zone située entre la plaine de Tafrata et la boutonnière de Debdou dans
la zone Centro-oriental du Maroc. Un levé pédologique a été réalisé dans la région d’étude, notamment l’étude d’aménagement
de la forêt de Debdou et Lamkam au sein de laquelle cette étude a adopté un échantillonnage stratifié fondé sur un total de 210
profils (0,50 m x 0,50 m de profondeur jusqu’à la roche mère). On a tenu compte du rapport avec les unités cartographiques
identifiées avec les variables environnementales; élévation, pente, géologie et surfaces géomorphologiques et une classification
supervisée par le SIG a été implémentée. La carte finale est une base fondamentale pour visualiser les indicateurs de désertification
attachés à la qualité de sol, selon les critères de la méthodologie de quantification de la désertification MEDALUS (Mediter-
ranean Desertification and Land Use). Le résultat obtenu est une carte thématique comportant des classes et des informations
semi détaillées des sols de la région.
Mots clés: carte du sol, levé pédologique, SIG, désertification, plaine de Tafrata, Debdou, Maroc Centro-oriental
Mapping of soils of Tafrata plain (Morocco) and evaluation of their sensitivity to desertification
Abstract
The main objective of this study was to obtain a digital soil map by exploiting information from satellite image classification
systems, geological databases and environmental variables that express soil-landscape relationships. This study was carried out on
an area located between the plain of Tafrata and the buttonhole of Debdou in the Centro-oriental zone of Morocco. A soil survey
was conducted in the study area, including the Debdou and Lamkam forest management study, in which this study adopted a
stratified sampling based on a total of 210 profiles (0.50 m x 0.50 m deep to the bedrock). We took into account the relationship
with the map units identified with the environmental variables; elevation, slope, geology and geomorphological surfaces and
a supervised classification by GIS has been implemented. The final map provided us with a fundamental basis for visualizing
the indicators of desertification attached to soil quality, according to criteria of the desertification quantification methodology of
MEDALUS (Mediterranean Desertification and Land Use). The obtained result is a thematic map with classes and semi-detailed
information of soils of the region.
Keywords: Soil map, Soil survey, GIS, Desertification, Tafrata plain, Debdou, Morocco Centro-oriental
1
Faculté des lettres et des sciences humaines Sais, Université Sidi Mohammed Ben Abdellah, Fès, Maroc
©
Revue Marocaine des Sciences Agronomiques et Vétérinaires [Link]
372 Seghir et al.: Cartographie des sols de la plaine de Tafrata au Maroc
La région objet de cette étude présente des caractères pa- Pendant le Trias et le Jurassique inférieur, les sillons
léogéographies bien distingués. Elle est située entre deux intra-cratoniques des Hauts et du Moyen Atlas ont été
formations géologiques qui sont le bassin de Guercif au individualisés entraînant la formation d’un système de
N-NO et la chaîne de Debdou au S-SE culminant à 1670m horsts et de grabens et/ou de demi-grabens dans ces zones.
d’altitude. Cette chaîne perchée dans une zone enclavée et Les dépôts mésozoïques et cénozoïques ont été assujettis
escarpée, comprend une structure géologique composée à deux événements compressifs consécutifs: un événe-
d’un socle rigide remontant au primaire, représenté par des ment atlasique (Eocène-Crétacé) qui était à l’origine de
schistes violacés ou saumâtres et des bancs de quartzite la formation des rides et des dépocentres et un événement
affectés d’accidents WNW-ESE. alpin, post-miocène, responsable des plissements affectant
les séries mésozoïques et cénozoïques (Michard, 1976).
Le bassin de Guercif est situé dans le couloir Taza-Oujda.
Il est représenté comme une vaste dépression dont la Les conditions climatiques sont celles du climat méditer-
superficie est de 6230 km2, remblayé par d’épaisses séries ranéen, étage bioclimatique semi-aride à variante tempé-
néogènes atteignant jusqu’à 1800 m d’épaisseur (Wernli, rée avec une tendance continentale et sèche en s’avançant
1988), masquées le plus souvent par des formations conti- vers l’intérieur du pays.
nentales [Link] est constitué de sédiments du
Les précipitations moyennes annuelles de l’ordre de 300
Trias, du Jurassique du Néogène et du Quaternaire.
mm/an sur le massif de Debdou, et 240 mm dans la plaine
Les études géologiques montrent que le bassin de Guer- de Tafrata. Le régime pluviométrique est du type AHPE
cif correspond à une dépression, développée durant le avec une période sèche s’étalant sur 6 mois, du mois de
Jurassique et le Tertiaire. Il correspond à la continuité mai à celui d’octobre. La végétation naturelle se caracté-
vers le Nord-Est du sillon moyen-atlasique avec un épais- rise par une faible densité presque dans toute la plaine et
sissement important de sédiments (calcaires, dolomites, un recouvrement moyen dans le massif de Debdou, qui
marnes et grès) dans sa partie occidentale au moment où sa est dominée par le Steppe d’Alfa et Romarin et le chêne
partie orientale correspondait à une plate-forme soumise à vert aussi le Thuya et Genévrier oxycédre.
d’importants dépôts de calcaires. (Bernini M et al., 1999).
De point de vue litho-stratigraphique on trouve: Matériel et Méthode
• Le Paléozoïque forme le socle rigide de la Gada de À l’exception de quelques travaux localisés et assez som-
Debdou. Il est représenté par des schistes violacés ou maires; la région ne connaît pas d’études pédologiques
brunâtres et des bancs de quartzites, affectés d’accidents abondantes; sauf les études de Beaudet (1971); Ruellan
WNW-ESE. Il affleure au sud de la plaine de Tafrata et à (1970) et Laouina (1990) qui ont donné un aperçu sur
l’Est de la plaine de Mahrouf; un petit batholite de gra- les principales caractéristiques des sols dans la zone de
nite hercynien est venu se mettre en place dans ce massif Debdou.
primaire. Par ailleurs, le Sud du Jbel Mazgout, à l’W de Les caractéristiques géologiques et géomorphologiques
la plaine de Sangal, est formé également de schistes pri- de la zone étudiée ont permis aussi de distinguer des zones
maires injectés de granites. assez bien déterminées sur le plan pédologique. Les sols
• Les dépôts Triasiques sont constitués de conglomérats, en zone à déficit hydrique, sont les zones qui bénéficient
de grès et d’argiles, recouverts par le sel et les basaltes. d’un relief marqué par une vaste pénéplaine reposée sur
une cuvette, dont les sédiments sont formés par les couches
• Les formations jurassiques sont recouvertes, en discor- Néogène et Plio-quaternaire, notamment la plaine de
dance angulaire, par des faciès du Tertiaire, composés de Tafrata. Les sols situés en zone ayant un profit hydrique
conglomérats, de grès, de marnes, de gypse, de lignite et se reposent sur des reliefs tabulaires (Gaâda), surtout sur
de calcaires lacustres. (Michard, 1976). les bordures et le piedmont de la Gaâda dont les sédiments
sont formés par des couches schisteuses du paléozoïque.
La couverture pédologique du bassin a été étudiée à Les indices de qualité ont été estimés en utilisant les para-
diverses échelles en se basant premièrement sur l’utili- mètres et les formules suivantes:
sation des données de la télédétection pour une connais-
IQS (Indice de la qualité des sols) = (matériaux Parentaux
sance spatialisée de certains facteurs liés à la propriété de
* Texture * Profondeur du sol * Pente de terrain) ¼
l’état de surface, la topographie, la brillance du sol et la
couverture végétale. Ensuite, sur l’utilisation du système IQC (Indice de Qualité du Climat) = (précipitations
d’information géographique (SIG) pour des opérations annuelles * orientation du terrain* indice d’aridité bio-
d’analyse d’autres facteurs liés à la classification et la climatique) 1/3
détermination des types pédologiques basée sur la litho- IQV (Indice de qualité de la végétation) (risque d’incendie
logie, la carte d’occupation du sol et la carte des pentes. * Résistance à la sécheresse *protection de l’érosion *
L’utilisation des images satellitaires constitue une aide couvert végétal) ¼
précieuse, en raison de la masse d’informations qu’elles
apportent dans le travail d’interprétation des types des
sols. Trois méthodes de traitement ont été appliquées à
l’image, il s’agit notamment de la méthode de:
• Correction géométrique de l’image satellitaire Landsat
ETM 2000 en composition colorée 3.2.7 avec une résolu-
tion de 28,5 m à l’aide d’une carte topographique.
• Correction radiométrique pour régler les valeurs des
pixels.
• Classification supervisée à partir d’une composition
colorée, on l’appelle aussi une segmentation en télédé-
tection, qui se fait selon des paramètres clés: la forme, la
texture et finalement le calcul des indices comme celui
d’indice de luminance et l’indice de végétation.
Le système de classification utilisé est le système français
(CPCS: Commission de Pédologie et de Cartographie du
Sol, 1967) et celui de l’école américaine (Soil Survey
Staff, 1975).
Troisièmement, nous nous sommes appuyés sur les résul-
tats de l’analyse des caractéristiques physico-chimiques
des sols qui a été réalisé dans la région d’étude, notamment
l’étude d’aménagement de la forêt de Debdou et Lamkam,
au sein de laquelle cette étude a adopté un échantillonnage
stratifié fondé sur un total de 210 profils (0,50 m x 0,50 m
de profondeur jusqu’à la roche mère).
En ce qui concerne l’évaluation de la qualité des sols, nous
avons basé notre approche sur la méthodologie de MEDA-
LUS comme une référence de base adoptée pour identifier
les zones sujettes à la désertification. Les régions sensibles
à désertification sont identifiées par la combinaison de 4
index de qualité (Sol; climat; végétation; aménagements
existants). Les trois premiers indices de qualité fournissent Figure 2: Carte des formations géologiques du bassin étudié
une image de la situation de l’environnement tandis que
le dernier exprime une évaluation de la pression résultant
d’activités anthropiques.
La méthodologie est basée sur la classification de chaque
indice de la qualité obtenue à partir de la moyenne géo-
métrique des cinq indices de qualité des facteurs environ-
nementaux et anthropiques disponibles retenus.
Les paramètres disponibles sont quantifiés par rapport à
leur influence sur le processus de désertification. Cepen-
dant, ce modèle ne précise pas le nombre et le type de
classes, laissant la flexibilité nécessaire pour s’adapter à
la disponibilité des données.
Les notes attribuées aux différents paramètres varient
entre 1 (meilleure valeur) et 2 (plus faible valeur). Figure 3: Coupe géologique traversant la partie Sud du bas-
sin de Tafrata-Guercif de direction WNW-ESE (Medioni R. 1997)
374 Seghir et al.: Cartographie des sols de la plaine de Tafrata au Maroc
IQG (Indice de la qualité de gestion) (intensité d’utili- Ces divers types des sols de la zone étudiée en général et
sation des terres * qualité d’aménagement des zones) ½ dans la forêt de Debdou en particulier sont soumis à une
forte dégradation sous l’effet des facteurs divers (naturels
La moyenne géométrique de ces paramètres, se référant à
et anthropiques); en raison des conditions climatiques
chacune des quatre indices mentionnés ci-dessus, et clas-
défavorables, la nature du relief, les facteurs géomorpho-
sifiée selon des intervalles réguliers (c’est-à-dire la même
logiques et l’état de la végétation ainsi que la pression
gamme) présentés au tableau 1, qui montre les classes de
exercée sur les terres par le système d’exploitation, le taux
qualité sur une échelle vers le bas.
de matière organique contenue dans le sol, a notablement
Tableau 1: Classe de qualité baissé en produisant un tassement superficiel. Cette forme
Classe Rang de dégradation facilite l’action des facteurs physiques tels
Fort 1,00 - 1,33 que l’érosion.
Moyen 1,34 - 1,66 Dans la région étudiée, l’érosion agit de manière beaucoup
Faible 1,67 - 2,00 plus rapide et efficace sur les terres dénudées, particuliè-
rement celles qui ont connu un manque de couverture
végétale, surtout dans la plaine de Tafrata. En effet, La
L’indice global final de l’ESA est obtenu comme une
présence de roches plus ou moins dures avec un relief
moyenne géométrique des indices de qualité.
relativement plat dans son ensemble a bien atténué l’éro-
ESA = (IQS * IQC * IQV * IQM)1\4 sion différentielle, laquelle commence à agir, par contre,
Mais en termes de notre étude, nous allons concentrer sur les versants au-dessus de 1000 m et sur des pentes
uniquement sur l’extraction des indicateurs les plus impor- dépassant 25º, soumis à des processus de reptation et de
tants de la qualité des sols. décapage des sols par ruissellement diffus au cours des
périodes pluviales.
Résultat et discussion L’érosion éolienne aussi est un facteur non négligeable eu
égard à l’importance de la succession des années sèches.
Suite à ces constats et conformément à ce qui procède, les Cependant, comme il a été déjà expliqué il faut considérer
différents types de sols qui existent relèvent des classes de l’impact du vent Chergui qui souffle le long de l’année
sols peu évolués, des sols minéraux bruts, des rendzines avec des vitesses dramatiques sur des terres dépourvues
et des sols calcimagnésiques (Figure 4). de l’essentiel de la protection végétale. Cette érosion est
responsable de l’emportement des particules fines et des
Classe des sols peu évolués substances nutritives nécessaires à la bonne fertilité des
sols. Cette forme d’érosion, due à ce processus de vannage,
Il s’agit des sols développés sur les cônes de déjection et est aussi responsable de la couleur blanchâtre des sols au
sur des substrats dégradés suite au défrichement ou à relief niveau des zones proches des bas-fonds, située sur les
très pentu, ils apparaissent dans des positions où les roches terres de Tafrata.
sont difficilement altérables sur les pentes ou le long des
vallées présentant une texture hétérogène. Leurs horizons
plus évolués avec une texture argileuse à limono-argileuse
ont une structure particulaire à grumeleuse. Ils s’étendent
sur les vallées de plaine de Tafrata et les dépressions de
Debdou. On distingue alors plusieurs classes des sols:
Classe des sols calcimagnésiques
Ce sont des sols caractéristiques des affleurements rocheux
riches en couches marneuses tendres. Ils s’étendent sur
une grande partie de la région de Debdou (Lala mimouna,
Ain Serraq, Al Ateuf). Ils résultent de la nature de la roche
mère qui est principalement calcaire.
Des sols minéraux bruts (lithosols)
Ces sols couvrent les zones à substrat schisteux ou gréseux
de la zone. Ils se trouvent surtout au niveau des sommets
des collines, des endroits qui ont subi un fort impact érosif.
Des rendzines (groupe de la classe des calcimagné-
siques)
Les rendzines correspondent à des sols calcimagnésiques
peu développés qui reposent directement sur les Mollasses
et les substrats calcaires. Ces formations, essentiellement
rencontrées sur des pentes sont dotées d’une faible capacité
d’infiltration.
Figure 4: carte des différents types des sols dans le bassin
Rev. Mar. Sci. Agron. Vét. (2019) 7 (2):371-378 375
Les cartes d’indices de particularité des sols L’étude géologique et lithologique de la zone permet de
Les sols du bassin sont caractérisés par une grande variété, dégager les principaux faciès des formations lithologiques
allant de types pédologiques pauvres aux plus développés, du bassin selon les trois classes.
cela est dû aux différentes formations géo-lithologiques Les résultats de la classification montrent une dominance
et aux conditions climatiques allant des zones arides aux de la classe des matériaux cohérents qui couvre presque
zones humides. En outre, les caractéristiques naturelles la partie Est et Sud-est du bassin (Figure 5).
du sol sont probablement influencées par la culture inten- Les classes des matériaux moyennement cohérents et
sive pratiquée par les populations qui habitaient depuis tendres se localisent dans la majeure partie de la plaine
longtemps. de Tafrata.
Afin de mettre en place l’indice de qualité des sols (IQS), Indice de la texture des sols
les paramètres suivants ont été utilisés, en corrélation
La texture des sols est définie par la taille des particules
avec la texture, les matériaux parentaux, la profondeur et
qui le composent (graviers, sable, limons, argiles) ayant
la pente du sol.
des effets sur la rétention de l’eau. En effet, la dissociation
Cartes et indice d’homogénéité et de matériaux est élevée dans les sols sableux que dans les sols argileux
parentaux et ce d’autant plus qu’ils sont moins humifères. La stabi-
lité structurale des argiles et les sols riches en humus leur
Les différentes formations utilisées pour classer cet indice confèrent une bonne résistance à l’érosion. L’aptitude au
sont regroupées dans le tableau 2. transport dépend essentiellement de la densité du matériau
Tableau 2: Normes des classes de l’indice des maté- (Latrille, 1979).
riaux parentaux Dans ce sens, les classes de texture utilisées selon l’ap-
Valeur proche préconisée par MEDALUS sont montrés dans le
Classe Description tableau 3.
indice
Matériaux Calcaires, dolomies, grés, dalle Tableau 3: les normes de classification de l’indice de
calcaire, granite, quartzite, 1 texture
cohérents basalte, conglomérats Valeur de
Classe Description
Matériaux Grés et basalte altérés, flyschs l’indice
Moyennement armés, croûte calcaire, calcaires 1,5 Texture équilibrée Bonne 1,00
cohérents granulaire compact Texture fine Moyenne 1,33
Matériaux Argile, marne, sable, formation Texture très fine Faible 1,66
superficielle, non ou peu 2 Grossière à très grossière Très faible 2,00
tendres encroutées, flysch, gypse, schistes
Cependant, en plus de la carte géologique, nous disposons
aussi des travaux d’aménagement de la forêt de Debdou
sur la texture pédologique du terroir; la représentation
cartographique de cet indice, selon la carte de la figure 6
a fait ressortir une dominance des sols à texture argileuse
à limono-argileuse avec des cailloux et graviers moyen-
nement abondants à abondants; ces sol ont en général une
structure particulière à grumeleuse fine, à très fine.
Indice de la pente du terrain
L’indice de pente a été calculé à partir d’un modèle numé-
rique d’élévation (DEM) des grilles de 10 m. La carte de la
pente du bassin versant représentée sur la figure 6, donne
la topographie générale du bassin. La classification adop-
tée suit les normes du modèle décrit par dans le tableau 4.
La lecture de la carte nous renseigne que 67,1 % du bassin
versant est occupé par des faibles pentes, ce qui présente
pour cette partie des risques quasi nuls en termes de perte
en sol par ravinement. Les pentes les plus raides (>30%)
Tableau 4: Normes des classes de l’indice de pente
Classe Pente Description Valeur
1 <5% Pente faible 1,00
2 5_15 % Pente moyenne 1,33
3 15_30 % Pente forte 1,66
4 > 30 % Pente très forte 2,00
Conclusion
D’après ce qui procède et suite à cette brève analyse, nous
pouvons distinguer dans notre zone d’étude cinq classes ou
types de sols (sols peu évolués, des sols minéraux bruts,
des rendzines et des sols calcimagnésiques), mais une
différence significative dans l’une des propriétés étudiées
peut être la base pour identifier les différentes classes de
sensibilité à la désertification.
Les informations contenues dans ces cartes semi-détail-
lées du sol peuvent être utilisées dans tous les domaines où
la composante du sol joue un rôle important dans la quan-
tification de l’érosion, la pollution du sol, la vocation pour
les cultures spécifiques, la vulnérabilité des aquifères, etc
de sorte qu’ils puissent être reconnus, protégés et promus.
Figure 8 : carte d’indice de profondeur du sol
Cependant, malgré la variété lithologique du territoire,
les sols dominants ont des caractéristiques clés assez
similaires en ce qui concerne la qualité et la capacité de
soutenir la croissance et l’entretien de la végétation, cela
est dû à l’absence de pluie persistante et l’ondulation du
terrain qui laissant un horizon de sol en particulier pauvre
et d’autre part à l’absence séculaire d’une tapée végétale
capable de stocker une quantité en matière organique.
Dans l’ensemble, et après l’observation des cartes réa-
lisées, on peut dire que la dégradation physique et bio-
logique dans cette zone aride provoque des impacts très
importants sur la société et sur la situation économique des
populations, ce qui nécessite l’intensification des efforts
pour s’adapter à ces conditions environnementales des Figure 9: sol dégradé dans la plaine de Maârouf (4 Mai 2017)
plus difficiles.
Bibliographie
Beaudet G. (1971). Le quaternaire marocain: état des
études. Rev. Géogr. Maroc, 20: 3-55.
Bernini M. Boccalettim., EL Mokhtari J., Gelati R., Morat-
tig., Papanig. (1999). Tectonics and sedimentation in the
Taza-Guercif basin, northern Morocco: implications for
the Neogene evolution of the Rif-Middle Atlas orogenic
system. Mar. Petr. Geol., 22: 115-128
El Harradji A. (1996). Structuration géologique et évo-
lution géomorphologique du Massif de Debdou (Maroc
Oriental). Rev. Géogr. Maroc, vol. XVI: 127-148.
Laaouina A. (1990). Le Maroc nord-oriental: reliefs,
modelés et dynamique du calcaire. Publ. Rectorat Univ.
Mohammed Ier, Oujda, 1. p: 605.
Latrille E. (1979). Cours sur l’érosion et les techniques de
lutte anti-érosives, 58 p.