Renaud ECUYER

Elève ingénieur de l’E.N.S.E.E.I.H.T., I.N.P.T.

Hydraulique – Mécanique des fluides

APPLICATION DE TOPMODEL
A DIFFERENTS BASSINS VERSANTS

Rapport de stage de recherche

préparé au sein du
Laboratoire de Mécanismes de Transferts en Géologie
(UMR 5563 CNRS-OMP-UPS-IRD)

Encadré par
Mr. D. LABAT et Mr. D. DARTUS
 REMERCIEMENTS

Je tiens à remercier Mme. V. Estupina-Borel et Mr. J. Chorda de l’Institut de

Mécanique des Fluides de Toulouse, pour leur collaboration relative à l’utilisation de
ARCVIEW.
Je souhaite également exprimer ma reconnaissance à Mr. G. Cochonneau et Mr. J-L.
Guyot pour m’avoir fait accéder aux données hydrologiques du projet Hibam.
Enfin, je tiens à remercier le Laboratoire de Mécanismes de Transferts en Géologie
pour son accueil, et tout particulièrement Mr. D. Labat pour sa disponibilité, sa qualité
d’écoute et ses conseils avisés.
 RESUME
TopModel, le modèle hydrologique distribué fondé sur le concept d'aires contributives
variables de Beven et Kirkby (1979), a été testé sur trois bassins de taille et de nature
différentes : le bassin du Strengbach dans les Vosges, le bassin du Nsimi dans le Sud-
Cameroun et six sous-bassins du Rio Negro au Brésil. Les paramètres du modèle ont été
calibrés à l’aide du critère de Nash-Sutcliffe, afin de mesurer la qualité de la reconstitutions
des débits observés. Cette étude s’est attachée à donner un aperçu le plus complet possible de
la base conceptuelle du modèle. Des valeurs d’efficience encourageantes ont été obtenues, en
particulier pour les sous-bassins du Rio Negro, incitant à utiliser ce modèle pour des grands
bassins.

MOTS CLES
Modélisation hydrologique, TOPMODEL, petits bassins, grands bassins, modèle numérique
de terrain, aires contributives saturées.

ABSTRACT
TopModel, the topography-based variable contributing area model of Beven and
Kirkby (1979), was tested on three different catchments: the Strengbach catchment located in
France, the Nsimi catchment located in South Cameroon and six subcatchments of the Rio
Negro located in Brazil. Parameters were calibrated using the Nash-Sutcliffe efficiency
criteria. This approach provides useful insights into the structure of the model. Acceptable
numerical fits were reached, especially for the subcatchments of the Rio Negro, leading to a
promising exportation of TOPMODEL to large catchment.

KEY WORDS
Hydrological modeling, TOPMODEL, small catchment, large catchment, digital terrain
model, saturated contributing areas.
 TABLE DES MATIERES

INTRODUCTION 2

1. PRESENTATION DE TOPMODEL 3
1.1. Les aires contributives saturées 3
1.2. Base conceptuelle de TOPMODEL 3
1.2.1. Estimation des lames d’au entrante et sortante 4
1.2.2. L’indice topographique 5
1.3. Première application de GRIDATB sur des bassins versants synthétiques 7
1.3.1. Etude de l’effet de la topographie du bassin versant 7
1.3.2. Etude de l’effet de la précision du M.N.T. 8
1.4. Algorithme de TOPMODEL 9
1.4.1. Modélisation de la zone non saturée 9
1.4.2. Récapitulatif et algorithme 9
1.5. Paramètres de TOPMODEL 10

2. MISE EN ŒUVRE ET EVALUATION DE TOPMODEL POUR DES 11

PETITS BASSINS VERSANTS
2.1. Etude d’un petit bassin versant de moyenne montagne : Le Strengbach 11
2.1.1. Le bassin versant du Strengbach 11
2.1.2. Modélisation du comportement hydrologique du Strengbach 11
2.1.3. Etude de la variation des paramètres pour le Strengbach 13
2.1.4 Extension des surfaces saturées 14
2.2. Etude d’un petit bassin versant tropical : Le Nsimi 14
2.2.1. Le bassin versant du Nsimi 14
2.2.2. Modélisation du comportement hydrologique du Nsimi 16
2.2.3 Etude de la variation des paramètres pour le Nsimi 17
2.2.4. Extension des surfaces saturées 18

3.- MISE EN ŒUVRE ET EVALUATION DE TOPMODEL POUR DES 19

GRANDS BASSINS VERSANTS
3.1 Etude d’un grand bassin versant tropical : le Rio Negro 19
3.2. Modélisation du comportement hydrologiques des six sous-bassins 20
3.3. Etude de la variation des paramètres pour le Rio Negro 23
3.4 Extension des surfaces saturées 26

CONCLUSION 27

BIBLIOGRAPHIE 29

-1-
INTRODUCTION

Ce stage s’est déroulé au sein du Laboratoire de Mécanismes de Transferts en

Géologie (L.M.T.G.), dans l’équipe « Eau Sol Erosion », sous la direction de D. Labat et en
collaboration avec D. Dartus de l’Institut de Mécanique des Fluides de Toulouse (I.M.F.T.).

L’objet scientifique de l’hydrologie est l’étude du cycle de l’eau. L’eau circule sur la
terre, sous la terre et dans l’atmosphère selon un cycle qui est bien connu.
Les nuages apportent des précipitations sous forme de pluie, de neige ou de grêle.
L’eau ruisselle sur la terre. Une partie est captée par la végétation. Le restant va
alimenter les rivières ou va s’infiltrer dans le sous-sol pour atteindre les nappes
souterraines.
L’eau des fleuves, des lacs et des océans s’évapore sous l’action des rayons du soleil et
se retrouve à l’état gazeux dans l’atmosphère.
La vapeur d’eau de l’atmosphère se condense en gouttelettes au contact des masses d’air
froid, ce qui entraîne la formation de nuages.

La modélisation hydrologique s’attache à représenter schématiquement les différents

phénomènes physiques liés à l’eau et à son comportement afin de mieux les étudier et
d’analyser l’influence qu’ils exercent. L’unité géographique sur laquelle se base l’analyse du
cycle hydrologique et de ses effets est le bassin versant. Le bassin versant est une surface
élémentaire en théorie hydrologiquement close, c’est-à-dire qu’aucun écoulement n’y pénètre
de l’extérieur et que tous les excédents de précipitations s’évaporent ou s’écoulent par une
seule section à l’exutoire.

Dans cette optique, nous nous sommes focalisés sur un modèle en particulier,
TOPMODEL, qui permet de prédire le comportement hydrologique des bassins versants en
fonction de leur topographie. Il sera appliqué à trois bassins versants de taille distincte, situés
sous différents climats. Initialement développé par K.J. Beven et M.J. Kirby vers la fin des
années 1970 à l’Université de Lancaster, TOPMODEL est articulé autour de deux postulats
centraux : le ruissellement se produit sur des surfaces contributives d’aire variable, et la
topographie influence la manière dont se produit le ruissellement. Il est devenu très populaire
ces dernières années avec la généralisation de l’utilisation des Modèles Numériques de
Terrain (M.N.T.), et permet de plus, de spatialiser la réponse hydrologique.

Cette étude s’inscrit dans une démarche exploratoire. Idéalement conçu pour des
bassins de moyenne montagne humide, les limites de TOPMODEL sont ici allégrement
franchies vues la nature et la taille de certains bassins versants étudiés.

Après avoir présenté la méthode d’investigation retenue et les équations de base du

modèle, nous l’appliquerons tout d’abord sur un petit bassin versant de moyenne montagne
(cas d’école d’application de TOPMODEL), puis sur un petit bassin versant tropical et enfin
sur un grand bassin versant tropical, afin d’étudier ses capacités d’exportation et ses limites.

-2-
1. PRESENTATION DE TOPMODEL
TOPMODEL (Beven and Kirkby 1979 ; Beven et al. 1995 ; Beven 1997 et Kirkby
1997) est un modèle hydrologique distribué fondé sur le concept d'aires contributives
variables, initialement proposé par Cappus (1960). Ce modèle prédit, à chaque pas de temps,
la distribution spatiale du contenu en eau au sein des différentes mailles du modèle numérique
de terrain (M.N.T.) du bassin versant considéré. Ce contenu en eau est calculé en fonction de
la distribution d'un indice de similarité hydrologique et du déficit en eau moyen sur le bassin
versant. Le déficit moyen est calculé à partir d'un bilan d'eau à chaque pas de temps (Beven et
al. 1995). Nous allons étudier à présent son principe de fonctionnement.

1.1. Les aires contributives saturées

On distingue deux concepts de génération de ruissellement.

