Atelier régional sur l’ échange de méthodologie pour

« le développement des bassins versants cultivés

sur les collines et terrains montagneux »

Ruissellement et Érosion
sur les terres cultivées en pente
Cas des petits bassins versants
de la Dorsale Tunisienne

Hamadi Habaïeb Taoufik Hermassi

Professeur Attaché de recheches
À l’I.N.A.T. À l’I.N.R.G.R.E.F.

Tunis 28/02 au 01/03/2008

 Plan
1. Introduction

2. Bilan hydrologique

3. Ruissellement : Mesure et Formulation

4. Érosion : Mesure et Formulation

1- Introduction
* Politique en matière d’aménagement
rural :
- maîtrise complète des eaux de
surface,
- protection des espaces vulnérables
et des infrastructures hydrauliques
aval.

* Lacs collinaires : pierre maîtresse de

la stratégie nationale en matière de
conservation des eaux et des sols.
Les objectifs de la stratégie :
- réduire les pertes en terres agricoles
(estimées couvrir 10 000 ha par an),
- réduire l’envasement des barrages (25 Mm3
par an),
- augmenter la recharge des nappes,
- mobiliser la plus grande partie possible des
500 Mm3 actuellement perdus dans la mer ou
dans les Sebkhas,
- créer des points de développement des
cultures irriguées.
- Un lac collinaire est une Mille lacs
retenue créée par un petit
barrage en terre.
collinaires
en Tunisie
- Il contient de 10.000 à
300.000 de m3 d'eau recueillie
sur des bassins versants de
quelques ha à quelques km².
- Il s'intègre de façon naturelle
dans le paysage ne créant pas
de nuisance particulière.
• Protection des La construction d'un
infrastructures en aval lac collinaire vise
• Mise à disposition d'une les objectifs suivants
ressource en eau de
manière disséminée dans
le paysage
• Captage du
ruissellement et recharge
des nappes phréatiques
• Amélioration de
l'environnement
 Localisation de la stratégie

IRD
Avec la banque des données hydrologiques sur les
lacs collinaires de la Dorsale tunisienne, la DG ACTA
dispose d’un ensemble d’informations uniques au
monde sur l’hydrologie des petits barrages et de leurs
bassins versants (0,5 à 5 km²) dans un
environnement semi-aride à sub-humide,
caractéristique du monde méditerranéen.

Cette information a été abondamment utilisée par les

étudiants tunisiens et européens. Elle continuera à
alimenter encore longtemps la réflexion des
chercheurs et des ingénieurs chargés de l’intégration
de ces aménagements dans les paysages ruraux
tunisiens et méditerranéens.
Objectifs des études hydrologiques
des lacs collinaires

Caractériser le régime des apports superficiels et des

apports en sédiments dans les petites retenues,

Evaluer les pertes par érosion et connaître la durée de

vie des petites retenues,

Évaluer les pertes par infiltration pour la recharge des

nappes d’eaux souterraines et les divers usages des
eaux superficielles ainsi stockées.
 Méthodes
Utiliser l’équation du bilan d’eau des retenues
pour déterminer :
les apports par ruissellement,
les volumes utilisés,
les pertes par infiltration.

Mesurer les autres termes du bilan.

Estimation de l’érosion et de l’envasement des

retenues collinaires
2- Bilan hydrologique
 Schématisation du bilan hydrologique
d’une retenue (DCES/IRD, 1999)

Vu Vsup

V p lu

Vd Vev

Vecs
V vi
Vi

∆V= (Vr + Vecs+ Vp+Vf) – (Vev + Vd+Vvi+Vi+Vu)

IRD
 Les équations utilisées
Conservation des volumes d’eau
∆V = (Vr + Vecs + Vp + Vf) - (Vev + Vd + Vvi + Vi + Vu)
Dérivée simplifiée de cette équation pour le calcul
des débits entrants dans la retenue
Qe = -Qp + Qdev + Qvid + dV/dt
Qe : débit entrant
Qp : débit dérivé de la pluie tombée sur le lac
Qdev : débit déversé
Qvid : débit de vidange
dV/dt : débit dérivé de la variation de stock d’eau
 Exemple de crue reconstituée
