la Modélisation Spatiale Appliquée à l’Environnement
I. Application du modèle EPM
Le modèle EPM
Le modèle EPM "Erosion Potential Method" de Gavrilovic a été conçu dans les années 50 par Gavrilovic pour
des bassins versants d'Ex-Yougoslavie. La méthode a été développée pour la prédiction des taux d'érosion
annuel des sols d'un ensemble de type d'érosion (érosion en nappes, érosion en rigoles et en ravins et sapement
de berges) pour la gestion des pratiques de lutte contre l’érosion. II se base sur la cartographie et la combinaison
de six paramètres qui sont : La pente, La sensibilité des sols a l'érosion, L’état érosif, la protection des sols,
La température et les précipitations
L'application du modèle de Gavrilovic nécessite la cartographie et l’intégration dans un SIG de tous les
paramètres nécessaires à l'application du modèle empirique et l'utilisation des techniques d'analyse spatiale pour
évaluer les pertes en sols et estimer le poids de chaque facteur et de leurs effets combines d'une part et de
démêler leur interdépendance d'autre part. Ces paramètres sont :
A- Les facteurs permanents de l'érosion hydrique sont :
- La sensibilité à l'érosion,
- La pente et l'état érosif.
B- Les facteurs dynamiques sont :
- Les précipitations,
- Les températures.
- L’utilisation des sols.
C- Les processus et formes de l'érosion réelle et leur hiérarchisation en fonction du degré du risque
La méthodologie
Le volume annuel moyen des sédiments érodés désigné par la lettre (W) est déterminé de la manière suivante :
Où :
W : volume annuel moyen des sédiments érodés (m3/km2/an)
T : Indice de température
Avec : 𝑡0 : Température moyenne annuelle en (C°)
H : Indice de Précipitation moyenne annuelle en (mm)
Z : Indice d’érosion (Erosion potentille)
𝑋𝑎 : Indice de protection du sol (Xa) (0.01=<Xa=<1)
Y : Indice de sensibilité des sols à l'érosion (0.1=<Y=<1)
𝜑 : Indice d’état erosive (0.1=< Φ =<1)
𝐽𝑎 : Indice de Pente en (%)
Les Indices du modèle EPM
Symboles Facteurs Source de données Nature de données Méthodes
Indice de Précipitation Stations météo (fichiers tabulaires Excel) Points tabulaires Interpolation spatiale (IDW, Kriging…)
H moyenne annuelle en
(mm) Images satilitaires (Worldclim, CHIRPS...) données Raster Interpolation spatiale (IDW, Kriging…)
Ja Indice de Pente Modéle numérique de terrain (MNT) données Raster Surface (Slope en %)
Indice de sensibilité des Cartes pédologiques Cartes+ Analyses physico-chimiques Calcule de l'Indice de sensibilité des sols
Y sols ou des roches à
l'érosion Cartes géologiques Cartes+ Analyses physico-chimiques Calcule de l'Indice de sensibilité des roches
Photos aériennes/images de google earth données Raster Dégitalisation des formes d'érosion
Phi Indice d’état érosif
Images satilitaires (Landsat, Sentinel…) données Raster Calcule de l'indice a partir des formules mathématiques
Photos aériennes/images de google earth données Raster Dégitalisation des types d'occupations des sols
Indice de protection du
Xa
sol
Images satilitaires (Landsat, Sentinel…) données Raster / RED, NIR… NDVI/Classification supervisée/SVM,RF…
Stations météo (fichiers tabulaires Excel) Points tabulaires Interpolation spatiale (IDW, Kriging…)
T Indice de température
Images satilitaires (Landsat, Sentinel…) données Raster / Bande thérmique DN → radiance → température en Kelvin/°C
1. Indice de Précipitation moyenne annuelle (H) en (mm)
le facteur (H) moyenne annuelle des précipitation est généralisé par le model d’analyse spatiale IDW ( Inverse
distance weighted) à partir des données des précipitations disponible (8stations) entre 1970 et 2015.
Implantation des stations pluviométriques
Spatial analyst tools – interpolation – IDW
2. Indice de Pente (Ja) en (%)
Les pentes constituent un paramètre très important dans le modelé EPM. L’augmentation de la vitesse
d’écoulement sous l’effet des pentes provoque fortement l’érosion.
Les pentes sont calculées à partir du fichier MNT.
Spatial analyst tools – surface –slope
Input: dem
Output: slope
Outup musearemnt : purcent
3. Indice de Sensibilité de sol à l’érosion (Y)
L’indice de sensibilité des sols à l’érosion est la réponse moyenne du sol à la capacité érosive des pluies
(Heusch, 1990). C'est un paramètre qui représente l’intégration de la valeur moyenne annuelle de la réaction de
l’horizon superficiel du sol aux processus érosifs, dus respectivement à l'impact des gouttes de pluies et au
ruissellement (Wischmeier & Smith, 1978). C’est également la facilité avec laquelle un sol se détache sous
l’énergie des gouttes de pluies ou du ruissellement, durant un événement pluvieux (Renard et al., 1997)
Plusieurs méthodes d’évaluation de la valeur du facteur d’érodibilité Y ont été mises au point. Cependant
l’approximation de la valeur du Y par l’équation de Wischmeier et Smith (1978) demeure la plus répandue à
travers le monde.
Il se détermine expérimentalement pour un échantillon de sol donné par la relation suivante :
M : (% sable fin + limon) (100 - % argile)
OM : % matière organique
S : code de structure de l’échantillon du sol
P : code de perméabilité de l’échantillon du sol
4. Indice de l’état érosif (𝜑)
L’indice (φ) de l’état érosif est calculé à partir des images satellitaires Landsat selon la formule de
(Milevski, 2008), en divisant la racine carrée de la bande rouge (RED) par la valeur maximale du
rayonnement (Qmax).
Ou
OLI4: La bande rouge (4) pour les images de landsat 8
Qmax: la valeur maximale du rayonnement pour la bande rouge à partir des métadonnées
(QUANTIZE_CAL_MAX_BAND_x, où x est le numéro de la bande)
5. Indice de protection du sol (Xa)
Pour bien évaluer le rôle de cet indice dans le risque d‘érosion, on s’est appuyé sur des images satellitaires
Landsat, ainsi que sur des critères identifiés par Gavrilovic en 1954 qui sont modifiés plus tard par un groupe de
chercheurs,
Pour ajuster les valeurs de l'indice de protection du sol, le facteur de couverture végétale NDVI a été effectué en
premier, à partir de l'image satellite Landsat, après les résultats obtenues ont été modifiés pour qu’ils soient
adaptés aux critères fixés par ZornM et Komac,B (2005), en limitant les valeurs positives de l'indice de
couverture végétale à 0,6 comme étant la valeur maximale, et qui représente les zones à forte densité de
couverture, et les valeurs négatives à -0,19 qui correspond aux Bad-lands et les terrains nus. Après avoir ajusté
le facteur de couverture végétal, le Xa a été calculé selon l'équation suivante
Xa = (XaNDVI - 0,61) * (-1,25)
Où, Xa : Indice de protection du sol
XaNDVI : Coefficient de couverture végétale modifié pour adaptation aux critères de l'indice de protection du sol.
Tableau 2. Type d’occupation et leurs valeurs selon Gavrilovic
Les valeurs du facteur Xa de la zone d’étude se situent
entre 0.1 pour les zones à forte densité de végétation
et 1 pour les zones des Badlands.
L’indice de végétation par différence normalisée (NDVI)
Permet de mieux différencier entre les composantes (Végétation, Sols nus, plan d’eau…)
Les valeurs du NDVI sont comprises entre -1 et +1, les valeurs négatives correspondant aux
surfaces d’eau. Pour les sols nus, le NDVI présente des valeurs proches de 0. Les formations
végétales quant à elles, ont des valeurs positives, généralement comprises entre 0,1 et 0,7 - les
valeurs les plus élevées correspondant aux couverts les plus denses. Pour calculer l’indice NDVI,
on utilise la bande rouge et proche infrarouge selon la formule suivante :
𝑁𝐼𝑅 − 𝑅𝐸𝐷
NDVI =
𝑁𝐼𝑅 + 𝑅𝐸𝐷
Avec : NIR : bande proche infrarouge et RED : bande rouge
𝑃𝐼𝑅 𝑏𝑎𝑛𝑑𝑒 5 −𝑅𝑂𝑈𝐺𝐸(𝑏𝑎𝑛𝑑𝑒4)
NDVI =
𝑃𝐼𝑅 𝑏𝑎𝑛𝑑𝑒 5 + 𝑅𝑂𝑈𝐺𝐸 𝑏𝑎𝑛𝑑𝑒 4
Avec : NIR : bande proche infrarougeet RED : bande rouge
7. Indice de température (T)
Dans cette étape, nous expliquons les étapes du processus d'estimation de la température de surface. Nous avons
besoins de la bande 10 de landsat 8 et le fichier texte MTL, les étapes de calcule sont comme suivant :
La température de surface terrestre (LST) récupérées par les capteurs infrarouges thermiques des satellites, le
processus de calcule est synthétisé en 2 étapes ci-dessous pour le capteur landsat 8:
a) Calcul de la Radiance spectrale à partir de la band 10
Radiance spectrale (L) = (ML * Qcal) + AL
D’où : ML est le facteur de redimensionnement multiplicatif spécifique à la bande à partir des métadonnées
(radiance_mult_band_x, où x est le numéro de la bande).
Qcal est correspond à la bande 10 ou/et 11, et AL le facteur de redimensionnement additif spécifique à la bande
à partir des métadonnées (RADIANCE_ADD_BAND_x, où x est le numéro de bande).
b) Conversion de la Radiance spectrale (L) en Température de brillance (BT)
BT en °C = (K2 / (ln ((K1 / L) + 1))) - 273,15
D’où : K1 est la constante de conversion thermique spécifique à la bande à partir des métadonnées
(K1_CONSTANT_BAND_x, où x est le numéro de la bande thermique). K2 est la constante de conversion
thermique spécifique à la bande à partir des métadonnées (K2_CONSTANT_BAND_x, où x est le numéro de la
bande thermique). L = TOA Par conséquent, pour obtenir les résultats en degrés Celsius, la température de
rayonnement est ajustée en ajoutant le zéro absolu (environ -273,15 ° C).
T= racine carré((BT/10)+0,1)
7. Erosion potentielle (Z)
Le facteur Z est l’un des indicateurs le plus important dans le module de Gavrilovic car un ensemble de
variables est utilisé pour calculer ce paramètre, il permet de suivre dans le temps et dans l’espace l’évolution de
l’érosion et les changements au niveau de l’occupation du sol.
Le Facteur Z réunit les facteurs ayant une incidence sur l’importance de l’érosion et se présente comme suit:
Les stades de l’érosion selon Z sont classés par (Gavrilovic Z. et al.2006) en cinq catégories (Tableau.1).
Tableau 1. Classe des états érosifs potentiels.
8. la perte en sol annuelle (W) en m3/km2/an
Le produit de la phase d’intégration des facteurs (H, T et Z) dans un SIG selon la formule de Gavrilovic est la
carte finale des volumes annuels des sédiments érodés,
T * H * 3.1415 * (SquareRoot(Power(Z",3)))
Les volumes annuels des sédiments érodés sont classés par (Gavrilovic Z. et al.2006) en six classes (Tableau.1).
Tableau 1. Classe des états érosifs potentiels.