Le premier processus proposé pour expliquer la formation du ruissellement est celui de
Horton (1945). Ce mécanisme considère que le ruissellement est contrôlé par la capacité
d’infiltration des sols. Le ruissellement peut ainsi être généré sur toute la surface de bassin
lorsque les intensités de pluie dépassent la capacité d’infiltration des sols. Il y a donc
ruissellement par refus à l’infiltration.
Cependant dans certains cas, ce concept ne permet pas d’expliquer la formation du
ruissellement. Une autre idée est alors apparue, celle des aires contributives, qui fait jouer un
rôle essentiel aux zones saturées. Cappus (1960) est probablement le premier, a avoir proposé
ce mécanisme, observant que sur certains bassins « le ruissellement est produit par la totalité de
la pluie tombant sur une zone imperméable ou saturée dont l’étendue varie avec le niveau
moyen de la nappe phréatique ». Ces surfaces imperméabilisées parce que la nappe affleure en
ces lieux la surface du sol, constituent ce que l’on nomme les zones contributives. C’est ce type
de mécanisme de génération de ruissellement que TOPMODEL s’attache particulièrement à
bien représenter.

Refus à l’infiltration Surfaces contributives

Figure 1 – Les deux mécanismes de ruissellement

1.2. Base conceptuelle de TOPMODEL

L’échelle de TOPMODEL est celle du versant. Les bilans sont essentiellement

effectués à cette échelle. Par versant, nous entendons un « tube de courant » bordé :
- à l’amont par la ligne de crête (condition de flux nul)
- latéralement par deux lignes de courant passant par les points A et B de la
figure suivante, qui à l’aval définissent un segment curviligne le long d’une ligne de niveau,
orthogonale aux lignes de courant (Figure 2).

-3-
 Figure 2 - Schématisation d’un versant

Dans le cas d’un élément situé en milieu de versant, sans ruissellement de surface en
régime permanent, pour tout temps t, toute la quantité d’eau alimentant la partie de la zone
saturée, située à l’amont de l’élément de versant considéré, est égale à la quantité d’eau
s’évacuant vers l’aval au droit de ce même élément de contour. A l’échelle d’un versant dont
la largeur du pied de versant est unitaire, on a ainsi égalité entre la lame d’eau entrante et la
lame d’eau sortante à tout temps t.

1.2.1. Estimation des lames d’au entrante et sortante

Si on suppose que l’écoulement subsuperficiel suit une loi de Darcy, la lame d’eau
sortante est proportionnelle à la transmissivité du sol, notée T, sur laquelle se transmet le flux
d’eau et au gradient de charge hydraulique, noté grad(h) (pente de la surface de la nappe en
aquifère libre).

De plus, dans TOPMODEL, la transmissivité est exprimée en terme de déficit local en

eau. Le déficit local en eau, noté d, définit la lame d’eau qu’il faudrait infiltrer pour faire
affleurer la zone saturée, en supposant que la zone non saturée soit initialement vide. On
supposera donc que la transmissivité latérale est fonction de ce déficit.
De plus si on admet que, localement, le gradient hydraulique ne varie pas au cours du
temps, celui-ci peut s’identifier à la pente locale, notée tg b. Remarquons que cette
simplification ne signifie pas que la zone saturée soit globalement parallèle à la surface du sol.
Cela signifie qu’elle l’est localement. Et de ce fait, sur un profil de versant, les zones saturée
ont l’allure illustrée par la Figure (3).

Figure 3 – Pente de la zone saturée

-4-
 On peut donc écrire pur tout instant t et pur tout versant i,

Lame sortante i ,t = T(d i,t ) . tg b i (1)

avec T(d i,t ) la transmissivité de l’épaisseur de sol sur laquelle se transmet le flux d’eau
tg b i la pente de la surface du sol

Historiquement (Beven et Kirby 1979, Beven et Germann 1982), un profil exponentiel

de la transmissivité a été proposé, où l’on suppose les sols homogènes :

T(d i,t ) = T0 exp  − d i, t  (2)

 M 

avec T0 la transmissivité de l’ensemble du profil

M le paramètre de décroissance de cette transmissivité

Beven et Germann (1982) ont justifié ce profil théorique d’après des mesures de profil
exponentiel de décroissance de la conductivité à saturation avec la profondeur. La
généralisation à d’autres profils théoriques (linéaire ou puissance, par exemple) est tout à fait
possible (Ambroise et al. 1996, Iorgulescu et Musy, 1997), mais n’a pas été considérée dans
cette étude.

La lame d’eau entrante s’exprime comme l’intégrale recharges élémentaires sur le

versant considéré par percolation verticale, en un point du bassin, de la zone non saturée vers
la zone saturée. On a donc :
Lame entrante i,t = ai R t (3)
avec ai l’aire amont de l’élément de versant, d’indice i et de largeur unitaire
Rt la moyenne des recharges élémentaires, sur la surface amont du pied de versant

1.2.2. L’indice topographique

L’équation de continuité implique la relation suivante :

Lame entrante i ,t = Lame sortante i,t ⇔ T(d i,t ) . tg b i = a i R t (4)

Cette relation permet de connaître l’état du déficit local en eau, en chaque point du bassin et
pour chaque pas de temps. Cette équation formalise l’idée maîtresse de TOPMODEL : celle
du rôle de la topographie dans la variabilité des états hydriques des sols. En effet, cette
relation montre que tous les points du bassin présentant le même rapport a auront à
tg.b
chaque pas de temps la même valeur de déficit local

La connaissance du profil de transmissivité nous permet de déterminer l’expression

reliant le déficit local au déficit moyen, noté Dt :

-5-
   
Dt - d i,t = - M  IT −Ln a i   (5)
  tgb.i 

avec IT la moyenne sur l’ensemble du bassin des valeurs de l’indice

 
topographique ITi =Ln a i 
 tgb.i 
M le paramètre de décroissance du profil de transmissivité
1
D t = ∫ d i , t . da déficit moyen
AA
avec A la surface du bassin
da une surface élémentaire autour du point i

A ce niveau, une interprétation physique de l’indice topographique doit être donnée.

L’indice topographique Ln ai  est une variable caractéristique de TOPMODEL qui traduit
tgb.i
 
le potentiel d’un point du bassin à se saturer plus ou moins facilement.
Ainsi l’indice topographique peut être perçu comme le rapport entre le flux entrant au point
considéré, donc sa capacité de collecte (l’aire drainée ai ), et le flux sortant du profil de ce
point, donc sa capacité d’évacuation (le gradient gravitaire tg bi).

Ce rapport traduit ainsi la propension des points du bassin à s’engorger, à se saturer

par extensions des zones saturées :
- Les fortes valeurs de l’indice topographique correspondent aux points dont le
flux entrant est bien plus grand que le flux maximum qu’il lui est possible de faire transiter.
Ces points seront donc rapidement saturés.
- Les faibles valeurs de l’indice topographique correspondent aux points dont la
transitance est grande en comparaison de ce qu’ils vont être susceptibles de devoir faire
transiter. Ces points ne seront probablement que rarement saturés.
A chaque pas de temps, il existe donc une valeur limite de l’indice topographique telle
que, à cet instant, tous les points du bassin ayant une valeur d’indice topographique supérieure
à cette valeur limite seront saturés (déficit nul).
La relation (5) permet de constater qu’à chaque instant, tous les points du bassin ayant
la même valeur d’indice topographique auront la même valeur de déficit local.
Il n’est donc pas nécessaire de simuler chacun de ces points de manière individuelle. Il suffit
de simuler un seul de ces points, et de considérer que les résultats obtenus sont généralisables
à tous les points de même valeur d’indice topographique : c’est ce que fait l’algorithme de
TOPMODEL. Cela permet également de supposer que toute surface saturée est contributive
au débit de ruissellement véhiculé par le réseau.
Donc pratiquement, l’algorithme de TOPMODEL effectue des simulations sur des
classes de valeurs d’indice topographique. Puis les réponses totales du bassin sont calculées
comme :
- la somme de chacune des réponses simulées pour chaque classe d’indice
- pondérée par la représentativité de chacune des classes d’indice sur le bassin étudié.

Nous allons maintenant présenter la méthode qui nous permet d’obtenir les classes de
valeurs d’indice topographique et leur représentativité sur le bassin.

-6-
 Un modèle numérique de terrain (M.N.T.) est une représentation numérique de la
topographie d’un bassin. C’est une grille régulière, que l’on superpose virtuellement sur le
bassin. On attribue à chaque nœud de la grille, la valeur de l’altitude du point matériel
correspondant sur le bassin. On procède donc à un échantillonnage à deux dimensions et à pas
constant de l’altitude du bassin.
GRIDATB est un programme qui permet d’obtenir les classes de valeurs d’indice
topographique et leur représentativité sur le bassin, à partir du M.N.T. du bassin étudié, en
utilisant un schéma multidirectionnel pour la propagation de l’aire drainée (Quinn et al. 1991,
Wendling (1992), Laroche (1994)). Dans ce schéma, l’aire drainée par unité de longueur de
chaque point du bassin, est calculée comme étant la surface du bassin représentée par tous les
points s’écoulant au travers du point considéré. La représentativité sur le bassin des classes
d’indice topographique étant fonction de la résolution du M.N.T., il est conseillé de fixer
celle-ci à une valeur inférieure ou égale à 50 m.