Kamech : crue du 27 au 30/11/1999

80 0

70 20
Qentrant
60 40
Qsortant
50 Intensités 60

I (mm/h)
Q (m3/s)

40 80
28-nov 30-nov
30 100

20 120

10 140

0 160
10:00 22:00 10:00 22:00 10:00 22:00 10:00 22:00 10:00
Exemple de données enregistrées

Es Sénéga

1600 0
1400 20
Digue
1200 40
Hauteur (cm)

Pluie (mm)
1000 Seuil 60
800 80
600 100
400 Vase 120
200 140
0 160
sept-93 janv-95 mai-96 oct-97 févr-99

Hauteur d’eau, pluie et vase

Les données recueillies et traitées
Données géographiques, cartographiques
et caractéristiques des aménagements
Données hydro-pluviométriques
Données d’envasement, d’érosion
Données de qualité des eaux
Données agro-climatologiques
Données agro-économiques
3- Ruissellement :
Mesure et Formulation
 Introduction
Objectif

Lr = f(P,IK,IL….)
 Méthodologie -1-

Analyser les relations entre la pluie Lr

et la lame ruisselée correspondante

P
Hypothèse : linéarité des relations
Lr = f (P)
Lr = A (P-PL)
Avec Lr = 0 pour P < PL
 Méthodologie -2-

Analyser l’effet de l’humectation du sol sur P

le paramètre de rétention PL

Ik
Déterminer la pluie limite de ruissellement
à partir des averses qui ne ruissellent pas
PL = P0*e-α*IK
 Méthodologie -3-

Comprendre le rôle de l’intensité de la IL30

pluie dans la genèse du ruissellement

Ik
Déterminer l’intensité de pluie limite de
ruissellement au pas de temps de 30 mn IL30 = IL0*e-β*IK
 Etude du ruissellement
Bassin versant d’El Mouidhi
Relation entre la pluie et la lame ruisselée

Lr = 0,82(P-11,5)

Lr = 0,14(P-20)

SAIDI, I-2007
 Bassin versant de Dekikira

Relation entre la pluie et la lame ruisselée

Lr = 0,32 (P- 22)

Lr = 0,48(P- 4) après 1998
avant 1998

Lr = 0,3(P- 12)
avant 1998

SAIDI, I-2007
 Tableau 3 : Variation de la lame ruisselé en fonction de la pluie

Equation du Equation du
Bassin versant
ruissellement max ruissellement min
Lr = 0,82 (P-11,5) Lr = 0,14 (P-20)
El Mouidhi
IP>60 IP<20
Sadine1 Lr=0,78 (P-4) Lr=0,77 (P-15)

Fidh ben naceur Lr=0,56(P-7) Lr=0,19(P-7) P courte

Sadine2 Lr = 0,5(P - 2) Lr=0,45(P -17)

Lr= 0,48 (P-4)
Dekikira Lr=0,3(P- 12)
Lr = 32 (P- 22)

Es Senega Lr = 0, 4(P- 5) Lr= 0,1(P- 5)

SAIDI, I-2007
Etude du ruissellement
Relation entre Lr et P à l’échelle annuelle
 Bassin versant de Dekikira

Relation entre la pluie limite de ruissellement et l’indice d’humectation

PL= 12 e-0,05 Ik

SAIDI, I-2007
 Variation de La pluie limite de ruissellement
en fonction de la valeur de l’indice d’humectation

Bassin versant Equation de PL Equation de PL

Sadine1 Pour Ik< 20; PL=15 e-0,06 IK Pour Ik>20 ; PL= 3 mm

Sadine2 Pour Ik<22; PL= 17e-0,06 IK Pour Ik>22 ; PL=2mm

Fidh ben naceur PL= 7 mm

Es Senega PL= 5mm

