1°) Choix de la zone d'étude et cartographie des occupations du sol
Bassin versant choisi : Andromba (région Analamanga).
Données : MNT ASTER (30 m), image Google Earth géoréférencée (ou via SASPlanet).
Méthode :
1. Capturer l'image de la zone (Google Earth Pro) et exporter au format GeoTIFF.
2. Charger l'image et les orthophotos dans QGIS.
3. Effectuer une classification supervisée (maximum de vraisemblance) pour distinguer
forêts, rizières, zones bâties, eaux, savanes, sols nus.
4. Générer un shapefile d'occupation du sol en 6 classes principales.
2°) Délimitation du bassin versant et analyse hydrophysique
Délimitation :
1. Utilisez l'outil « D8 Flow Direction » sous QGIS pour calculer la direction d'écoulement à
partir du MNT.
2. Appliquer « Basin » pour extraire le bassin versant de la cellule de rejet principale.
Caractéristiques physiques : superficie, altitude min/max, pente moyenne.
Combinaison avec l'occupation du sol :
o Calculateur pour chaque classe le pourcentage de couverture (ex. 30 % rizières, 25 %
forêts, etc.).
o Inférer la capacité d'infiltration (forêt élevée, rizière modérée, sol nu faible).
Informations tirées : disponibilité potentielle de l'eau, zones à risque d'érosion, corridors
hydriques.
3°) Évaluation socio-économique
Réseau hydrographique : longueur totale de cours d'eau, présence de plans d'eau.
Superficie cultivée : extraites du shapefile d'occupation du sol (ex. 150 ha de rizières).
Démographie :
o Nombre de ménages fixé à 5 (consigne).
o Population estimée = 5 ménages × 5 personnes/foyer = 25 hab.
Activités économiques : agriculture (rizières, cultures de rente), pêche en plans d'eau, petits
élevages.
4°) Évaluation de l'érosion hydrique (RUSLE)
Choix de l'année d'observation : 1985 (pluviométrie moyenne annuelle P = 1 200 mm hypothétique ).
a) Facteur R (érosivité de la pluie)
Formule Nguyen (1996) :
R=0,548257×P−59,9=0,548257×1200−59,9≈603 MJ mm/hun h unn R = 0{,}548257 \fois P - 59{,}9 =
0{,}548257 \fois 1200 - 59{,}9 \environ 603\,\mathrm{MJ\,mm/ha\;h\;an}R=0 , 548257×P−59 , 9=0 ,
548257×1200−59 , 9≈603MJmm/hahun
Formule Kassam (1992) (si P < 2000 mm) :
R=117,6×1 00105P R = 117{,}6 \fois 1{,}00105^PR=117 , 6×1 , 0010 5P
b) Facteur K (érodibilité du sol)
Extraire le shapefile pédologique et associer à chaque texture la valeur K correspondante (base
de données pédologiques locales).
c) Facteur topographique LS
1. Carte de pente en degré (θ) via l'outil « Slope » sur le MNT.
2. Calcul de β et m pour chaque cellule :
β=péchéθ0,0896+3(péchéθ)0,8+0,56,m=β1+β \beta = \frac{\sin\theta}{0{,}0896 + 3(\sin\theta)^{0.8} +
0{,}56}, \quad m = \frac{\beta}{1+\beta}β=0 , 0896+3 ( péchéθ )0,8+0 , 56péchéθ,m=1+ββ
3. Accumulation de flux (outil « Flow accumulation ») → carte A(i,j).
4. Calcul de L pour chaque cellule :
L=(UN+D2)m+1−UNm+1m Dm+2 22,13m L = \frac{(A + D^2)^{m+1} - A^{m+1}}
{m\,D^{m+2}\,22{,}13^m}L=mDm + 222 , 1 3m( UN+D2)m + 1−UNm + 1
5. Carte de tan θ puis S via
S={10,8 péchéθ+0,03,tannéθ<9%(péchéθ/péché5 143°)0,6,tannéθ≥9% S = \begin{cases} 10{,}8\,\sin\
theta + 0{,}03, & \tan\theta < 9\%\\ (\sin\theta / \sin5{,}143°)^{0{,}6}, & \tan\theta \ge 9\% \
end{cases}S={10 , 8péchéθ+0 , 03 ,( péchéθ /péché5 , 143° )0 , 6,tannéθ<9%tannéθ≥9%
d) Facteur C (occupation du sol)
Reclasser le shapefile d'occupation du sol en valeurs C d'après le tableau fourni (ex. forêt = 0,005;
rizière = 0,05; sols nus = 0,75; etc.).
e) Facteur P (conservation de l'eau)
1. Calcul de la pente en pourcentage (%) via l'outil « Slope (%) ».
2. Reclassement selon les classes de pente et d'occupation du sol : intersection des deux couches →
attribuer P-factors du tableau.
f) Carte finale de la perte en sol (A)
Appliquer la formule RUSLE :
UN=R×K×L×S×C×P A = R \times K \times L \times S \times C \times PUN=R×K×L×S×C×P
Générer une carte raster de A (t/ha·an), classée en 5 classes par méthode Natural Breaks (Jenks).
Remarque : toutes les cartes intermédiaires et finales doivent être exportées en haute résolution et
intégrées dans un rapport Word, avec mention des auteurs, méthodes et sources de données.