1.3. Première application de GRIDATB sur des bassins versants synthétiques

Afin d’étudier le comportement de GRIDATB, nous nous sommes d’abord intéressés à

des bassins synthétiques. Nous avons donc créé un programme qui permet de générer des
bassins versants élémentaires de type xn sur la première partie aval puis en forme de cuvette
dans la deuxième partie amont à l’aide du logiciel MATLAB (Figure 4).

20 400
8000
350
7000
15 300
6000
250
5000
10
200
4000

150 3000
5
100 2000

50 1000
0
40 0
40 0 40
30 30 40 30 40
20 20 30 40 20 30
10 20 30 20
10 10 20 10 10
0 0 10 0
0 0 0

15
14 15

12

10
10
10
8

4
5 5
2

0
40
30 0
40 40
20 30 0 40
40 30
20 30 20 30
10 20 30 40 20
10 10 20 10 10
0 0 0 10 0 0
0

a b c

Figure 4 – Topographie et distribution de l’indice topographique IT

pour trois bassins synthétiques (a : n=1 , b : n=2, c : n=3)

1.3.1. Etude de l’effet de la topographie du bassin versant

Pour vérifier l’effet de la topographie du bassin sur la distribution de l’indice

topographique, on trace sur un même graphique cette distribution pour les trois bassins. Ces

-7-
distributions ont été tracées avec une précision de 50 m (dx=50) et une extension spatiale de
40 mailles suivant x et suivant y (nx=ny=40) (Figure 5).

0,3

0,25
bassin 1
pourcentage
0,2
bassin 2
0,15
bassin 3
0,1

0,05

0
0 2 4 6 8 10 12 14
indice topographique

Figure 5 - Distribution de l'indice topographique pour les trois bassins pour nx=ny=40,dx=50

Rappelons l’expression de l’indice topographique :

ITi =Ln a i 
 tgb.i 

Supposons que pour un même point l’aire drainée par unité de longueur soit
sensiblement la même pour les trois bassins. Si pour ce même point on augmente la concavité
du bassin, ceci a pour effet de diminuer localement la valeur de la pente, donc d’augmenter la
valeur de l’indice topographique et à fortiori d’étaler sa distribution.
On en conclut que plus la forme du bassin est incurvée, plus l’étalement de la
distribution de l’indice topographique est importante.

1.3.2. Etude de l’effet de la précision du M.N.T.

Afin d’étudier l’effet de la précision du M.N.T. sur la distribution de l’indice

topographique, on trace sur un même graphique cette distribution pour n=3, en gardant les
mêmes proportions, pour diverses précisions dx (Figure 6).

0,3

0,25 nx=ny=40 dx=50

nx=ny=50 dx=40
pourcentage

0,2
nx=ny=20 dx=100
0,15 nx=nx=32 dx=62.5

0,1

0,05

0
0 2 4 6 8 10 12 14
indice topographique

Figure 6 - Distribution de l'indice topographique

pour n=3 pour différentes résolutions du MNT

-8-
 On vérifie que la qualité de la distribution de l’indice topographique dépend
essentiellement de la précision du M.N.T. : plus la précision est importante, meilleure est la
qualité de la distribution. On remarque également que les distributions pour dx=40 et dx=50
sont sensiblement les mêmes. En effet, ces deux distributions possèdent les mêmes valeurs de
classe d’indice topographique à 10-2 prés. Ceci confirme le fait qu’il est conseillé de travailler
avec une précision dont la valeur est inférieure ou égale à 50 m.

1.4. Algorithme de TOPMODEL

1.4.1. Modélisation de la zone non saturée

A l’intérieur de la zone non saturée, TOPMODEL distingue un compartiment

correspondant à la zone racinaire. Lorsque la pluie arrive au sol, deux processus peuvent se
produire :
- si la pluie est inférieure à la conductivité à saturation à la surface du sol K 0 = T0 , elle
M
s’infiltre complètement
- sinon, une partie seulement égale à K0 s’infiltre, et l’excédent ruisselle. En pratique, vu les
valeurs de K0 , ce type de ruissellement ne sera jamais activé dans nos simulations de débits.
L’eau infiltrée tombe alors dans le compartiment racinaire : de profondeur SRMax
(Soil Root Maximum level). Tant que ce compartiment n’est pas rempli, l’eau ne pénètre pas
dans la zone non saturée. Ce réservoir a donc pour fonction de simuler l’eau de pluie qui est
stockée dans la partie superficielle du sol en début d’épisode pluvieux.Une partie de l’eau de
ce réservoir peut également être évaporée : on calcule ainsi à chaque pas de temps un flux
sortant incluant à la fois l’évaporation, mais aussi l’interception et la transpiration des plantes.
L’eau, issue du réservoir racinaire, est alors transférée dans le réservoir de la zone non
saturée. L’eau est transférée en un pas de temps, au fond de la zone non saturée : cette eau
peut alors alimenter la zone saturée (ou nappe).
Dans le cas où le point est saturé, la zone non saturée est inexistante : toute la pluie qui
tombe sur ce point ruisselle et est transférée directement à l’exutoire par une fonction de
transfert linéaire caractérisée par la vitesse ChVel (Channel Velocity). Il s’agit de la
composante rapide de la formation des crues.

1.4.2. Récapitulatif et algorithme

Nous avons désormais à notre disposition tous les éléments pour mettre en œuvre
TOPMODEL. Néanmoins, dans un souci de clarté, les étapes essentielles de la modélisation
sont énumérées ci-après :

Détermination du premier déficit moyen

Calcul de la répartition des déficits locaux autour du déficit moyen, en fonction de
la topographie
Gestion de la zone non saturée
déficit strictement positif :
remplissage éventuel de la zone racinaire
pompage de l’évapotranspiration
infiltration dans la zone non saturée
alimentation de la nappe

-9-
 déficit nul :
génération de ruissellement dans le cas où la pluie est non nulle
Calcul de la recharge globale de la nappe
Calcul du débit d’exfiltration de la nappe
Mise à jour du déficit moyen
Fin du pas de temps (retour au point pour le début du pas de temps suivant)

A la fin de chaque pas de temps, on dispose d’une lame d’eau correspondant à

l’exfiltration de la nappe et éventuellement d’une lame d’eau ruisselée. Pour calculer le débit
à l’exutoire, ces deux composantes sont ensuite transférées à l’exutoire.

1.5. Paramètres de TOPMODEL

Tout au long du développement des concepts et de la théorie de TOPMODEL, nous

avons fait apparaître divers paramètres (Tableau 7) :

Paramètre Description unité

M paramètre de décroissance dans le profil exponentiel de la transmissivité m
Ln(To) logarithme de la transmissivité de l’ensemble du profil m²/h
SRMax profondeur du réservoir racinaire m
SRinit valeur initiale du déficit local de la zone racinaire m
ChVel vitesse caractéristique de la fonction de transfert m/h

Tableau 7 – Liste des différents paramètres de calage de TOPMODEL

Des tests de sensibilité de TOPMODEL à ces différents paramètres ont été développés
par Franchini et al. (1996). Pour caler ces paramètres, il faut disposer d’une mesure de la
qualité des reconstitutions en débits observés. Classiquement, c’est le critère de Nash, ou
efficience (Nash et Sutcliffe, 1970), qui est retenu pour évaluer les simulations de débits. Ce
critère se définit de la manière suivante :

n
∑ (Q −Qsim ,t )
2
obs ,t

Efficience =1− t =1
(6)
∑ (Q )
n 2
obs ,t −Qobs
t =1

avec Qobs,t le débit observé à l’instant t

Qsim,t le débit simulé par le modèle à l’instant t
Qobs la moyenne des débits observés
n le nombre de pas de temps de l’épisode

Pour que l’efficience soit proche de 1, c’est-à-dire pour que la reconstitution soit jugée
bonne, il faut non seulement que la forme de la crue corresponde à la forme réelle, mais aussi
que les volumes générés soient proches de ceux observés.