Dekikira Pour Ik< 15 ; PL = 12e-0,05 IK Pour Ik>15 ; PL= 4 mm

Pour IP<20 ; PL= 20e-0,14 Ik Pour IP>60 ; PL= 11,5e-0,10 Ik
El Mouidhi
PL= 2 pour IK>15 mm PL= 2 pour IK>15 mm

SAIDI, I-2007
 Bassin versant de Dekikira
Relation entre l’intensité limite de ruissellement et l’indice d’humectation

IL30 = 12 e-0,05 Ik

SAIDI, I-2007
 Variation de l’intensité limite de ruissellement
en fonction de la valeur de l’indice d’humectation

Bassin versant Equation de IL30 Equation de IL30

Sadine1 Pour Ik< 20 ; IL30= 20 e-0,05 IK Pour Ik>20 ; IL30= 9mm/h

Sadine2 Pour Ik<20 ; IL30= 20 e-0,05 IK Pour Ik>20 ; IL30= 9mm/h

Fidh ben Naceur IL30=10mm/h

Es Senega IL30=15mm/h

Dekikira Pour Ik< 25 ; IL30 = 12 e-0,05 IK Pour Ik> 25 ; IL30 = 4mm/h

El Mouidhi Pour Ik<10 ; IL30 = 11,5 e-0,05 Ik Pour Ik>10 ; IL30= 5mm/h

SAIDI, I-2007
Bassin versant d’El Mouidhi
Très faibles travaux CES ce qui diminue la rétention des sols;
un fort pourcentage de parcours couplé à un problème de surpâturage
ce qui encourage l’érosion et favorise le ruissellement.
Bassin versant de Sadine 1
Des formations géologiques à prépondérance argileuse
(imperméabilité);
Des caractéristiques morphologiques favorables au ruissellement :
une forte pente et un relief bien développé ;
Une pluviométrie moyenne interannuelle supérieure à 400 mm ;
Une faible densité de drainage (0,06)
 4- Erosion:
Mesure et Formulation
Etude de l’envasement des lacs collinaires

Les objectifs d ’une mesure de l ’envasement sont

de deux ordres :
- connaître le volume d ’eau disponible et la
vitesse de sédimentation d ’un réservoir,
- actualiser les courbes hauteur-volume et hauteur-
surface d ’un lac pour reconstituer les apports
liquides et calculer les bilans hydrologiques du
bassin.
 Le canot pneumatique et son équipement

I.R.D.
 Le théodolite électronique

I.R.D
 L’envasement des lacs collinaires
A
250 m F id h A li

V o lu m e à la c o n s tr u c tio n ( 1 9 9 1 ) = 1 3 4 7 0 0 m 3
200 m V o lu m e a u 2 7 /0 6 /9 7 = 9 0 1 5 0 m 3

150 m

100 m
A

50 m

0 m

-5 0 m

-1 0 0 m

-4 0 0 m -3 5 0 m -3 0 0 m -2 5 0 m -2 0 0 m -1 5 0 m -1 0 0 m -5 0 m 0 m

IRD
 Evolution de l’envasement du lac collinaire de Fidh Ali
Profil en long

12

Coup2 1991
10 Coup2 1993
Coup2 1996

8 Coup2 1997
Hauteur à l'échelle (m)

0
0 50 100 150 200 250 300 350 400 450 500
Distance cummulée (m)

IRD
 Evolution de l’envasement dans le lac de Kamech
La première campagne de nivellement sur le site de Kamech a été effectuée le
01/01/1991 c. a. d. avant la mise en eau du lac (1994). Depuis cette mise en eau, 7
nivellements ont eu lieu (1994, 1995, 1996, 1998, 1999, 2002, 2004).

220000
70000

200000

60000 2004 180000

2004
1999 160000
50000 1999
1994 140000
1998
Surfaces en m²

40000
1991

volumes en m3
120000
1995
100000
30000
1995
1994
80000
1998
20000
1991
60000
2002
40000 2002
10000

20000