- 10 -
2. MISE EN ŒUVRE ET EVALUATION DE TOPMODEL POUR DES
PETITS BASSINS VERSANTS
Nous allons, à présent, appliquer TOPMODEL sur deux petits bassins versants situés sous
deux climats différents : la bassin du Strengbach et le bassin du Nsimi. Ceci permet une
première approche des limites à l’utilisation de TOPMODEL. Notons que TOPMODEL à
déjà été utilisé dans le cadre de petits bassins : dans le Massif des Maures (Wendling 1992 et
Saulnier 1998), en Lozère (Durand et al. 1992), en Angleterre (Robson et al. 1993), en Suisse
(Iorgulescu et al. 1994), en Slovaquie (Holko et al. 1997) et en Ardèche (Datin 1998).

2.1. Etude d’un petit bassin versant de moyenne montagne : Le Strengbach

2.1.1. Le bassin versant du Strengbach

Le bassin-versant du Strengbach est situé sur le côté alsacien des Vosges, au N-W de
la commune d'Aubure et à une soixantaine de kilomètres au S-O de Strasbourg. D'une
superficie de 80 ha, son altitude varie entre 1146 m à son point culminant et 883 m à son
exutoire. La pente moyenne est de 15° (Figure 8).

Figure 8 – Modèle Numérique de Terrain du bassin versant du Strengbach

Le bassin versant est boisé à 97%, dont une majorité de résineux (épicéas
essentiellement) mais aussi 30% environ de hêtres. Sur le plan climatique, les précipitations
s'échelonnent de 1000 à 1500 mm par an, le relief influant considérablement sur leur
variabilité spatiale. Le climat est de type océanique montagnard et les moyennes de
température mensuelles varient de -1,2°C à 14°C. Le bassin versant du Strengbach est drainé
par un cours d'eau principal qu'alimentent quatre affluents. Le module de débit annuel et de 20
l/s. Les mesures se font à l'aide 48 pluviomètres de 400 cm2 répartis sur l'ensemble du site
(Latron 1990, Brendle 1992 et Ladouche 1997).

2.1.2. Modélisation du comportement hydrologique du Strengbach

Les seules données d’entrée requises par TOPMODEL sont les mesures de pluie, de
débit, et d’évapotranspiration potentielle (combinaison de l’évaporation physique et de la
transpiration végétale). Pour le Strengbach, on dispose de données au pas de temps journalier
pour une période s’étalant de 1989 à 1995.

- 11 -
 Afin d’obtenir le M.N.T. du bassin à partir du relevé des points topographiques avec
une précision de 29 m, nous avons utilisé le module hydrologique du logiciel ARCVIEW-
GIS. Suite à ce traitement topographique, nous avons appliquer GRIDATB pour obtenir la
distribution et la représentativité des classes de valeur d’indice topographique sur le bassin
(Figure 9). Enfin, nous avons utilisé la version Windows initiale de TOPMODEL, distribuée
par l’Université de Lancaster, afin de caler les différents paramètres. Ce calage est effectué à
l'aide de l'analyse de Monte Carlo et du logiciel GLUE, disponible avec TOPMODEL, en
utilisant le critère de Nash-Sutcliffe comme critère de qualité.

0,21

Pourcentage
0,14

0,07

0
0 5 10 15
Indice topographique

Figure 9 – Evaluation de l’indice topographique sur le bassin du Strengbach

(à gauche : distribution spatiale, à droite : fonction de distribution)

On obtient le résultat suivant :

1 365 729 1093 1457 1821 2185

0 0,05

0,04
0,1 débits observés
débits simulés Débit ( m / j )
0,03
Pluie ( m )

0,2
0,02

0,3
0,01

0,4 0

Figure 10 –Comparaison entre débits observés et simulés à

l’exutoire du Strengbach sur un période de 2191 jours

Le modèle permet de reproduire à la fois les débits de pointe (t = 138, 813, 1149,
1543) mais aussi les débits d’étiage. A l’échelle d’une crue, il rend aussi bien compte du
temps de monté que de la forme de la courbe de récession. Il ne génère pas de crue factice et
répond aux impulsions de pluie. Ceci conduit à une valeur d’efficience élevée de l’ordre de
0.8 (calculée à l’aide de l’équation (6)). Idéalement conçu pour des bassins de moyenne
montagne humide, il est encourageant de trouver une valeur d’efficience aussi élevée pour le
Strengbach. Notons que des valeurs d’efficience similaires, ont été obtenues avec
TOPMODEL sur le bassin versant du Ringelbach, qui se situe à quelques kilomètres du
Strengbach (Ambroise et al. 1996).

- 12 -
 2.1.3. Etude de la variation des paramètres pour le Strengbach

Dans une deuxième partie, afin d’étudier une possible variation des paramètres de
calage au cours du temps, nous allons successivement appliquer TOPMODEL, pour des
périodes plus courtes : pour environ deux années (Figure 11) et à l’échelle d’une crue (Figure
12). Les nouvelles simulations et valeurs obtenues sont regroupées dans le Tableau 13 :

1 101 201 301 401 501 601

0 0,07
0,02
0,06
0,04
0,06 débits sim ulés 0,05

Débits ( m / j )
débits observés
Pluie ( m )

0,08 0,04
0,1
0,12 0,03

0,14 0,02
0,16
0,01
0,18
0,2 0

Figure 11 – Comparaison entre débits observés et simulés à

l’exutoire du Strengbach sur un période de 681 jours

1 31 61 91
0 0,06

0,04

ef f ic ienc e = 0,97 débits obs erv és 0,04

Débits (m / j)
Pluie ( m )

0,08
débits s imulés
0,12
0,02

0,16

0,2 0

Figure 12 – Comparaison entre débits observés et simulés à

l’exutoire du Strengbach sur une période de 117 jours
Paramètre Valeur
Période d’étude 2191 j 681 j 117 j
M 0.033 m 0.031 m 0.033 m
Ln(To) 5 m²/h 5 m²/h 9 m²/h
SRMax 5.2 m 5.2 m 5.2 m
SRinit 0.012 m 0.035 m 0.03 m
ChVel 950 m/h 640 m/h 852 m/h
Efficience 0.79 0.86 0.97

Tableau 13 – Variation des paramètres de calage pour

le Strengbach pour des périodes d’étude de plus en plus courtes

- 13 -
 A l’échelle de deux années et d’une crue, le modèle reproduit tout à fait
convenablement les débits observés avec des efficience respectives de l’ordre de 0.85 et 0.95.
Dans le cas du Strengbach, on remarque que tous les paramètres sont rigoureusement
du même ordre de grandeur. Les petites variations observés, notamment pour ChVel, ne sont
pas significatives pour le modèle à cet ordre. Les valeurs de SRinit sont différentes, car
rappelons le, celui-ci est sensé prendre en compte l’état hydrique du sol en début d’épisode
pluvieux, et ce dernier est différent pour les trois simulations. Enfin, les valeurs d’efficience
augmentent inversement à la période d’étude, vu que le nombre de débits à simuler diminue.

2.1.4. Extension des surfaces saturées

Pour étudier l’extension des surfaces saturées, nous avons utilisé la version initiale de
TOPMODEL, en code FORTRAN, fournie également par l’Université de Lancaster, celle-ci
nous donne pour chaque pas de temps les valeurs des différentes lames et du déficit moyen.
Puis, nous avons créé un programme en langage FORTRAN qui utilise les résultats fournis
par cette version, afin de calculer les valeurs des déficits locaux et donc de visualiser ces
surfaces (Figure 14).

12

10
pourcentage

0
0 500 1000 1500 2000

Figure 14 - Évolution du pourcentage de surfaces saturées

en fonction du temps (valeur moyenne de 1%)

Ladouche (1997) a montré que l’extension maximale cartographiée représente 2,8 %

de la superficie totale du bassin tandis que l’extension minimale mesurée ne couvre que 0,3 %
du bassin versant. Notre simulation conduit à un pourcentage moyen de surface saturée de
l’ordre de 1 %, celle-ci semble donc conforme aux données expérimentales.
Dans le cas du Strengbach, TOPMODEL rend donc convenablement compte de
l’extension spatiale des surfaces saturées.

2.2. Etude d’un petit bassin versant tropical : Le Nsimi

2.2.1. Le bassin versant du Nsimi

Le bassin versant du Mengong appartient au bassin du Nyong (Figure 15). Il est situé à
120 km au SSE de Yaoundé et à 13 km au SW de Zoetelé dans le village de Nsimi.

- 14 -
 Figure 15 – Localisation géographique du bassin versant du Nsimi

D’une superficie de 60 ha, le bassin versant du Mengong présente un profil

dissymétrique caractérisé par des versants convexo-concaves raccordés par un important bas-
fond marécageux représentant 15 à 20% selon la saison de sa superficie totale. Les altitudes
vont de 670 m à 705 m. (Figure 16).

Source 1

Figure 16 – Topographie, réseau hydrographique et layons du bassin versant du Nsimi

Le cheminement du Mengong dans la vallée marécageuse peut être découpé en trois

domaines de pentes caractéristiques :
• une pente de 4% qui part des sources jusqu’à 150 m environ en aval,
• une pente presque nulle dans le cours moyen du ruisseau entre 150 m et 450 m
responsable du mauvais drainage
• une pente de 1% entre 450 m et 860 m (exutoire).

Le bassin a été quadrillé par une dizaine de layons rectilignes (chemins

entretenus) orienté perpendiculairement à l’axe du ruisseau et le long desquels sont installés
les appareils de mesure et de suivi hydrochimiques et hydrologiques : fosses pédologiques
atteignant le plus souvent la nappe, sondages profonds, source et exutoire aménagé, réseaux
de pluviométres et de piézomètres, limnigraphes et pluviographes (Ndam Ngoupayou 1997 et
Bodin 1998).

- 15 -
 2.2.2. Modélisation du comportement hydrologique du Nsimi

Pour le Nsimi, on dispose des mesure journalières de pluie, débit, et évapotranspiration

potentielle la période allant de 1994 à 1998. Néanmoins cette étude se bornera à une seule
année hydrologique.
L’obtention du M.N.T. à partir du module hydrologique de ARCVIEW-GIS n’a pas été
concluante, vu la disposition des points topographiques. En effet, pour le Nsimi, le relevé de
ces points s’est fait le long des layons, ce qui pose des problèmes lors de l’interpolation
effectuée par ARCVIEW-GIS. Afin de trouver l’interpolation la plus adaptée, nous avons
utilisé le logiciel WINSURF. Finalement, nous avons choisi l’interpolation Minimum
Curvature qui est largement employée en Sciences de la Terre, avec une précision de 10 m.
Après avoir déterminé la distribution de l’indice topographique, nous avons calé les
paramètres de façon analogue au bassin précédent (Figure 17).

0,25

0,2
Pourcentage

0,15

0,1

0,05

0
0 4 8 12 16
Indice topographique

Figure 17 – Evaluation de l’indice topographique sur le bassin du Nsimi

(à gauche : distribution spatiale, à droite : fonction de distribution)

On obtient le résultat suivant :

1 31 61 91 121 151 181 211 241 271

0 0 ,0 1 2

0 ,0 1
0 ,0 4

d é b its o b s e rvé s 0 ,0 0 8
d é b its s im u lé s
Débits ( m / j )

0 ,0 8
Pluie ( m )

0 ,0 0 6

0 ,1 2
0 ,0 0 4

0 ,1 6
0 ,0 0 2

0 ,2 0

Figure 18 – Comparaison entre débits observés et simulés à

l’exutoire du Strengbach sur une période de 277 jours

TOPMODEL permet de reproduire de manière convenable les débits de pointe,

néanmoins certaines crues demeurent sous-estimées (t=91, 110, 241). Les étiages et les
montés de crues sont quant à eux correctement reproduits. Par contre, les courbes de décrues

- 16 -
ne sont pas reproduites aussi convenablement que pour le Strengbach. Toutefois, comme pour
le Strengbach, le modèle ne génère pas de crue factice et répond aux impulsions de pluie. Ceci
conduit à une efficience moyenne de l’ordre de 0.65, qui indique que TOPMODEL peut être
tout à fait utilisé pour ce type de bassin.
Les différences peuvent s’expliquer de différentes manières. Les pentes sont très
faibles, voir quasi nulles sur certaines parties du bassin versant. Or la donnée fondamentale
pour TOPMODEL est l’indice topographique qui est justement calculé à partir de ces pentes.
Il est donc naturel que TOPMODEL éprouve des difficultés pour modéliser ces débits. De
plus une étude sur un bassin versant tropical du Nigeria, pays voisin du Cameroun, a montré
que la pente de la zone saturée pouvait ne pas être localement parallèle à la surface du sol, une
hypothèse fondatrice du modèle peut donc ne pas être validé par ce bassin (Campling et al.
2002).

2.2.3. Etude de la variation des paramètres pour le Nsimi

Dans une deuxième partie, afin d’étudier une possible variation des paramètres de calage
au cours du temps, nous allons appliquer TOPMODEL, pour une période plus courte : environ
un mois de l’année hydrologique (Figure 19). Les nouvelles valeurs et simulation obtenues
sont regroupées dans le Tableau 20:

1 11 21 31
0 0 ,0 1 5

0 ,0 4 0 ,0 1 2

Débits ( m / j )
Pluie ( m )

0 ,0 8 d é b its o b s e rvé s 0 ,0 0 9
d é b its s im u lé s
0 ,1 2 0 ,0 0 6

0 ,1 6 0 ,0 0 3

0 ,2 0

Figure 19 – Comparaison entre débits observés et simulés à

l’exutoire du Nsimi sur une période de pour 33 jours

Paramètre Valeur
Période d’étude 277 j 33 j
M 0.94 m 0.3 m
Ln(To) 0.29 m²/h 0.28 m²/h
SRMax 0.003 m 0.0001 m
SRinit 0.13 m 0m
ChVel 365 m/h 480 m/h
Efficience 0.63 0.78
Tableau 20 – Variation des paramètres de calage
pour le Nsimi pour deux périodes d’étude

- 17 -
 A l’échelle du mois, le modèle reproduit tout à fait convenablement les débits observés
avec une efficience de l’ordre de 0.8. La principale différence se trouve dans le retard à la
première crue, sans doute imputable à l’initialisation du modèle.
Dans le cas du Nsimi, deux paramètres, SRMax et SRinit, ont varié d’un ordre de
grandeur. La variation de SRMax peut être du au fait que l’hypothèse de parallélisme de la
pente de la zone non saturée n’est pas respectée ainsi que des faibles valeurs de pentes.

2.2.4. Extension des surfaces saturées

De même que pour le bassin précédent, nous avons essayé de visualiser les surfaces
saturées (Figure 21) :

8,5

8
Pourcentage

7,5

6,5

6
0 50 100 150 200 250

Figure 21 - Évolution du pourcentage de surfaces saturées en fonction du temps

On voit que les résultats fournies par TOPMODEL, bien que permettant de réaliser
une bonne simulation hydrologique ne nous autorisent pas de prévoir l’évolution des zone
saturées dans le cas du Nsimi. Ceci vient du fait que les pentes sont quasi nulles et de la
présence de la zone marécageuse qui représente, rappelons le, 15 à 20 % du bassin versant.

- 18 -
3. MISE EN ŒUVRE ET EVALUATION DE TOPMODEL POUR DES
GRANDS BASSINS VERSANTS
Nous allons à présent appliquer TOPMODEL sur de très grands bassins versants tropicaux
de tailles et de propriétés différentes. Les données sont émises de la banque du projet HIBAM
(ANA-IRD-HIBAM). Notons que des applications de TOPMODEL, ou de modèles dérivés
de TOPMODEL, à de grands bassins versants ont déjà été proposées pour des bassins
versants du Japon et de Chine (Takeuchi et al. 1999), de la Corée (Deg-Hyo et al. 2000), du
Mekong (Nawarathna et al. 2001 et 2002), du Japon et des USA (Tianqui et al. 2002).

3.1. Etude d’un grand bassin versant tropical : le Rio Negro

Affluent de rive gauche de l'Amazone d’une longueur de 2 200 km, le Rio Negro né
en Colombie, sous le nom de Guainía, il détermine la frontière avec le Venezuela, avant
d'arroser le nord du Brésil et de confluer en aval de Manaus. Le bassin versant du Rio Negro,
qui se situe au Brésil, s’étend sur 720 114 km2 avec une densité de population voisine de 2
habitants au km2. Sa végétation se compose principalement de forêts, ainsi que de quelques
prairies et cultures agricoles. Ce bassin possède 38 stations de mesures qui définissent donc
autant de sous-bassins versants (Figure 22).

Figure 22 – Extension des sous bassins versants et M.N.T. du Rio Negro

Pour quatorze de ces sous-bassins, on dispose des mesures de pluie, de débit, et

d’évapotranspiration potentielle au pas de temps journalier, pour des périodes s’étalant sur
plusieurs années. La choix des bassins d’étude s’est effectuée selon un critère temporel, en
retenant le plus grand nombre de sous-bassins qui avaient des périodes de mesures
communes, soit six sous-bassins versants (Tableau 23).

Stations Surface Période commune

Code Nom Rio Date Début Date Fin
14110000 Cucui Rio Negro 71061 km² 01/07/1983 30/05/1990
14220000 Louro Poco Rio Aiari 6174 km² 01/07/1983 30/05/1990
14350000 Jusante Da Cacho Rio Curicuriari 11449 km² 01/07/1983 30/05/1990
14430000 Pari Cachoeira Rio Tiquie 1401 km² 01/07/1983 30/05/1990
14440000 Caracarai Rio Branco 19189 km² 01/07/1983 30/05/1990
14495000 Fazenda Caju Pir Rio Uraricoera 22504 km² 01/08/1983 31/04/1987