0 0
0 100 200 300 400 500 600 700 800 900 1000 1100 0 100 200 300 400 500 600 700 800 900 1000 1100
-20000
-10000

Cotes en cm cotes en cm

Evolution du barème côte/surface Evolution du barème côtes/volumes

IRD
 Evolution du volume envasé et de l’ablation cumulée
dans Kamech de 1991 à 2004
Volume Volume Ablation
Année utile en envasé en cumulée en
m3 m3 mm
1991 142000 0 0
1994 140920 1080 0,43
1995 139440 2560 1,04
1996 130715 11285 4,59
1998 127247 14753 6,01
1999 121000 21000 8,55
2002 112920 29080 11,84
2004 108000 34000 13,85

IRD
 Volumes des vases dans le lac Kamech de 1991 à 2004
40000

35000
Volume envasé en m3

30000

25000

20000

15000

10000

5000

0
1990 1992 1994 1996 1998 2000 2002 2004

Année

IRD
 Estimation de l ’envasement
-Les formules empiriques :
* Formule de Tixeront
Ts(t/km²/période) = a * Lr0,15(mm/période)

* Formule de Sogreah
Ts(t/km²/période) = b * Lr0,15(mm/période)
b dépend de perméabilité du bassin

* Formule de Fournier
Ts(t/km²/période) = (1/36) * (Pm/Pan)(mm) * (Dm2(m) /S0,46(km2) )

* Formule de Fersi :
Ts(t/km²/période) = d * Lr1,25(mm/période)
* Formule de Frigui :
Ts(t/km²/période) = c * Lr0,63(mm/période) * S-0,26(ha)

* Formule de la FAO :
Ts(m3/ha/période) = Fm * C1 * C2 * C3
Fm (indice de Fournier : ΣPim2/Pan),
C1 (coef. de texture : nature pédologique du sol),
C2 (coef. de topographie : classes des pentes),
C3 (coef. d’exploitation du sol : occupation des sols)
 Les pertes en sol calculées par les formules empiriques
(bassin versant Sbaihia)

Formules Perte en sol calculée Perte en sol

(m3/ha/an) observée (m3/ha/an)
Tixeront 3,6
Sogreah 9,55
Fournier 11,55 10,24
Fersi 35,14
Frigui 9,03
FAO 17,31
Mejri, T-2005
 ’érosion
’é
Etude de l’é ravinaire

Nivellement topographique Typologie-Morphométrie

Comparaison MNT2 – MNT1

Bilans abrasions-dépôts Recherche de relations

Analyse extensive

Khallili, A-2005
 Le nivellement topographique

Levés de terrain au
tachéomètre laser

I.R.D. 2001
 La Typologie
2 grands types des états de couverture végétale du sol

Khallili, A 2005
Vif
H 0 – 25 % H 25 – 50%
Vif1 Vif2

Enherbée
H 50 – 75 % H 75 – 100 %
H2 H1

H: Taux d’enherbement Ravine 2F

 Le bilan quantitatif Ravine

Dépôt

Ablation

Khallili, A-2005
 Bilan global des abrasions et des dépôts

Bassin versant de Kamech

Octobre 2006 – janvier 2007
A l’échelle de la ravine (0,172 ha)

abrasion en m3 59
dépôts en m3 26

Bilan global exporté en m3 33
 Étude de l’érosion
Relation entre Ab et Lr à l’échelle annuelle

Guermazi, L-2005
 Ruissellement et envasement (B.V. Kamech)
Période Pluie V ruisselé L.Ruisselée Ablation
(mm) (m3) (mm) (mm)
94-95 495 59400 24 1,04
95-96 1036 1091112 444 3,55
96-98 1105 192777 78 1,42
98-99 700 330498 182 2,54

4
3,5
3
Ablation en mm

2,5 Ablation (mm) = f(Lame ruisselée)

2
1,5
1 Ab= 0,89 Ln(Lr) - 2,03
0,5
0
0 100 200 300 400 500
lam m e ruisselée en m m

IRD
 Merci
pour votre attention