Tableau 23 – Caractéristiques des six sous-bassins versants du Rio Negro sélectionnés

- 19 -
 3.2. Modélisation du comportement hydrologiques des six sous-bassins

Pour chacun de ces bassins, nous avons déterminé le M.N.T. avec une précision de
2,28 km à l’aide de ARCVIEW-GIS, ainsi que les distributions des indices topographiques
(Figure 24).

a b c

d e f

Figure 24 – Distribution spatiale de l’indice topographique surles six sous-

bassins du Rio Negro (a : Rio Negro, b : Rio Aiari, c : Rio Curicuriari,
d : Rio Tiquie,e : Rio Branco, f : Rio Uraricoera)

Les bassins du Rio Branco et du Rio Uraricoera présentent des distributions spatiales
de l’indice topographique relativement classiques avec des faibles valeurs sur les versants et
des valeurs plus élevées le long du cours d’eau principal.
Les bassins du RioNegro et du Rio Curicuriari sont des cas intermédiaires avec une
proportion plus importante d’indice topographique élevés, c’est-à-dire, une tendance à une
saturation rapide marquée.
Enfin, les bassins du Rio Aiari et du Rio Tiquie sont composés de plaines, révélées par
des valeurs d’indice topographique élevées, qui recouvrent quasiment l’ensemble du bassin.

Puis, nous avons effectué les calages des différents paramètres.

On obtient les résultats suivants (Figures 25 à 30) :

- 20 -
 1 365 729 1093 1457 1821 2185
0 0,02

0,05 0,016
débits observés

Débits (m/ j )
0,1
Pluie ( m )

débits simulés 0,012

0,15
0,008
0,2

0,25 0,004

0,3 0

Figure 25 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Negro sur une période de 2500 jours
1 365 729 1093 1457 1821 2185
0 0,012

0,1 débits observés 0,009

Débits (m / j )
débits simulés
Pluie ( m )

0,2
0,006
0,3
0,003
0,4

0,5 0
.

Figure 26 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Aiari sur une période de 2500 jours
1 365 729 1093 1457 1821 2185
0 0,025

0,02
0,1 débits observés

Débits ( m / j )
Pluie ( m )

débits simulés 0,015

0,2
0,01
0,3
0,005

0,4 0

Figure 27 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Curicuriari sur une période de 2500 jours
1 365 729 1093 1457 1821 2185
0 0,04

0,1 0,03
Débits ( m / j )

débits observés
Pluie ( m )

débits simulés
0,2 0,02

0,3 0,01

0,4 0

Figure 28 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Tiquie sur une période de 2500 jours

- 21 -
 1 365 729 1093 1457 1821 2185
0

0,1 0,012
débits observés

Débit ( m / j )
Pluie ( m )

0,2 débits simulés

0,008
0,3
0,004
0,4

0,5 0

Figure 29 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Branco sur une période de 2500 jours
1 365 729 1093
0 0,025

0,02
0,05

Débits ( m / j )
débits observés
Pluie ( m )

0,015
0,1 débits simulés
0,01
0,15
0,005

0,2 0

Figure 30 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Uraricoera sur une période de 1400 jours

Pour le Rio Negro (Figure 25), TOPMODEL permet de reproduire les débits de
pointe, néanmoins certaines crues demeurent sous-estimées (t=1093, 1821, 2185). Les montés
de crues et les courbes de récession sont quant à elles correctement reproduites. Cependant,
les étiages ne sont pas aussi convenablement modélisés que pour les deux petits bassins
versants précédemment étudiés.
Pour le Rio Aiari (Figure 26), excepté pour t=1821, les montés de crues et les débits
de crêtes sont correctement reproduits (t=41, 365, 729, 1093, 1457, 2185). Toutefois, les
courbes de décrues et les étiages ne sont pas reproduits convenablement.
Pour le Rio Curucuriari (Figure 27), TOPMODEL permet de reproduire les débits de
pointe, néanmoins certaines crues demeurent sous-estimées (t=1662, 1821, 2185). Les montés
de crues, les courbes de récession et les étiages sont quant à eux correctement reproduits.
Pour le Rio Tiquie (Figure 28), excepté pour t=1457, les débits de crêtes sont
systématiquement sous-estimés (t=365, 729, 1093, 1821). Néanmoins, les temps de montés,
les courbes de décrues et les débits d’étiages sont conformes aux valeurs observées.
Pour le Rio Branco (Figure 29), excepté pour t=365, les débits de pointes sont
correctement reproduits (t=729, 1093, 1457, 1821, 2185). Les montés de crues et les étiages
sont simulés de façon convenable. Cependant, les formes des courbes de récession simulées
ne correspondent pas à celles observées.
Enfin, pour le Rio Uraricoera (Figure 30), les montées, les récessions, les étiages et les
débits de pointes sont correctement reproduits (t=365, 729, 1093).

Toutefois, comme pour le Strengbach et le Nsimi, le modèle ne génère pas de crue

factice et répond aux impulsions de pluie. Ceci conduit à une efficience moyenne de l’ordre
de 0.65 pour les Rio Negro, Curucuriari, Tiquie et Branco, de l’ordre de 0.55 pour le Rio
Aiari et de l’ordre de 0.8 pour le Rio Uraricoera, ce qui indique que TOPMODEL peut être

- 22 -
tout à fait utilisé pour ce type de bassins. La basse valeur d’efficience obtenue pour le Rio
Aiari (eff=0.53), s’explique par la présence de pentes faibles. La distribution de l’indice
topographique pour ce bassin (Figure 24) révèle en effet que ce bassin est quasiment plat.
Néanmoins, les différentes valeurs d’efficience obtenues sont très encourageantes, vu
que nous avons largement dépassé la précision minimale recommandé, à savoir 50 m pour les
M.N.T., et que nous avons supposé les pluies uniformément réparties sur le bassin, ce qui
n’est pas le cas puisque nous avons une seule station de mesure par exutoire de bassin.

3.3. Etude de la variation des paramètres pour le Rio Negro

Dans une deuxième partie, afin d’étudier une possible variation des paramètres de
calage au cours du temps, nous allons successivement appliquer TOPMODEL, pour des
périodes plus courtes : pour environ deux années (Figure 31) et à l’échelle d’une crue (Figure
32) pour le Rio Negro, à l’échelle d’une crue pour les Rio Aiari (Figure 33), Branco (Figure
37) et Uraricoera (Figure 38), pour une année (Figure 34) et à l’échelle d’une crue (Figure 35)
pour le Rio Curucuriari et pour environ un an et demi (Figure 36) pour le Rio Tiquie. Les
nouvelles simulations et valeurs obtenues sont regroupées dans le Tableau 39 :

1 101 201 301 401 501 601

0 0,02

0,04 0,016

Débits ( m / j )
Pluie ( m )

0,08 débits observés 0,012

débits simulés
0,12 0,008

0,16 0,004

0,2 0

Figure 31 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Negro sur une période de 665 jours

1 31 61 91 121
0 0,02
0,02 0,018
0,04 0,016
débits observés
Débits ( m / j )

0,06 0,014
Pluie ( m )

0,08 débits simulés 0,012

0,1 0,01
0,12 0,008
0,14 0,006
0,16 0,004
0,18 0,002
0,2 0

Figure 32 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Negro sur une période de 149 jours

- 23 -
 1 31 61 91
0 0,009

0,08 0,008
débits observés

Débits ( m / j )
débits simulés 0,007
Pluie ( m )

0,16
0,006
0,24
0,005
0,32 0,004

0,4 0,003

Figure 33 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Aiari sur une période de 104 jours
1 101 201 301
0 0,016

0,012

Débi ( m / j )
0,1
Pluie ( m )

0,008

0,2
0,004
débits observés
débits simulés
0,3 0

Figure 34 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Curicuriari sur une période de 360 jours
1 31 61 91 121
0 0,018

0,04
0,014

Débits ( m / j )
Pluie ( m )

0,08
0,01
0,12
0,006
0,16 débits observés
débits simulés
0,2 0,002

Figure 35 – Comparaison entre débits observés et simulés à

l’exutoire du Curicuriari sur une période de 135 jours
1 101 201 301 401 501
0 0,03

débits observés
Débits ( m / j )

0,1 debits simulés 0,02

Pluie ( m )

0,2 0,01

0,3 0

Figure 36 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Tiquie sur une période de 226 jours

- 24 -
 1 31 61 91 121
0 0,01

0,04 débits observés 0,008

débits simulés

Débit ( m / j )
Pluie ( m )

0,08 0,006

0,12 0,004

0,16 0,002

0,2 0

Figure 37 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Branco sur une période de 141 jours
1 31 61 91 121 151 181 211
0 0,016

0,05
0,012
débits observés

Débit ( m / j )
Pluie ( m )

0,1 débits simulés

0,008
0,15
0,004
0,2

0,25 0

Figure 38 – Comparaison entre débits observés et simulés à

l’exutoire du Rio Uraricoera sur une période de 226 jours
Paramètre Rio Negro Rio Curicuriari
Période 2500 j 665 j 149 j 2500 j 360 j 135 j
M 0.186 m 0.142 m 0.02 m 0.18 m 0.14 m 0.14 m
Ln(To) 12.6 m²/h 12 m²/h 9.46 m²/h 35 m²/h 58 m²/h 61 m²/h
SRMax 3.10-3 m 3.10-3 m 10-3 m 2.2 10-3 m 3.10-3 m 10-3 m
SRinit 0.126 m 0m - 0.13 m 0m - 0.028 m - 0.065 m
ChVel 1770 m/h 1500 m/h 494 m/h 2500 m/h 10000 m/h 9000 m/h
Efficience 0.66 0.75 0.93 0.63 0.8 0.86
Paramètre Rio Aiari Rio Tiquie Rio Branco Rio Uraricoera
Période 2500 j 104 j 2500 j 566 j 2500 j 141 j 1400 j 226 j
M 0.455 m 0.483 m 0.195 m 0.22 m 0.12 m 0.3 m 0.12 m 0.1 m
Ln(To) 11.6 m²/h 12.2 m²/h 13.6 m²/h 12.7 m²/h 8.8 m²/h 6.36 m²/h 10.7 m²/h 10.5 m²/h
SRMax 0.068 m 6.10-3 m 10-4 m 10-6 m 0.075 m 0.17 m 2.10-3 m 2.7 10-3 m
SRinit 8.514 m 0.5 m 0m - 0.25 m 0.5 m - 0.051 m 0m 0.08 m
ChVel 1570 m/h 3700 m/h 452 m/h 198 m/h 2300 m/h 1275 m/h 3400 m/h 3500 m/h
Efficience 0.53 0.79 0.65 0.70 0.64 0.83 0.81 0.88

Tableau 39 – Variation des paramètres de calage pour les six sous-bassins versants
du Rio Negro pour des périodes d’étude de plus en plus courtes

De façon générale, pour des périodes d’étude plus courtes, la qualité des simulations
augmentent de manière significative, il en résulte des valeurs d’efficience plus élevées de
l’ordre de : 0.75 (Figure 31) et 0.95 (Figure 32) pour le Rio Negro, 0.8 (Figure 33) pour le Rio
Aiari, 0.8 (Figure 34) et 0.85 (Figure 35) pour le Rio Curucuriari, 0.7 (Figure 36) pour le Rio
Tiquie, 0.85 (Figure 37) pour le Rio Branco et 0.9 (Figure 38) pour le Rio Uraricoera.

- 25 -
 Globalement, les valeurs de M et Ln(To) restent du même ordre de grandeur sur
chaque bassin, notons toutefois que les valeurs de Ln(To) sont très élevées pour le Rio
Curucuriari. Néanmoins, les valeurs de SRMax varient entre un et deux ordres de grandeur,
pour les bassins les Rio Aiari, Tiquie et Branco, et celles de ChVel d’un ordre de grandeur
pour le Rio Curucuriari. De plus, les valeurs négatives de SRinit, probablement sensées
modéliser le ruissellement en début d’épisode pluvieux, sont peu physiques. Cependant, ces
valeurs négatives n’apparaissent que pour de courts épisodes et non pas pour les périodes
d’études globales. Ces variations sont vraisemblablement dues à la faible précision des
M.N.T. et à la présence d’une seule station de mesure par bassin.

3.4. Extension des surfaces saturées

Nous allons à présent, nous intéresser à l’ extension des surfaces saturées pour ces
différents bassins (Figure 40) :

4,5 12
4 10
Pourcentage
Pourcentage

3,5 8
3
6
2,5
2 4
1,5 2
1 0
0 500 1000 1500 2000 2500 0 500 1000 1500 2000 2500

a b
6 21
5
Pourcentage

Pourcentage

4 14
3
2 7
1
0 0
0 500 1000 1500 2000 2500 0 500 1000 1500 2000 2500

c d
14
12
Figure 40 - Évolution du pourcentage de
surfaces saturées en fonction du temps (a : Rio
Pourcentage

10
8 Negro (2.4% en moyenne), b : Rio Aiari (4.6%
6 en moyenne), c : Rio Tiquie (2.7% en
4 moyenne), d : Rio Branco (11.4% en
2 moyenne), e : Rio Uraricoera (6.4% en
0 350 700 1050 1400 moyenne))
e

La modélisation de l’évolution des surfaces saturées est satisfaisante pour cinq des six
bassins. En effet, pour le Rio Curucuriari, la lame d’eau ruisselée, et à fortiori le pourcentage
de surfaces saturées est nulle à chaque pas de temps. Ce résultat s’explique par les très fortes
valeurs de Ln(To) pour ce bassin (Tableau 35). Rappelons que To est proportionnelle à la
conductivité à la surface du sol, autrement dit lorsque les valeurs de To deviennent trop
importantes toute la pluie s’infiltre et il ne peut y avoir de ruissellement. Enfin, le fort
pourcentage moyen trouvé pour le Rio Branco (11,4%) peut être à l’origine des difficultés
rencontrées par TOPMODEL pour modéliser les différentes récessions pour ce bassin (Figure
29).

- 26 -
CONCLUSION

Cette étude exploratoire se proposait de vérifier si TOPMODEL permettait de prédire

le comportement hydrologique de trois bassins versants de nature et de tailles différentes,
ainsi que d’illustrer ses principales caractéristiques et limitations.

La première partie nous a donné l’occasion de rassembler puis de mettre en place tous
les outils nécessaires à la chaîne de prévision. Après avoir soigneusement étudié le modèle
hydrologique TOPMODEL afin de modéliser la transformation des pluies en débits, nous
avons présenté le programme GRIDATB et étudié son comportement sur quelques bassins
versants synthétiques.
Les deux parties suivantes, plus exploratoires, ont été l’occasion de mettre en œuvre et
d’évaluer la méthode d’étude choisie. Nous avons appliqué TOPMODEL tout d’abord sur un
petit bassin versant de moyenne montagne, puis sur un petit bassin versant tropical et enfin
sur différents grands bassins versants tropicaux.

Pour le bassin versant du Strengbach, toutes les simulations ont été d’un très bon
niveau. Les valeurs des paramètres sont toujours restés physiques et du même ordre de
grandeur, le pourcentage moyen de surface saturée trouvé, conforme aux données
expérimentales. Ce premier bassin fut l’occasion de s’assurer que TOPMODEL donnait de
bon résultats sur un bassin de nature et de taille en conformirté avec les conditions
d’application définies par Beven (1995 et 1997).
TOPMODEL a correctement modéliser le comportement hydrologique du Nsimi, mais
a présenter quelques difficultés pour certains débits d’étiages. Les paramètres sont restés
physiques et seul la valeur de la profondeur de la zone racinaire (SRMax) a présenté des
variations significatives. Néanmoins, nous n’avons pas obtenu une bonne estimation de
l’évolution des surfaces saturées. Nous avons montré que ce problème venait de la quasi
nullité des pentes, de la présence d’une zone marécageuse et de l’éventuel non respect de
l’hypothèse de parallélisme la pente de la zone saturée à la surface du sol.
Enfin, pour le Rio Negro, nous avons également obtenu des simulations satisfaisantes,
toutefois le modèle a montrer ses limites, notamment pour certaines récessions. En ce qui
concerne les paramètres, les valeurs de la profondeur de la zone racinaire (SRMax) et de la
vitesse caractéristique de la fonction de transfert (ChVel) ont présenté des variations. La
transmissivité à saturation (Ln(To)), bien qu’ayant resté constante, a présenté des valeurs très
élevées pour un des bassins. L’étude de l’évolution des surfaces saturées a révélé qu’un fort
pourcentage moyen pouvait être la cause des problèmes rencontrés par le modèle pour
reproduire les récessions, et qu’un forte valeur de To, pouvait être à l’origine d’une lame d’eau
ruisselée et donc d’un pourcentage, nul à chaque pas de temps. Néanmoins, les différentes
valeurs d’efficience obtenues sont très encourageantes, vu que nous avons largement dépassé
la précision minimale recommandé, à savoir 50 m pour les M.N.T., et que nous avons supposé
les pluies uniformément réparties sur le bassin, ce qui n’était pas le cas puisque nous avions
qu’une seule station de mesure par exutoire de bassin.

Notre étude nous a permis de montrer l’universalité et la grande robustesse de

TOPMODEL. Initialement conçu pour des bassins de moyenne montagne humide, les limites

- 27 -
de TOPMODEL ont été allégrement franchies vu la nature et la taille des bassins versants
étudiés. En effet, les données disponibles ne nous permettaient pour certains bassins, ni de
satisfaire la valeur minimale de précision du M.N.T. , ni de respecter le pas de temps horaire.
Cependant, nous avons réussi à modéliser le comportement hydrologique des trois bassins
versants de manière tout à fait satisfaisante.

Différentes perspectives sont envisageables pour la suite de cette étude.

Tout d’abord, il serait avisé de disposer de mesures des surfaces saturées au pas de
temps journalier. Ces mesures permettraient de vérifier les modélisations de l’évolution des
surfaces saturées, et surtout de valider des relations entre ces surfaces et les débits mesurés.
De même, il serait intéressant de vérifier si d’autres modèles permettent de simuler le
comportement hydrologique des trois bassins versants étudiés. Le modèle S.H.E. (Système
Hydrologique Européen) qui constitue sans doute le modèle à base physique actuellement le
plus complet pourrait être testé. Ce qui permettrait de comparer les différentes reconstitutions
et de choisir le modèle le plus adapté.
Enfin, il serait judicieux de modifier les codes FORTRAN originaux de GRIDATB et
de TOPMODEL afin de prendre en compte les caractéristiques des différents bassins.
Notamment l’introduction de certains paramètres pourrait permettre de prendre en compte le
non respect de l’hypothèse de parallélisme de la pente de la zone saturée à la surface du sol
ainsi que les crues dues aux fontes de neiges.

- 28 -
BIBLIOGRAPHIE

Ambroise B., Beven K & Freer J., Application of a generalized TOPMODEL to the small Ringelbach
catchment, Vosges, France, Water Resources Research, Vol. 32, No. 7, Pages 2147-2159 (1996).

Beven K., TOPMODEL: a critique, Hydrological Processes, Vol. 11, 1069-1085 (1997).

Beven K., DTM9501 and GRIDATB : Digital Terrain Analysis

A users guide to the distribution versions (95.01), Centre for Research on Environmental Systems and
Statistics-Institute of Environmental and Biological Sciences, Lancaster University (1995).

Beven K. & Kirkby M.J., A physically based variable contributing area model of basin hydrology, Hydrol. Sci.
Bull., 24, 43-69 (1979).

Beven K. & Germann P., Macropores and water flow in soils, W.R.R., 18, 1311-1325 (1982).

Beven K., Quinn P., Romanowicz R., Freer J., Fisher J. & Lamb R., TOPMODEL
A users guide to the distribution version for DOS (95.02), Centre for Research on Environmental
Systems and Statistics-Institute of Environmental and Biological Sciences, Lancaster University (1995).

Beven K., Quinn P., Romanowicz R., Freer J. & Lamb R., TOPMODEL, in Singh, V.P. (Ed), Computer
Models of Watershed Hydrology, Water Resources Publications, Highlands Ranch, Colo., 627-668
(1995).

Bodin J., Fonctionnement hydrodynamique d’une nappe de bassin versant en zone tropicale humide :
modélisation du site expérimental de Nsimi au Sud – Cameroun, Institut Français de Recherche
Scientifique pour le Développement en Coopération (1998).

Brendle T., Caratérisation hydrodynamique des formations superficielles du basin versant du Strengbach
(Aubure, Haut-Rhin), Mémoire de maîtrise, Université Strasbourg I (1992).

Campling P., Gobin A., Beven K. & Feyen J., Rainfall-runoff modelling of a humid tropical catchment: the
TOPMODEL approach, Hydrological processes, Vol.16, 231-253 (20002).

Cappus P., Etude des lois de l’écoulement, application au calcul et à la prévision des débits, La HouilleBlanche
(1960).

Cosandey C. & Robinson M., Hydrologie continentale, Collection U-Géographie, Armand Colin (2000).

Datin R., Outils opérationnels pour la prévision des crues rapides : traitement des incertitudes et intégration des
prévisions météorologiques. Développements de TOPMODEL pour la prise en compte de la variabilité
spatiale de la pluie. Application au bassin versant de l’Ardèche, Thèse de Doctorat, INP Grenoble
(1998).

Deg-Hyo B., Jin-Hoon K. & Miller N. L., An extended TOPMODEL for flow forecasting in large drainage
basins, 4th International Conference on Hydro-Science and Engineering (ICHE), Seoul (2000).

Durand P., Robson A. & Neal C., Modelling the hydrology of submediterranean mountain catchments (Mont-
Lozère) using TOPMODEL: initial results, Journal of Hydrology, 1-14 (1992).

Franchini M., Wendling J, Obled Ch. & Todini E., Physical interpretation and sensitivity analysis of the
TOPMODEL, Journal of Hydrology, 175, 293 – 338 (1996).

Holko L. & Lepistö A., Modelling the hydrological behaviour of a mountain catchment using TOPMODEL,
Journal of Hydrology, 361-377 (1997).

Horton J.H., Erosional development of streams and their drainage basins, .Bull. Geol. Soc. Amer., 56, 275 –370
(1945).

- 29 -
Iorgulescu I. & Jordan J.-P., Validation of TOPMODEL on a small Swiss catchment, Journal of Hydrology,
225 - 273 (1994).

Iorgulescu I. & Musy A., Generalization of TOPMODEL for a power law transmissivity profile, Hydrological
Processes, 11, 1353 - 1355 (1997).

Kirkby M. J., TOPMODEL : a personal view, Hydrological Processes, Vol. 11, 1087 – 1097 (1997).

Ladouche B., Etude des flux hydriques par le traçage isotopique naturel à l’échelle d’un bassin versant forestier
(Strengbach, Vosges), Thèse de Doctorat, Université Paris VI (1997).

Laroche Ch., Comparaison de deux petits bassins de la zone méditerranéenne. Test de sensibilité du pas
d’espace du MNT sur la conductivité hydraulique à saturation, Mémoire de DEA, UJF Grenoble (1994).

Latron J., Caractérisation géomorphologique et hydrologique du bassin versant du Strengbach (Aubure),

Mémoire de maîtrise, Université Strasbourg I (1990).

Mendicino G. & Sole A., The information content theory for the estimation of the topographic index
distribution used in TOPMODEL, Hydrological Processes, Vol. 11, 1099 - 1114 (1997).

Musy A., Hydrologie appliquée, Editions HGA, Bucarest (1998).

Nash J. E. & Sutcliffe J. V., River flow forecasting through conceptual models. Part I – A discussion of
principles, Journal of Hydrology, 10, 282 – 290 (1970).

Nawarathna B., Ao T., Kazama S. Sawamoto M. & Takeuchi K., Influence of human activities on the
BTO¨MC model runoff simulations in large-scale watersheds, XXIX IAHR Congress Proceedings,
Theme A, 93 – 99 (2001).

Nwarathna B., Kazama S. & Sawamoto M., Improvement of calibration procedure of the block wise
TOPMODEL with Muskingum-Cunge routing method using sub basins simulated results : application
to part of the lower Mekong River basin, Proceeding 13th Congress ADP/IAHR, Vol.1, 540 – 545
(2002).

Ndam Ngoupayou J.R., Bilans hydrogéochimiques sous forêts tropicale humide en Afrique : du bassin
expérimental de Nsimi – Zoételé aux réseaux hydrographiques du Nyong et de la Sanaga (Sud –
Cameroun), Thèse de Doctorat, Université Paris VI (1997).

Quinn P.F., Beven K.J. & Lamb R., The ln(a/tanβ) index; how to calculate it and how to use it within the
TOPMODEL framework, Centre for Research on Environmental Systems and Statistics, Lancaster
University (1991).

Robson A.J., Whitehead P.G. & Johson R.C., An application of a physically based semi-distributed model to
the Balquhidder catchments, Journal of Hydrology, 357-370 (1993).

Saulnier G.M., Information pédologique spatialisée et traitement topographiques améliorés dans la modélisation
hydrologique par TOPMODEL, Thèse de Doctorat, INP Grenoble (1996).

Takeuchi K., Ao T. & Ishidaira H., Introduction of block-wise use of TOPMODEL and Muskingum-Cunge
method for the hydro-environmental simulation of a large ungauged basin, Hydrological Sciences, 44
(1999).

Tianqi A., Yoshitani J., Takeuchi K., Fukami K., Mutsuura T. & Ishidaira H., Toward the application of the
physically based distributed hydrological model BTOPMC to ungauged basins, Meeting – IAHS –
University of Brasilia, (20 – 22 November 2002).

Wendling J., Modélisation pluie-débit : comparaison d’approches conceptuelles/physico-déterministes,

globales/semi-distribuées. Essai de prise en compte de la variabilité spatiale des pluies (Applications au
bassin versant du Réal Collobrier), Thèse de Doctorat, INP Grenoble (1992).

- 30 -