Cours D'hydrologie Générale
Cours D'hydrologie Générale
COURS D’HYDROLOGIE
GENERALE
Professeur André Musy
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU
Informations générales
Recommandations
3. Répondre soi-même aux questions posées, vérifier les réponses qui se trouvent
toutes dans le texte présenté sur Internet.
5. Ne pas hésiter à utiliser le forum. Tenter d'être clair, précis et succinct dans les
questions posées tout comme dans les réponses aux commentaires que vous pouvez
formuler.
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU II
SOMMAIRE
CHAPITRE 1 : CYCLE ET BILAN HYDROLOGIQUES ...................................... 9
1.1 Introduction .....................................................................................................................9
1.2 L'eau, généralités............................................................................................................9
1.3 Définition et composantes du cycle hydrologique .................................................10
1.3.1 Définition .....................................................................................................................10
1.3.2 Les précipitations......................................................................................................11
1.3.3 L'évaporation/l'évapotranspiration.....................................................................12
1.3.4 L'interception et le stockage dans les dépressions ........................................12
1.3.5 L'infiltration et la percolation................................................................................13
1.3.6 Les écoulements ........................................................................................................13
1.4 La répartition des eaux...................................................................................................13
1.4.1 A l'échelle du globe ...................................................................................................14
1.4.2 A l'échelle des continents .......................................................................................17
1.4.3 A l'échelle d'une zone géographique : la Suisse..............................................17
1.5 Le bilan hydrique..............................................................................................................18
CHAPITRE 2 : LE BASSIN VERSANT ET SON COMPLEXE............................ 21
2.1 Définition du bassin versant.........................................................................................21
2.2 Comportement hydrologique.........................................................................................23
2.2.1 Le temps de concentration.....................................................................................24
2.2.2 Les courbes isochrones...........................................................................................24
2.3. Caractéristiques physiques et leurs influences sur l'écoulement des eaux.25
2.3.1 Les caractéristiques géométriques ......................................................................25
[Link] La surface ............................................................................................................25
[Link] La forme ...............................................................................................................26
[Link] Le relief .................................................................................................................27
2.3.2 Le réseau hydrographique .....................................................................................30
[Link] La topologie : structure du réseau et ordre des cours d'eau ..............31
[Link] Les longueurs et les pentes caractéristiques du réseau ......................32
[Link] Le Degré de développement du réseau ......................................................34
[Link] L'endoréisme.......................................................................................................37
2.3.3 Les caractéristiques agro-pédo-géologiques.....................................................37
[Link] La couverture du sol ........................................................................................37
[Link] La nature du sol ................................................................................................39
[Link] La géologie du substratum ............................................................................41
2.4 Informations digitales et modèles numériques.......................................................42
2.4.1 Généralités sur les MNA et MNT..........................................................................43
2.4.2 En Suisse.....................................................................................................................43
CHAPITRE 3 LES PRECIPITATIONS ........................................................... 44
3.1 Principes météorologiques .............................................................................................44
3.1.1 Définition des précipitations .................................................................................44
3.1.2 Les nuages ..................................................................................................................44
3.1.3 Mécanismes de formation des précipitations ..................................................45
3.1.4 Types de précipitations ...........................................................................................46
3.1.5 Régime des précipitations ......................................................................................47
3.2 Mesures des précipitations............................................................................................48
3.2.1 Mesures de la hauteur d'eau précipitée ............................................................48
3.2.2 Réseau d'observation et publication des données .........................................49
[Link] Le réseau d'observation ..................................................................................49
[Link] Publication des données pluviométriques ................................................49
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU III
3.3 Analyse ponctuelle ...........................................................................................................50
3.3.1 Notion d'averses et d'intensités............................................................................50
3.3.2 Statistique descriptive des séries chronologiques..........................................52
3.3.3 Notion de temps de retour .....................................................................................53
3.3.4 Les courbes IDF (intensité – durée – fréquence).............................................54
[Link] Lois de pluviosité...............................................................................................54
[Link] Utilisation des courbes IDF ...........................................................................54
[Link] Construction de courbes IDF ........................................................................54
3.3.5 La structure des pluies ...........................................................................................57
3.4 Evaluation régionale des précipitations ....................................................................57
3.4.1 Passage des pluies ponctuelles aux pluies moyennes sur une surface..58
[Link] Calcul de la moyenne arithmétique ............................................................58
[Link] Calcul de la moyenne pondérée - méthode des polygones de
Thiessen ..............................................................................................................................58
[Link] La méthode des isohyètes (isovaleurs) .......................................................59
3.4.2 Le hyétogramme moyen ..........................................................................................60
3.4.3 Notion d'abattement des pluies............................................................................61
CHAPITRE 4 : EVAPORATION ET INTERCEPTION ....................................... 63
4.1 Introduction........................................................................................................................63
4.1.1 L'interception ...............................................................................................................63
4.1.2 L'évaporation et la transpiration..........................................................................64
4.2 L'évaporation - L'évapotranspiration..........................................................................65
4.2.1 Processus physique de l'évaporation..................................................................65
[Link] Description et formulation du processus physique...............................65
[Link] Facteurs météorologiques intervenant dans le processus
d'évaporation .....................................................................................................................66
[Link] Facteurs physiques du milieu intervenant dans le processus
d'évaporation .....................................................................................................................71
[Link] Estimation de l'évaporation des nappes d'eau libre ..............................73
4.2.2 Evapotranspiration d'un sol couvert par de la végétation ..........................74
[Link] Rappel sur les processus physiques de la transpiration des
végétaux ..............................................................................................................................74
[Link] Notions d'évapotranspiration de référence, maximale et réelle .........75
[Link] Facteurs intervenant dans le processus d'évapotranspiration ..........76
[Link] Evaluation de l'évapotranspiration .............................................................77
4.3 L'interception .....................................................................................................................80
4.3.1 Facteurs météorologiques intervenant dans le processus d'interception
....................................................................................................................................................81
4.3.2 Facteurs végétatifs intervenant dans le processus d'interception ..............81
4.3.3 Quelques limitations et ordre de grandeur du processus de l'interception
....................................................................................................................................................82
CHAPITRE 5 : L'INFILTRATION ET LES ECOULEMENTS ............................. 85
5.1 Introduction........................................................................................................................85
5.2 L'infiltration ........................................................................................................................85
5.2.1 Définitions et paramètres descriptifs de l'infiltration ...................................85
5.2.2 Facteurs influençant l'infiltration........................................................................87
5.2.3 Variation du taux d'infiltration au cours d'une averse ................................88
5.2.4 Modélisation du processus d'infiltration...........................................................90
[Link] Relations empiriques .......................................................................................90
[Link] Modèles à base physique ................................................................................91
5.3 Les écoulements................................................................................................................93
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU IV
5.3.1 Généralités ..................................................................................................................93
5.3.2 L'écoulement de surface .........................................................................................95
5.3.3 L'écoulement de subsurface ..................................................................................96
5.3.4 L'écoulement souterrain .........................................................................................96
5.3.5 Ecoulement dû à la fonte des neiges..................................................................97
5.3.6 Bilan annuel des écoulements .............................................................................97
5.3.7 Introduction au transport solide (dans les cours d'eau)..............................98
[Link] Transport solides dans les cours d'eau .....................................................98
[Link] Notions de transport spécifique et érosion mécanique sur un bassin
versant .................................................................................................................................99
[Link] Distribution du transport spécifique à travers le monde ..................100
CHAPITRE 6 : LE STOCKAGE ET SES VARIATION .................................... 101
6.1 Introduction......................................................................................................................101
6.2 Les stocks d'eau de surface.........................................................................................101
6.2.1 Quelques définitions ..............................................................................................101
6.2.2 Introduction à la limnologie et caractéristiques générales des lacs ......102
6.3 Les stocks d'eau souterraine ......................................................................................104
6.3.1 Distinction zone saturée et zone non saturée ...............................................104
6.3.2 L'eau du sol...............................................................................................................105
[Link] Caractéristiques de la phase liquide du sol ...........................................105
[Link] Etat énergétique de l'eau dans le sol ........................................................106
[Link] Comportement dynamique : la loi de Darcy...........................................107
[Link] Calcul du stock d'eau ....................................................................................108
6.3.3 L'eau du sous-sol ou souterraine ......................................................................109
[Link] Définitions : aquifères et type de nappes ................................................109
[Link] Définitions et mesure de la surface piézométrique..............................111
[Link] Caractéristiques principales de l'aquifère...............................................111
[Link] Vitesses d'écoulement réelle et fictive, débit d'une nappe.................112
[Link] Calcul du stock d'eau ....................................................................................113
[Link] Concept de tarissement des nappes .........................................................113
6.4 Les stocks d'eau sous forme solide...........................................................................114
6.4.1 La couverture neigeuse .........................................................................................114
[Link] Evaluation du stock neigeux.......................................................................115
[Link] Ecoulement de l'eau à l'intérieur du massif de neige .........................115
[Link] La fonte de la neige.........................................................................................116
6.4.2 La couverture glaciaire..........................................................................................117
[Link] Les glaciers........................................................................................................117
[Link] La glace recouvrant les lacs et les rivières..............................................118
CHAPITRE 7 : LA MESURE HYDROLOGIQUE............................................. 119
7.1 Introduction à la mesure des précipitations..........................................................119
7.1.1 Les pluviomètres .....................................................................................................119
7.1.2 Les pluviographes ...................................................................................................120
[Link] Le pluviographe à siphon .............................................................................121
[Link] Le pluviographe à augets basculeurs.......................................................121
7.1.3 Le radar ......................................................................................................................122
7.1.4 Les erreurs de mesure...........................................................................................122
7.2 La mesure de l'évaporation, de la transpiration et de l'évapotranspiration 123
7.2.1 Facteurs influençant la mesure de l'évaporation .........................................123
[Link] Le rayonnement solaire et la durée d'insolation...................................124
[Link] La température ................................................................................................124
[Link] L'humidité de l'air ...........................................................................................125
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU V
[Link] La pression atmosphérique..........................................................................125
[Link] Le vent ................................................................................................................125
7.2.2 Mesure de l'évaporation des nappes d'eau libre...........................................126
[Link] Les évaporimètres ...........................................................................................126
[Link] Les balances d'évaporation ..........................................................................126
[Link] Les bacs d'évaporation ..................................................................................126
7.2.3 Mesure de l'évaporation à partir des sols nus ..............................................127
[Link] Les verrières ou châssis vitrés....................................................................127
[Link] Le lysimètre .......................................................................................................127
7.2.4 Mesure de l'évapotranspiration..........................................................................127
7.3 La mesure des débits.....................................................................................................128
7.3.1 La mesure des hauteurs d'eau ...........................................................................129
[Link] Le limnimètre....................................................................................................129
[Link] Le limnigraphe à flotteur ..............................................................................130
[Link] Le limnigraphe "bulle à bulle".....................................................................130
[Link] Autres capteurs pour la mesure des hauteurs d'eau..........................130
7.3.2 La mesure des débits.............................................................................................131
[Link] Le jaugeage par exploration du champ de vitesse ...............................131
[Link] La détermination du débit à l'aide d'ouvrages calibrés ......................135
[Link] Les jaugeages par dilution ...........................................................................136
7.4 La mesure du transport solide dans les cours d'eau..........................................138
7.4.1 La mesure du transport en suspension ..........................................................139
[Link] Matériel de prélèvements..............................................................................139
[Link] Calcul du débit solide en suspension ......................................................140
[Link] Mesures en continu........................................................................................141
7.4.2 La mesure du transport de fond ........................................................................141
7.5 La mesure de l'infiltration............................................................................................142
7.6 La mesure de l'humidité du sol..................................................................................143
7.6.1 La sonde à neutrons ..............................................................................................144
7.6.2 La technique TDR (Time Domain Reflectometry) .........................................144
CHAPITRE 8 : L'ORGANISATION ET LE CONTROLE DES DONNEES ........... 146
8.1 L'organisation des données .........................................................................................146
8.1.1 Acquisition des données .......................................................................................146
8.1.2 Traitement primaire des données......................................................................146
8.1.3 Contrôle des données ............................................................................................147
8.1.4 Organisation des données....................................................................................147
8.1.5 Diffusion et publication ........................................................................................147
8.2 Le contrôle des données ...............................................................................................148
8.3 Recherche des erreurs et corrections des mesures.............................................149
Hypothèses de l'analyse statistique :...........................................................................149
8.4 Introduction aux tests statistiques...........................................................................150
8.4.1 Démarche d'un test : exemple introductif - Les faiseurs de pluie ..........150
8.4.2 Les grandes catégories de tests..........................................................................153
[Link] Tests selon leurs propriétés mathématiques .........................................153
[Link] Tests selon leur objet .....................................................................................153
[Link] Tests selon la nature des informations....................................................155
8.5 Contrôle des données : application au débit de la Viege...................................155
8.5.1 1 paramètre – échelle locale ................................................................................157
[Link] Tests paramétriques.......................................................................................157
[Link] Tests non paramétriques..............................................................................162
8.5.2 1 paramètre – échelle régionale..........................................................................164
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU VI
[Link] Méthode du double cumul ...........................................................................164
[Link] Méthode du cumul des résidus ..................................................................165
8.5.3 Plusieurs paramètres – échelle locale ..............................................................167
8.5.4 Plusieurs paramètres – échelle régionale........................................................167
[Link] Méthode des bilans.........................................................................................167
[Link] Méthode des débits spécifiques maximaux ............................................168
[Link] Méthodes statistiques multivariées...........................................................168
8.5.5 Estimation des données manquantes, correction des données ..............169
CHAPITRE 9 : LES REGIMES HYDROLOGIQUES ....................................... 171
9.1 Introduction aux régimes hydrologiques ................................................................171
9.2 Classification des régimes hydrologiques...............................................................173
9.2.1 Le régime simple .....................................................................................................174
[Link] Le régime glaciaire ..........................................................................................174
[Link] Le régime nival .................................................................................................175
[Link] Le régime pluvial pur (ou océanique) .......................................................176
[Link] Le régime pluvial tropical .............................................................................176
9.2.2 Le régime mixte .......................................................................................................177
[Link] Le régime nivo-glaciaire ................................................................................177
[Link] Le régime nivo-pluvial....................................................................................177
[Link] Le régime pluvio-nival....................................................................................178
9.2.3 Le régime complexe ................................................................................................179
9.3 Les régimes hydrologiques en Suisse ......................................................................179
CHAPITRE 10 : LES PROCESSUS HYDROLOGIQUES ................................. 181
10.1 Introduction ...................................................................................................................181
10.2 Les différents types d'écoulements.........................................................................181
10.3 Les précipitations directes à la surface libre du cours d'eau ........................184
10.4 Les écoulements de surface......................................................................................184
10.4.1 Ecoulement par dépassement de la capacité d'infiltration.....................184
10.4.2 Ecoulement sur surfaces saturées .................................................................185
10.4.3 Localisation des surfaces saturées.................................................................186
10.5 Les écoulements de subsurface...............................................................................186
10.5.1 Effet piston .............................................................................................................187
10.5.2 Ecoulement par macropores.............................................................................187
[Link] Définition d'un macropore.........................................................................188
[Link] Types et origines des macropores ...........................................................188
[Link] Détermination expérimentale et dynamique .......................................189
[Link] Relation entre macroporosité et infiltration.........................................189
[Link] Relations entre macropores et écoulement de subsurface .............191
[Link] Ecoulement par "tubes" éphémères et pérennes................................191
10.5.3 Intumescence de la nappe.................................................................................191
10.5.4 Ecoulement de retour .........................................................................................192
10.6 Les écoulements souterrains....................................................................................193
10.6.1 Processus général.................................................................................................193
10.6.2 Introduction à l'utilisation des traceurs isotopiques................................194
[Link] Conditions d'utilisation ..............................................................................195
[Link] Variations de la composition isotopique moyenne des eaux..........195
[Link] Exemple d'application : modèle à deux composantes ......................196
[Link] Autres types de traceurs ............................................................................197
[Link] Modèle à trois composantes, la méthode EMMA ...............................198
10.7 Essai de synthèse relative à la notion de processus de génération de
l'écoulement .............................................................................................................................199
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU VII
CHAPITRE 11 : LA REPONSE HYDROLOGIQUE ......................................... 202
11.1 Introduction ...................................................................................................................202
11.1.1 Réponses hydrologiques.....................................................................................202
11.1.2 Transformation de la pluie en hydrogramme de crue..............................203
11.2 Analyse des événements pluies-débits..................................................................203
11.3 Genèse des crues..........................................................................................................205
11.4 Facteurs d'influence de la réponse hydrologique ..............................................205
11.4.1 Facteurs liés à la pluviosité ..............................................................................206
[Link] Influence de la durée de l'averse .............................................................207
[Link] Influence de la distribution spatiale.......................................................208
[Link] Influence des variations d'intensité en fonction du temps .............208
11.4.2 Importance des conditions antécédentes d'humidité ..................................209
11.5 Conclusion.......................................................................................................................210
ANNEXE 1 : L'ANALYSE FREQUENTIELLE................................................ 211
A.1 Généralités........................................................................................................................211
A.1.1 Définition de l'analyse fréquentielle .................................................................211
A.2 Choix du modèle fréquentiel.......................................................................................212
A.2.1 Considérations théoriques...................................................................................212
A.2.1.1 Loi normale .......................................................................................................212
A.2.1.2 Loi log-normale ................................................................................................212
A.2.1.3 Loi de Gumbel ..................................................................................................212
A.2.2 Comportement asymptotique .............................................................................213
A.2. 3 L'expérience et la coutume.................................................................................213
A.2.4 Utilisation des tests d'adéquation.....................................................................213
A.2.5 Utilisation de divers diagrammes......................................................................214
A.2.5.1 Le diagramme des moments .......................................................................214
A.2.5.2 Le diagramme des L-moments ...................................................................214
A.3 Ajustement du modèle fréquentiel............................................................................215
A.3.1 Présentation de la loi de Gumbel ......................................................................215
A.3.2 Techniques d'ajustement .....................................................................................216
A.3.2.1 Méthode graphique.........................................................................................216
A.3.2.2 Méthode des moments ..................................................................................217
A.3.2.3 Méthode des L-moments ..............................................................................218
A.3.2.4 Méthode des moindres rectangles .............................................................218
A.3.2.5 Méthode du maximum de vraisemblance...............................................218
A.4 Contrôle de l'ajustement ..............................................................................................219
A.4.1 Examen visuel de l'ajustement ..........................................................................219
A.4.2 Le test chi-carré de K. Pearson ..........................................................................219
A.4.3 Le test de Kolmogorov-Smirnov .........................................................................221
A.4.4 Test d'Anderson–Darling ......................................................................................222
A.5 Analyse des incertitudes ..............................................................................................222
A.5.1 L'intervalle de confiance .......................................................................................223
A.5.1.1 Erreur-type d'un quantile ............................................................................223
ANNEXE 2 : REFERENCES BIBLIOGRAPHIQUES....................................... 227
Cours du Professeur André Musy et mis sous ce format par Lawani MOUNIROU VIII
Chapitre 1 : Cycle de l’eau et bilan hydrologique
La question de la disponibilité et d'accès à l'eau est sans aucun doute un des problèmes
majeurs auquel devra faire face l'humanité durant le siècle à venir. Aujourd'hui on estime en
effet qu'un habitant sur cinq de la planète n'a pas accès à l'eau en suffisance et un sur trois
a une eau de qualité. Dans ce contexte, il peut être utile de rappeler que "la mesure
quantitative et qualitative des éléments du cycle hydrologique et la mesure des autres
caractéristiques de l'environnement qui influent sur l'eau constituent une base essentielle
pour une gestion efficace de l'eau". (Déclaration de Dublin, 1992). De fait, la compréhension
et l'analyse du cycle de l'eau est la base de toute étude et réflexion au sujet de la gestion des
eaux.
L'eau est la source principale et originelle de toute vie. Elle se présente, dans la nature, sous
trois états :
♣ Solide : neige et glace.
♣ Liquide : eau chimiquement pure ou chargée en solutés.
♣ Gazeux : à différents degrés de pression et de saturation.
Le changement de phase de l'eau dépend essentiellement de la température et de la pression
mais aussi du degré de pollution de l'atmosphère. La figure suivante donne les différentes
conditions de pression et de température pour les trois états de l'eau, ainsi que les
transformations de phase.
Les mécanismes des mouvements de l'eau dans la nature sont déterminés par l'énergie
thermique solaire, la gravité, l'attraction solaire, l'attraction lunaire, la pression
atmosphérique, les forces intermoléculaires, les réactions chimiques, nucléaires et les
activités biologiques, et enfin les activités humaines. L'énergie thermique du soleil produit
une circulation de l'air dans l'atmosphère, en raison du fait que la surface terrestre est
réchauffée de façon inégale. La force de gravité est responsable des phénomènes de
précipitations, de ruissellement, d'infiltration et de courant de convection. L'attraction solaire
et lunaire est à l'origine des marées et des courants marins. Les différences de pression
atmosphérique occasionnent les déplacements horizontaux de l'air. Les vents sont eux-
mêmes responsables du mouvement des couches superficielles dans les lacs et les océans.
Les forces intermoléculaires dans le sol provoquent les phénomènes capillaires ainsi que la
viscosité et influencent donc la vitesse d'écoulement. L'eau est une des composantes de
plusieurs réactions chimiques organiques ou inorganiques. Un autre type de transformation
de l'eau est le processus physiologique qui se produit dans l'organisme animal. Finalement,
l'homme intervient directement sur les processus de mouvement et de transformation de
l'eau. Son action peut conduire à une meilleure gestion de sa plus précieuse ressource
naturelle, mais elle peut aussi causer de nombreux problèmes, notamment en perturbant le
cycle hydrologique, tant au niveau quantitatif que qualitatif.
Sous l'effet du rayonnement solaire, l'eau évaporée à partir du sol, des océans et des autres
surfaces d'eau, entre dans l'atmosphère. L'élévation d'une masse d'air humide permet le
refroidissement général nécessaire pour l'amener à saturation et provoquer la condensation
de la vapeur d'eau sous forme de gouttelettes constituant les nuages, en présence de noyaux
de condensation. Puis la vapeur d'eau, transportée et temporairement emmagasinée dans les
nuages, est restituée par le biais des précipitations aux océans et aux continents. Une partie
de la pluie qui tombe peut être interceptée par les végétaux puis être partiellement restituée
sous forme de vapeur à l'atmosphère. La pluie non interceptée atteint le sol. Suivant les
conditions données, elle peut alors s'évaporer directement du sol, s'écouler en surface
jusqu'aux cours d'eau (ruissellement de surface) ou encore s'infiltrer dans le sol. Il peut aussi
y avoir emmagasinement temporaire de l'eau infiltrée sous forme d'humidité dans le sol, que
peuvent utiliser les plantes. Il peut y avoir percolation vers les zones plus profondes pour
contribuer au renouvellement des réserves de la nappe souterraine. Un écoulement à partir
de cette dernière peut rejoindre la surface au niveau des sources ou des cours d'eau.
L'évaporation à partir du sol, des cours d'eau, et la transpiration des plantes complètent ainsi
le cycle.
Le cycle de l'eau est donc sujet à des processus complexes et variés parmi lesquels nous
citerons les précipitations, l'évaporation, la transpiration (des végétaux), l'interception, le
ruissellement, l'infiltration, la percolation, l'emmagasinement et les écoulements souterrains
qui constituent les principaux chapitres de l'hydrologie. Ces divers mécanismes sont rendus
possibles par un élément moteur, le soleil, organe vital du cycle hydrologique.
Sont dénommées précipitations toutes les eaux météoriques qui tombent sur la surface de la
terre, tant sous forme liquide (bruine, pluie, averse) que sous forme solide (neige, grésil, grêle)
et les précipitations déposées ou occultes (rosée, gelée blanche, givre,...). Elles sont
provoquées par un changement de température ou de pression. La vapeur d'eau de
l'atmosphère se transforme en liquide lorsqu'elle atteint le point de rosée par refroidissement
ou augmentation de pression. Pour produire la condensation, il faut également la présence de
certains noyaux microscopiques, autour desquels se forment des gouttes d'eau condensées.
La source de ces noyaux peut être océanique (chlorides, en particulier NaCl produit par
l'évaporation de la mer), continentale (poussière, fumée et autres particules entraînées par
des courants d'air ascendants) ou cosmiques (poussières météoriques). Le déclenchement des
précipitations est favorisé par la coalescence des gouttes d'eau. L'accroissement de poids leur
confère une force de gravité suffisante pour vaincre les courants ascendants et la turbulence
de l'air, et atteindre le sol. Enfin, le parcours des gouttes d'eau ou des flocons de neige doit
être assez court pour éviter l'évaporation totale de la masse. Les précipitations sont
exprimées en intensité (mm/h) ou en lame d'eau précipitée (mm) (rapport de la quantité d'eau
précipitée uniformément répartie sur une surface).
1.3.3 L'évaporation/l'évapotranspiration
♣ l'évapotranspiration réelle (ETR) : somme des quantités de vapeur d'eau évaporées par le
sol et par les plantes quand le sol est à une certaine humidité et les plantes à un stade
de développement physiologique et sanitaire spécifique.
L'évaporation est une des composantes fondamentales du cycle hydrologique et son étude est
essentielle pour connaître le potentiel hydrique d'une région ou d'un bassin versant. En
général, des analyses spécifiques d'évaporation devront être faites pour des études de bilan et
de gestion de l'eau par les plantes. Cependant, ces analyses approfondies sont moins
nécessaires pour les études de projets d'aménagement où l'eau est plutôt considérée sous un
aspect d'agent dynamique.
La pluie (ou dans certains cas la neige) peut être retenue par la végétation, puis redistribuée
en une partie qui parvient au sol et une autre qui s'évapore. La partie n'atteignant jamais le
sol forme l'interception. Son importance est difficile à évaluer et souvent marginale sous nos
climats, donc souvent négligée dans la pratique. Le stockage dans les dépressions est, tout
comme l'interception, souvent associé aux pertes. On définit l'eau de stockage comme l'eau
retenue dans les creux et les dépressions du sol pendant et après une averse.
L'infiltration désigne le mouvement de l'eau pénétrant dans les couches superficielles du sol
et l'écoulement de cette eau dans le sol et le sous-sol, sous l'action de la gravité et des effets
de pression. La percolation représente plutôt l'infiltration profonde dans le sol, en direction
de la nappe phréatique. Le taux d'infiltration est donné par la tranche ou le volume d'eau qui
s'infiltre par unité de temps (mm/h ou m3/s). La capacité d'infiltration ou l'infiltrabilité est la
tranche d'eau maximale qui peut s'infiltrer par unité de temps dans le sol et dans des
conditions données. L'infiltration est nécessaire pour renouveler le stock d'eau du sol,
alimenter les eaux souterraines et reconstituer les réserves aquifères. De plus, en absorbant
une partie des eaux de précipitation, l'infiltration peut réduire les débits de ruissellement.
De par la diversité de ses formes, on ne peut plus aujourd'hui parler d'un seul type
d'écoulement mais bien des écoulements. On peut distinguer en premier lieu les écoulements
rapides des écoulements souterrains plus lents. Les écoulements qui gagnent rapidement les
exutoires pour constituer les crues se subdivisent en écoulement de surface (mouvement de
l'eau sur la surface du sol) et écoulement de subsurface (mouvement de l'eau dans les
premiers horizons du sol). L'écoulement souterrain désigne le mouvement de l'eau dans le
sol. On peut encore ajouter à cette distinction les écoulements en canaux ou rivières qui font
appel à des notions plus hydrauliques qu'hydrologiques (à l'exception des méthodes de
mesures comme nous le verrons ultérieurement).
Au-delà de cette distinction simpliste – ces notions seront réexaminées plus en détail au
chapitre 9 consacré à l'étude des processus de génération des crues – on remarquera que les
écoulements peuvent aussi se signaler par leur domaine d'application. L'écoulement de
surface caractérise un écoulement sur une surface et s'exprime généralement par un rapport
volume / surface / temps [L3/L2/T]. Il est ainsi souvent exprimé en millimètre par année
hydrologique dans les études de bilans ou encore en litres par secondes et par hectares dans
le cadre de projet d'aménagement des terres et des eaux (drainage ou irrigation). Les
écoulements souterrains et en rivière font explicitement référence à la notion de débit, à
savoir à un volume d'eau traversant une section par unité de temps [L3/T].
Nous pouvons concevoir la répartition des eaux sur la terre selon différents points de vue :
♣ Une répartition spatiale du bilan de l'eau sur les continents et à l'échelle d'une zone
géographique.
La terre, vue de l'espace, apparaît comme une planète recouverte en grande partie d'eau
(planète bleue). Les océans occupent en effet une superficie à peu près égale à 70% de la
surface du globe et représentent 97% de la masse totale d'eau dans la biosphère. Le tableau
1.1 donne quelques grandeurs indicatives tandis que la figure 1.3 présente la disponibilité
mondiale d'eau.
On peut encore remarquer que la superficie des terres émergées de l'hémisphère Nord est
deux fois supérieure à celle de l'hémisphère sud. De plus la distribution spatiale des aires
continentales et océaniques à la surface du globe est inhomogène. La distribution
quantitative des eaux sur terre fait apparaître que les eaux dites douces ne représentent
qu'environ 3% du volume total des eaux du globe. Elles se retrouvent à 99% dans les calottes
polaires, les glaciers et les eaux souterraines de grandes profondeurs qui représentent des
réserves d'eau douce difficilement accessibles. Toutefois, dans certaines régions
montagneuses (Andes, Rocheuses, Alpes), les eaux de fonte alimentent la plupart des cours
d'eau et le débit des fleuves est étroitement lié au taux de fonte des glaciers.
Les eaux souterraines occupent le 2ème rang des réserves mondiales en eau douce après les
eaux contenues dans les glaciers. Elles devancent largement les eaux continentales de
surface. Leur apport est d'autant plus important que, dans certaines parties du globe, les
populations s'alimentent presque exclusivement en eau souterraine par l'intermédiaire de
puits, comme c'est le cas dans la majorité des zones semi-arides et arides. En Suisse, l'eau
potable a pour origine principale l'eau souterraine (70 - 80%) et secondaire l'eau de surface
(20 - 30%). On doit cependant garder à l'esprit que plus de la moitié de l'eau souterraine se
trouve à plus de 800 mètres de profondeur et que son captage demeure en conséquence
difficile. En outre, son exploitation abusive entraîne souvent un abaissement irréversible des
nappes phréatiques et parfois leur remplacement graduel par de l'eau salée (problème
rencontré en zone maritime telle qu'en Libye, Sénégal, Egypte, etc.).
Les eaux continentales de surface (lacs d'eau douce, rivières, fleuves, etc.) sont, à l'inverse
des eaux souterraines, très accessibles. Par contre, elles sont quantitativement infimes et
sont susceptibles d'être plus facilement polluées malgré l'effort fait depuis une dizaine
d'années pour en améliorer la qualité. Le Canada possède à lui seul 30 % des réserves
mondiales d'eau douce et 6% du ruissellement terrestre.
Quant aux eaux météoriques, elles peuvent paraître quantitativement très modestes, du
moins dans certaines régions. Néanmoins, elles constituent une étape essentielle du cycle de
l'eau. Le pourcentage d'eau disponible pour l'homme est certes très faible, mais suffisant
grâce à la circulation ou au recyclage de cette eau.
Dans chacun des ces grands réservoirs terrestres, l'eau se renouvelle au fil des ans. La
vitesse de renouvellement des eaux dans les réservoirs est mesurée par un flux : le temps de
séjour moyen ou temps de résidence est obtenu en divisant la taille du réservoir par le flux
d'entrée (somme de tous les flux entrants) ou de sortie (somme de tous les flux sortants)
(tableau 1.2).
Tableau 1.2 - Temps de renouvellement de l'eau dans les principaux réservoirs (Tiré de Gleick
(1993), Jacques (1996))
égales, par la distribution des altitudes. Cependant, cette représentation du cycle de l'eau
reste quand même approximative et les pourcentages attribués aux divers mécanismes de
transport de l'eau peuvent être quelque peu différents suivant les auteurs. Les trois
processus principaux, à savoir les précipitations, l'évaporation et le ruissellement,
décroissent de l'équateur vers les pôles.
Sur un même parallèle, l'intensité de l'évaporation sur les continents est pratiquement
uniforme. En général, la quantité totale de précipitations en un point est inversement
proportionnelle à sa distance à l'océan. Pour une même position géographique, les quantités
totales de précipitations et de ruissellement sont directement proportionnelles à l'élévation
moyenne du bassin versant jusqu'à une certaine altitude (optimum pluviométrique). Parmi
les composantes du cycle hydrologique, l'évaporation est la moins sensible aux changements
d'environnement géographique, suivie des précipitations et du ruissellement.
A l'échelle continentale, les principaux éléments de la répartition des eaux sont donnés par le
tableau 1.3 ci-après. Le pourcentage des précipitations qui ruisselle est plus important dans
l'hémisphère Nord (~40%) que dans l'hémisphère sud (Australie : ~35%, Afrique : ~20% et
Amérique du sud : ~10%).
Tableau 1.3 - Principaux éléments de la répartition des eaux à l'échelle du globe
Pour conclure sur le cycle hydrologique, on peut dire qu'il est caractérisé par
l'interdépendance de ses composantes, par sa stabilité et son équilibre dynamique. Si un
processus est perturbé, tous les autres (cycle de l'azote, cycle du phosphore, etc.) s'en
ressentent ! En particulier, le cycle hydrologique peut être influencé à des degrés divers par
les activités humaines. En effet, l'homme agit directement sur le processus de transformation
de l'eau, et cela de plusieurs façons : la construction de réservoirs, le transport de l'eau pour
des besoins industriels, le captage des eaux phréatiques, l'irrigation, le drainage, la
correction des cours d'eau, l'utilisation agricole des sols, l'urbanisation, les pluies
provoquées, etc., sont des exemples de l'intervention humaine.
♣ les précipitations,
♣ le ruissellement de surface et l'écoulement souterrain,
♣ l'évaporation.
Il est intéressant de noter que dans chacune des phases on retrouve respectivement un
transport d'eau, un emmagasinement temporaire et parfois un changement d'état. Il s'ensuit
que l'estimation des quantités d'eau passant par chacune des étapes du cycle hydrologique
peut se faire à l'aide d'une équation appelée "hydrologique" qui est le bilan des quantités
d'eau entrant et sortant d'un système défini dans l'espace et dans le temps. Le temporel
introduit la notion de l'année hydrologique. En principe, cette période d'une année est choisie
en fonction des conditions climatiques. Ainsi en fonction de la situation météorologique des
régions, l'année hydrologique peut débuter à des dates différentes de celle du calendrier
ordinaire. Au niveau de l'espace, il est d'usage de travailler à l'échelle d'un bassin versant
(entité structurelle définie en détails au chapitre 2) mais il est possible de raisonner à un
autre niveau (zone administrative, entité régionale, etc.).
L'équation du bilan hydrique se fonde sur l'équation de continuité et peut s'exprimer comme
suit, pour une période et un bassin donnés :
(1.1)
Avec :
P : précipitations (liquide et solide) [mm],
S : ressources (accumulation) de la période précédente (eaux souterraines,
humidité du sol, neige, glace) [mm],
R : ruissellement de surface et écoulements souterrains [mm],
E : évaporation (y compris évapotranspiration) [mm],
S + DS : ressources accumulées à la fin de la période [mm].
On exprime généralement les termes du bilan hydrique en hauteur d'eau (mm par exemple),
on parle alors de lame d'eau (précipitée, écoulée, évaporée, stockée, etc.). Cette équation
exprime simplement que la différence entre le débit d'eau entrant et le débit d'eau sortant
d'un volume donné (par exemple un bassin versant) au cours d'une période déterminée est
égale à la variation du volume d'eau emmagasinée au cours de la dite période. Elle peut
s'écrire encore sous la forme simplifiée suivante :
(1.2)
Avec :
E : évaporation [mm] ou [m3],
I : volume entrant [mm] ou [m3],
O : volume sortant [mm] ou [m3],
DS : variation de stockage [mm] ou [m3].
Si le bassin versant naturel est relativement imperméable, la variation de stock sur une
période donnée peut être considérée comme nulle (DS=0). Dès lors, on peut introduire le
déficit d'écoulement D dans l'équation qui s'écrit :
(1.3)
Ce déficit d'écoulement représente essentiellement les pertes dues à l'évaporation. Il peut être
estimé à l'aide de mesures ou de méthodes de calcul. A titre illustratif, les formules de Turc
et Coutagne sont les suivantes :
1. Formule de Turc
(1.4)
Avec :
D : déficit d'écoulement [mm],
P : pluie annuelle
T : température moyenne annuelle [°C].
L = 300 + 25 T + 0.05 T3.
2. Formule de Coutagne
(1.5)
Avec :
D : déficit d'écoulement [mm],
P : pluie annuelle [mm],
m= 1/(0.8 + 0.16 T) : coefficient régional (m=0.42 pour la
France).
L'application de la méthode du bilan hydrique est limitée par la difficulté de quantifier les
variables. Effectivement, les processus hydrologiques sont difficiles à observer directement
sur le terrain et donc à mesurer. Notons que les erreurs de mesure éventuelles des termes
qu'on retrouve dans l'équation hydrologique simplifiée se répercutent directement sur les
valeurs calculées de l'évaporation. Devant ces imprécisions, on suggère l'emploi de cette
méthode dans le cas d'un avant-projet par exemple, pour vérifier l'état du système et surtout
la validité (la fiabilité) des mesures qui le décrit.
Le bassin versant représente, en principe, l'unité géographique sur laquelle se base l'analyse
du cycle hydrologique et de ses effets.
Plus précisément, le bassin versant qui peut être considéré comme un " système " est une
surface élémentaire 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.
Le bassin versant en une section droite d'un cours d'eau, est donc défini comme la totalité de
la surface topographique drainée par ce cours d'eau et ses affluents à l'amont de cette
section. Il est entièrement caractérisé par son exutoire, à partir duquel nous pouvons tracer
le point de départ et d'arrivée de la ligne de partage des eaux qui le délimite.
Généralement, la ligne de partage des eaux correspond à la ligne de crête. On parle alors de
bassin versant topographique.
Fig. 2.2 - Distinction entre bassin versant réel et bassin versant topographique
D'après Roche - Hydrologie de surface, Ed. Gauthier-Villars, Paris 1963.
Cette différence entre bassins réel et topographique est tout particulièrement importante en
région karstique. Lorsque l'on s'intéresse au ruissellement, la délimitation du bassin versant
doit aussi tenir compte des barrières artificielles (routes, chemins de fer, etc.). En effet,
l'hydrologie du bassin versant, et notamment la surface drainée, peuvent être modifiées par
la présence d'apports latéraux artificiels (réseaux d'eaux usées ou potables, drainages,
routes, pompages ou dérivations artificielles modifiant le bilan hydrologique).
Fig. 2.3 - Exemples de modifications de la délimitation du bassin versant suite à la mise en place
d'un réservoir et la construction d'une route
Il convient donc également de définir, en plus des délimitations topographiques, les limites
souterraines de ce système. De plus, il est aussi nécessaire de tenir compte des effets
anthropiques relatifs aux eaux du système.
La réaction hydrologique d'un bassin versant à une sollicitation particulière (Fig. 2.4) est
caractérisée par sa vitesse (temps de montée tm, défini comme le temps qui s'écoule entre
l'arrivée de la crue et le maximum de l'hydrogramme) et son intensité (débit de pointe Qmax,
volume maximum Vmax, etc.). Ces deux caractéristiques sont fonction du type et de l'intensité
de la précipitation qui le sollicite mais aussi d'une variable caractérisant l'état du bassin
versant : le temps de concentration des eaux sur le bassin.
Fig. 2.4 - Principes d'analyse du comportement hydrologique du bassin versant et hydrogramme résultant.
La figure 2.5 fourni un exemple d'hydrogramme de crue résultant d'un hyétogramme donné.
Le hyétogramme est la courbe représentant l'intensité de la pluie en fonction du temps.
Fig. 2.5 - Exemple de réaction hydrologique pour le bassin versant de Bois-Vuacoz (Haute-Mentue)
Le temps de concentration tc des eaux sur un bassin versant se définit comme le maximum
de durée nécessaire à une goutte d'eau pour parcourir le chemin hydrologique entre un point
du bassin et l'exutoire de ce dernier. Il est composé de trois termes différents :
♣ th : Temps d'humectation. Temps nécessaire à l'imbibition du sol par l'eau qui tombe
avant qu'elle ne ruisselle.
♣ tr : Temps de ruissellement ou d'écoulement. Temps qui correspond à la durée
d'écoulement de l'eau à la surface ou dans les premiers horizons de sol jusqu'à un
système de collecte (cours d'eau naturel, collecteur).
♣ ta : Temps d'acheminement. Temps mis par l'eau pour se déplacer dans le système de
collecte jusqu'à l'exutoire.
Le temps de concentration tc est donc égal au maximum de la somme de ces trois termes,
soit :
(2.1)
Théoriquement on estime que tc est la durée comprise entre la fin de la pluie nette et la fin du
ruissellement (cf. chapitre 11). Pratiquement le temps de concentration peut être déduit de
mesures sur le terrain ou s'estimer à l'aide de formules le plus souvent empiriques.
Les courbes isochrones représentent les courbes d'égal temps de concentration des eaux sur
le bassin versant. Ainsi, l'isochrone la plus éloignée de l'exutoire représente le temps mis
pour que toute la surface du bassin versant contribue à l'écoulement à l'exutoire après une
averse uniforme (Fig. 2.6). Le tracé du réseau des isochrones permet donc de comprendre en
Fig. 2.6 - Représentation d'un bassin avec ses lignes isochrones et diagramme surface-temps de
concentration du bassin par élément de surface. On remarquera la forme des isochrones au
voisinage des éléments constitutifs du réseau hydrographique.
[Link] La surface
Le bassin versant étant l'aire de réception des précipitations et d'alimentation des cours
d'eau, les débits vont être en partie reliés à sa surface.
La surface du bassin versant peut être mesurée par superposition d'une grille dessinée sur
papier transparent, par l'utilisation d'un planimètre ou, mieux, par des techniques de
digitalisation.
[Link] La forme
(2.2)
Avec :
KG est l'indice de compacité de Gravélius,
A : surface du bassin versant [km2],
P : périmètre du bassin [km].
circulaire et supérieur à 1 lorsque le bassin est de forme allongée, tel qu'illustré par la figure
2.8.
1. La courbe hypsométrique
Ajoutons que lorsqu'on désire caractériser des bassins versants de haute montagne,
on a l'habitude de tracer des courbes hypsométriques glaciaires, en planimétrant les
surfaces recouvertes de glace.
Les courbes hypsométriques demeurent un outil pratique pour comparer plusieurs
bassins entre eux ou les diverses sections d'un seul bassin. Elles peuvent en outre
servir à la détermination de la pluie moyenne sur un bassin versant et donnent des
indications quant au comportement hydrologique et hydraulique du bassin et de son
système de drainage.
b. L'altitude moyenne
(2.3)
Avec :
Hmoy : altitude moyenne du bassin [m] ;
Ai : aire comprise entre deux courbes de niveau [km2] ;
hi : altitude moyenne entre deux courbes de niveau [m] ;
A : superficie totale du bassin versant [km2].
L'altitude moyenne est peu représentative de la réalité. Toutefois, elle est
parfois utilisée dans l'évaluation de certains paramètres hydrométéorologiques
ou dans la mise en œuvre de modèles hydrologiques.
c. L'altitude médiane
La pente moyenne est une caractéristique importante qui renseigne sur la topographie
du bassin. Elle est considérée comme une variable indépendante. Elle donne une
bonne indication sur le temps de parcours du ruissellement direct - donc sur le temps
de concentration tc - et influence directement le débit de pointe lors d'une averse.
Plusieurs méthodes ont été développées pour estimer la pente moyenne d'un bassin.
Toutes se basent sur une lecture d'une carte topographique réelle ou approximative.
La méthode proposée par Carlier et Leclerc (1964) consiste à calculer la moyenne
pondérée des pentes de toutes les surfaces élémentaires comprises entre deux
altitudes données. Une valeur approchée de la pente moyenne est alors donnée par la
relation suivante :
(2.4)
Où :
im : pente moyenne[m/km ou 0/00],
L : longueur totale de courbes de niveau [km],
D : équidistance entre deux courbes de niveau [m],
A : surface du bassin versant [km2].
Cette méthode de calcul donne de bons résultats dans le cas d'un relief modéré et
pour des courbes de niveau simples et uniformément espacées. Dans les autres cas, il
convient de styliser les courbes de niveau pour que leur longueur totale ait un sens
réel vis-à-vis de la pente.
Le calcul de la pente moyenne tout comme celui de leur exposition (orientation des
pentes) peut-être assez facilement automatisée en se basant sur des données
numériques représentant la topographie des bassins versants (Modèle Numérique
d'Altitude). Le recours à ces données et méthodes et vivement encouragé. La dernière
section de ce chapitre est consacrée aux informations digitales et aux modèles
numériques.
4. L'indice de pente ip
Cet indice se calcule à partir du rectangle équivalent. Il est égal à la somme des
racines carrées des pentes moyennes de chacun des éléments pondérés par la surface
intéressée, soit :
(2.5)
où :
ip: indice de pente [%],
L: longueur du rectangle [m],
xi: distance qui sépare deux courbes sur la rectangle [m] (la largeur du rectangle étant
constante, cette distance est égale au facteur de pondération),
d: distance entre 2 courbes de niveau successives (peut être variable) [m],
d/xi: pente moyenne d'un élément [%].
coefficient de compacité : .
En combinant ces trois relations, on obtient :
(2.6)
Le réseau hydrographique se définit comme l'ensemble des cours d'eau naturels ou artificiels,
permanents ou temporaires, qui participent à l'écoulement. Le réseau hydrographique est
sans doute une des caractéristiques les plus importantes du bassin. Le réseau
hydrographique peut prendre une multitude de formes. La différenciation du réseau
hydrographique d'un bassin est due à quatre facteurs principaux.
♣ Le climat : le réseau hydrographique est dense dans les régions montagneuses très
humides et tend à disparaître dans les régions désertiques.
Afin de caractériser le réseau hydrographique, il est souvent utile de reporter son tracé en
plan sur une carte à une échelle adéquate. L'utilisation de photographies analogiques ou
numériques est utile à cette identification. Divers paramètres descriptifs sont utilisés pour
définir le réseau hydrographique.
La classification est facilitée par un système de numérotation des tronçons de cours d'eau
(rivière principale et affluents). L'ordre des cours d'eau est donc une classification qui reflète
la ramification du cours d'eau. La codification des cours d'eau est également utilisée pour la
codification des stations de mesures, permettant ainsi un traitement automatisé des
données. Il existe plusieurs types de classifications des tronçons des cours d'eau, dont la
classification de Strahler (1957) qui est la plus utilisée.
Un bassin versant se caractérise principalement par les deux longueurs suivantes, illustrées
sur la figure ci-dessous.
La longueur d'un bassin versant (LCA) est la distance curviligne mesurée le long du cours
d'eau principal depuis l'exutoire jusqu'à un point représentant la projection du centre de
gravité du bassin sur un plan (Snyder, 1938).
• La longueur du cours d'eau principal (L) est la distance curviligne depuis l'exutoire
jusqu'à la ligne de partage des eaux, en suivant toujours le segment d'ordre le plus
élevé lorsqu'il y a un embranchement et par extension du dernier jusqu'à la limite
topographique du bassin versant. Si les deux segments à l'embranchement sont de
même ordre, on suit celui qui draine la plus grande surface.
Fig. 2.11 - Longueurs caractéristiques d'un bassin versant, LCA : longueur du bassin versant ; L :
longueur du cours d'eau principal
Fig. 2.12 - Profil en long de la Broye avec représentation de ses affluents (D'après Parriaux :
Contribution à l'étude des ressources en eau du bassin de la Broye)
La pente moyenne du cours d'eau détermine la vitesse avec laquelle l'eau se rend à l'exutoire
du bassin donc le temps de concentration. Cette variable influence donc le débit maximal
observé. Une pente abrupte favorise et accélère l'écoulement superficiel, tandis qu'une pente
douce ou nulle donne à l'eau le temps de s'infiltrer, entièrement ou en partie, dans le sol.
Le calcul des pentes moyennes et partielles de cours d'eau s'effectue à partir du profil
longitudinal du cours d'eau principal et de ses affluents. La méthode la plus fréquemment
utilisée pour calculer la pente longitudinale du cours d'eau consiste à diviser la différence
d'altitude entre les points extrêmes du profil par la longueur totale du cours d'eau.
(2.7)
Où :
Pmoy : pente moyenne du cours d'eau [m/km] ;
DHmax : dénivellation maximale de la rivière [m] (différence d'altitude entre le point le plus
éloigné et l'émissaire) ;
L : longueur du cours d'eau principal [km].
On préférera parfois utiliser d'autres méthodes plus représentatives : par exemple celle qui
consiste à assimiler la pente moyenne à la pente de la droite tracée entre les points situés à
15% et 90% de distance à partir de l'exutoire, suivant le cours d'eau principal (Benson,
1959) ; ou encore, comme le préconise Linsley (1982), on prendra la pente de la ligne, tracée
depuis l'exutoire, dont la surface délimitée est identique à la surface sous le profil en long
(Fig. 2.13).
Fig. 2.13 - Calcul de la pente moyenne du cours d'eau selon Linsley (1982)
♣ Courbe aire-distance
A partir de données sur un bon nombre de bassins versants (Hack, 1957), une relation a pu
être établie entre la longueur L [km] de la rivière et l'aire A [km2] du bassin versant :
(2.8)
On peut aussi définir la courbe aire-distance, qui met en relation la longueur moyenne des
cours d'eau d'ordre u donné et l'aire tributaire moyenne des cours d'eau du même ordre u, et
ceci ordre par ordre. Cette courbe permet de visualiser la répartition des superficies du
bassin par rapport à l'exutoire ou par rapport au point de mesure du débit. Cette répartition affecte
en effet la concentration du ruissellement et donc influence la réponse hydrologique du bassin versant.
(2.9)
Avec :
♣ La densité hydrographique
(2.10)
Où :
F : densité hydrographique [km-2] ;
Ni : nombre de cours d'eau ;
A : superficie du bassin [km2].
Il existe une relation assez stable entre la densité de drainage Dd et la densité hydrographique
F, de la forme :
(2.11)
♣ Le rapport de confluence
Sur la base de la classification des cours d'eau, Horton (1932) et Schumm (1956) ont établi
différentes lois :
(2.12)
Loi des nombres :
(2.13)
Loi des longueurs :
(2.14)
Loi des aires :
Avec :
RB : rapport de confluence des cours d'eau ("bifurcation ratio") ;
RL : rapport des longueurs des cours d'eau ; RA : rapport des aires des cours d'eau ;
u : ordre d'un cours d'eau u varie entre 1 et w (w est l'ordre du cours d'eau principal,
classification selon Strahler) ;
Nu : nombre des cours d'eau d'ordre u ; Nu+1 : nombre des cours d'eau d'ordre suivant ;
Lu : longueur moyenne des cours d'eau d'ordre u ;
Au : aire tributaire moyenne des cours d'eau d'ordre u.
[Link] L'endoréisme
L'endoréisme est un phénomène rencontré dans certains bassins versants pour lesquels le
réseau hydrographique n'est relié à aucun autre réseau. L'eau est alors acheminée et
concentrée en un point du bassin qui peut être un lac, une mare ou une accumulation
souterraine. Ce phénomène est généralement observé en zones arides (ex : mare d'Oursi au
Burkina Faso, lac Tchad, mer Morte, etc.).
♣ La couverture végétale
L'activité végétative et le type de sol sont intimement liés et leurs actions combinées
influencent singulièrement l'écoulement en surface. Le couvert végétal retient, selon sa
densité, sa nature et l'importance de la précipitation, une proportion variable de l'eau
atmosphérique. Cette eau d'interception est en partie soustraite à l'écoulement.
La forêt, par exemple, intercepte une partie de l'averse par sa frondaison. Elle exerce une
action limitatrice importante sur le ruissellement superficiel. La forêt régularise le débit des
cours d'eau et amortit les crues de faibles et moyennes amplitudes. Par contre, son action
sur les débits extrêmes causés par des crues catastrophiques est réduite.
A l'inverse, le sol nu, de faible capacité de rétention favorise un ruissellement très rapide.
L'érosion de la terre va généralement de paire avec l'absence de couverture végétale.
Etant donné l'importance du rôle joué par la forêt, on traduit parfois sa présence par un
indice de couverture forestière K :
(2.15)
On peut calculer ce type d'indice avec d'autres couvertures végétales telle que les cultures.
On soulignera encore que la surface du cours d'eau constitue aussi un plan d'eau et que le
canal d'une rivière permet aussi de laminer une crue.
Fig. 2.15 - Illustration de l'effet de laminage par un plan d'eau. Le cas du Léman sur le débit du Rhône
Il est toujours possible de calculer un indice analogue à celui de la couverture forestière pour
les surfaces enneigées et celles des glaciers.
Les surfaces imperméables jouent un très grand rôle en hydrologie urbaine. Elles augmentent
l'écoulement de surface, réduisent les infiltrations et la recharge des nappes, et diminuent le
temps de concentration. On calcule souvent un taux d'imperméabilité qui est le rapport entre
les surfaces imperméables et la surface totale.
♣ Le coefficient de ruissellement
Pour caractériser la capacité d'un bassin versant à ruisseler un indice est très souvent utilisé
en hydrologie de surface : le coefficient de ruissellement (Cr). Son calcul et son emploi sont
simples, mais notons qu'il peut conduire à commettre de grossières erreurs. Ce coefficient est
défini comme suit :
(2.16)
Ce coefficient est fortement influencé par la couverture du sol comme le montre le tableau
suivant dans lequel les quelques valeurs de ce coefficient issues des normes suisses SNV
sont présentées. Ces valeurs reflètent la capacité des sols à ruisseler en fonction uniquement
de la couverture du sol. On remarque notamment le très fort taux du coefficient de
ruissellement donné pour les routes et toitures. Comme on l'a vu, cela s'explique par le fait
que ces surfaces sont pratiquement imperméables.
Bois 0,1
Rochers 0,7
La nature du sol intervient sur la rapidité de montée des crues et sur leur volume. En effet, le
taux d'infiltration, le taux d'humidité, la capacité de rétention, les pertes initiales, le
coefficient de ruissellement (Cr) sont fonction du type de sol et de son épaisseur.
Pour étudier ce type de réactions, on peut comparer le coefficient de ruissellement sur
différentes natures de sol (intérêt d'une carte pédologique détaillée dans les études de
prédétermination des crues). La littérature fournit des valeurs du coefficient de ruissellement
pour chaque type de sol et, très souvent, en rapport avec d'autres facteurs tels que la
couverture végétale, la pente du terrain ou l'utilisation du sol. Un exemple est donné dans le
tableau 2.2 pour la Suisse, et en secteur rural.
Tableau 2.2 Différentes valeurs de coefficient de ruissellement pour les cas suisses. Cr est
une fonction de la pente et de la couverture du sol. (Tiré de Sautier, Guide du Service Fédéral
des Améliorations foncières)
Couverture du sol
On peut introduire, dès à présent, une caractéristique du sol importante : l'état d'humidité du
sol qui est un des facteurs principaux conditionnant les temps de concentration. Cet état est
cependant très difficile à mesurer car très variable dans l'espace et le temps. On a souvent
recours à d'autres paramètres qui reflètent l'humidité du sol et qui sont plus faciles à
obtenir. En hydrologie, on fait souvent appel à des indices caractérisant les conditions
d'humidité antécédentes à une pluie. Il en existe de nombreux qui sont pour la plupart basés
sur les précipitations tombées au cours d'une certaine période précédant un événement. Ils
sont généralement notés IPA, c'est-à-dire Indices de Précipitations Antécédentes (API en
anglais).
La forme la plus classique de cet indice repose sur le principe de décroissance logarithmique
avec le temps du taux d'humidité du sol, au cours des périodes sans précipitations :
(2.17)
Avec :
IPA0 : valeur initiale de l'indice des précipitations antécédentes [mm] ;
IPAt : valeur de cet indice t jours plus tard [mm] ;
K : facteur de récession, K< 1. Il est variable d'un bassin à l'autre, ainsi que d'une saison à
l'autre pour un même bassin ;
t : temps [jour].
L'Institut d'Aménagement des Terres et des Eaux de l'EPFL (IATE/HYDRAM), après différents
travaux de recherche sur parcelles expérimentales, a adopté un indice de la forme suivante :
(2.18)
Où :
IPAi : indice de précipitations antérieures au jour i [mm] ;
IPAi-1 : indice de pluies antécédentes au jour i-1 [mm] ;
Pi-1 : précipitations tombées au jour i-1 [mm] ;
K : coefficient inférieur à 1, en général compris entre 0,8 et 0,9.
La figure 2.16 illustre le calcul de l'IPA au cours d'une année à la station de Payerne (VD).
Fig. 2.16 - Variation de l'indice IPA en fonction du temps à Payerne (VD) en 1991 (K = 0.9).
La connaissance de la géologie d'un bassin versant s'avère importante pour cerner l'influence
des caractéristiques physiographiques. La géologie du substratum influe non seulement sur
l'écoulement de l'eau souterraine mais également sur le ruissellement de surface. Dans ce
dernier cas, les caractères géologiques principaux à considérer sont la lithologie (nature de la
roche mère) et la structure tectonique du substratum. L'étude géologique d'un bassin versant
dans le cadre d'un projet hydrologique a surtout pour objet de déterminer la perméabilité du
substratum. Celle-ci intervient sur la vitesse de montée des crues, sur leur volume et sur le
soutien apporté aux débits d'étiage par les nappes souterraines. Un bassin à substratum
imperméable présente une crue plus rapide et plus violente qu'un bassin à substratum
perméable, soumis à une même averse. Ce dernier retient l'eau plus aisément, et en période
de sécheresse, un débit de base sera ainsi assuré plus longtemps. Néanmoins, le substratum
peut absorber une certaine quantité d'eau dans les fissures et diaclases des roches
naturellement imperméables ou dans les formations rocheuses altérées.
La demande de données spatiales s'est accrue ces dernières années car l'on sait désormais
qu'il est essentiel de connaître la distribution spatiale de la réponse hydrologique pour bien
comprendre les processus sous-jacents de la génération de l'écoulement. De plus, la
représentation et la connaissance du terrain sont essentielles pour comprendre les processus
d'érosion, de sédimentation, de salinisation et de pollution via des cartes de risque.
Nous n'aborderons ici que les modèles numériques d'altitude (MNA) et de terrain (MNT).
A partir de ces trois approches, il est possible de déterminer plusieurs attributs du modèle
numérique d'altitude tels que des attributs topographiques (élévation, orientation, pente,
surface, courbure) qui influencent diverses grandeurs intervenant directement dans les
processus d'écoulement.
2.4.2 En Suisse
En Suisse, le nouveau modèle numérique du terrain MNT25 est disponible pour toute la
superficie du pays depuis fin 1996 ([Link] Ce
modèle est établi à partir de la digitalisation des courbes de niveaux des feuilles
topographiques à l'échelle 1:25'000. Dans une seconde étape, le modèle matriciel du MNT25
est interpolé avec une maille de 25 m. Ce jeu de données est uniquement destiné à l'emploi
numérique. Il répond aux exigences demandées pour des applications d'une très grande
précision. La précision altimétrique du MNT25 est d'environ 1,5 m sur le Plateau, entre 5 et 8
m dans les Alpes.
Sont dénommées précipitations, toutes les eaux météoriques qui tombent sur la surface de la
terre, tant sous forme liquide (bruine, pluie, averse) que sous forme solide (neige, grésil, grêle)
et les précipitations déposées ou occultes (rosée, gelée blanche, givre,...). Elles sont
provoquées par un changement de température ou de pression. Les précipitations
constituent l’unique « entrée » des principaux systèmes hydrologiques continentaux que sont
les bassins versants.
Les processus responsables de la formation des nuages sont décrits dans les manuels de
climatologie et leur exposé détaillé sort du cadre de ce cours. Signalons toutefois que la
forme, l'ampleur, le développement des nuages dépendent de l'importance et de l'étendue
horizontale des mouvements verticaux ascendants qui leur donnent naissance. Quant aux
types de nuages, on distingue deux morphologies de base : les nuages stratiformes et
cumuliformes. On classe généralement les nuages aussi en fonction de leur altitude : nuages
supérieurs, nuages moyens, nuages inférieurs et nuages à développement vertical.
Fig. 3.1 - Exemple de nuages supérieurs (à gauche) et de nuages moyens (à droite) : le cirrus et l'altocumulus
([Link]
La saturation n'est cependant pas une condition suffisante à la condensation ; cette dernière
requiert également la présence de noyaux de condensation (impuretés en suspension dans
l'atmosphère d'origines variées - suie volcanique, cristaux de sable, cristaux de sel marin,
combustions industrielles, pollution) autour desquels les gouttes ou les cristaux se forment.
Lorsque les deux conditions sont réunies, la condensation intervient sur les noyaux ; il y a
alors apparition de gouttelettes microscopiques qui grossissent à mesure que se poursuit
l'ascendance, celle-ci étant le plus souvent la cause génératrice de la saturation. Les noyaux
de condensation jouent en faite un rôle de catalyseur pour la formation de gouttelettes d’eau.
Pour qu’il y ait précipitations il faut encore que les gouttelettes ou les cristaux composant les
nuages (les hydrométéores) se transforment en gouttes de pluie. Ce phénomène est lié à
l'accroissement de ces éléments dont la masse devient suffisante pour vaincre les forces
d'agitation. Ce grossissement peut s'expliquer par les deux processus suivant :
à température positive suffisamment épaisse (souvent à partir de 300 m dans les nuages
stables) et si la durée de chute le permet, il fond et donne lieu à de la pluie. Le même
processus de grossissement a lieu entre deux gouttelettes à des températures différentes
(la plus froide grossit au détriment de la plus chaude).
♣ Les précipitations convectives. Elles résultent d'une ascension rapide des masses d'air
dans l'atmosphère. Elles sont associées aux cumulus et cumulo-nimbus, à
développement vertical important, et sont donc générées par le processus de Bergeron.
Les précipitations résultantes de ce processus sont en général orageuses, de courte
durée (moins d'une heure), de forte intensité et de faible extension spatiale.
♣ Les précipitations orographiques. Comme son nom l'indique (du grec oros, montagne),
ce type de précipitations résulte de la rencontre entre une masse d’air chaude et humide
et une barrière topographique particulière. Par conséquent, ce type de précipitations
n’est pas « spatialement mobile » et se produit souvent au niveau des massifs
montagneux. Les caractéristiques des précipitations orographiques dépendent de
l'altitude, de la pente et de son orientation, mais aussi de la distance séparant l'origine
♣ Les précipitations frontales ou de type cyclonique. Elles sont associées aux surfaces de
contact entre deux masses d'air de température, de gradient thermique vertical,
d'humidité et de vitesse de déplacement différents, que l'on nomme « fronts ». Les fronts
froids (une masse d’air froide pénètre dans une région chaude) créent des précipitations
brèves, peu étendues et intenses. Du fait d’une faible pente du front, les fronts chauds
(une masse d’air chaude pénètre dans une région occupée par une masse d’air plus
froide) génèrent des précipitations longues, étendues, mais peu intenses.
Nom Caractéristiques
Finalement, les précipitations sont un des processus hydrologiques les plus variables. D'une
part, elles sont caractérisées par une grande variabilité dans l'espace et ceci quelle que soit
l'échelle spatiale prise en compte (régionale, locale, etc.). D'autre part, elles sont caractérisées
par une grande variabilité dans le temps, aussi bien à l'échelle annuelle qu'à celle d'un
événement pluvieux
Les différents instruments permettant la mesure des précipitations sont décrits dans le
chapitre 7 "mesures hydrologiques". Citons toutefois les deux appareils de mesures
fondamentaux que sont :
Pour un bassin versant donné ou une région donnée, les stations pluviométriques forment un
réseau d'observations. Elles fournissent des mesures ponctuelles.
Les données relatives aux stations sont d'une haute importance pour les statistiques
climatiques, la planification et la gestion des ressources et les projets de construction ; la
nature et la densité des réseaux doivent donc tenir compte du phénomène observé, du but
des observations, de la précision désirée, de la topographie, de facteurs économiques ou
d'autres encore.
La représentativité des précipitations par les mesures est fonction du réseau d'observation.
Plus celui-ci est dense, meilleure est l'information et plus l'ensemble des mesures est
représentatif de la lame d'eau tombée sur une surface donnée. Cependant le réseau est le
résultat d'un compromis entre la précision désirée et les possibilités ou charges
d'exploitation. Le réseau devra donc être planifié. Il existe plusieurs théories sur la
planification optimale d'un réseau, mais elles donnent des résultats approximatifs, qui
doivent toujours être adaptées aux contraintes locales et financières.
L'hydrologue devra donc faire appel à son expérience de terrain pour planifier un réseau. Il
tiendra compte du relief et du type de précipitations (frontales, orographiques, convectives). Il
s'assurera également des facilités d'accès, de contrôle et de transmission des informations
(par l'homme ou par télétransmission : téléphone, satellite, etc.).
La publication des données pluviométriques est du ressort des services publics (en Suisse et
à l'échelle nationale, l'Institut Suisse de Météorologie) qui le font généralement sous forme
d'annuaires. En Suisse, la publication de référence s'intitule "Ergebnisse der täglichen
Niederschlagen" (résultats des mesures de précipitations journalières). Les annuaires
pluviométriques regroupent, pour chacune des stations de mesure, les résultats suivants :
Un certain nombre de ces grandeurs est accessible en temps réel par le biais du site Internet
de l'Institut Suisse de Météorologie ([Link] Il est alors possible de
consulter l'évolution et la répartition spatiale de plusieurs paramètres hydroclimatiques.
Certaines de ces valeurs peuvent être régionalisées et présentées sous forme de cartes
d'isohyètes (cartes d'équivaleurs de précipitations). Il existe d'autres ouvrages de synthèse qui
sont davantage dirigés vers une analyse synthétique des précipitations (exemple de l'«Atlas
hydrologique de la Suisse» [Link]
Les mesures ponctuelles acquises au niveau des pluviomètres ou des pluviographes sont
analysées et soumises à différents traitements statistiques.
(3.1)
Où :
im : intensité moyenne de la pluie [mm/h, mm/min] ou ramenée à la surface [l/[Link]],
h : hauteur de pluie de l'averse [mm],
t : durée de l'averse [h ou min].
L'intensité des précipitations varie à chaque instant au cours d'une même averse suivant les
caractéristiques météorologiques de celle-ci. Plutôt que de considérer l'averse entière et son
intensité moyenne, on peut s'intéresser aux intensités observées sur des intervalles de temps
au cours desquels on aura enregistré la plus grande hauteur de pluie. On parle alors
d'intensité maximale.
La courbe des hauteurs de pluie cumulées représente en ordonnée, pour chaque instant t,
l'intégrale de la hauteur de pluie tombée depuis le début de l'averse.
Le critère de continuité d'un épisode pluvieux varie selon le bassin versant. Généralement,
deux averses sont considérées comme distinctes : (1) si la précipitation DH tombant durant
l'intervalle de temps Dt qui les sépare est inférieure à un certain seuil et (2) si cet intervalle
de temps est lui-même supérieur à une certaine valeur définie compte tenu du type de
problème étudié. En représentant les averses sous forme de hyétogrammes, la problématique
de la séparation des averses se résume comme suit (figure 3.4) :
Fig. 3.4 - Conditions pour la distinction de deux averses consécutives (1) DH durant Dt < seuil (par
exemple 2 mm) et (2) Dt > durée choisie en fonction du problème (par exemple 1 heure)
Cette notion d'averse est très importante en milieu urbain et de petits bassins versants car
elle s'avère déterminante pour l'estimation des débits de crue.
L'ensemble des données d'une station de mesures pluviométriques constitue une information
considérable qu'il est souhaitable de condenser à l'aide de caractéristiques bien choisies. On
applique ainsi les lois et d'autres techniques de la statistique aux relevés pluviométriques
pour en tirer des informations utiles aux études et travaux envisagés. On détermine de la
sorte :
L'ensemble de ces valeurs ponctuelles, condensées sous forme statistique, est utilisé pour
déterminer la fréquence et les caractéristiques d'un événement pluvieux isolé ou encore pour
étudier la variabilité de la pluviométrie dans l'espace.
Les projets d'aménagements hydrauliques ou hydrologiques sont souvent définis par rapport
à une averse type associée aux fréquences probables d'apparition.
Lorsque l'on étudie des grandeurs comme les précipitations (caractérisées à la fois par leur
hauteur et leur durée) ou les débits de crue d'un point de vue statistique, on cherche donc et,
en règle générale, à déterminer par exemple la probabilité pour qu'une intensité i ne soit pas
atteinte ou dépassée (i.e. soit inférieure ou égale à une valeur xi).
Cette probabilité est donnée, si i représente une variable aléatoire, par la relation suivante :
(3.2)
On définit alors le temps de retour T d'un événement comme étant l'inverse de la fréquence
d'apparition de l'événement. Soit :
(3.3)
Ainsi, l'intensité d'une pluie de temps de retour T est l'intensité qui sera dépassé en moyenne
toutes les T années.
Une pluie peut être caractérisée par plusieurs paramètres qui peuvent avoir, au sein de la
même pluie, des temps de retour très différents. Citons notamment :
L'analyse des pluies a permis de définir deux lois générales de pluviosité qui peuvent
s'exprimer de la manière suivante :
♣ Pour une même fréquence d'apparition - donc un même temps de retour - l'intensité
d'une pluie est d'autant plus forte que sa durée est courte.
♣ Ou encore, en corollaire, à durée de pluie égale, une précipitation sera d'autant plus
intense que sa fréquence d'apparition sera petite (donc que son temps de retour sera
grand).
Ces lois permettant d'établir les relations entre les intensités, la durée et la fréquence
d'apparition des pluies peuvent être représentées selon des courbes caractéristiques : on
parle généralement de courbes Intensité-Durée-Fréquence (IDF) (Fig. 3.5). La notion de
fréquence est en faite exprimée par la notion de temps de retour.
Les courbes IDF ne sont pas une fin en soi, mais sont construites dans un but bien précis.
Elles permettent d'une part de synthétiser l'information pluviométrique au droit d'une station
donnée et, d'autre part de calculer succinctement des débits de projet et d'estimer des débits
de crue ainsi que de déterminer des pluies de projet utilisées en modélisation hydrologique.
Les courbes IDF sont établies sur la base de l'analyse d'averses enregistrées à une station au
cours d'une longue période. Les courbes obtenues peuvent donc être construites de manière
analytique ou statistique.
1. Représentation analytique
Différentes formules sont proposées pour représenter l'intensité critique d'une pluie en
fonction de sa durée.
La forme la plus générale (avec T variable) est la suivante :
(3.4)
Avec :
i : intensité totale [mm/h], [mm/min] ou intensité spécifique [l/[Link]],
T : période de retour en années,
t : durée de référence [h ] ou [min],
k, a, b, c : paramètres d'ajustement.
Montana suggère une formulation plus simple :
(3.5)
Avec :
i: intensité maximale de la pluie [mm/h],
t: durée de la pluie [minutes ou heures],
T; intervalle de récurrence (ou temps de retour) [années],
a,b: constantes locales, dépendant généralement du lieu (0.3<0.8).
Pour une fréquence de dépassement donnée, cette formule de Montana a été adaptée pour la
Suisse et a abouti à la formulation suivante (Bürki et Ziegler, 1878) :
(3.6)
où a est une constante définie localement et qui varie selon le lieu géographique de
l'observation.
La figure 3.6 (a, b et c) représente les courbes IDF, calculées pour les différentes
régions de la Suisse qui figurent dans les normes suisses pour la construction routière
(Norme Suisse SNV 640-350 p 5). Les courbes sont présentées selon un découpage
suisse en zones d’égale intensité, et un tableau dans lequel figurent les valeurs de K
en fonction du temps de retour, ainsi que les valeurs de B pour ces différentes zones.
Elles s'expriment par :
(3.7)
Où :
r : intensité spécifique moyenne d'une pluie d'une durée de t minutes atteinte ou dépassée en
moyenne une fois toutes les T années [l/s/ha],
K : coefficient fonction du lieu et du temps de retour et B une constante de lieu [min].
Fig. 3.6.a - Courbes IDF pour les différentes régions de la Suisse, déterminées par l'EAWAG
Fig. 3.6 b. - Délimitation des zones d’intensité égale pour la Suisse (SNV 640-350)
Fig. 3.6.c. - Tableau des coefficients K et B utilisées dans les calculs des courbes IDF selon
la norme SNV
2. Représentation statistique
Les courbes IDF sont établies sur la base de l'analyse d'averses enregistrées à une station au
cours d'une longue période. L'analyse fréquentielle peut s'appliquer si on ne présuppose pas
une loi connue (de type Montana, etc.) et si on s'intéresse à des événements rares, donc
extrêmes. Les données recueillies sont alors ajustées, à un pas de temps choisi, à une loi
statistique qui doit décrire relativement bien la répartition des extrêmes. La loi de Gumbel est
la plus utilisée. Si l'opération est répétée sur plusieurs pas de temps, on obtient la variation
de l'intensité avec la durée de la pluie pour différents temps de retour, c'est à dire des
courbes IDF de la station considérée sur la période analysée.
La structure d'une averse est définie comme la distribution de la hauteur de pluie dans le
temps. Cette distribution influence de manière notoire le comportement hydrologique du
bassin versant.
Fig. 3.7 Exemple de hyétogramme et de structure correspondante pour une précipitation enregistrée
au nord de Lausanne du 13 novembre 1991 à 12h00 au 14 novembre à 12h00
3.4 Evaluation régionale des précipitations
Le passage des mesures ponctuelles des précipitations à une estimation spatiale de celles-ci,
souvent nécessaire en hydrologie, est délicat. Les méthodes les plus simples et les plus
couramment utilisées sont les méthodes de calcul de moyennes ou les méthodes
d'interpolation des données pluviométriques collectées localement. Ces méthodes permettent
notamment le calcul des lames d'eau moyennes à l'échelle du bassin, la cartographie des
précipitations, et le calcul de hyétogrammes moyens. Des méthodes faisant appel à la notion
d'abattement des pluies existent également.
3.4.1 Passage des pluies ponctuelles aux pluies moyennes sur une surface
Parmi les méthodes généralement proposées pour calculer la moyenne des pluies à partir de
l'ensemble des mesures ponctuelles obtenues à plusieurs stations pluviométriques sur le
bassin ou à proximité, on distingue la méthode de la moyenne arithmétique, la méthode des
polygones de Thiessen ou l'utilisation d'isohyètes. Le choix de la méthode dépendra
notamment de la longueur de la série de données dont on dispose, la densité du réseau de
mesure, et la variation du champ pluviométrique.
La méthode la plus simple qui consiste à calculer la moyenne arithmétique des valeurs
obtenues aux stations étudiées, s'applique uniquement si les stations sont bien réparties et
si le relief du bassin est homogène.
Cette méthode est souvent peu recommandée car peu représentative. Il faut lui préférer des
méthodes graphiques (tracé d'isohyètes) ou statistiques qui permettent de donner un poids
différent à chacun des points de mesures (moyennes pondérées).
La méthode du polygone de Thiessen est la plus couramment utilisée, parce que son
application est aisée et qu'elle donne en général de bons résultats. Elle convient notamment
quand le réseau pluviométrique n'est pas homogène spatialement (pluviomètres distribués
irrégulièrement).
Cette méthode permet d'estimer des valeurs pondérées en prenant en considération chaque
station pluviométrique. Elle affecte à chaque pluviomètre une zone d'influence dont l'aire,
exprimée en %, représente le facteur de pondération de la valeur locale. Les différentes zones
d'influence sont déterminées par découpage géométrique du bassin sur une carte
topographique 1 (voir figure 3.8). La précipitation moyenne pondérée Pmoy pour le bassin, se
calcule alors en effectuant la somme des précipitations Pi de chaque station, multipliées par
leur facteur de pondération (aire Ai), le tout divisé par la surface totale A du bassin. La
précipitation moyenne sur le bassin s'écrit :
(3.8)
Avec :
Pmoy : précipitation moyenne sur le bassin,
A : aire totale du bassin (=å Ai),Pi : précipitation enregistrée à la station i,
Ai : superficie du polygone associée à la station i.
Fig. 3.8 - Exemple de détermination des polygones de Thiessen (vous pouvez changer les valeurs de
pluie ponctuelle et calculer la pluie moyenne).
1Les stations disponibles étant reportées sur une carte géographique, on trace une série de
segments de droites reliant les stations adjacentes. On élève des perpendiculaires au centre
de chacune des droites (médiatrices); les intersections de ces perpendiculaires déterminent
des polygones. Dans chaque polygone, la hauteur de précipitation choisie est celle relevée à
la station située à l'intérieur de celui-ci. Les côtés des polygones et/ou la ligne de partage des
eaux représentent les limites de l'aire (et du poids) accordée à chaque station. L'aire de
chaque polygone Ai est déterminée par planimétrie ou numériquement. D'autres critères pour
la détermination des valeurs de pondération peuvent être adoptés. Ceux-ci peuvent être
fonction de l'averse, du relief, de la position géographique, etc.
La méthode la plus rigoureuse mais qui présente l'inconvénient de demeurer lourde en dépit
des moyens actuels, est fondée sur l'utilisation des isohyètes.
Les isohyètes sont des lignes de même pluviosité (isovaleurs de pluies annuelles, mensuelles,
journalières, etc.). Grâce aux valeurs pluviométriques acquises aux stations du bassin et aux
autres stations avoisinantes, on peut tracer le réseau d'isohyètes. Le tracé des isohyètes n'est
pas unique comme celui des courbes de niveau. Il doit être dessiné avec le maximum de
vraisemblance compte tenu de la région, du réseau, de la qualité de la mesure, etc. Il existe
aujourd'hui des méthodes automatiques qui effectuent le tracé d'isovaleurs par des moyens
statistiques élaborés (technique de krigeage).
Lorsque les courbes isohyètes sont tracées, la pluie moyenne peut être calculée de la manière
suivante :
(3.9)
Avec :
Pmoy : précipitation moyenne sur le bassin,
A : surface totale du bassin,
Ai : surface entre deux isohyètes i et i+1,
K : nombre total d'isohyètes,
Pi : moyenne des hauteurs h de précipitations entre deux isohyètes i et i+1.
Fig. 3.10 - Exemple de coefficient de réduction des hauteurs de précipitations régionales en fonction
de la durée des précipitations et de la surface considérée (d'après les données de l'Atlas
hydrologique de la Suisse, zone 1).
On peut donc définir le coefficient d'abattement dit " probabiliste " comme le rapport de la
pluie moyenne de fréquence donnée à la pluie ponctuelle de même fréquence :
(3.10)
Avec :
K : coefficient d'abattement,
Pm : pluie moyenne sur la surface, de fréquence donnée,
P : pluie ponctuelle de même probabilité.
Cette définition implique qu'en chaque point, la pluie suit une même loi de probabilité. Cette
condition d'isotropisme de la pluie sur la surface est assez bien respectée pour une région
homogène et peut s'appliquer dans le cas de petits bassins versants.
La figure 4.1 représente schématiquement les différents éléments intervenant dans les
processus d'interception et d'évapotranspiration, qui font l'objet de ce chapitre.
Fig. 4.1 - Principaux éléments intervenant dans les concepts d'interception et d'évapotranspiration.
4.1.1 L'interception
Parmi les éléments de perte qui interviennent lors de l'estimation d'un bilan hydrologique
d'un bassin versant, il faut mentionner qu'une partie non négligeable de l'eau des
précipitations n'atteint pas le sol. Cette eau peut être en effet interceptée par des obstacles
au cours du trajet vertical mais aussi horizontal de l'eau. On sait aujourd'hui qu'il existe
ainsi un mécanisme d'interception horizontal des brouillards ou des rosées qui prend toute
son importance dans certaines régions du globe (e.g. les forêts situées à proximité de la côte
chilienne).
Dans ce chapitre, nous abordons exclusivement les notions liées à l'interception verticale des
précipitations, définie comme la fraction de l'eau qui n'atteint jamais le sol. On comprend
déjà ici que l'interception telle qu'elle est définie par les hydrologues est l'interception
évaporée. C'est pourquoi les auteurs anglo-saxons parlent le plus souvent de « interception
losses » soit littéralement « pertes par interception ». De façon analytique, les pertes par
interception s'expriment par la relation suivante :
Où :
I : interception (pluie n'atteignant jamais le sol) [mm],
Pi : pluie incidente [mm],
Ps : pluie atteignant le sol drainé au travers du couvert végétal (canopée) [mm],
Pt : pluie atteignant le sol par écoulement le long des branches et des troncs [mm].
Dans la troposphère, soit la couche de l'atmosphère au voisinage du sol (son épaisseur est de
2 à 3 kilomètres environ), l'air ambiant n'est jamais sec mais contient une part plus ou moins
importante d'eau sous forme gazeuse (vapeur d'eau) qui est fournie par :
♣ L'évaporation physique au-dessus des surfaces d'eau libre (océans, mers, lacs et
cours d'eau), des sols dépourvus de végétation et des surfaces couvertes par de la
neige ou de la glace.
♣ La transpiration des végétaux qui permet à la vapeur d'eau de s'échapper des plantes
vers l'atmosphère.
Fig. 4.2 - Importance relative (en %) de l'évapotranspiration (ET) par rapport à la précipitation
incidente (P) à différentes échelles spatiales.
C'est par le mouvement des molécules d'eau que débute l'évaporation. A l'intérieur d'une
masse d'eau liquide, les molécules vibrent et circulent de manière désordonnée et ce
mouvement est lié à la température : plus elle est élevée, plus le mouvement est amplifié et
plus l'énergie associée est suffisante pour permettre à certaines molécules de s'échapper et
d'entrer dans l'atmosphère. Dalton (1802) a établi, suite à des travaux sur le sujet, une loi
qui exprime le taux d'évaporation d'un plan d'eau en fonction du déficit de saturation de l'air
(quantité d'eau es-ea que l'air peut stocker) et de la vitesse du vent u. Cette loi est formulée
selon la relation suivante :
(4.2)
Avec :
E : taux d'évaporation (ou flux d'évaporation ou vitesse d'évaporation),
ea : pression effective ou actuelle de vapeur d'eau dans l'air,
es : pression de vapeur d'eau à saturation à la température de la surface évaporante,
f(u) : constante de proportionnalité (avec vitesse du vent u).
Cette relation exprime aussi que, en théorie et dans des conditions de pression et de
température données, le processus d'évaporation est possible jusqu'à ce que la pression de
vapeur effective atteigne une limite supérieure qui n'est autre que la pression de vapeur
saturante (l'évaporation cesse dès que es = ea). Ainsi, pour qu'il y ait évaporation, il faut que
le gradient de pression due à la vapeur d'eau soit positif.
(4.3)
La quantité d'eau pouvant être évaporée à partir d'une surface dépend de la quantité de
chaleur provenant du soleil. Cette quantité de chaleur varie, d'une part, selon les conditions
géographiques (gradient de latitude), et d'autre part, selon l'élévation de la surface liquide par
rapport au niveau de la mer (gradient altimétrique). Les échanges de chaleur entre
l'atmosphère, la surface du sol et la surface des lacs et des océans qui sont les agents de
l'évaporation, s'effectuent par convection et conduction. Cette énergie échangée est, en tous
points, compensée par un transfert d'eau qui s'évapore à un endroit pour se condenser à un
autre et retomber sous forme de précipitations. Ces échanges de chaleur entretiennent le
cycle de l'eau.
Les mouvements horizontaux et verticaux qui brassent l'atmosphère mettent en jeu des
échanges et des transformations d'énergie. L'une des causes fondamentales de cette agitation
réside dans la distribution des températures à la surface terrestre ainsi qu'au sein de
l'atmosphère elle-même. L'évaporation est donc fonction des rapports énergétiques entre
l'atmosphère et le plan d'eau évaporant.
réflexion diffuse, cette proportion pouvant atteindre 80% lorsque le ciel est couvert. Le
rayonnement solaire global atteignant la surface du sol comporte ainsi deux
composantes, d'une part le rayonnement solaire incident transmis par l'atmosphère,
et d'autre part le rayonnement solaire diffus réfléchi par l'atmosphère en direction du
sol (Fig. 4.3). Cette énergie solaire arrivant sur terre est essentiellement constituée de
rayonnements de courtes longueurs d'ondes (0,1 à 10 mm) (hautes fréquences).
Le rayonnement global est partiellement réfléchi par la surface du sol, selon la nature,
la couleur, l'inclinaison ou encore la rugosité de celui-ci. On définit l'albédo comme le
pourcentage de lumière solaire réfléchie à la surface terrestre pour une zone irradiée.
L'albédo varie considérablement suivant divers composantes terrestres (tableau 4.1)
ou atmosphériques et climatiques (nuages, angle d'incidence du rayonnement solaire,
saison et moment de la journée, etc.). De même, certaines composantes de
l'atmosphère telles que les poussières modifient l'albédo du globe.
Tableau 4.1. - Valeurs de l'albédo pur différentes surfaces.
Etant donné les températures régnant dans les différentes couches atmosphériques,
ces dernières émettent un rayonnement de grandes longueurs d'onde, comprises entre
5 et 100 m (infrarouge). Cette émission est due essentiellement à la vapeur d'eau, au
gaz carbonique ainsi qu'à l'ozone.
Les aérosols, poussières, cristaux, etc., en suspension interviennent également dans
les transferts radiatifs atmosphériques tandis que la présence de nuages accroît
sensiblement l'importance de l'émission atmosphérique.
(4.4)
Avec :
RS : rayonnement solaire direct et diffus atteignant le sol [Wm-2],
RA : rayonnement atmosphérique dirigé vers le sol [Wm-2],
RT : rayonnement terrestre [Wm-2],
a : albédo de la surface.
(4.5)
Avec :
RN : rayonnement net à la surface de la terre.
C : flux de chaleur dans le sol rendant respectivement compte du transport de chaleur par
conduction dans le sol ainsi que du transfert d'eau sous la forme de vapeur.
S : flux de chaleur sensible dû à la convection thermique au voisinage de la surface
évaporante. Ce flux de chaleur n'affecte donc qu'une hauteur restreinte de l'atmosphère.
L : flux de chaleur latente ou flux de vapeur à l'interface sol-atmosphère dû à la vaporisation
lorsqu'il n'y a pas d'accumulation de vapeur au sein du couvert végétal ou au sein de la
tranche d'eau qui recouvre le sol.
Cette équation traduit simplement le fait que le flux d'énergie perdu par la surface de la terre
au cours de l'évaporation soit égal au flux apporté par rayonnement, diminué du flux
d'énergie perdue par convection dans l'air et dans le sol.
(4.6)
Avec :
ea : pression de vapeur d'eau effective ou actuelle,
es : pression de vapeur d'eau à saturation.
L'humidité relative est donc le rapport entre la quantité d'eau contenue dans une masse d'air
et la quantité maximale d'eau que peut contenir cette masse d'air. Ainsi, lorsqu'une masse
d'air se refroidit, elle garde la même quantité d'eau. Par contre, la valeur de sa quantité
maximale diminue avec la température. Cette diminution implique qu'à un certain moment,
l'air devient saturé car Hr = 100%. On nomme la température pour laquelle la pression de
vapeur saturante est = à la pression de vapeur actuelle la température du point de rosée. On
exprime parfois l'humidité de l'air en kg d'eau par kg d'air humide (humidité spécifique) ou
encore en gramme d'eau par m3 d'air humide (humidité absolue). La figure 4.5 ci-après donne
une illustration des relations entre pression de vapeur, température et humidité relative.
4) Pression atmosphérique
La pression atmosphérique représente le poids d'une colonne d'air par unité de surface
considérée. Elle constitue un indicateur de la variation des types de masse d'air passant au-
dessus d'un point donné et intervient dans le calcul des humidités spécifique et absolue.
Plus la pression totale au-dessus d'un liquide est élevée, plus grande est sa tension de
vapeur ; mais cet effet reste négligeable pour des pressions totales inférieures à 106 Pa (ou
10 bars). Par contre, certains auteurs considèrent que le taux d'évaporation augmente
lorsque la pression atmosphérique diminue. Cette relation inverse n'est pas encore
clairement démontrée, car la variation de la pression barométrique est généralement suivie
d'autres variations, comme celles de la température et du régime du vent.
5) Le vent
Le vent joue un rôle essentiel sur les processus d'évaporation car c'est lui qui permet, par le
mélange de l'air ambiant, de remplacer au voisinage de la surface évaporante, l'air saturé par
de l'air plus sec. En effet, l'air au voisinage de la surface évaporante va se saturer plus ou
moins rapidement et par conséquent stopper le processus d'évaporation. Un verre d'eau placé
dans une enceinte fermée à l'abri de tout mouvement de l'air ne pourrait évaporer son
contenu bien longtemps même dans une atmosphère extrêmement sèche. Le vent, par le bais
de sa vitesse mais aussi de sa structure verticale et de ses turbulences, joue un rôle
prépondérant dans le processus d'évaporation. Les turbulences permettent entre autre
l'ascension de l'air humide, tandis que l'air sec descend et se charge d'humidité.
Les facteurs physiques qui affectent l'évaporation d'une surface dépendent étroitement des
propriétés de cette surface et sont donc variables selon qu'il s'agit de l'évaporation à partir
d'une surface d'eau libre, d'un sol nu ou d'une surface recouverte de neige ou de glace.
L'évaporation d'une surface d'eau libre dépend non seulement de propriétés physiques et
géométriques de cette surface (profondeur, étendue) mais aussi des propriétés physiques de
l'eau (outre la température déjà évoquée ci-dessus, on peut citer la salinité).
de par son inertie thermique, présentera une réponse évaporative nettement différente.
Cependant, les volumes totaux évaporés peuvent être sensiblement les mêmes dans les
deux cas.
♣ Etendue - L'étendue de la surface d'eau libre joue un rôle important sur les quantités
évaporées puisque l'évaporation, à vitesse du vent égale, est proportionnelle à la surface
évaporante ainsi qu'à l'humidité relative.
L'évaporation d'un sol nu est conditionnée par les mêmes facteurs météorologiques que ceux
intervenant dans l'évaporation d'une surface d'eau libre. Toutefois, si la quantité d'eau à
disposition n'était pas un facteur limitant dans le cas de l'évaporation à partir d'une surface
d'eau libre, elle le devient dans la situation d'un sol nu. En résumé, l'évaporation d'un sol nu
est donc influencée d'une part par la demande évaporative mais aussi par la capacité du sol à
répondre à cette demande et sa capacité à transmettre de l'eau vers la surface, fonction de
diverses caractéristiques.
♣ Capillarité - Dans le cas où le sol est relativement peu humide et dans la situation
d'un sol nu en l'absence de nappe, le régime d'évaporation est fixé par la plus petite
des contraintes entre les conditions météorologiques et la capacité du sol à
transmettre de l'eau vers sa surface. Dans ce cas, les remontées capillaires
permettent d'amener de l'eau jusqu'au front d'évaporation.
♣ Couleur du sol et albédo - Les sols de couleur claire présentant des valeurs
d'albédo élevées vont absorber moins de rayonnement que des sols foncés. Toutefois,
dans le cas où la quantité d'eau n'est pas un facteur limitant, les écarts entre
l'évaporation d'un sol clair et celui d'un sol foncé ne sont généralement que de l'ordre
de quelque pour cent, l'avantage étant donné au sol foncé.
3) Evaporation de la neige
Dans le cas de la neige et de la glace, le processus mis en œuvre est celui de la sublimation.
On estime toutefois qu'en règle générale, les quantités évaporées à partir d'une surface
couverte de neige sont assez faibles puisque la neige fond à zéro degré et qu'à cette
température, la pression de vapeur saturante est faible. L'évaporation doit donc cesser
lorsque le point de rosée atteint cette température et la fonte de la neige prend le relais sur le
processus d'évaporation.
Les différentes méthodes pour évaluer le taux d'évaporation sont soit directes (bacs
évaporants, etc.) ou indirectes (méthodes faisant appel au bilan d'énergie, d'eau ou au
transfert de masse) soit elles utilisent des formules empiriques. Les méthodes de type direct
et indirect sont abordées dans le chapitre 7 « mesures ». La plupart des formules empiriques
reposent sur des relations entre l'évaporation à un endroit donné et les facteurs
atmosphériques responsables de celle-ci. Elles sont pour la plus part établies d'après
l'équation de Dalton présentée au début du chapitre. Elles permettent toutes d'évaluer
l'évaporation et ne prennent donc pas en considération les effets dus à la présence de
végétation. Nous ne présenterons ici que trois formulations empiriques et pseudo-empiriques
du taux d'évaporation :
(4.7)
Avec :
E : évaporation physique d'un grand réservoir [mm],
Hr : l'humidité relative [%],
N : durée d'insolation effective pendant la période de calcul [h],
nj : le nombre total de jour de la période considérée.
♣ La formule de Rohwer :
(4.8)
Où :
E : pouvoir évaporant de l'air [mm],
u : vitesse du vent [m/s],
es : pression de vapeur saturant [kPa],
ea : pression de vapeur actuelle de l'air [kPa].
♣ La formule de Penman :
(4.9)
avec
Où :
E : évaporation physique d'un grand réservoir [mm],
Les formulations et les valeurs (tables) des différentes constantes météorologiques citées ci-
dessus peuvent être consultées sur le site de la FAO (Food and agriculture Organization of
the United Nation) aux adresses suivantes :
♣ formules :[Link]
♣ tables des valeurs :[Link]
La transpiration peut se définir comme l'émission ou l'exhalation de vapeur d'eau par les
plantes vivantes. La plante prélève l'eau du sol par l'intermédiaire de ses racines munies de
cellules épidermiques. Le développement du système radiculaire est lié à la quantité d'eau
disponible dans le sol ; les racines peuvent atteindre des profondeurs très variables, d'une
dizaine de centimètres à plusieurs mètres. L'absorption de l'eau est réalisée par osmose ou
par imbibition. L'eau circule à l'intérieur des canaux du système vasculaire de la plante pour
atteindre les feuilles. Le siège de l'évaporation se situe alors essentiellement au niveau des
parois internes des stomates. Une certaine évaporation peut se produire directement au
travers de la cuticule des feuilles (Fig. 4.6).
.
Fig. 4.6 – Représentation schématique du chemin de l'eau à travers la plante.
Pour la culture de référence, en l'occurrence le gazon, on a donc : ETR <= ETM <= ET0.
Pour tous les autres végétaux, seule la relation ETR<=ETM est toujours valable tout au long
de l'année.
On peut également noter qu'il existe deux résistances aux flux évaporatoires à partir d'un
couvert végétal, d'une part une résistance aérodynamique et d'autre part une résistance de
surface, toutes deux dues à la présence de la végétation (Fig. 4.7).
(4.10)
Avec :
ra : résistance aérodynamique [s/m] ; : constante de von Karman ( =0.41) ;
: vitesse du vent [m/s] ; z0 : hauteur de frottement [m],
z : hauteur de l'anémomètre (= h +2 où h est la hauteur de la végétation en m) [m],
d0 : translation du plan origine de la relation logarithmique entre la vitesse du vent et la
hauteur [m].
Un second élément qui prend aussi la forme d'une résistance est la résistance de surface (rs)
ou résistance de la canopée (Fig. 4.7). Elle représente la contrainte physiologique imposée par
la végétation au mouvement de l'eau à travers de ses stomates. L'importance de la résistance
de la couverture végétale a été montrée dans l'étude des processus d'évaporation d'un couvert
végétal humide. En effet, lorsque le couvert végétal est humide, les pertes par évaporation
sont essentiellement contrôlées par la canopée et non plus par le rayonnement car la
végétation agit comme un puits pour le transfert d'énergie par advection. L'évaporation de
l'eau entraîne un gradient thermique entre l'air ambiant et le végétal suffisant pour fournir
un flux de chaleur. Ceci est aussi confirmé par les quantités non négligeables qui sont
évaporées durant la nuit. Ces divers processus sont conditionnés par les valeurs de la
résistance de la canopée.
L'évapotranspiration d'un sol couvert par de la végétation est difficile à estimer. Pour faciliter
la tâche et dans un souci d'homogénéisation des modèles, les chercheurs sont arrivés à
déterminer les besoins en eau des cultures, équivalent à l'ETM, par la correction de
l'évapotranspiration potentielle (ET0) d'une culture de référence, qui est normalement le
gazon, par un coef appelé "coefficient cultural" (kc) en utilisant la formule suivante (Fig. 4.8) :
L'échelle de temps sur laquelle les besoins sont calculés peut être l'heure, la journée, la
décade, le mois ou la phase de croissance, selon l'objectif poursuivi et la disponibilité de
données. La valeur du coefficient kc est largement affectée par la nature de la culture, sa
hauteur, sa durée de cycle, et son taux de croissance, mais aussi par la fréquence des pluies
ou de l'irrigation au début du cycle de la culture. kc est toujours établi expérimentalement au
début, pour une région et une culture données, puis ensuite confiné dans des tables pour
une utilisation ultérieure dans la même région ou dans une région similaire. Les valeurs du
coefficient kc sont théoriquement comprises entre 0 et 1, selon le stade de la culture.
Fig. 4.8 - Besoin en eau des cultures (ETM) et évapotranspiration de référence (ET0).
(d'après FAO, 1998 et modifié)
La détermination de l'ET0 peut être faite : soit directement à l'aide des lysimètres (cf. chapitre
7 « mesures ») ; soit indirectement à l'aide de formules empiriques et théoriques (ou à bases
physiques) qui combinent des variables climatiques.
La formule de Turc (1961) est en revanche une relation qui peut être appliquée dans les
régions tempérées pour estimer l'évapotranspiration de référence. Elle s'écrit dans son
expression mensuelle ou décadaire : ¨
(4.12)
(pas de temps mensuel)
(4.13)
(pas de temps décadaire)
Avec :
t : température moyenne de la période considérée t en [°C],
ET0 : évapotranspiration de référence mensuelle ou décadaire [mm],
RG : rayonnement global mensuel ou décadaire [cal/cm2/jour].
Cette formulation est très simple d'emploi mais ne permet pas de prendre en compte les effets
du vent. De plus, elle n'est pas applicable à des échelles de temps réduites (pas de temps
horaire ou journalier) qui sont justement celle qui intéresse l'ingénieur lors de projets
d'irrigation.
(4.14)
Où :
ET0 : évapotranspiration de référence calculée par la relation de Penman
[mm/s],
Rn : rayonnement net [W/m2],
Pour l'exécution pratique des calculs, certaines grandeurs définies ci-dessus sont considérées
comme constantes et certaines sont à calculer sur la base des données météorologiques
disponibles (en règle générale : la température, la vitesse du vent, la pression, le
rayonnement global, l'humidité et l'albédo). Les valeurs des différentes constantes
météorologiques citées ci-dessus peuvent être consultées dans des tables sur le site de la
FAO à l'adresse suivante :[Link]
Une fois ces valeurs précisées, on peut déterminer la résistance aérodynamique (comme
effective dans l'air (en kPa et avec la température en degrés Celsius) etc. On a :
(4.15)
(4.16)
avec en [kPa/C°2]
(4.17)
En conséquence, on retiendra ici qu'il est possible d'estimer l'évaporation ainsi que
l'évapotranspiration de référence par le biais de formulations plus ou moins complexes qui
requièrent toutes la connaissance d'un certain nombre de paramètres climatiques. Ce sont,
en fin de compte, la disponibilité en données météorologiques qui conditionnera le choix
d'une formulation au détriment d'une autre ainsi que ses possibilités d'application pour la
région d'étude concernée.
4.3 L'interception
Les pertes importantes, qui résultent du phénomène d'interception par le couvert végétale
mais aussi par les surfaces plus ou moins perméables comme les constructions, les routes (à
ne pas négliger en hydrologie urbaine), vont dépendre, comme pour l'évaporation, de facteurs
météorologiques et de la nature de la couverture du sol.
Le processus d'interception en hydrologie étant lié aux pertes par évaporation, on retrouve les
mêmes facteurs météorologiques intervenant dans ce processus que ceux que l'on a évoqués
au sujet de l'évaporation. La structure de l'épisode pluvieux va toutefois jouer un rôle essentiel
sur le processus d'interception. En effet, on admet aujourd'hui que même durant la pluie,
une fraction de l'eau interceptée peut s'évaporer. Ainsi, la durée de la précipitation va
influencer directement les volumes interceptés : si une averse est fractionnée, une plus
grande partie de l'eau interceptée mécaniquement par le feuillage peut s'évaporer par rapport
à la situation d'une précipitation constante. La figure suivante (Fig. 4.9) illustre ce
phénomène. De la même manière, brumes et brouillards concèdent une plus grande part
d'eau à l'interception que les averses orageuses.
Fig. 4.9 - Evolution du stockage sur la canopée pour deux structures de précipitations différentes. (a)
cas d'une précipitation observée non uniforme, (b) cas d'une précipitation uniforme de même durée et
de même volume total. On note ainsi l'importance de la structure des précipitations.
Les facteurs essentiels qui conditionnent les quantités d'eau pouvant être interceptées par un
couvert végétal sont les suivants :
plus, il est à noter que la végétation basse telle que des fougères ont des capacités de
stockages qui sont loin d'être négligeables puisque celles-ci peuvent être du même
ordre de grandeur que celle des feuillus.
Finalement, il est délicat de faire des comparaisons pertinentes entre les valeurs proposées
dans la littérature pour différents types de végétation, étant donné d'une part de la
complexité des processus d'interception, et d'autre part des relations entre les facteurs liés à
la végétation elle-même et ceux liés aux conditions météorologiques.
Le tableau 4.2 suivant donne toutefois quelques chiffres sur l'ordre de grandeur de
l'interception pour différentes composantes du milieu et en relation avec le type de climat.
Re-précisons les quelques définitions utilisées dans ce tableau :
♣ La précipitation au sol nette ou précipitation nette est la quantité d'eau de pluie qui
atteint effectivement la surface du sol. Elle est égale à la précipitation brute moins la
précipitation interceptée ou encore à la somme des quantités d'eau tombées
directement au sol et celles provenant de l'écoulement supercortical.
Tableau 4.2. Redistribution des précipitations brutes en différents composants, pour les
arbres individuels ou pour la végétation prise dans son ensemble, en relation avec le type
de climat, l'intensité des pluies et les espèces ligneuses
Végétation, Evénement
Pertes Précipitations Ecoulement Précipitations
pays pluvieux +
d'interception au sol supercortical au sol nettes
et source spécifications
Climat tempéré
jeunes
épicéas,
2,5mm 64%
Angleterre
17,8 mm 21%
(Jackson
1975)
bois durs, 2.5 mm
40% 60% 0 60%
Est des USA 20 mm
10% 86% 4% 90%
(Helvey et à base
13% 87%
Patric 1965 annuelle
à base
peuplements annuelle:
22% 65% 13% 78%
forestiers, - hêtre
21% 74% 6% 79%
Allemagne - chêne
23% 72% 6% 78%
(Lunt 1934) - érable
59% 40% 1% 41%
-épicéa
Climat tropical à pluviosité élevée
forêts 2.5mm 48% 52%
Les écoulements représentent une partie essentielle du cycle hydrologique. On a déjà vu que
l'eau précipitée sur un bassin versant va se répartir en eau interceptée, évaporée, infiltrée et
écoulée. La quantité d'eau collectée puis transportée par la rivière résultera des précipitations
directes à la surface même du cours d'eau et des écoulements de surface et souterrains
parvenant à son exutoire. La proportion entre ces deux types d'écoulements est définie par la
quantité d'eau infiltrée dans le sol. Les différents processus d'infiltration et d'écoulements
participant à la génération de crue sont représentés de manière schématique dans la figure
5.1. L'analyse des écoulements et la compréhension des processus générateurs font l'objet du
chapitre 10 et 11 et, par conséquent, seront traités succinctement dans ce chapitre.
5.2 L'infiltration
L'infiltration qualifie le transfert de l'eau à travers les couches superficielles du sol, lorsque
celui-ci reçoit une averse ou s'il est exposé à une submersion. L'eau d'infiltration remplit en
premier lieu les interstices du sol en surface et pénètre par la suite dans le sol sous l'action
de la gravité et des forces de succion. L'infiltration influence de nombreux aspects de
l'hydrologie, du génie rural ou de l'hydrogéologie. Afin d'appréhender le processus
d'infiltration, on peut définir :
♣ Le régime d'infiltration i(t), nommé aussi taux d'infiltration, qui désigne le flux
d'eau pénétrant dans le sol en surface. Il est généralement exprimé en mm/h. Le
régime d'infiltration dépend avant tout du régime d'alimentation (irrigation, pluie), de
l'état d'humidité et des propriétés du sol.
♣ L'infiltration cumulative, notée I(t), est le volume total d'eau infiltrée pendant une
période donnée. Elle est égale à l'intégrale dans le temps du régime d'infiltration (Fig.
5.2).
(5.1)
Avec :
I(t) : infiltration cumulative au temps t [mm],
i (t) : régime ou taux d'infiltration au temps t [mm/h].
Fig. 5.2 - Evolution générale du régime d'infiltration et de l'infiltration cumulative au cours du temps
(Ks = conductivité hydraulique à saturation)
♣ La pluie nette représente la quantité de pluie qui ruisselle strictement sur la surface
du terrain lors d'une averse. La pluie nette est déduite de la pluie totale, diminuée
des fractions interceptées par la végétation et stockée dans les dépressions du
terrain. La séparation entre la pluie infiltrée et la pluie écoulée en surface s'appelle
fonction de production. Ce concept est développé dans le chapitre 11 « la réponse
hydrologique ».
Fig. 5.3 - Régime d'infiltration en fonction du temps pour différents types de sol (d'après Musy, Soutter, 1991).
♣ La couverture du sol - La végétation influence positivement l'infiltration en
ralentissant l'écoulement de l'eau à la surface, lui donnant ainsi plus de temps pour
pénétrer dans le sol. D'autre part, le système radiculaire améliore la perméabilité du
sol. Enfin, le feuillage protège le sol de l'impact de la pluie et diminue par voie de
conséquence le phénomène de battance.
Finalement, les facteurs les plus influents, pour une même topographie, sont le type de sol,
sa couverture et son taux initial d'humidité.
La variabilité spatiale et temporelle de la teneur en eau dans le sol est décrite par des profils
d'infiltration, ou plus généralement profils hydriques, successifs, représentant la distribution
verticale des teneurs en eau dans le sol, à différents instants donnés. Dans un sol homogène
et lorsque la surface du sol est submergée, le profil hydrique du sol présente : une zone de
saturation, située immédiatement sous la surface du sol ; une zone proche de la saturation
appelée zone de transmission, qui présente une teneur en eau proche de la saturation et en
apparence uniforme ; et finalement une zone d'humidification qui se caractérise par une
teneur en eau fortement décroissante avec la profondeur selon un fort gradient d'humidité
appelé front d'humidification qui délimite le sol humide du sol sec sous-jacent. (Fig. 5.4) :
Fig. 5.4 - Caractéristiques du profil hydrique au cours d'une infiltration (avec (qo) teneur initiale en
eau et (qf) teneur finale) (Tiré de Musy, Soutter 1991) .
Finalement la pluie qui arrive à la surface du sol y pénètre assez régulièrement selon un front
d'humectation qui progresse en fonction des apports, selon le jeu des forces de gravité et de
succion. La figure 5.5 montre comment au cours d'une infiltration, la zone de transmission
s'allonge progressivement tandis que la zone et le front d'humidification se déplacent en
profondeur, la pente de ce dernier augmentant avec le temps.
Fig. 5. 5 - Evolution du profil hydrique au cours de l'infiltration (Tiré de Musy et Soutter 1991).
Au cours d'une averse, la capacité d'infiltration du sol décroît d'une valeur initiale jusqu'à
une valeur limite qui exprime le potentiel d'infiltration à saturation. En fait, elle diminue très
rapidement au début de l'infiltration mais par la suite, la décroissance est plus progressive et
tend en règle générale vers un régime constant, proche de la valeur de la conductivité
hydraulique à saturation. Cette décroissance, due essentiellement à la diminution du
gradient de pression, peut être renforcée entre autre par le colmatage partiel des pores et la
formation d'une croûte superficielle suite à la dégradation de la structure du sol provoquant
la migration de particules.
Si l'on compare l'intensité de la pluie et la capacité d'infiltration d'un sol, il existe deux
possibilités :
♣ Tant que l'intensité de la pluie est inférieure à la capacité d'infiltration, l'eau s'infiltre
aussi vite qu'elle est fournie. Le régime d'infiltration est dans ce cas déterminé par le
régime d'alimentation. C'est le cas au début du processus. Le temps nécessaire pour
égaler la capacité d'infiltration est variable. Il dépend principalement des conditions
antécédentes d'humidité du sol et de l'averse. Le temps requis est d'autant plus long que
le sol est sec et que le régime d'alimentation est voisin de la conductivité hydraulique à
saturation Ks.
Fig. 5.6 - Régime d'infiltration et capacité d'infiltration d'un sol (Tiré de Musy et Soutter, 1991).
Parmi les nombreux modèles existants, on peut retenir deux grandes approches, à savoir :
♣ une approche basée sur des relations empiriques, à 2, 3 ou 4 paramètres,
♣ une approche à base physique.
(5.2)
(3 paramètres)
Avec :
i(t) : capacité d'infiltration au temps t [mm/h],
io :capacité d'infiltration respectivement initiale dépendant surtout du type de sol
[mm/h],
if : capacité d'infiltration finale [mm/h],
t : temps écoulé depuis le début de l'averse [h],
g : constante empirique, fonction de la nature du sol [min-1].
L'utilisation de ce type d'équation, quoique répandue, reste limitée, car la détermination des
paramètres, i0, if, et g présente certaines difficultés pratiques.
(5.3)
Avec :i(t) :
capacité d'infiltration au temps t [mm/h],
if : capacité d'infiltration finale [mm/h],
a et b : coefficients d'ajustement.
D'autres formules peuvent être utilisées pour déterminer le régime d'infiltration de l'eau du
sol (cf. tableau 5.1). Elles font toutes appel à des coefficients empiriques à évaluer en fonction
du type de sol rencontré.
Ces modèles décrivent d'une manière simplifiée le mouvement de l'eau dans le sol, en
particulier au niveau du front d'humidification et en fonction de certains paramètres
physiques. Parmi les modèles présentés dans le tableau 5.1, les deux modèles suivants sont
les plus connus :
♣ Le modèle de Green et Ampt - Un autre modèle tout aussi connu que le précédent
est celui de Green et Ampt (tableau 5.1). Ce modèle repose sur des hypothèses
simplificatrices qui impliquent une schématisation du processus d'infiltration (Fig.
5.7).
Fig. 5.7 - Schématisation du processus de l'infiltration selon Green et Ampt (Tiré de Musy et Soutter,
1991).
Il est basé sur la loi de Darcy (cf. chapitre 6) et inclut les paramètres
hydrodynamiques du sol tels que les charges hydrauliques totales, au niveau du front
d'humidification (Hf est la somme de la hauteur d'eau infiltrée depuis le début de
l'alimentation - Zf - et de la charge de pression au front d'humidification - hf ) et en
surface (H0 = ho = charge de pression en surface).Une des hypothèses du modèle de
Green et Ampt stipule que la teneur en eau de la zone de transmission est uniforme.
L'infiltration cumulative I(t) résulte alors du produit de la variation de teneur en eau et
de la profondeur du front d'humidification. Ce modèle s'avère satisfaisant dans le cas
de son application à un sol dont la texture est grossière. Cette méthode reste
cependant empirique puisqu'elle nécessite la détermination expérimentale de la valeur
de la charge de pression au front d'humidification.
s : sorptivité [cm.s-0,5]
A : composante gravitaire fonction de la
Philip
conductivité hydraulique à saturation [cm/s]
a : constante
Dooge Fmax : capacité de rétention maximale
Ft : teneur en eau au temps t
Ks : conductivité hydraulique à saturation
[mm/h]
h0 : charge de pression en surface [mm]
hf : charge de pression au front
Green&Ampt
d'humidification [mm]
zf : profondeur atteinte par le front
d'humidification [mm]
De par la diversité de ses formes, on ne peut plus aujourd'hui parler d'un seul type
d'écoulement mais bien des écoulements. On distingue dans un premier temps deux grands
types d'écoulements, à savoir : les écoulements « rapides » et par opposition, les
écoulements souterrains qualifiés de « lents » qui représentent la part infiltrée de l'eau de
pluie transitant lentement dans les nappes vers les exutoires. Les écoulements qui gagnent
rapidement les exutoires pour constituer les crues se subdivisent en écoulement de surface
et écoulement de subsurface :
A cet ensemble de processus peut encore s'ajouter l'écoulement dû à la fonte des neiges. Les
différentes composantes de l'écoulement dans le cas simple d'une averse uniforme dans le
temps et dans l'espace, sont également représentées schématiquement dans la figure 5.9
suivante.
Ces processus qui se produisent à des vitesses très différentes, mobilisent des eaux d'âge,
d'origine et de cheminement très distincts, et permettent d'expliquer la plupart des
comportements hydrologiques rencontrés sur les bassins versants, depuis les crues de
« ruissellement pur » jusqu'aux crues où la contribution à l'écoulement final est
essentiellement hypodermique ou phréatique.
Les éléments les plus importants dans la génération des crues sont finalement les
écoulements de surface et de subsurface et les précipitations directes à la surface du cours
d'eau, l'écoulement souterrain n'entrant que pour une faible part dans la composition du
débit de crue (Fig. 5. 10).
♣ en eau qui s'infiltre et qui contribue, par un écoulement plus lent à travers les
couches de sol, à la recharge de la nappe et au débit de base,
Cependant des crues sont fréquemment observées pour des pluies d'intensité inférieure à la
capacité d'infiltration des sols. Dans ce cas, d'autres processus tel que l'écoulement sur des
surfaces saturées en eau, permettent d'expliquer la formation des écoulements. Des zones
de sol peuvent être saturées soit par contribution de l'eau de subsurface restituée par
exfiltration (d'une nappe perchée par exemple), soit par contribution directe des
précipitations tombant sur ces surfaces saturées.
Il existe ainsi deux modes principaux d'écoulement de surface qui peuvent se combiner (cf.
chapitre 10) :
Une partie des précipitations infiltrée chemine quasi horizontalement dans les couches
supérieures du sol pour réapparaître à l'air libre, à la rencontre d'un chenal d'écoulement.
Cette eau qui peut contribuer rapidement au gonflement de la crue est désignée sous le
terme d'écoulement de subsurface (aussi appelé, dans le passé, écoulement hypodermique
ou retardé). L'importance de la fraction du débit total qui emprunte la voie subsuperficielle
dépend essentiellement de la structure du sol. La présence d'une couche relativement
imperméable à faible profondeur favorise ce genre d'écoulement. Les caractéristiques du sol
déterminent l'importance de l'écoulement hypodermique qui peut être important. Cet
écoulement tend à ralentir le cheminement de l'eau et à allonger la durée de l'hydrogramme.
Lorsque la zone d'aération du sol contient une humidité suffisante pour permettre la
percolation profonde de l'eau, une fraction des précipitations atteint la nappe phréatique.
L'importance de cet apport dépend de la structure et de la géologie du sous-sol ainsi que du
volume d'eau précipité. L'eau va transiter à travers l'aquifère à une vitesse de quelques
mètres par jour à quelques millimètres par an avant de rejoindre le cours d'eau. Cet
écoulement, en provenance de la nappe phréatique, est appelé écoulement de base ou
écoulement souterrain. A cause des faibles vitesses de l'eau dans le sous-sol, l'écoulement
de base n'intervient que pour une faible part dans l'écoulement de crue. De plus, il ne peut
pas être toujours relié au même événement pluvieux que l'écoulement de surface et provient
généralement des pluies antécédentes. L'écoulement de base assure en générale le débit des
rivières en l'absence de précipitations et soutient les débits d'étiage (l'écoulement souterrain
des régions karstiques fait exception à cette règle).
L'écoulement par fonte de neige ou de glace domine en règle générale l'hydrologie des régions
de montagne ainsi que celles des glaciers ou celles des climats tempérés froids. Le processus
de fonte des neiges provoque la remontée des nappes ainsi que la saturation du sol. Selon les
cas, il peut contribuer de manière significative à l'écoulement des eaux de surface. Une crue
provoquée par la fonte des neiges dépendra : de l'équivalent en eau de la couverture
neigeuse ; du taux et du régime de fonte et finalement des caractéristiques de la neige.
L'écoulement total Et représente la quantité d'eau qui s'écoule chaque année à l'exutoire d'un
bassin versant considéré. L'écoulement est la somme des différents termes : écoulement
superficiel Es, écoulement hypodermique Eh et écoulement de base (ou écoulement
souterrain) Eb qui résulte de la vidange des nappes. L'écoulement total s'exprime ainsi :
(5.4)
Le bilan hydrologique d'un bassin versant est également caractérisé par trois coefficients
essentiels :
♣ le coefficient d'écoulement total Cet, défini par le rapport entre les quantités d'eau
écoulées et les quantités d'eau précipitées P :
(5.5)
(5.6)
(5.7)
Pour de fortes précipitations, Es >> Eh. Par ailleurs, il n'est pas toujours évident de distinguer
quantitativement sur le terrain Es et Eh. Par conséquent on adopte souvent Cr » Ces. Cr varie
en général entre 0 et 1 (voir chapitre 2) mais peut être supérieur à 1 dans le cas où des
échanges entre bassins, via le système géologique, sont supposés exister (exemple des
milieux karstiques).
Les écoulements de surface transportent avec eux les produits de la désagrégations des
roches des régions hautes vers les zones basses et en définitive vers la mer. Cette section est
une introduction brève aux problématiques du transport solide dont l'étude est devenue
essentielle dans de nombreux domaines, de l'étude des processus d'érosion et de
sédimentation (dans les retenues par exemple) aux études sur la pollution des cours d'eau.
Le transport solide est par définition la quantité de sédiment (ou débit solide) transportée par
un cours d'eau. Ce phénomène est limité par la quantité de matériaux susceptible d'être
transportée (c'est à dire la fourniture sédimentaire). Il est principalement réglé par deux
propriétés du cours d'eau :
♣ Sa compétence - Elle est mesurée par le diamètre maximum des débris rocheux que
peut transporter le cours d'eau. Cette caractéristique est essentiellement fonction de
la vitesse de l'eau. Les variations de la compétence en fonction de la vitesse et la
granulométrie du substrat ont été étudiées par Hjulstrom (Fig. 5.11).
Ces deux propriétés du cours d'eau ne sont pas directement liées. Ainsi dans un fleuve, la
compétence décroît vers l'aval, ce qui n'est pas le cas de la capacité.Le transport des
sédiments par les cours d'eau est donc déterminé par les caractéristiques des particules
(taille, forme, concentration, vitesse de chutes et densité des particules). Ce qui permet de
distinguer :
♣ la charge de fond (bed load), formée de matériaux trop gros pour être mis en
suspension compte tenu de leur densité et de la vitesse du courant. Ces particules
roulent sur le fond ou se déplacent par saltation. Le transport par saltation
correspond à un déplacement par bonds successifs.
Fig. 5.12 - Classification des différentes "couches" de transport solide (d'après Wen Shen, & Julien, 1992).
Les principales méthodes utilisées pour évaluer ces deux charges sont décrites dans le
chapitre 7. On calcule en générale un flux de matières transportées par unité de temps que
l'on peut ramener à la surface du basin versant (transport spécifique).
Les notions d'érosion mécanique sur un bassin versant (ou prédictions des pertes en sols) et
de transport spécifique dans les fleuves (flux annuel de MES rapporté à la superficie du
bassin versant) regroupent deux processus différents. Ces deux notions permettent de
distinguer d'une part les processus de détachement et de transport de matériaux du sol
avant leur entrée dans le système "rivière" et d'autre part leur transport dans la rivière elle
même. Pour le premier point, on peut parler des agents de l'érosion qui sont principalement
les pluies, les ruissellements qui en découlent et le vent, ainsi que des facteurs qui vont
conditionner les quantités de particules arrachées : caractéristiques des pluies, des sols, de
la végétation, de la topographie et enfin les activités humaines. Les taux de particules
transportées vont à leur tour être régis par de nombreux facteurs dont la vitesse de l'eau, les
caractéristiques du lit, la granulométrie des particules... Le matériel particulaire ainsi
transporté par le cours d'eau ne reflétera qu'en partie les phénomènes d'érosion sur les
versants puisqu'une partie des sédiments arrachés au bassin pourra se déposer
(éventuellement temporairement) entre les sources d'érosion et l'exutoire du bassin de
drainage. D'autre part, l'érosion des berges pourra contribuer à la charge en suspension
mesurée dans le cours d'eau tandis que la présence de lacs, réservoirs entraînent une
sédimentation des particules. Pour ces différentes raisons, il est donc généralement admis
que le transport spécifique de matières particulaires calculé dans les fleuves ne peut être
assimilé à un taux de dénudation mécanique des versants.
Pour l'ensemble du monde, tous continents réunis la quantité totale des sédiments évacués
en suspension est aujourd'hui de 13,505 106 tonnes par an sur une aire de drainage externe
de 88,6 106 kilomètres carrés (contre 148,9 106 km2 pour l'ensemble des continents), ce qui
correspond à une transport spécifique de 152 tonnes par kilomètre carré et par an (Meade,
1983). Cependant, la distribution est très variable d'un point à l'autre. Sur les grandes îles
du Pacifique (Indonésie), le transport spécifique est près de six fois supérieure à la moyenne
globale (1 000 t . km-2 an-1). En Asie du Sud-Est, l'érosion mécanique des régions
himalayennes est également très forte (380 t . km-2 an-1 pour l'ensemble du continent).
L'Europe (50 t . km-2 . an-1), l'Afrique (35 t . km-2 . an-1), l'Australie (28 t . km-2 an-1) ne
contribuent que pour une faible part (6%) au bilan global.
(6.1)
Avec :
E : évaporation [mm] ou [m3],
I : volume entrant [mm] ou [m3],
O : volume sortant [mm] ou [m3],
DS : variation de stockage [mm].
Le stockage d'eau se présente sous différentes formes. On peut distinguer trois grands types
de réservoirs :
♣ Les dépressions de la surface du sol dans lesquelles l'eau peut s'accumuler. C'est le
stock d'eau de surface.
♣ Le sol et le sous-sol dans lesquelles l'eau est emmagasinée. C'est le stock d'eau
souterraine.
♣ Les couvertures neigeuses et glaciaires qui constituent le stock d'eau sous forme
solide.
La rétention de surface comprend toute l'eau accumulée sur, ou au-dessus du sol. Elle
comprend l'eau interceptée par le couvert végétal, l'évaporation durant les précipitations et le
stockage dans les dépressions du sol qui est le volume d’eau emmagasiné dans les petites
dépressions du sol jusqu'à leur niveau de déversement. Elle ne comprend pas la rétention
superficielle qui est la partie de la pluie qui demeure à la surface du sol durant la
précipitation et qui ruisselle ou s'infiltre quand la pluie a cessé.
Toute l'eau captée dans les dépressions de surface, des plus petites, dues à la rugosité du
sol, aux plus grandes plaines inondées, lacs, marais, étangs, etc., est désignée comme le
stock d'eau de surface (Fig. 6.1).
♣ Les petites dépressions de surface qui se remplissent dès que l'intensité des
précipitations est supérieure à la capacité d'absorption du sol. Lors d'averses
suffisamment importantes, ces dépressions sont comblées et le surplus prend part au
ruissellement de surface. Le volume total d'eau pouvant être retenu dans ces
dépressions de surface est appelé capacité de rétention de surface. Après l'averse, l'eau
emmagasinée dans ces dépressions s'infiltre dans le sol, ou est utilisée par les végétaux
ou encore s'évapore directement. Ces dépressions ne sont que de petits réservoirs
temporaires, qui peuvent cependant agir comme tampons durant une averse sur un
bassin versant.
♣ Les lacs, les étangs ou les plaines inondées sont des réservoirs d'eau de surface,
naturels ou artificiels, de volume et superficie pouvant être très importants. Ils
interviennent directement dans le bilan hydrologique par les échanges d'eau avec le sol
(relations eau de surface-nappe), en favorisant l'évaporation à leur surface ou encore, en
retardant l'écoulement en rivière par laminage. L'étude de ce type de réservoir fait appel
à la limnologie.
L'étude de l'état ou du comportement d'un lac nécessite la connaissance d'un certain nombre
de ses caractéristiques physiques dont :
Les variations de niveaux du plan d'eau sont un facteur important. Toute surface lacustre est
soumise à des variations de niveau, du fait des apports d'eau, de l'évaporation, des pompages
et des écoulements à l'émissaire. Le vent agit aussi fortement sur le fonctionnement et la
morphologie des lacs. Il crée notamment un déplacement général des eaux superficielles vers
le coté du lac sous le vent. L'amplitude de la dénivellation produite est fonction de la forme et
de la profondeur du lac. Elle est plus forte dans les lacs peu profonds et allongé. Le lac Érié
au Canada est un bon exemple; des vents de 30 nœuds soufflant dans l'axe du lac
produisent une dénivellation de 1 mètre (3 pieds). Le régime d'un lac est finalement défini par
la connaissance du niveau moyen sur quelques années et des niveaux maximum et minimum
enregistrés durant ces années.
Les apports à un lac varient généralement suivant les saisons. Ces mouvements saisonniers
sont principalement dus aux variations saisonnières du climat. Sous nos latitudes, la fonte
de la neige ou des glaciers provoque en générale une augmentation du niveau des lacs. Le lac
Léman par exemple, a un maximum en été causé par l'apport de la fonte des neiges et des
glaciers et un minimum en février-mars. Actuellement, le niveau de nombreux lacs est
toutefois régulé par des vannes à la sortie des réservoirs et les mouvements saisonniers sont
très fortement atténués, voire supprimés. Rappelons que la présence de plans d'eau influence
fortement le comportement hydrologique d'un bassin versant, notamment par leur capacité
de stockage qui a pour effet de laminer les crues (cf. chapitre 2).
Cette section s'intéresse à l'eau qui pénètre dans le sol et y séjourne, un court instant ou de
longues années (phase souterraine du cycle de l'eau). Les contraintes qui régissent la
circulation de l'eau dans toute l'épaisseur du sol et du sous-sol amène à distinguer l'eau du
sol et l'eau des réservoirs souterrains. Ces deux « compartiments » sont étudiés séparément.
Au-dessous de la surface du sol, deux zones peuvent être identifiées de haut en bas (Fig.
6.2) :
♣ la zone non saturée, système à trois phases (solide, liquide, gaz) ou seule une partie
des espaces lacunaires sont remplis d'eau, le reste étant occupé par l'air du sol,
♣ la zone saturée, système à deux phases (solide, liquide) où tous les pores sont
remplis d'eau.
♣ L'eau du sol, assimilée à celle se trouvant dans la zone non saturée. La zone de l'eau
du sol est le siège des racines des végétaux et constitue surtout une limite
supérieure importante des nappes (alimentation, évaporation) ; elle est également le
lieu de transit de matières et de substances. Ces processus font partie du continuum
sol-plante-atmosphère.
Le sol dans sa partie non saturée apparaît comme un complexe dynamique à trois phases :
liquide, solide et gazeuse. La variabilité temporelle et spatiale de la phase liquide d'un sol se
manifeste aussi bien sur le plan quantitatif que qualitatif. L'évolution de la quantité (volume)
et de la qualité (composition de l'eau) découle d'une dynamique de transferts liée aux
propriétés même de l'eau et aux caractéristiques du sol.
La teneur en eau des éléments minéraux varie généralement entre 5 et 40%. La présence de
matière organique augmente cette valeur qui peut dépasser 100% (par exemple les tourbes
où la teneur en eau pondérale peut atteindre 800%).
La variabilité spatiale et temporelle de la teneur en eau dans le sol est décrite par des profils
hydriques successifs, représentant la distribution verticale des teneurs en eau dans le sol, à
différents instants donnés. La surface comprise entre deux profils successifs, aux temps t1 et
t2, représente le volume d'eau par unité de surface stocké ou perdu dans l'intervalle de temps
(Fig. 6.3).
Fig. 6.3 - Exemples de profils hydriques au temps t1 et t2 (d'après Musy et Soutter, 1991).
La dynamique de l'eau résulte de l'action de différents champs de forces auxquelles elle est
soumise : force de gravité, de capillarité, d'adsorption, etc. On parle ainsi eau gravitaire
lorsque l'effet de la gravité est prépondérant, d'eau capillaire lorsque l'effet des forces de
capillarité prédomine, ou encore d'eau hygroscopique pour signaler la supériorité des forces
d'adsorption. Signalons cependant que cette description crée une discontinuité arbitraire
entre les diverses fractions de la phase liquide. Il est donc préférable de décrire le
comportement dynamique de la phase liquide en se basant sur les principes généraux de la
thermodynamique et donc sur une quantification de l'état énergétique de la phase liquide en
un point du sol et à un instant donné.
L'état énergétique de la phase liquide dans le sol est ainsi caractérisé par la somme de son
énergie interne (mise en jeu à l'échelle atomique), de son énergie cinétique et de son énergie
potentielle. L'énergie cinétique pouvant être négligée en raison des très faibles vitesses
d'écoulement, on ne tient compte que de l'énergie potentielle.
H=h+z (6.2)
Avec :
H : charge hydraulique [m], c'est-à-dire la pression exprimée en hauteur d'eau équivalente,
soit la pression exercée par une colonne d'eau verticale de même hauteur ;
h : charge de pression [m], c'est-à-dire la pression effective de l'eau du sol, en hauteur d'eau,
par rapport à la pression atmosphérique ;
z : charge de gravité [m], c'est-à-dire la hauteur de l'eau au-dessus du plan de référence.
La distribution des potentiels de pression, de gravité et du potentiel total dans le sol le long
d'une verticale est représentée graphiquement par des profils de charge de pression, de
gravité et de charge totale (Fig. 6.4).
Fig. 6.4 - Profils de charge de pression, de gravité et de charge totale d'un système en équilibre
hydrostatique (tiré de Musy, Soutter, 1991).
Les mouvements d'eau dans le sol, leur direction et leur importance sont régis par les
différences d'énergie potentielle totale de l'eau, celle-ci se déplaçant d'un point à énergie
élevée vers un point de plus basse énergie, pour tendre vers un équilibre.
La loi de comportement dynamique de la phase liquide d'un sol traduit l'existence d'une
relation entre les forces auxquelles est soumis le fluide et sa vitesse d'écoulement. Cette loi,
appelée, loi de Darcy propose de calculer le flux d'eau total comme le produit d'une
constante de proportionnalité (la conductivité hydraulique à saturation) et d'un gradient,
celui de la charge hydraulique en fonction de la profondeur. La loi de Darcy s'exprime comme
suit :
(6.3)
Avec :
q : flux transitant [mm/h]
H : charge hydraulique totale [m]
Deux cas sont alors à distinguer selon que l'on se situe en milieu saturé ou non. Dans le cas
d'un milieu non saturé, la conductivité hydraulique n'est plus constante ; elle varie avec la
teneur en eau q tout comme la pression effective de l'eau du sol qui est négative. Au
contraire, en milieu saturé, la pression effective de l'eau du sol est positive ; elle correspond à
la profondeur de submersion en dessous de la surface d'eau libre.
La quantification des flux se fait à l'aide de profils hydriques et repose sur l'application de
l'équation de continuité. La loi de continuité exprime que la variation de la teneur en eau
dans le temps est égale aux variations spatiales du flux :
ou encore (6.4)
Avec :
Dq : variation de la teneur en eau [m3/m3] º .100[%], valeur positive ou négative suivant que
le sol perd ou stocke de l'eau ;
Dq : variation du flux transitant [mm/h] ;
Dz : variation de la profondeur [mm] ;
Dt : variation du temps [h].
Soient deux profils hydriques mesurés respectivement aux temps t1 et t2, la variation de
stock DS entre les cotes altimétriques z1 et z2 durant l'intervalle de temps Dt = t2 - t1 est
représentée par la surface de profondeur unitaire comprise entre ces deux profondeurs et les
deux profils hydriques correspondants (Fig. 6.5). On a alors les équations suivantes :
(6.5)
(6.6)
(6.7)
Où :
qz1 et qz2 : flux d'eau moyen entre t1 et t2à travers les sections de cote respectives z1 et z2,
Dt : intervalle de temps compris entre t1 et t2,
DS z2 - z1 : surface comprise entre les deux profils hydriques et les profondeurs z1 et z2.
Fig. 6. 5 - Calcul des variations du stock d'eau du sol ( d'après Musy et Soutter, 1991 ).
La discipline des sciences hydrologiques qui s'occupent des eaux souterraines est
l'hydrogéologie. Celle-ci a pour objet d'une part la connaissance des conditions géologiques et
hydrologiques et des lois physiques qui régissent l'origine, la présence, les mouvements et les
propriétés des eaux souterraines, et d'autre part l'application de ces connaissances à la
prospection d'eaux souterraines, le captage, l'exploitation, la protection et la gestion des eaux
souterraines. L'hydrogéologie met aussi l'accent sur la relation entre les eaux souterraines et
l'environnement géologique, c'est-à-dire la chimie, les modes de migration des substances
chimiques, l'accumulation de l'eau, etc.
L'hydrogéologie se base sur l'analyse de deux entités essentielles, l'aquifère et la nappe d'eau
souterraine :
♣ Un aquifère est une formation géologique perméable (sol ou roche) dont les pores ou
fissures communiquent et sont suffisamment larges pour que l'eau puisse y circuler
librement sous l'effet de la gravité (exemples : sables, graviers, craie fissurée, grès, etc.).
L'aquifère constitue ainsi un réservoir des nappes d'eau souterraines.
♣ La nappe d'eau souterraine est constituée par l'ensemble des eaux comprises dans la
zone saturée de l'aquifère dont toutes les parties sont en continuité hydraulique.
♣ Une nappe libre est une nappe dont la limite supérieure dans la formation poreuse est à
surface libre, sans contraintes physiques. On appelle nappes phréatiques, les
premières nappes libres rencontrées. La pression exercée sur le toit de cette nappe est
égale à la pression atmosphérique.
♣ Une nappe captive est une nappe d'eau souterraine emprisonnée dans une formation
géologique perméable, entre deux formations imperméables (Fig. 6.6). L'eau contenue
dans la nappe captive est donc soumise à une pression supérieure à la pression
atmosphérique. La surface fictive de cette nappe correspondant à la surface
piézométrique est située au-dessus de la limite supérieure de l'aquifère confiné. Lorsque
la charge hydraulique est supérieure au niveau du sol, l'eau jaillit naturellement (cf.
puits artésien dans Fig. 6.6). Ce phénomène est appelé l'artésianisme et on appelle alors
ce type de nappe captive, nappe arté[Link] qu'une nappe captive présente
également une surface libre, par où l'eau peut s'infiltrer. Cette zone d'alimentation
s'appelle la surface de captage.
♣ Une nappe semi-captive appartient à un aquifère dont le toit ou/et le substratum est
constitué par une formation semi-perméable. Les échanges d'eau avec cette formation
semi-perméable superposée ou sous-jacente, réalisés dans certaines conditions
hydrodynamiques favorables (différences de charge), sont appelés drainance.
♣ Une nappe perchée est une nappe libre, permanente ou temporaire, formée dans une
zone non saturée, et qui surmonte une nappe libre de plus grande extension.
La surface d'une nappe ou surface piézométrique est la surface de la zone saturée d'un
aquifère à nappe libre, mais peut aussi correspondre au toit d'un aquifère à nappe captive.
C'est une donnée dimensionnelle importante. Sa forme permet d'étudier les caractéristiques
de l'écoulement des eaux souterraines et la réserve de la nappe. Dans un aquifère à nappe
libre, elle ne doit pas être confondue avec la surface libre, dont elle diffère dès que la frange
capillaire saturée n'est plus négligeable.
La surface libre d'une nappe correspond au lieu des points d'une nappe où la pression de
l'eau est égale à la pression atmosphérique. Celle-ci est un cas particulier de surface
piézométrique (surface d'équipression)
On distingue ainsi :
♣ La transmissivité est le débit d'eau qui s'écoule d'un aquifère, par unité de largeur,
sous l'effet d'une unité de gradient hydraulique. Elle est égale au produit de la
conductivité hydraulique à saturation et de la puissance (hauteur) de la nappe.
Rappelons que l'écoulement de l'eau à travers les formations perméables, en milieu saturé,
est régi par la loi de Darcy. La vitesse d'écoulement de l'eau est en fait une vitesse fictive de
l'eau à travers la section totale d'écoulement. La figure 6.7 montre bien que, compte tenu du
fait que la section d'écoulement n'est pas du tout celle de l'ensemble du massif sol, l'eau
devra circuler beaucoup plus rapidement dans les cheminements disponibles (effet de
tortuosité).
Le débit d'une nappe Q est le volume d'eau par unité de temps, traversant une section
transversale d'aquifère sous l'effet d'un gradient hydraulique déterminé.
Le débit d'une nappe souterraine, à travers une section de sol, peut s'exprimer par
l'équation :
Q = Ks . i . A ,
Q = Ks . i . H . l (6.8)
Q=T.i.l
Où :
Q : débit d'une nappe souterraine [m3/s] ;
Ks : conductivité hydraulique [m/s] ;
i : gradient de charge hydraulique [m/m] ;
A : section de sol [m2], A = H . l ;
H : épaisseur de l'aquifère [m] ;
l : largeur moyenne de la section d'écoulement [m] ;
T : transmissivité [m2/s].
Pour évaluer le volume des eaux souterraines, on procède soit par estimation du niveau
imperméable par une étude géologique appropriée, soit par détermination du coefficient
d'emmagasinement de la roche ou encore par des mesures des niveaux piézométriques.
La réserve exploitable d'eau souterraine d'une nappe libre ou captive est donnée par la
différence du niveau piézométrique actuel avec le niveau auquel on accepte de rabattre la
nappe, multiplié ensuite par sa surface moyenne et son coefficient d'emmagasinement.
On appelle "tarissement simple" tout tarissement de nappe, de source, de cours d'eau qui se
déroule en conditions semblables à la décharge, en régime non influencé (dû à l'apport de
pluie par exemple, pendant la période de tarissement), d'une nappe captive ou d'une nappe
libre, profonde ou phréatique. Le tarissement simple peut être décrit par différentes lois.
Nous ne développerons ici que la "loi exponentielle simple" qui est l'une des lois les plus
appliquées. Celle-ci s'exprime par la relation suivante avec le temps t en seconde :
(6.9)
Où :
Q : débit d'étiage au temps t [m3/s] ;
: coefficient de tarissement ;
(6.10)
Où : V : volume d'eau disponible contenu dans les réserves d'un bassin versant.
Dans le cas particulier d'une loi décroissante exponentielle, et en prenant t = 0, on obtient
(Fig. 6.8) :
(6.11)
Le calcul du volume d'eau disponible permet d'évaluer la possibilité du soutien à l'étiage (plus
petit débit observé dans un cours d'eau) en période sèche d'une région donnée.
La couverture neigeuse est une composante essentielle du stockage dans les régions
montagneuses. La neige accumulée sur un bassin versant constitue une réserve
potentiellement utilisable pour l'alimentation en eau d'une région et le remplissage de
réservoirs.
Sur les bassins montagneux, l'écoulement en rivière est pour une grande part composé de la
fonte de la neige. Celle-ci influence le ruissellement de surface en modifiant la surface
d'écoulement.
L'épaisseur et l'étendue du manteau neigeux peuvent être évalués par différentes méthodes :
♣ Les photographies prises par satellites (analogiques ou digitales) peuvent également être
utilisées dans la détermination générale de l'étendue du manteau nival, aussi bien en
montagne qu'en plaine.
♣ la prospection in situ reste cependant le système le plus employé pour estimer les
variations de l'épaisseur de neige.
Les mesures du manteau nival sur de grandes surfaces, combinées avec les valeurs de
densité de la neige estimées localement, permettent une évaluation de l'équivalent en eau
pour toute une région. L'équivalent en eau moyen du stock neigeux sur l'ensemble du bassin
versant peut être déduit à partir des mesures de l'équivalent en neige obtenues aux diverses
stations ou zones témoins, en appliquant par exemple une méthode de pondération de type
polygones de Thiessen.
La connaissance du volume d'eau emmagasiné sous forme de neige ne suffit souvent pas à
l'hydrologue ; il doit également estimer le temps de fonte et d'écoulement du stock neigeux.
Pendant la période de fonte, le couvert de neige est formé de deux parties distinctes, à
savoir : la partie supérieure, non saturée, qui peut tout de même contenir une certaine
quantité d'eau (l'eau s'y écoule verticalement, par percolation) et la partie sous-jacente, en
contact avec le sol, qui est constituée par de la neige saturée en eau (Fig. 6.9). Cette dernière
fournit le ruissellement superficiel qui alimente les rivières et les lacs. L'écoulement se fait
parallèlement au terrain suivant la loi de Darcy.
Fig. 6.9 - Illustration des processus d'écoulement au sein d'une couche neigeuse.
La vitesse à laquelle l'eau accumulée sous forme de neige apparaît dans les rivières n'est pas
seulement déterminée par le taux de fonte de neige, mais aussi par le temps pris par l'eau
pour atteindre ces rivières. Le couvert neigeux traversé par cette eau contrôle le type
d'écoulement et sa vitesse.
La fonte de la neige résulte d'un transfert de chaleur à la couverture neigeuse et dépend des
éléments suivants :
Le calcul de taux de fonte du manteau neigeux est un problème délicat qui nécessite de poser
différentes hypothèses simplificatrices. On admet par exemple, que la chaleur latente de la
glace est de 80 cal/g, que la neige est de la glace pure et que la température de la neige est
de zéro degré. Or, durant les mois d'hiver, il n'est pas rare de constater que cette dernière
hypothèse n'est pas respectée et que la température de la neige est négative. De plus, durant
la période de fonte, la couverture neigeuse n'est pas isotherme puisqu'une partie d'eau
liquide peut se trouver occluse dans la neige. Ce constat a conduit les scientifiques à
proposer, par analogie avec les notions de teneur en eau et de capacité de rétention du sol,
une teneur en eau de la neige ainsi qu'une valeur limite de rétention ("absence de fonte")
nommée capacité au champ de la neige. La figure 6.10 ci-après illustre ces principes en
relation avec la répartition altimétrique de la neige sur le bassin versant.
Fig. 6.10 - Distribution des conditions de fonte des neiges dans un bassin versant montagneux
(D'après Ward et Robinson, 1989).
Une méthode de calcul de fonte de neige relativement simple, originaire des Etats-Unis, est la
méthode d'indice de température ou de la méthode des degrés-jour qui relie le phénomène de
fonte à la température de l'air. Elle présente l'avantage d'utiliser des données météorologiques
généralement accessibles.
La hauteur d'eau de fonte provenant de la fonte, sur i jours, est calculée par la formule
suivante :
(6.12)
Où :
hf i jours : hauteur d'eau de fonte en i jours [cm],
k : coefficient exprimant l'influence des conditions naturelles et climatiques du bassin
(excepté la température) sur la fonte de la neige [cm/°C],
Ti : température moyenne journalière de l'air, au-dessus de zéro [°C] pour le jour j,
déterminée pour l'altitude moyenne du bassin,
To : température de référence, généralement admise comme égale à la température de
congélation [°C].
On peut distinguer deux types de couvertures glaciaires : les glaciers permanents et la glace
qui se forme au-dessus des plans d'eau (lacs et rivières).
Un glacier est défini comme une masse de glace à la surface du sol (l'hydrologue englobe
dans la notion de glacier, toutes glaces et neige pérennes), constituée de la recristallisation de
la neige ou d'autres précipitations, se déplaçant lentement vers l'aval.
Rappelons que l'équivalent en eau des glaciers ne représente que 2 % de la totalité de l'eau
du globe, mais 77 % des ressources en eau douce. La glace du globe se retrouve
essentiellement en Antarctique (13,9 106 km2 et 90% de la glace totale) et au Groenland
(1,8 106 km2 et 9% de la glace totale). Seulement 1% de glace se retrouve dans les autres
régions du globe. Toutefois, cela peut représenter une quantité de glace importante à une
échelle locale. Par exemple, le volume total des glaciers suisses actuels pourrait recouvrir
tout ce pays d'une couche de glace de 150 cm d'épaisseur, ce qui correspond à peu près aux
précipitations moyennes annuelles de la Suisse.
Le bilan annuel d'un glacier est, en général, calculé par des méthodes indirectes. Les études
glaciologiques étant très complexes et très coûteuses, on se contente, pour de nombreux
glaciers, d'observer la fluctuation de leur front. Ce calcul du bilan peut se faire selon trois
possibilités, soit par bilan d'énergie, soit par bilan hydrologique ou encore par bilan
géodésique.
Les quantités de glace recouvrant les cours d'eau, les lacs et les réservoirs peuvent causer
divers problèmes, entre autre gêner la navigation, endommager certains ouvrages ou former
des embâcles. Ces dernières peuvent par la suite générer des débâcles brutales pouvant
provoquer de sérieuses inondations.
Le régime caractérisant la formation de glace recouvrant des lacs et des rivières peut-être
estimé par les éléments suivants :
L'épaisseur de la glace est le seul élément qui peut être déterminé par des mesures, au
moyen d'une tarière de sondage et à la règle, à des endroits représentatifs de la rivière, lac ou
réservoir. Les autres caractéristiques sont évaluées visuellement.
Sur les cours d'eau et lacs importants, les observations aériennes sur la formation de la glace
ou la débâcle sont d'une grande valeur. Les données de télédétection (infrarouges), fournies
par les satellites permettent également une estimation des caractéristiques de la glace sur les
lacs et les réservoirs.
La mesure des précipitations est l'une des plus complexes en météorologie car on observe
une forte variation spatiale selon le déplacement de la perturbation, le lieu de l'averse, la
topographie et les obstacles géographiques locaux gênant sa captation.
On exprime généralement les précipitations en hauteur ou lame d'eau précipitée par unité de
surface horizontale (mm). Si on rapporte cette hauteur d'eau à l'unité de temps, il s'agit d'une
intensité (mm/h). Rappelons que :
1 mm = 1 l / m2 = 10 m3 / ha = 1 000 m3 / Km²
La précision de la mesure est au mieux de l'ordre de 0,1 mm. En Suisse, toute précipitation
supérieure à 0,5 mm est considérée comme pluie effective.
L'accumulation de la pluie dans un réservoir cylindrique est enregistrée par l'élévation d'un
flotteur. Lorsque le cylindre est plein, un siphon s'amorce et le vide rapidement. Les
mouvements du flotteur sont enregistrés par un tambour rotatif à vitesse constante, entouré
d'un papier, et déterminent le tracé du pluviogramme.
Cet appareil comporte, en dessous de son entonnoir de collecte de l'eau, une pièce pivotante
dont les deux compartiments peuvent recevoir l'eau tour à tour (augets basculeurs). Quand
un poids d'eau déterminé (correspondant en général à 0,1 ou 0,2 mm de pluie) s'est
accumulé dans un des compartiments, la bascule change de position : le premier auget se
vide et le deuxième commence à se remplir (Fig. 7.2). Les basculements sont comptés soit
mécaniquement avec enregistrement sur papier enroulé autour d'un tambour rotatif, soit
électriquement par comptage d'impulsions (par exemple système MADD) : appareil
permettant l'acquisition d'événements en temps réel, développé par l'HYDRAM en 1983. Les
pluviographes à augets basculeurs sont actuellement les plus précis et les plus utilisés (Fig.
7.3).
7.1.3 Le radar
Les erreurs instrumentales sont multiples ; elles ont presque toutes pour conséquence de
sous-estimer les quantités précipitées. On distingue :
♣ Les pertes par mouillage (environ 0,5 %) : déficit équivalent à l'eau qui humecte les
parois intérieures du pluviomètre.
Les erreurs d'observation sont en principe systématiques mais ne sont pas trop graves du
moment que l'on ne change pas d'observateur (possibilité de corrections).
♣ Les erreurs de positionnement de l'appareil (on peut avoir une bonne mesure mais
de quelque chose de "faux").
Les facteurs qui conditionnent l'évaporation sont les suivants : les rayonnements solaires et
atmosphériques, la température de l'eau et de l'air, l'humidité de l'air, la pression
atmosphérique, le vent, la profondeur et la dimension de la nappe d'eau, la qualité de l'eau et
les caractéristiques du bassin (exposition des versants au soleil, au vent, pentes, sol,...).
Certains de ces paramètres (facteurs météorologiques) sont facilement mesurables. La figure
7.4 montre une station météo équipée de l'ensemble des instruments de mesures de ces
paramètres.
On mesure couramment le rayonnement solaire parvenant au sol. Les mesures portent d'une
part sur l'intensité du rayonnement direct, et d'autre part sur le rayonnement global tant
sous forme de rayonnement diffus que sous forme de rayonnement direct. Les instruments
utilisés sont désignés sous le nom général d'actinomètres. Pour la mesure du rayonnement
net, on utilise des pyranomètres à thermopiles, à lames ou plus rarement à distillation.
Il existe plusieurs appareils, nommés héliographes, qui évaluent chaque jour la durée totale
de l'insolation pour une station. Ils déterminent la somme des intervalles de temps au cours
desquels l'intensité du rayonnement solaire direct a dépassé un certain seuil.
[Link] La température
On mesure le taux d'humidité dans l'air avec un instrument appelé hygromètre. Les plus
simples sont les hygromètres organiques. Ils sont basés sur la propriété des substances
organiques de se contracter ou se dilater selon l'humidité. Le cheveu humain, dégraissé,
s'allonge de 2,5 % lorsque l'humidité relative passe de 0 à 100 %. La lecture peut être
aisément faite sur un tambour ou sur un cadran qu'on étalonne en fonction de l'humidité
relative. L'appareil relié à un système d'enregistrement constitue un hygrographe.
[Link] Le vent
Les instruments de mesure du vent sont de deux types ; certains évaluent la vitesse, d'autres
la direction. En surface, les anémomètres mesurent la vitesse du vent. Ils sont installés à 10
mètres au-dessus du sol, à un endroit dégagé de tout obstacle (bâtiment, arbre,...). Les plus
fréquemment utilisés sont les anémomètres totalisateurs, constitués de trois ou quatre
branches terminées respectivement par une coupelle hémisphérique. Le système se rattache
aussi à un dispositif d'enregistrement pour former un ensemble appelé anémographe. Pour la
mesure en altitude troposphèrique, on se sert d'un ballon rempli d'hydrogène qui s'élève dans
l'atmosphère. Connaissant sa vitesse d'ascension et son déplacement horizontal en fonction
du temps, on calcule aisément la vitesse du vent qui l'entraîne. La direction du vent est,
quant à elle, déterminée à l'aide d'une girouette ou d'une manche à air. La direction du vent
est donnée selon les points cardinaux (cf. Fig. 7.4).
Les verrières sont constituées d'un cadre métallique sans fond, de 1 m2 de section, posé sur
le sol. Une vitre inclinée recouvre ce châssis. L'eau du sol s'évapore et la vapeur se condense
sur la paroi froide de la vitre. L'eau condensée est alors collectée par une gouttière et
recueillie dans un récipient. Ce type de mesure doit cependant subir des corrections pour
tenir compte des effets du vent et de la température à l'air libre.
[Link] Le lysimètre
Le lysimètre est une cuve étanche enterrée, à parois verticales, ouverte en surface et remplie
par une portion de terrain d'une épaisseur de 0,5 à 2 mètres. La végétation et les conditions
à chaque niveau, surtout la teneur en eau, sont maintenues sensiblement identiques à celles
du terrain en place. Les variations de stock d'eau peuvent alors être mesurées avec précision.
Le lysimètre est pourvu à sa base d'un dispositif recueillant l'eau de drainage. On peut
déduire l'évaporation à la surface du terrain de ces variations de stock par pesée, ou encore
des mesures de l'eau du sol et de drainage et des données de précipitations indiquées par un
pluviomètre à proximité. L'aire horizontale de la portion de terrain isolé doit être
suffisamment grande pour obtenir une bonne précision de la hauteur d'eau évaporée, en
théorie à 0,01 mm près.
Fig. 7.7 - Schéma d'une cuve lysimétrique (Tiré de Musy et Soutter, 1991).
La mesure de l'évapotranspiration est une mesure complexe. A l'inverse des autres termes du
bilan hydrique, elle est le plus souvent indirecte (en procédant par bilan hydrologique sur
une parcelle expérimentale ou sur un bassin versant). Cependant, la mesure de
l'évapotranspiration réelle (Etr) peut être effectuée de façon ponctuelle et directe en se basant
par exemple, sur les pertes en eau d'une case lysimétrique portant de la végétation.
Le niveau d'eau dans un canal est facilement observable, mais n'est représentatif que de la
section d'observation et peut être soumis à des modifications dans le temps. Seule la variable
débit reflète physiquement le comportement du bassin versant, et peut être interprétée dans
le temps et l'espace. Généralement, on ne dispose pas d'une mesure directe et continue des
débits mais d'un enregistrement des variations de la hauteur d'eau en une section donnée
(station hydrométrique). On passe alors de la courbe des hauteurs d'eau en fonction du
temps H=f(t) (appelée limnigramme ) à celle des débits Q=f(t) (appelée hydrogramme ) par
l'établissement d'une courbe de tarage Q=f(H) (Fig. 7.8).
Fig. 7.8 - Passage d'un limnigramme à un hydrogramme par l'intermédiaire de la courbe de tarage.
appelle jaugeage l'ensemble des opérations destinées à mesurer le débit d'une rivière. Vous
pouvez voir ici un film présentant la méthode de jaugeage (RealMedia, 2.1 Mo).
Il est nécessaire de procéder régulièrement à des vérifications de la courbe de tarage au cours
du temps, pour tenir compte d'éventuelles déficiences de l'appareil de mesure ou
modifications de la section du cours d'eau (voir Fig 7.8 bis).
Fig. 7.8 bis - Courbe de tarage pour différentes sections d'un même cours d'eau.
La mesure des hauteurs d'eau (la limnimétrie) ou de la variation d'un plan d'eau s'effectue
généralement de manière discontinue par la lecture d'une règle graduée (échelle
limnimétrique) fixée sur un support. Pour connaître en continu les variations d'un plan
d'eau, on utilise des limnigraphes qui fournissent sur un support un enregistrement continu
des variations du niveau d'eau dans la rivière en fonction du temps (enregistrement
graphique sur bande papier, enregistrement magnétique sur cassette, etc.).
[Link] Le limnimètre
mesurer des différences de hauteur d'eau lorsqu'on l'immerge verticalement dans le cours
d'eau. La pression de l'eau est transmise par l'intermédiaire d'une membrane solidaire de la
partie mobile du condensateur.
Pour mesurer le débit d'un écoulement naturel (cours d'eau, canal, dérivation...), il existe
quatre grandes catégories de méthodes.
♣ Les méthodes "hydrauliques" tiennent compte des forces qui régissent l'écoulement
(pesanteur, inertie, viscosité...). Ces méthodes obéissent aux lois de l'hydraulique.
Rappelons que la vitesse d'écoulement n'est jamais uniforme dans la section transversale
d'un cours d'eau. Le principe de cette méthode consiste donc à calculer le débit à partir du
champ de vitesse déterminé dans une section transversale du cours d'eau (en un certain
1. Le jaugeage au moulinet
La vitesse d'écoulement est mesurée en chacun des points à partir de la vitesse de rotation de
l'hélice située à l'avant du moulinet (nombre de tours n par unité de temps). La fonction v = f
(n) est établie par une opération d'étalonnage (courbe de tarage du moulinet). Suivant le mode
opératoire adopté pour le jaugeage, le moulinet peut être monté sur une perche rigide ou sur
un lest profilé appelé "saumon" (Fig. 7.13).
Dans le cas du montage sur perche, le moulinet peut être manœuvré de deux manières :
Les différents modes opératoires du jaugeage au moulinet monté sur un lest sont présentés
dans le tableau 7.1.
Tableau 7.1. – Méthodes et limites des différents modes opératoires du jaugeage au
moulinet monté sur un lest.
• Mesure à l'aide d'un canot (Fig. 7.14) Profondeur < 10 m et vitesse < 2 m/s
• Mesures à partir d'un bateau mobile Lorsque la rivière est large (> 200 m),
uniforme et sans présence de hauts-fonds
afin d'y manœuvrer facilement.
(7.1)
2. Le jaugeage au flotteur
Différents principes de mesure peuvent être mis en œuvre basés sur le développement récent
des instruments utilisant des sondes électromagnétiques. On peut citer :
La construction d'un déversoir ou d'un canal calibré (Fig. 2.15) pour la détermination des
débits d'un cours d'eau a pour but l'obtention d'une relation entre le niveau de l'eau H et le
débit Q aussi stable que possible, et en principe sans jaugeage sur le terrain. Le débit est
alors obtenu par des formules hydrauliques et par étalonnage sur modèles. Les canaux
jaugeurs et les déversoirs calibrés sont notamment utilisés dans le cas de petits cours d'eau
aux lits étroits, instables, encombrés de blocs et à faible tirant d'eau, pour lesquels
l'installation de stations à échelles limnimétriques et l'exécution de jaugeages au moulinet ne
sont pas recommandés. Leur fonctionnement obéit aux lois de l'hydraulique classique.
Cette méthode de jaugeages par dilution s'applique à des torrents ou des rivières en forte
pente où l'écoulement est turbulent ou pour lesquels on ne trouve pas de section se prêtant à
des jaugeages au moulinet.
Le principe général consiste à injecter dans la rivière une solution concentrée d'un traceur
(sel, colorant,...) et à rechercher dans quelle proportion cette solution a été diluée par la
rivière, par prélèvements d'échantillons d'eau à l'aval du point d'injection (Fig. 7.16). Cette
dilution est notamment fonction du débit, supposé constant le long du tronçon, concerné
pendant la durée de la mesure. On a la relation suivante dans laquelle le rapport C1 / C2
représente la dilution :
(7.2)
Où :
Q : débit du cours d'eau [l/s] ;
C1 : concentration de la solution injectée dans le cours d'eau [g/l] ;
C2 : concentration de la solution restante dans des échantillons prélevés à l'aval du point
d'injection dans le cours d'eau [g/l] ;
k : coefficient caractéristique du procédé et du matériel utilisé.
Les conditions suivantes sont nécessaires pour que les méthodes par intégration ou dilution
puissent être appliquée :
♣ à la hauteur des prélèvements, le mélange doit être tel qu'en chaque point de la
section du cours d'eau, doit passer la même quantité de traceur.
(7.4)
valeur moyenne .
(7.3)
Avec :
Q : débit du cours d'eau [l/s ou m3/s] ;
M : masse de traceur injecté [g] ; M = V . C1 ;
V : volume de la solution lâchée dans le cours d'eau [l ou m3] ;
C1 : concentration de la solution lâchée dans le cours d'eau [g/l] ;
: concentration moyenne du traceur dans les échantillons, obtenue par intégration [g/l] ;
C2(t) : concentration de l'échantillon prélevé au temps t [g/l];
T : durée du prélèvement [s].
Dans ce cas, on injecte en un point du cours d'eau une masse connue de sel (NaCl) diluée
dans un volume d'eau de la rivière. On place une sonde conductimétrique en aval de
l'injection, à une distance suffisamment longue pour que le mélange soit bon. La sonde
mesure la conductivité électrique de l'eau au cours du passage du nuage de sel. On peut
alors tracer la courbe conductivité en fonction du temps.
Une relation linéaire existe entre la conductivité de l'eau et sa concentration en sel dissous.
On peut donc en déduire la courbe concentration en fonction du temps. Le débit est alors
obtenu par intégration de la concentration au cours du temps.
La quantité de sédiments (ou, flux solide, charge solide 1, débit solide 2) transportée par un
cours d'eau à une section donnée pendant un temps Dt (Dt=1 jour, 1 mois, 1 année) est
composée de la charge en suspension (suspended load) et du transport de fond (glissement ou
roulement sur le fond et saltation).
Différentes méthodes de mesures sont possibles :
On s'intéresse ici particulièrement aux mesures sur les cours d'eau. Signalons que la
question sempiternelle dans tous les programmes de surveillance du transport solide est de
savoir comment peut-on estimer celui-ci avec un coût non prohibitif, sachant que le bilan
exact des matériaux en suspension transportés demeure inaccessible. Outre les erreurs
analytiques produites, la majeure source d'erreur dans la mesure de la charge solide d'un
cours d'eau est en relation avec la variabilité des concentrations en sédiment à travers le
temps et la possibilité du programme d'échantillonnage de caractériser précisément cette
variabilité. Ce dernier point peut être déterminé dans une large mesure par la fréquence
d'échantillonnage adoptée.
1. Le terme charge solide est utilisé pour une période déterminée (e.g. charge annuelle).
2. Chez les hydrologues, on parle de débit solide, qui correspond au poids total des
matériaux transportés par les cours d'eau, d'une manière ou d'une autre, passant à
travers une section par unité de temps. On l'exprime généralement en kg.s-1. On
distingue ensuite le débit solide en suspension et le débit de charriage associés aux
deux modes de transport des matériaux.
Une large gamme d'options est aujourd'hui disponible pour mesurer la quantité de sédiments
en suspension transportée par un cours d'eau. La méthode la plus rigoureuse pour obtenir
une estimation de la charge solide en suspension consiste à procéder, comme pour la mesure
du débit liquide, à une intégration de différentes concentrations et des vitesses sur plusieurs
verticales. Cette technique nécessite un matériel de prélèvement adapté aux caractéristiques
de la section de mesure. Le contrôle en continue de la charge solide est possible grâce aux
programmes d'échantillonnages intensifs avec des pompes automatiques ou, de manière
indirecte, avec l'installation de turbidimètres.
Outre les prélèvements manuels réalisés dans des récipients généralement en polypropylène,
il existe du matériel de prélèvements plus ou moins automatisé qui peut être classé en trois
catégories principales :
♣ Les préleveurs instantanés - Ils sont constitués d'un récipient largement ouvert qui
peut se refermer, de façon quasi instantanée, au moyen d'une commande appropriée.
♣ Les préleveurs à pompe - Un embout formé d'un tube métallique faiblement coudé,
fixé sur un lest ou une perche, permet d'effectuer, à l'aide d'un tuyau flexible et d'une
pompe, des prélèvements en divers points de la section de mesures.
♣ Les préleveurs par intégration - Ils prélèvent des échantillons durant un intervalle
de temps suffisamment long pour atténuer les fluctuations de concentrations. Ils
peuvent selon les modèles, fonctionner point par point ou par intégration le long
d'une verticale. Dans ce dernier cas l'échantillon recueilli permet de mesurer la
concentration moyenne pondérée par les débits sur toute une verticale. Le plus
simple d'entre eux est constitué d'une bouteille à large col fixée à une perche. A
travers le bouchon, passe un ajustage d'admission qui doit être dirigé face au courant
tandis qu'un deuxième conduit, dirigé vers l'aval, permet à l'air de s'échapper. Des
systèmes obéissant au même principe peuvent être installés sur des saumons de
lestage.
Fig. 7.18 – Bouteille de prélèvements fixée sur une perche avec système d'ajustage (d'après Nouvelot, 1993).
En considérant la section S d'un cours d'eau de largeur L, chaque verticale V peut être
définie par son abscisse l (distance à l'une des 2 rives), et sa profondeur totale P. Si en un
point d'une verticale V, situé à la profondeur p, sont mesurées à la fois la vitesse du courant
v et la concentration c de matériaux en suspension, le débit solide sur la surface dS de la
(7.5)
♣ Les nasses constituées d'une poche de grillage montée sur un cadre métallique qui
laisse passer les matières en suspension, mais retient les matériaux grossiers.
♣ Les sondeurs à ultrasons permettent de suivre le déplacement des dunes dans les
fonds sableux à faibles pentes.
♣ Hors des parcelles et des petits bassins versants dont les exutoires peuvent être
équipés de pièges ou de fosses à sédiments, la mesure du transport de fond reste
imprécise. Les dispositifs communément utilisés perturbent en effet de manière non
négligeable le régime du transport de fond.
principe de l'infiltration à charge variable. Après remplissage des deux cylindres, les
variations du niveau d'eau dans le cylindre central sont mesurées au cours du temps.
Cette méthode permet donc d'évaluer l'infiltration verticale de l'eau dans le sol.
L'humidité du sol peut être déterminée de plusieurs façons soit par méthode directe, qui
consiste à peser les échantillons avant et après étuvage, soit par des méthodes indirectes, qui
sont établies sur des relations entre les propriétés physiques (conductivité électrique,
température) ou chimiques des sols et leur teneur en eau. Afin de suivre dans le temps
l'évolution de l'humidité du sol, il est nécessaire de recourir à des méthodes indirectes qui
sont non destructives, telles que les mesures neutroniques, les mesures de conductivité
électrique ou de la constance diélectrique dans le sol.
La mesure neutronique de la teneur en eau du sol repose sur les propriétés de réflexion que
possèdent les molécules d'eau à l'égard d'un flux de neutrons. Rappelons que parmi les
divers éléments que l'on trouve dans le sol, ce sont les atomes d'hydrogène qui possèdent le
noyau dont la masse est la plus proche de celle du neutron. Les deux parties essentielles
d'une sonde à neutrons, isolées l'une de l'autre, sont l'émetteur et le détecteur de neutrons.
Elles sont fixées à un câble qui transmet les impulsions électriques émises par le détecteur à
un compteur. Le blindage (fig. 7.19) sert à neutraliser la source radioactive lors de son
transport.
Lorsque la sonde est en place dans le sol, des neutrons rapides sont émis par la source
(mélange de americium et de beryllium) dans toutes les directions. Ils se heurtent au noyau
des divers atomes qui se trouvent sur leur trajectoire et voient ainsi leur énergie cinétique et
leur vitesse diminuer progressivement. Si le sol présente une concentration d'atomes
d'hydrogène suffisante, le ralentissement des neutrons émis par la source se produit alors
qu'ils se trouvent encore à proximité de celle-ci. Les neutrons ralentis par collisions
successives se propagent dans des directions aléatoires, si bien qu'il se forme un nuage
neutronique dont la densité est plus ou moins constante. Une partie de ces neutrons, qui
dépendent de la concentration en atomes d'hydrogènes, sont renvoyés directement en
direction du détecteur en créant des impulsions. Le nombre d'impulsions pendant un
intervalle de temps est enregistré par un compteur. La conversion de la valeur enregistrée par
le compteur en une teneur en eau se fait par le biais d'une courbe d'étalonnage.
Cette technique a l'avantage de permettre des mesures rapides et répétées sur un site sans
perturbation du sol et avec une bonne précision.
La définition de la constante diélectrique relative (er) d'un matériau est le rapport entre le
potentiel mesuré entre deux électrodes dans le vide Vo et le potentiel mesuré entre ces deux
électrodes identiquement chargées et espacées, immergées dans un matériau diélectrique V.
A titre d'information, les composantes diélectriques des matériaux constituants les sols sont
indiqués dans le tableau 7.2. La constante diélectrique de l'eau est nettement plus élevée que
celle des autres constituants du sol. Par conséquent, les constantes diélectriques des sols
sont étroitement dépendantes de leur humidité.
(7.6)
Cette méthode a l'avantage d'être non destructive, facile à mettre en œuvre et nécessite que
peu d'information sur le milieu sondé. Son principal défaut est le faible volume échantillonné
par les sondes.
La compréhension des processus intervenant dans le cycle de l'eau ainsi que l'étude de leurs
variations spatiales et temporelles nécessitent de disposer de données. Celles-ci sont
essentielles et constituent un préalable à toute analyse hydrologique, que ce soit dans le but
de procéder à une étude du cycle de l'eau, d'impacts environnementaux ou pour procéder au
dimensionnement d'ouvrages hydrauliques.
De façon générale, pour permettre le passage de l'acquisition des données à leur utilisation
effective dans le cadre d'une analyse hydrologique on distingue les étapes suivantes :
acquisition, traitement, contrôle et validation, organisation, diffusion et publication. Une
collaboration entre l'Institut d'Aménagement des Terres et des Eaux de l'EPFL et des bureaux
privés a permis par exemple la création d'un logiciel complet de traitement des données
hydrologiques nommé CODEAU qui comprend une partie de ces étapes.
Le traitement des données inclut aussi le contrôle primaire des données qui comprend les
contrôles de cohérence à l'exclusion de tous traitements statistiques. Il s'agit par exemple,
dans le cas d'une acquisition manuelle des données, de les convertir en fichiers
informatiques. Dans ce cas, on procède généralement à une double saisie des données puis
les fichiers sont comparés afin de déceler d'éventuelles erreurs de saisie. Dans la situation où
l'on procède à l'acquisition de données de précipitations et de débits, on vérifie encore la
cohérence temporelles des données acquises, à savoir par exemple qu'une crue est bien la
conséquence d'un épisode pluvieux.
Avant de pouvoir exploiter les données et bien qu'elles soient dans un format adéquat, il
importe de contrôler la fiabilité et la précision de ces dernières. Le contrôle permet de valider
les données avant leur organisation au sein d'une banque de données pour leur mise à
disposition à des fins opérationnelles. Lors de cette opération, on introduit des indices de
qualité de la donnée ainsi que des indices indiquant que celle-ci est reconstituée, calculée
voire manquante. Par exemple, le logiciel CODEAU utilise pour ce faire toute une série
d'indice ou flags permettant de qualifier des données présentant une rupture de continuité,
une ou plusieurs mauvaises valeurs, des valeurs manquantes ou à vérifier etc.
La diffusion des données en vue de leur utilisation opérationnelle peut-être faite selon deux
modes :
Des erreurs peuvent en effet être perpétrées lors de l'une ou de l'autre des quatre phases du
déroulement classique des opérations, à savoir : la mesure ; la transmission de l'information ;
le stockage de l'information ; le traitement de l'information (prétraitement et analyse). Il est
donc indispensable, avant d'utiliser des séries de données, de se préoccuper de leur qualité et
de leur représentativité en utilisant diverses techniques en général de type statistique ou
graphiques.
Une erreur de mesure est définie comme étant la différence entre la vraie valeur (qui est
l'idéal recherché, mais qui n'est en principe et malheureusement jamais connue) et la valeur
mesurée. Il est commode, tant pour les présenter que pour différencier la façon de les
aborder, de considérer deux types d'erreur : les erreurs aléatoires et les erreurs
systématiques.
♣ Les erreurs aléatoires (ou accidentelles) - Elles affectent la précision des données et
sont non corrélées. Ce type d'erreur est dû à des raisons nombreuses et variées,
généralement inconnues, affectant différemment chaque mesure individuelle.
Généralement on considère que ces erreurs sont les réalisations d'une variable
aléatoire normale centrée en 0 et de variance s 2. Ces erreurs étant inévitables, il faut
en estimer l'importance afin de pouvoir en tenir compte lors de l'évaluation de
l'incertitude finale. Dans la mesure du possible, la technique de mesure induisant les
erreurs aléatoires les plus faibles devrait être préférée.
♣ les erreurs systématiques - Elles affectent la fiabilité des données et sont totalement
corrélées. On parle aussi d'inconsistance. Supposons qu'aucune erreur aléatoire
n'affecte les mesures. La différence entre la vraie valeur et la valeur mesurée, si elle
existe, est alors due à une erreur systématique. L'origine des erreurs systématiques est
le plus souvent liée à la calibration de l'appareil de mesure qui n'est pas parfaite ou à
un phénomène extérieur qui perturbe la mesure (erreur d'appareillage, changement
d'observateur…).
Selon la nature des erreurs constatées ou supposées la recherche de ces dernières fait appel
à différentes techniques et méthodes
♣ « in situ » qui consiste à vérifier sur place la manière dont les données ont été
organisées, traitées et/ou transformées.
Les calculs statistiques sont basés sur un certain nombre d'hypothèses qui doivent en
principe être vérifiées. Parmi celles-ci, citons :
♣ Les données sont homogènes - Une série de données est réputée non homogène
lorsque:
♣ La série doit être suffisamment longue - La longueur de la série influe sur les
erreurs d'échantillonage, notamment sur le calcul des moments d'ordre
supérieurs donc sur les tests inhérents à leur fiabilité.
Cet exemple est adapté de G. Saporta, Probabilités, Analyse des données et statistique,
Technip, 1990. Des relevés effectués pendant de nombreuses années ont permis d'établir que
la hauteur annuelle des pluies dans la Beauce (France) [en mm] suit une loi normale
Tableau 8.1 - Hauteur annuelle des pluies dans la Beauce (France) [en mm] de 1951 et 1959.
Année 1951 1952 1953 1954 1955 1956 1957 1958 1959
Les agriculteurs hésitant à opter pour le procédé forcément onéreux des faiseurs de pluie
tenaient pour l'hypothèse H0 et il fallait donc que l'expérience puisse les convaincre ; c'est-à-
dire que les faits observés contredisent nettement la validité de l'hypothèse H0 dite
« hypothèse nulle » (H1 est l'hypothèse alternative).
Ils choisirent a=0,05 comme niveau de probabilité, c'est-à-dire qu'ils étaient prêts à accepter
H1 si le résultat obtenu faisait partie d'une éventualité improbable qui n'avait que 5 chances
sur 100 de se produire. Autrement dit, ils admettaient implicitement que des événements
rares ne sauraient se produire sans remettre en cause le bien-fondé de l'hypothèse de départ
H0 ; ce faisant, ils assumaient le risque de se tromper dans 5 cas sur 100, cas où
précisément les événements « rares » arrivent quand même.
moyenne des observations qui nous apporte le plus de renseignements sur . est
la « variable de décision ».
Si est vraie, comme l'expérience a porté sur ans, doit suivre une loi normale :
décision la suivante : si est trop grand, c'est-à-dire si est supérieur à un seuil qui n'a
que 5 chances sur 100 d'être dépassé, on optera pour H1 avec une probabilité 0,05 de se
tromper. Si on ne pourra pas rejeter H0 faute de preuves suivantes. est appelé valeur
critique.
Il est facile de calculer la valeur critique grâce aux tables de la loi normale et on trouve :
♣ Si , conserver H0.
relevées indiquent que =610,2 mm. La conclusion était donc de conserver H0 ;c'est-à-dire
que l'insémination était sans effet notable sur le niveau des pluies : les valeurs observées
pouvaient donc être dues au hasard en l'absence de toute influence de l'iodure d'argent.
Cependant, rien ne dit que ne pas rejeter H0 mette à l'abri de se tromper : en effet, les
faiseurs de pluie ont peut-être raison, mais on ne s'en est pas aperçu.
Il y avait en fait deux manières de se tromper : croire les faiseurs de pluie, alors qu'ils
n'étaient pour rien dans le résultat obtenu (probabilité ; ne pas croire les faiseurs
de pluies, alors que leur méthode est bonne et que seul le hasard (malencontreux pour eux),
dû au faible nombre d'observations, a donné des résultats insuffisants pour convaincre les
agriculteurs.
Supposons que les faiseurs de pluie ont raison, on commet alors une erreur chaque fois que
prend une valeur inférieure à 655 mm, c'est-à-dire avec une probabilité :
Ces erreurs correspondent à des risques différents en pratique ; ainsi dans l'exemple des
faiseurs de pluie le risque de première espèce consiste à acheter un procédé d'insémination
qui ne vaut rien ; le risque de deuxième espèce à laisser perdre une occasion d'augmenter le
niveau de pluie et peut-être de récoltes plus abondantes. Dans la pratique des tests
statistiques, il est de règle de se fixer comme donnée (les valeurs courantes sont par
exemple 0,05 ; 0,01 ou 0,1) de préférence en fonction du risque de première espèce couru, ce
qui fait jouer à H0 un rôle prééminent.
Sur la base de cet exemple, la démarche d'un test peut se résumer comme suit :
On peut classer les tests soit selon leur objet (but) ou selon leurs propriétés mathématiques.
Un test est dit paramétrique si son objet est de tester certaines hypothèses relatives à un ou
plusieurs paramètres d'une variable aléatoire de loi spécifiée. Dans la plupart des cas, ces
tests sont basés sur la considération de la loi normale et supposent donc explicitement
l'existence d'une variable aléatoire de référence suivant une loi normale. La question se
pose alors de savoir si les résultats restent encore valables lorsque n'est pas normale : si
les résultats sont valables on dit que le test en question est robuste. La robustesse d'un test
par rapport à un certain modèle est donc la qualité de rester relativement insensible à
certaines modifications du modèle. Un test est dit non paramétrique s'il ne fait pas appel à
des paramètres ou d'hypothèses précises concernant la distribution sous-jacente.
Les tests sont généralement classés en quatre groupes principaux qui contiennent la plupart
des tests statistiques généralement utilisés en hydrologie.
échantillons de taille et , peut-on admettre qu'ils ont été prélevés dans une
même population indépendamment l'un de l'autre ?
répartition .
On veut tester :
Le choix de est dicté par des considérations pratiques car est trop vague
pour obtenir une région critique. Dans la pratique on se contentera de vérifier l'égalité des
(8.1)
observations par :
(8.2)
Les divers tests exposés ci-dessous vont être appliqués à l'aide du jeu de données suivant
représentant les débits de pointe annuels en m3/s de la Viège à Viège de 1922 à 1996
(tableau 8.2). Cette série de débits possède la particularité qu'un changement de nature
anthropique a eu lieu en 1964. En effet, à cette date, un barrage a été construit en amont du
point de mesure (barrage du Mattmark).
Tableau 8.2 - Débits de pointe annuels en m3/s de la Viège à Viège de 1922 à 1996.
1996 49
La figure 8.3 présente nos données sous forme de série temporelle. Les moyennes avant et
après barrage sont représentées sous forme de droite.
Fig. 8.3- Débits de pointe annuels de la Viège à Viège [m3/s] de 1922 à 1996.
Les principales caractéristiques statistiques des deux sous séries avant et après barrage ainsi
que de la série complète sont résumées dans le tableau ci-dessous (tableau 8.3) :
Tableau 8.3 - Principales caractéristiques statistiques des deux sous séries (avant et après
barrage).
écart-
série longueur moyenne médiane asymétrie kurtosis
type
1922-
1963 42 190.26 176 48.52 1.47 2.88
1964-
1996 33 117.27 110 55.51 2.47 6.54
1922-
75 158.14 150 62.99 0.90 0.96
1996
Il s'agit par exemple de tester la série de débits de pointe ci-dessus (1 paramètre) à l'endroit
où ils ont été mesurés donc à Viège (échelle locale).
♣ Tests de conformité
Exemple : On aimerait tester si la moyenne des débits de pointe de la Viège à Viège pour la
période avant barrage est égale à 200 m3/s. Pour tous les tests ci-dessous le seuil de
Comme la variance n'est pas connue, nous utilisons le test t de Student. Le score est donné
par :
= .
Comme le test est bilatéral la valeur critique est donnée par le quantile à 97,5% d'une loi t de
. On ne peut donc pas rejeter l'hypothèse nulle que la moyenne des débits est
égale à 200 m3/s.
Pour la conformité de la variance on se base sur l'étude de la distribution d'échantillonnage
H0:
H1:
Fonction
discriminante
, suit une loi du chi-carré à degrés de liberté
Exemple : Un hydrologue vous affirme que la variance des débits de pointe de la Viège à
Viège de 1922 à 1963 est de (soit ). Votre expérience d'hydrologue vous incite
à penser que celle-ci est en fait supérieure à la valeur donnée par votre collègue.
♣ Tests d'homogénéité
Exemple : Motivée par le fait qu'une intervention humaine a eu lieu sur le bassin versant de
la Viège en 1964, la série de débits de pointes est séparée en deux échantillons :
avec :
♣ Tests d'adéquation
Exemple : on se pose la question de savoir si toute notre série de débits de pointe (i.e. 1922-
1996) suit une distribution normale. Posons Z la variable aléatoire modélisant les débits.
Les deux paramètres de la loi normale sont estimés respectivement par la moyenne et la
théoriques, on obtient : La valeur critique nous est donnée par le quantile d'une
Fig. 8.4 - Superposition de la densité normale à l'histogramme normalisé (12 classes) de toute la série des débits
22 degrés de liberté ont été enlevés à la loi chi-carré car les deux paramètres de la loi normale ont été estimés.
♣ Tests d'autocorrélation
Il faut tout d'abord remarquer que la façon la plus simple et la plus immédiate d'évaluer
l'indépendance sérielle, adoptée dans de nombreux ouvrages d'hydrologie statistique,
consiste à calculer le coefficient d'autocorrélation de rang un de la série, et d'appliquer
ensuite l'un des tests « classiques » paramétriques ou non-paramétriques proposés pour le
coefficient de corrélation rx,y « standard », à savoir :
Ces trois tests ont un comportement nominal tout à fait correct, lorsqu'ils sont appliqués à
une série bivariée « classique », mais s'avèrent inapplicables dans le cas du coefficient
d'autocorrélation, comme cela ressort de certains développements. (cf Meylan P. et Musy A.,
Hydrologie Fréquentielle, Edition HGA Bucarest, 1999).
Pour une taille n « assez grande » (Anderson fixe une limite de 75 valeurs!) le coefficient
d'autocorrélation suit une loi normale de moyenne et variance :
(8.3)
Pour des tailles inférieures d'échantillons la distribution est assez compliquée. Anderson
donne par conséquent des tables des valeurs critiques du coefficient de corrélation (tableau
8.5).
Exemple : Testons la nullité du coefficient d'autocorrélation de décalage 1 sur notre première
série de débits.
Tableau 8.5 - Table des valeurs critiques du coefficient d'autocorrélation d'Anderson pour les tests
unilatéraux. Les valeurs entre parenthèses ont été interpolées par Anderson.
N a= 5 % a= 10 % N a= 5 % a= 10 %
20 0,299 0,432
Il est rappelé que les tests non paramétriques ne font pas appel à des paramètres ou
d'hypothèses précises concernant la distribution sous-jacente.
♣ Tests de conformité
Le test de conformité de la moyenne classique non paramétrique est le test de Wilcoxon pour
un échantillon (cf. S. Morgenthaler, Introduction à la Statistique, PPUR, 1997, p. 255).
Nous allons reprendre les mêmes hypothèses que dans le cas paramétrique ci-dessus (test de
Student) :
(u)
Exemple : On obtient dans notre cas W+ = 497. Pour un échantillon de taille supérieur à 15,
l'approximation normale suivante est possible :
Et = 608
Comme 497<608, on ne peut pas rejeter l'hypothèse nulle que la moyenne de ces débits de
pointe est de 200 m3/s.
♣ Tests d'homogénéité
a) Test de Wilcoxon
Pour tester l'homogénéité de données issues de deux populations on utilise les deux
statistiques équivalentes de Mann-Whitney et Wilcoxon (cf. S. Morgenthaler,
Introduction à la Statistique, PPUR, 1997, p.251), ainsi que le test de la médiane.
Exemple : test de Wilcoxon pour deux échantillons. Comme on s'attend à avoir une
diminution significative des débits après 1964, on pose les hypothèses suivantes :
et =676.
La valeur critique est 1750. Comme , on rejette l'hypothèse nulle ce qui est
conforme à notre attente.
b) Test de la médiane
variante consiste à utiliser la moyenne ). Chaque observation est affectée d'un signe + si
elle est supérieure à la médiane, d'un signe - dans le cas contraire. Tout groupe de valeurs
« + » est une séquence positive (ou suite positive), et tout groupe de valeurs « – » est une
séquence négative (ou suite négative).Il s'agit de déterminer le nombre total de séquences
Pour un seuil de signification compris entre 91 et 95%, les conditions de vérification du test
sont les suivantes :
(8.5)
et
Exemple : On veut vérifier l'homogénéité de la série des débits de pointe de la Viège sur la
période totale d'observation.
débits 240 171 186 158 … 145 155 230 270 330 55 63 49
signe + + + + … - + + + + - - -
• Tests d'autocorrélation
A la suite des travaux d'Anderson, Wald et Wolfowitz ont développé un test non paramétrique
du coefficient d'autocorrélation.
(8.6)
Pour n «suffisamment grand», cette statistique suit une distribution normale de moyenne et
de variance :
(8.7)
(8.8)
avec: .
Ces méthodes s'appliquent par exemple pour contrôler des données pluviométriques à
plusieurs stations situées dans un bassin versant.
Les deux tests développés dans ce paragraphe sont essentiellement utilisés en hydrologie.
Leur but spécifique est de comparer un ou plusieurs échantillons, acquis à des stations
voisines, afin de déceler une éventuelle inhomogénéité (dont la cause la plus classique est
une modification de l'une des stations, par exemple, le déplacement du pluviomètre).
(8.9)
Exemple : Nous voulons tester graphiquement l'homogénéité des données des débits de
pointe de la Viège à Viège. Pour cela, nous allons utiliser comme station de référence, les
débits du Rhône à Brigue (après s'être assuré que cette série ne présente pas
d'inhomogénéité).La figure 8.5 ci-dessous présente l'application de la méthode des doubles
cumuls dans ce cadre. Elle nous montre une nette cassure de pente à la station que l'on
désire contrôler à partir de 1964. Cette méthode est capable de détecter une anomalie (la
construction du barrage).
La méthode du double cumul a l'avantage d'être simple, bien connue de tous et très
rapidement réalisable. Par contre, l'interprétation des graphes obtenus n'est pas toujours
aisée et, surtout, la méthode ne propose aucune graduation en probabilité des défauts
constatés: un test, au sens statistique, est donc exclu. Enfin, elle permet de détecter l'erreur
mais pas de la corriger, du moins pas directement. Après une analyse approfondie d'une
situation critique constatée, une correction peut être toutefois envisagée.
La méthode du cumul des résidus, due à Philippe Bois de l'Ecole nationale supérieure
d'hydraulique de Grenoble, est une extension de l'idée de la méthode du double cumul, à
laquelle elle ajoute un contenu statistique autorisant la pratique d'un véritable test
d'homogénéité : c'est donc un progrès décisif.
A nouveau pour une série double de valeurs xi (série de base) et yi (série à contrôler), l'idée de
base consiste à étudier, non pas directement les valeurs xi et yi (ou Sxi et Syi) mais les cumuls
des résidus e i de la régression linéaire de y en x :
(8.10)
ou encore : - . (8.11)
Fig 8.6 - Diagramme de dispersion des débits de pointe annuels du Rhône à Brigue et de la Viège à Viège (1922-1996)
De la théorie classique de la régression il découle que la somme des résidus est nulle par
construction et que leur distribution est normale, d'écart-type :
(8.12)
(8.13)
Le report graphique des résidus cumulés Ej (en ordonnée) en fonction des numéros d'ordre j
des valeurs (en abscisse, j = 0 à n, avec E0 = 0) devrait, pour une corrélation avérée entre X et
Y, donner une ligne partant de 0, oscillant aléatoirement autour de zéro entre j = 0 et j = n, et
aboutissant à 0 pour j = n. La présence d'une inhomogénéité se manifeste par des déviations
non aléatoires autour de la valeur nulle.
Bois a décrit et testé de nombreux types d'inhomogénéités. Il a en outre montré que, pour un
niveau de confiance 1 - a choisi, le graphe des Ej en fonction de j ( j = 0 à n ) doit être inscrit
dans une ellipse de grand axe n et de demi petit axe :
(8.14)
avec est le ème quantile de la loi normale centrée réduite.
croissants: l'anomalie se situe donc autour de 1996-32=1964, ce qui a déjà été constaté avec
la méthode du double cumul.
Fig. 8.7 - Résultat du test de Bois pour les débits de pointe annuels de la Viège à Viège avec comme série
de référence les débits du Rhône à Brigue, Ellipse de confiance à 95 %.
Cette méthode est basée sur les relations existant entre certaines grandeurs hydrologiques
telles que les précipitations et les débits ou la température et l'humidité relative de l'air. Si on
suppose cette relation connue, ou tout au moins approchée à l'aide de modèles divers
(physiques ou mathématiques), il est possible de détecter les données qui ne cadrent pas
avec cette relation.
Les méthodes qui prévalent à ce niveau sont d'ordre hydrologique ou statistique, selon les cas
envisagés.
(8.15)
Avec :
Un contrôle possible des données est de mesurer tous les paramètres du bilan et vérifier
Un autre exemple de contrôle hydrologique simple est celui qui utilise les débits spécifiques
maximaux. On sait que ces débits spécifiques varient de manière inversement proportionnelle
à la surface sur laquelle ils sont calculés. Ainsi, si on dispose de plusieurs stations
hydrométriques sur un même réseau hydrographique, on peut tracer la courbe décroissante
qmax = f(A). Si un des points, correspondant à une station n'est pas « hydrologiquement » bien
placé il y a lieu de craindre une erreur (figure 8.8). Cependant il faut être attentif au fait que
cette erreur ne provient pas forcément d'une mauvaise mesure du débit, elle peut provenir
d'une erreur d'estimation de la superficie des bassins versants !
Fig. 8.8 - Méthode des débits spécifiques maximaux : visiblement une erreur est présente dans le débit
spécifique de la station 4. Il s'agit alors de chercher à découvrir les sources de cette erreur.
On peut estimer les données manquantes ou erronées à une station à partir des valeurs
provenant des stations voisines soumises aux mêmes conditions climatiques et situées dans
la même zone géographique. Trois méthodes sont proposées pour les données
pluviométriques :
(8.16)
Où :
Pour reconstituer des données de hauteurs d'eau, on peut utiliser des critères de
proportionnalité analogues dans le cas de stations limnimétriques placées sur un même
cours d'eau (transposition géographique, conservation des volumes,…).
Exemple :
♣ Soit Y une variable hydrométéorologique dont on connaît valeurs : y1, y2, …, yn.
♣ Soit X une autre variable observée de façon concomitante, avec les observations
x1, x2, …, xn.
♣ mais on possède par ailleurs (sur une période d'observation antérieure par
(8.17)
(8.18)
Notons que cette méthode néglige la dispersion naturelle des autour de la droite de
régression ; on diminue ainsi la dispersion des valeurs reconstituées et on fausse la
distribution statistique de l'échantillon complet. Il est nécessaire d'ajouter à la valeur
moyenne conditionnelle yi donnée par la formule précédente, un écart aléatoire simulé
d'écart-type estimé s e, obtenu à partir de la loi de probabilité de e.
Les relevés des débits d'une rivière pendant une longue série d'années montrent des
variations saisonnières systématiques (position des hautes et basses eaux) en fonction des
principaux facteurs influençant l'écoulement : le régime des précipitations, la nature du
bassin versant, sa situation géographique, l'infiltration, etc. Le régime hydrologique d'un
cours d'eau résume l'ensemble de ses caractéristiques hydrologiques et son mode de
variation. Il se définit par les variations de son débit habituellement représentées par le
graphique de l'écoulement mensuel moyen (calculé sur un certain nombre d'années et aussi
appelé débit "inter-mensuel" ou module mensuel). La figure 9.1 représente des valeurs de
modules mensuels moyens de certain cours d'eau dans le monde.
On utilise aussi le coefficient mensuel de débits, qui est défini comme le rapport du débit
mensuel moyen au module inter-annuel (moyenne inter-annuelle calculée sur un certain
nombre d'années). Celui-ci permet de représenter la répartition, en pourcentage, des débits
mensuels au cours de l'année.
(9.1)
(9.2)
De même, les courbes des fréquences relatives des débits sur une longue série d'années,
définissent la variation saisonnière des quantiles de débits (Fig. 9.2). Les courbes cotées 10,
25,..., 90 % indiquent les valeurs des débits mensuels qui ont respectivement 10, 25,..., 90
chances sur 100 de ne pas être atteintes ou dépassées.
Fig. 9.2 - Exemple de courbes fréquentielles (fréquence de non-dépassement) des débits mensuels
.
9.2 Classification des régimes hydrologiques
Une des classifications des régimes hydrologiques des rivières les plus simples est celle de
Pardé (1933), qui distingue trois types de régimes :
Cette classification peut être éventuellement rectifiée en fonction des causes hydrologiques
provoquant les hautes eaux ; c'est le cas pour les phénomènes d'embâcle et de débâcle.
L'embâcle désigne une accumulation, due à un obstacle (présence d'un pont, d'un barrage,
d'un rétrécissement, d'un coude, etc.), de glaçons ou de bois dans un cours d'eau qui crée un
barrage. Lorsque ce barrage cède pour différentes raisons c'est la débâcle. Lorsqu'il s'agit
d'un démantèlement de couche de glace sur les cours d'eau, cela traduit l'effet du dégel. La
débâcle produit alors un charriage de glaçons de tailles plus ou moins grosses, pouvant à
leur tour occasionner, lorsqu'ils sont arrêtés par un autre accident hydrographique ou
autres, des barrages provisoires (embâcle) qui provoquent souvent des inondations.
Par ailleurs, la géologie peut modifier sensiblement les écoulements et par delà le régime
d'alimentation des cours d'eau. Ceci est particulièrement vrai dans les régions karstiques (ex.
dans le Jura).
Il est caractérisé par un seul maximum et un seul minimum annuel du coefficient mensuel
des débits et traduit la prépondérance d'un seul mode d'alimentation (régime glaciaire, nival
ou pluvial). Ce caractère peut cependant cacher la combinaison de plusieurs influences et
confère ainsi aux régimes des rivières concernées une simplicité apparente.
Le régime glaciaire se retrouve en général quand 15 à 20% du bassin est occupé par des
glaciers. Sous nos climats, le régime glaciaire se caractérise entre autres par :
♣ Amplitude des variations mensuelles des débits très grande (rapport entre les
coefficients mensuels extrêmes), due au rapport crue/étiage très élevé.
♣ Grande régularité d'une année à l'autre du régime car la température est de tous
les paramètres météorologiques le moins irrégulier.
Suivant l'altitude moyenne des bassins versants, ces caractéristiques seront plus ou moins
prononcées. Par exemple l'amplitude des variations mensuelles de débits est supérieure à 25
pour les bassins versants de haute altitude (altitude moyenne supérieure à 2500 mètres), et
varie de 12 à 35 pour les bassins de 2300 à 2600 mètres d'altitude moyenne. Le Rhône en
amont du lac Léman est caractérisé par un régime de type glaciaire (Fig. 9.1 et Fig. 9.3).
Le régime nival pur présente sous une forme atténuée certaines des caractéristiques du
régime glaciaire. Le maximum a lieu cependant plus tôt (juin). Il se subdivise en régime nival
de montagne et nival de plaine.
Le régime nival de montagne, se retrouve dans les zones montagneuses où la majorité des
précipitations arrive sous forme de neige. Il est caractérisé par :
♣ une fonte progressive de la neige, qui commence d'abord aux altitudes les plus
basses et provoque une crue en mai-juin (pour l'hémisphère Nord)
La rivière Fraser à Hope aux Canada (Fig. 9.1) est caractérisée par ce régime.
Le régime nival de plaine intéresse les régions continentales et maritimes à faible altitude du
nord de l'Europe. Ses caractéristiques sont les suivantes :
♣ Très grande variabilité au cours de l'année, due à des basses eaux d'été très
marquées (températures élevées et forte ET0)
♣ Ecoulement important.
Les fleuves sibériens, comme la Lena (Fig. 9.1) ont un régime nival de plaine.
On peut aussi distinguer le régime nival de transition que l'on rencontre sur les bassins
versants d'altitude moyenne comprise entre 1200 et 1600 mètres. Il se rapproche davantage
d'un type complexe dans ce sens qu'il présente quatre saisons hydrologiques. Ses
caractéristiques sont les suivantes :
♣ Courbe des coefficients mensuels des débits montrant deux maxima (fort en mai-
juin, et plus modéré en novembre-décembre) et deux minima.
Bien que le régime pluvial appartienne aux régimes simples, il présente des caractéristiques
différentes de celles des régimes précédents. Il se distingue par :
♣ Des Hautes eaux (avec un maximum plus ou moins marqué) en hiver et des
basses eaux en été. Bien qu'il soit fréquent que les pluies de la saison de basses
eaux soient égales ou supérieures à celles de la saison des hautes eaux, les
températures étant élevées, l'évaporation est importante).
♣ Ecoulement généralement assez faible (exemple la. Seine : 6 l/s/km2, cf. Fig 9.1)).
C'est le régime des cours d'eau de faible à moyenne altitude (500 - 1000 mètres). Il se
retrouve dans les régions tempérées sans neige.
On distingue le régime pluvial tropical dont l'allure des courbes de variation de Cm ressemble
au régime glaciaire. Il présente les caractéristiques suivantes :
♣ Une grande variabilité des débits au cours de l'année avec des minima pouvant
atteindre des valeurs très faibles. Par exemple à Koulikoro (fleuve Niger), le débit
instantané peut excéder 8000 m3/s en septembre mais rester inférieur à 100
m3/s à la fin du printemps (voir aussi Fig. 9.1).
♣ Une relative régularité d'une année à l'autre ; on observe cependant des années
marquées par un net déficit des pluies (cas des années 1971 et 1973) en région
subsaharienne.
Le fleuve Sénégal (Fig. 9.1) et le Niger amont (Fig. 9.1 et Fig. 9.4) sont des rivières
caractéristiques du régime tropical.
Il se caractérise par deux maxima et deux minima des coefficients mensuels au cours de
l'année hydrologique. Suivant les modes d'alimentation principaux on distingue le régime
nivo-glaciaire, glacio-nival, nivo-pluvial, pluvio-nival.
♣ Grandes variations d'une année à l'autre, mais cependant moindres que pour le
régime nival.
♣ Ecoulement important.
Il se caractérise par :
♣ Deux maxima nets, l'un assez prononcé vers avril-mai à la fonte des neiges, et
l'autre en automne (vers novembre) plus modéré. Ce second maximum, dépendant
des pluies tombées en automne, peut être faible (de coefficient inférieur à 1).
♣ L'amplitude (rapport entre les coefficients mensuels extrêmes) est comprise entre
2 et 5.
Fig. 9.5 - Exemple de courbe des coefficients mensuels de débit pour le régime nivo-pluvial préalpin
(Emme à Emmenmat)
[Link] Le régime pluvio-nival
La tendance pluviale est d'autant plus marquée que le bassin se situe à basse altitude (650 à
750 mètres). Le régime pluvio-nival est caractérisé par :
♣
Fig. 9.6 - Coefficient mensuel de débits pour le régime pluvio-nival (Orbe à Orbe).
La rivière Po en Italie présente un régime pluvio-nival (Fig. 9.1). En Suisse, l'Orbe à Orbe est
caractérisé par un régime de ce type même si le maximum pluvial est du même ordre de
grandeur que le maximum hivernal (Fig. 9.6).
Le régime complexe est généralement rencontré sur les grands fleuves, dont les affluents,
d'amont en aval, influencent de façon très diverse l'écoulement général. Le régime des grands
fleuves se présente comme une synthèse de ceux de leurs sous-bassins constitutifs, le plus
souvent très variés du point de vue altitude, climat, etc. Habituellement, ces influences
diverses tendent à atténuer les débits extrêmes et à accroître la régularité annuelle des débits
moyens mensuels, de l'amont vers l'aval (voir Fig. 9.1 pour le Rhin à Rees et le Rhône à l'aval
du lac Léman à Chancy).
♣ Pour le régime du sud des Alpes les régimes : nival méridional, pluvio-nival
méridional, nivo-pluvial méridional, pluvial méridional.
Fig. 9.7 - Variations des coefficients mensuels de débits pour les 16 types de régimes déterminés en Suisse suite à
l'étude de 95 bassins représentatifs (d'après les données de l'"Atlas hydrologique de la Suisse")
L'objectif de ce chapitre est de présenter les aspects à la fois généraux et modernes des
théories relatives aux mécanismes de génération de l'écoulement depuis les idées de Horton
(1933) - novatrices à l'époque -, jusqu'à celles de ce jour, basées sur l'idée de cheminements
préférentiels. On se gardera toutefois d'affirmer que la recherche en ce domaine est achevée
car le fonctionnement hydrologique, malgré l'importance des apports de ces dernières
années, reste encore bien mal compris. Le développement de nouvelles théories a même
donné lieu à toute une série de problèmes et soulevé de nouveaux questionnements.
Les notions abordées dans les chapitres précédents sont utilisées pour décrire plus
précisément les principaux processus de génération des écoulements. Nous ne détaillerons
pas les protocoles expérimentaux qui doivent être mis en place afin de les identifier. Seul le
traçage environnemental sera brièvement présenté à cette fin.
Les développements relatifs à l'étude des processus étant relativement récents, il règne
encore une certaine confusion concernant la terminologie utilisée pour leur dénomination. De
ce fait, lorsque nous distinguerons les différents processus, nous indiquerons les autres
terminologies connues.
Globalement, les processus régissant la génération des écoulements sont encore mal connus.
Ceci est dû en partie au fait que les réponses aux deux questions essentielles de l'hydrologie
sont multiples. On peut finalement tenter de répondre - du moins partiellement - aux deux
questions « quel est le devenir de l'eau des pluies ? » (Penmann, 1963) et « quelle est la
provenance de l'eau des rivières ? » (Hewlett, 1961), en distinguant dans un premier temps
quatre cheminements principaux de l'eau à la rivière (Ward et Robinson, 1990) :
Les figures 10.1a et 10.1b suivantes illustrent schématiquement ces deux processus
d'écoulement sur un versant et les différents types d'écoulements cités ci-dessus :
(a) Ecoulement par dépassement de la capacité d'infiltration Ra.
(b) Ecoulement par saturation Rb. (b') Ecoulement par saturation Rb.
Fig. 10.1 - Ecoulement par dépassement de la capacité d'infiltration (a) et écoulement par saturation (b).
A cet ensemble de processus peut encore s'ajouter un cinquième élément qui est le
ruissellement dû à la fonte des neiges ("snowmelt runoff") que nous avons traité au chapitre 6
"le stockage et ses variations".
Finalement les éléments les plus importants dans la génération des crues sont les
écoulements de surface et de subsurface.
La classification adoptée ici n'est pas la seule existante. Par exemple, Dunne (1978) propose
de classer les écoulements en cinq processus essentiels qui sont respectivement (Fig. 10.2) :
Fig. 10.2 - Chemins de l'écoulement de l'eau sur le versant (D'après Dunne, 1978).
Parmi les processus de génération de crue, le premier qui peut être souligné est bien
évidemment la contribution directe de la précipitation sur la surface des cours d'eau du
bassin versant. Cependant, ce type de processus est considéré comme marginal du fait que la
surface occupée par les cours d'eau de type pérenne ne représente qu'une très faible fraction
de la surface totale du versant. L'importance des précipitations directes croît toutefois suite à
des précipitations de longues durées ou avec le développement du réseau hydrologique qui
présente alors d'importantes zones lacustres ou marécageuses.
Fig. 10.3 - Taux d'infiltration et infiltration cumulée pour une pluie uniforme, définition du temps de submersion.
L'écoulement sur surfaces saturées se produit lorsque la capacité du sol à stocker l'eau est
épuisée et lorsque la capacité à transmettre latéralement le flux d'eau est dépassée1. Par
conséquent, l'eau ne pourra plus s'infiltrer et va s'écouler en surface.
La figure 10.4 résume les deux situations rencontrées jusqu'ici, à savoir l'écoulement sur
surfaces saturées ainsi que l'écoulement par dépassement de la capacité d'infiltration.
Fig. 10.4 - Processus de génération d'écoulement par dépassement de la capacité d'infiltration ainsi
que sur surfaces saturées.
1L'écoulement par percolation profonde est supposé négligeable.
La figure 10.5 montre un bassin versant et son réseau hydrographique ainsi que l'extension
successive des surfaces saturées au cours d'un événement de pluvieux (concept de "variable
source area").
Fig. 10.5 - Extension du réseau hydrographique et des surfaces saturées durant une crue pour 4
pas de temps (D'après Chorley ).
La condition essentielle pour voir apparaître des écoulements de subsurface est que la
conductivité hydraulique latérale du milieu doit être nettement supérieure à la conductivité
verticale. Dans ce cas l'eau, s'écoule latéralement en zone non-saturée, par un mécanisme
type « toit de chaume 2» ou dans le cas d'un écoulement en zone saturée, lorsqu'une nappe
perchée se forme dans la partie supérieure d'un horizon peu perméable.
2L'analogie avec un toit de chaume se comprend aisément. Le toit de chaume, par l'organisation
structurée de ses fibres permet un écoulement latéral tandis que le mouvement vertical de l'eau est
difficile. Dans le cas du sol, une conductivité hydraulique latérale élevée couplée à une conductivité
verticale faible permet ce type de transport de l'eau.
Parmi les analyses du processus des écoulements de subsurface, visant à expliquer la forte
proportion d'eau ancienne dans l'hydrogramme de crue, certains auteurs ont supposé
l'existence d'un mécanisme de transmission quasi instantanée d'une onde de pression. Ce
mécanisme, nommé "effet piston", suppose qu'une impulsion d'eau reçue par le versant est
transmise à l'aide d'une onde de pression vers l'aval, provoquant une exfiltration 3 immédiate
en bas de versant. Ce phénomène peut s'expliquer par analogie avec une colonne de sol
saturée sur laquelle on applique une charge d'eau. L'eau se déplace sous l'effet de la gravité
en chassant celle qui se trouve à l'autre extrémité de la colonne. Dans ce type d'explication, il
faut alors distinguer la vitesse "fictive" de l'eau dans le sol qui est relativement lente et qui
détermine le temps moyen de transit de l'eau sur le versant de la vitesse de propagation de
l'onde de pression qui peut exprimer la vitesse de réaction du bassin versant.
Malgré la simplicité de cette explication, "l'effet piston" est limité par le fait qu'une impulsion
d'une certaine quantité d'eau ne s'accompagne d'une exfiltration équivalente (ou presque) que
dans les cas où le sol présente une capacité de stockage très faible.
3L'exfiltration se produit lorsque l'apport par l'écoulement de surbsurface (ou même souterrain)
dépasse la capacité du sol à transmettre un flux d'eau transversal (rencontre avec un milieu peu
perméable par exemple). Dans ce cas l'écoulement de subsurface émerge à la surface du sol.
L'écoulement par macroporosité a vu son intérêt grandir auprès des hydrologues depuis une
quinzaine d'années. On lui attribue principalement le rôle d'accélérer la recharge de la nappe
4 Cette partie consacrée aux macropores s'inspire largement d'un article de Beven et Germann (1982).
La définition du macropore dépend d'une part d'arbitraire quant au choix d'une taille
effective ainsi que de l'expérience pratique que l'on a sur le terrain. D'autre part, l'analogie du
comportement de l'eau dans le sol avec celui d'un ensemble de capillaires devient discutable
dès lors que la porosité augmente. Ceci implique qu'il peut être pertinent de classer les types
de pores selon la conductivité hydraulique et non plus uniquement selon des critères
dimensionnels (par exemple le diamètre équivalent), bon nombre d'études ayant montré à ce
jour que le critère de la taille n'était pas suffisant pour effectuer un classement des types de
porosité du sol.
Pour terminer, on peut encore mentionner une définition relativement facile à mémoriser : un
macropore est un pore où les phénomènes de capillarité sont inexistants.
♣ Pores formés par la végétation - Ils sont formés lors de la pénétration des
racines dans le sol et deviennent libres pour le passage de l'eau lorsque la plante
meurt. La structure du réseau des macropores va dépendre du type de plante
ainsi que du stade de développement végétatif.
La présence de macropores connectés avec la surface du sol peut jouer un rôle important
dans le processus d'infiltration. Le comportement hydraulique met alors en défaut la loi de
Darcy car le sol ne présente plus les propriétés hydrauliques requises (Musy et Soutter,
1991). Les flux hydrauliques peuvent alors varier de plusieurs ordres de grandeur pour des
distances de quelques centimètres seulement. Il s'ensuit que l'idée d'une conductivité
hydraulique moyenne ayant un sens physique est vouée à l'échec. Des études ont mis en
évidence que le champ des pressions se développant durant l'infiltration était très irrégulier
et que le comportement de type darcien n'était guère représentatif du phénomène physique
sous-jacent. Celui-ci peut être décrit ainsi :
Soit une précipitation P(t) arrivant à la surface d'un sol présentant des macropores. Dans un
premier temps, l'eau va s'infiltrer dans la matrice du sol au taux I1(t) puis, lorsque la capacité
d'infiltration est dépassée par l'intensité de la pluie, l'infiltration s'effectue dans les
macropores au taux IM(t). Toutefois, le flux d'eau dans les macropores va être transféré dans
des pores de tailles plus faibles sous l'effet des forces de capillarité. Ceci s'effectue au régime
I2(t) et constitue un second type d'infiltration de l'eau dans la matrice. Finalement, une fois
que le sol est saturé, l'écoulement à travers les macropores EM(t) peut avoir lieu, tout comme
l'écoulement de surface ES(t). La Figure 10.6 ci-dessous illustre cette schématisation.
Un certain nombre d'observations montre donc qu'il est nécessaire de formuler de manière
nouvelle les processus d'infiltration et de redistribution dans le cas de sols présentant une
macroporosité. Pourtant, le rejet des idées de Darcy n'est pas pertinent. On conseille plutôt
de proposer un couplage entre les processus d'infiltration dans la matrice poreuse ou la loi de
Darcy est justifiée avec le comportement de l'eau dans les macropores. Cependant il convient
de préciser certains éléments :
♣ La nature des écoulements dans le domaine des macropores - On n'a que peu de
connaissances du transport de l'eau dans les macropores si ce n'est que ces
derniers peuvent véhiculer beaucoup d'eau sans pour autant être saturés.
On admet en règle générale que l'hydrogramme de crue est souvent contrôlé par les
écoulements de subsurface. Si ces écoulements transitent principalement par la présence de
macropores, la transmission de l'eau peut se faire à des vitesses du même ordre de grandeur
que celles obtenues pour l'écoulement de surface. Cependant, et comme nous l'avons
remarqué, les macropores ne sont pas les seules et uniques causes d'écoulement de
subsurface. Finalement et dans le cas où l'on considérerait uniquement la contribution due à la
macroporosité, la loi de Darcy ne peut pas être appliquée principalement dans deux cas de figure :
♣ Flux quasi turbulent ou turbulent dans les zones non saturées ou saturées.
De la même façon que les écoulements par macropores, les écoulements par tubes (pipeflow)
peuvent accélérer le drainage du sol et l'écoulement de l'eau. Il est cependant difficile d'établir
une distinction exacte entre les macropores et ces tubes. On considère généralement que les
tubes sont de très gros macropores. En plus de ce critère strictement géométrique, les tubes
présentent un degré de connectivité plus élevé que celui offert par un réseau de macropores
sans que l'on puisse pour autant affirmer que les tubes forment un réseau continu. Les tubes
conduisent le plus souvent l'eau dans un milieu non saturé et permettent la connexion de
sources éloignées d'écoulement avec le cours d'eau.
Si la nappe ou la frange capillaire est proche de la surface du sol, une petite quantité d'eau
suffit à saturer le profil. Si la capacité du sol à transmettre l'écoulement de subsurface
diminue (ceci peut se produire par exemple si le type de sol change), ce dernier revient en
surface et ruisselle. On retrouve alors dans l'écoulement de surface de l'eau antérieure à l'eau
de pluie. Cette contribution d'eau "ancienne" va ainsi augmenter le volume de l'écoulement
de surface. En analysant la répartition spatiale des zones où se produisent ce type de
phénomène, on met en évidence que celles-ci peuvent subir une extension rapide là où la
frange capillaire est proche de la surface du sol. La répartition spatiale des zones propices
aux écoulements de retour est donc liée à la topographie du bassin versant.
On doit toutefois insister sur le fait que ces surfaces propices à l'écoulement de retour
deviennent aussi favorables au développement du processus d'écoulement sur surfaces
saturées. L'hydrogramme résultant de ces deux types de processus sera donc composé d'eau
ancienne et d'eau nouvelle. Une fois encore, on peut noter que le processus de formation
d'une crue est la résultante de plusieurs phénomènes concomitants (Fig. 10.8)
L'eau souterraine joue un rôle aussi important dans la génération des écoulements et
notamment dans la composante "débit de base" de l'hydrogramme. Dans cette section, nous
considérons essentiellement les écoulements profonds dont l'étude complète appartient à
l'hydrogéologie (cf. chapitre 6) et présentons une méthode d'investigation expérimentale
devenue désormais classique : l'utilisation de traceurs pour déterminer notamment la
contribution de l'écoulement souterrain à l'écoulement total dans un cours d'eau.
En terme de processus, l'eau qui rejoint la nappe peut être considérée comme de l'eau
souterraine. Une partie de cette eau, après percolation, va transiter à travers l'aquifère avec
une vitesse de quelques mètres par jour à quelques millimètres par an avant de rejoindre le
cours d'eau souvent par le biais d'un phénomène de résurgence de la nappe. L'écoulement de
base assure ainsi le débit des rivières en l'absence de précipitations et soutient les débits d'étiage.
Fig. 10.9 - Deux situations distinctes où la nappe peut contribuer au débit de la rivière (a) et où la
nappe peut drainer le cours d'eau (b).
Les traceurs environnementaux sont apparus durant les années soixante et ont rapidement
vu leurs applications se diversifier. Parmi celles-ci, on peut noter l'analyse de la recharge des
nappes d'eau souterraines, leur datation ainsi que l'estimation de leur contribution au débit
des rivières et leurs interactions avec l'eau provenant d'autres sources. Ces méthodes ont
aussi donné la possibilité d'effectuer des études qualitatives de contamination des nappes
souterraines ainsi que de leur salinisation. Parmi les traceurs environnementaux utilisés, les
traceurs isotopiques ont été privilégiés. Les isotopes naturels les plus utilisés sont les
isotopes stables : hydrogène lourd ou deutérium (2H) et l'oxygène 18 (18O). Ces deux traceurs
sont des constituants de l'eau naturelle et sont de fait de bons éléments pour l'analyse du
cheminement de l'eau.
Principalement, les méthodes de traçages isotopiques se basent sur le fait que la composition
isotopique de l'eau contenue dans les sols est différente de celle de l'eau de pluie et de celle
de l'eau des rivières. On considère alors que l'eau contenue dans les sols est de l'eau
"ancienne" tandis que l'eau de pluie est dénommée eau "nouvelle". Il est alors possible de
déterminer les contributions d'eau ancienne et d'eau nouvelle d'un débit mesuré Q à l'aide
d'un système de deux équations à deux inconnues :
(10.1)
La première équation exprime simplement que le débit total du cours d'eau est la somme
du débit d'eau ancienne et du débit d'eau nouvelle tandis que la seconde équation
exprime le fait que le produit de la concentration en isotope ( ) avec le débit total est la
somme des produits des concentrations et des débits pour les deux sources d'eau ancienne
( ) et nouvelle ( ). En mesurant les différentes concentrations ainsi que le débit total, il est
alors possible de déterminer les débits dus à l'eau ancienne et à l'eau nouvelle :
(10.2)
(10.3)
D'un point de vue opérationnel, on détermine les teneurs en isotopes de la pluie ( ), de l'eau
encore le débit de la rivière . Ainsi, les seules inconnues sont bien les contributions de
l'eau ancienne et de l'eau nouvelle.
♣ La pluie est caractérisée par un signal unique. D'une façon générale, les
variations isotopiques de la pluie sont essentiellement temporelles. La teneur
isotopique varie fortement en fonction de la durée et du type de précipitations.
La mesure des teneurs en isotopes peut se faire à l'aide d'un spectromètre de masse. On
détermine ainsi les rapports 18O/16O et 2H/1H. On étudie ensuite la concentration en isotopes
stables par rapport à une référence qui est la teneur moyenne de l'isotope considéré dans les
eaux océaniques nommée SMOW (Standard Mean Ocean Water). On calcule ainsi le « delta »,
pour un élément donné, entre le rapport isotopique d'un échantillon et celui d'un standard
(SMOW), exprimé en pour mille (exemple pour 18O) :
(10.4)
Les valeurs négatives de concentration signifient que l'échantillon est plus pauvre que la
valeur moyenne des eaux océaniques tandis que les valeurs positives entraînent le constat
inverse. On utilise enfin très souvent la relation linéaire qui existe entre la concentration en
Oxygène 18 et la concentration en deutérium (Fig. 10.10). Pour les pluies, cette relation porte
le nom de droite des eaux météoriques et s'exprime comme suit pour les eaux météoriques
mondiales :
(10.5)
D'une façon générale, on constate que la valeur de la pente de cette droite est assez
constante tandis que l'ordonnée à l'origine qui marque l'excès de Deutérium peut dépasser la
valeur de 10. Ce dépassement se produit lorsque la vapeur d'eau d'origine océanique a été
enrichie de manière significative par une évaporation au niveau des continents ou de mers
fermées. Dans le cas du bassin de la Méditerranée par exemple, l'équation (10.5) s'écrit :
(10.6)
L'équation globale (10.5) peut aussi s'établir de manière locale (pluies locales)
On peut encore noter que l'analyse de la relation entre oxygène 18 et deutérium permet de
mettre en évidence des eaux qui ont subi un processus d'évaporation (par exemple pour un
cours d'eau et à différentes phases d'une crue). Le suivi des différentes valeurs de 18O et 2H
au cours du temps permet de tracer la droite d'évaporation qui présente une pente plus faible
que la droite des eaux météoriques ainsi qu'une valeur plus faible de l'ordonnée à l'origine.
L'intersection de cette droite avec la droite des eaux météoriques donne ainsi la possibilité de
déterminer la composition isotopique de l'eau avant son évaporation.
Fig. 10.10 - Relation entre Oxygène 18 et Deutérium dans les eaux naturelles (D'après Fontes)
L'exemple ci-dessous montre la séparation entre eau nouvelle et eau ancienne à l'aide de
l'Oxygène 18. L'événement analysé est une crue ayant eu lieu sur le bassin versant de Bois-
Vuacoz (24 ha) entre le 7 et le 8 septembre 1993 (les mesures pluviométriques sont effectuées
au pas de temps horaire). L'épisode pluvieux est considéré comme étant de volume important
(45 mm) mais d'intensité modérée (intensité maximale = 9 mm/h). La figure 10.11 ci-dessous
montre successivement le résultat de la décomposition isotopique en terme de débit,
l'évolution des teneurs en Oxygène 18 de l'eau de la pluie ainsi que de l'eau de la rivière et,
en dernier lieu, l'évolution relative de la quantité d'eau ancienne, eau nouvelle dans le débit
mesuré. On peut noter l'importance de l'eau ancienne au début de l'événement ainsi que
l'augmentation rapide de la fraction d'eau nouvelle dans le débit de la rivière.
Fig. 10.11 - Séparation isotopique des écoulements et variation des teneurs isotopiques pour
l'événement du 7-8 septembre 1993 (D'après Iorgulescu, 1997).
Pour déterminer le cheminement de l'eau, il est aussi possible de recourir à d'autres types de
traceurs naturels ou artificiels, chimiques et isotopiques. Les traceurs artificiels (rhodamine,
fluoréscéine) sont bien adaptés à des études locales sur de petites surfaces mais présentent
différents inconvénients à l'échelle du bassin versant. Leur répartition spatiale et temporelle
ainsi que leur coût et la difficulté d'obtenir un état d'équilibre avec leur concentration
limitent fortement leur usage à l'échelle du versant. Les traceurs chimiques ont été
relativement utilisés avant la généralisation des méthodes isotopiques. Les analyses portent
sur les différents éléments chimiques majeurs à savoir :
Les traceurs chimiques ont toutefois le très grand inconvénient de ne pas être conservatifs.
Ceci implique par exemple qu'une pluie va s'enrichir par pluvio-lessivage avant même qu'elle
n'atteigne le sol. Les isotopes naturels ne présentent pas ces inconvénients. On recourt ainsi
le plus souvent aux molécules contenant de l'oxygène 18 ainsi que celle contenant du
deutérium (cf. ci-dessus).
Par analogie avec ce que l'on a présenté au sujet de l'oxygène 18 et du deutérium, on peut
développer un modèle d'analyse de l'origine et du cheminement de l'eau en recourant à
plusieurs traceurs. En reprenant les équations du système (10.1), ce modèle de mélange, à j
traceurs et i composantes, peut s'exprimer comme suit :
(10.7)
Où :
Qt : débit total de la rivière (somme des débits due à chaque composante Qi),
dtj : concentration du j-ème traceur dans le débit total,
dij concentration du j-ème traceur dans la i-ème composante.
Dans le but de déterminer le comportement hydrologique des bassins versants, une méthode
originale (EMMA) a été proposée pour déterminer la contribution des différentes sources
d'écoulement. Cette méthode reposait initialement sur le fait que les sols dans les régions à
climat tempéré présentent une signature chimique verticale et/ou horizontale. La chimie de
la rivière est alors supposée être un mélange d'eau de subsurface de profondeur variable (ex :
eau de la nappe, eau du sol). Les composantes sont définies selon leur aptitude à expliquer la
variabilité des compositions chimiques observées dans la rivière. Ainsi le domaine défini dans
l'espace des traceurs par les concentrations des composantes doit contenir l'ensemble des
concentrations chimiques de l'eau de rivière.
L'évidence de la diversité des processus de génération de crues s'est posée tout au long de ce
chapitre. Ces mécanismes n'agissent pas de manière isolée mais forment un domaine continu
de processus. Il s'ensuit que sur un même bassin versant, plusieurs processus concomitants
peuvent intervenir durant un même événement pluvieux. De la même manière, il se peut que
le type de processus change selon le type d'événement. Les crues générées durant l'été ou
l'hiver ne font pas nécessairement appel au même type de processus. On comprend ainsi qu'il
est difficile de caractériser l'ensemble des processus de génération de crue. Dunne (1978) a
proposé une classification des processus dominant en fonction de trois critères qui sont
respectivement le temps de montée de la crue (peak lag time), le débit spécifique maximum
(peak runoff rate) et la surface du bassin versant. Les figures 10.13 et 10.14 présentent ces
classifications. On y distingue l'écoulement hortonien (Horton overland flow), l'écoulement
par saturation (saturation overland flow), l'écoulement de subsurface (Throughflow) ainsi que
deux types d'écoulements par tubes.
Les figures 10.13 et 10.14 permettent ainsi de noter que si l'on s'intéresse spécifiquement à
des problèmes liés aux crues - c'est à dire pour des valeurs de débits importantes, il
conviendra d'être attentif au fait qu'il s'agit essentiellement de processus de type hortonien et
ceci indépendamment de la taille du bassin versant. De plus, les crues sont caractérisées par
des temps de montée relativement faibles qui augmentent avec la taille du bassin versant. A
l'inverse, des crues de faible volume dont le temps de montée peut être de l'ordre de une à
plusieurs heures sont représentatives de processus d'écoulement de subsurface. On
soulignera encore que les figures présentées ci-dessous, bien qu'elles n'aient pas la
prétention de se révéler l'unique outil d'identification des processus, permettent néanmoins,
sur la base d'observations élémentaires, de fixer une gamme possible de processus de
génération des écoulements pour les situations étudiées.
Fig. 10.13 - Courbes enveloppes du temps de montée en fonction de la surface du bassin versant
(D'après Jones)
Fig. 10.14 - Courbes enveloppes du débit spécifique maximum en fonction de la surface du bassin
versant (D'après Jones)
Enfin, la figure 10.15 permet de replacer les différents types d'écoulement dans leur contexte
géomorphologique en faisant intervenir les facteurs sols, climat et topographiques.
La manière dont réagit le bassin lorsqu'il est soumis à une sollicitation se nomme réponse
hydrologique. Elle est schématiquement représentée à la figure 11.1 ci-dessous.
En fait, une averse tombant sur un bassin versant aura pour conséquence, en une station de
contrôle située sur le cours d'eau, une réponse pouvant être nulle (absence de modification
de l'écoulement ou absence de crue) ou positive (écoulement modifié ou crue). En fonction de
ce qui a été vu précédemment, cette réponse peut être :
♣ Rapide- La réponse rapide est imputable aux écoulements de surface ou, par
exemple, à un effet piston, ou encore à l'effet de la macroporosité du sol.
De plus, la réponse peut être différenciée selon que cette dernière est :
♣ Totale- Dans ce cas, la réponse hydrologique est composée à la fois par des
écoulements de surface et souterrains.
♣ Partielle - c'est à dire lorsque la réponse est la résultante d'un ou l'autre des
processus décrit précédemment
Le rôle de l'ingénieur est donc d'une part d'identifier les processus hydrologiques et leur part
respective intervenant dans la réponse du bassin versant et, d'autre part, les modalités du
passage de l'impulsion pluviométrique à la réponse hydrologique. La question qui se pose
alors est de comprendre et interpréter les mécanismes de transformation de la pluie à
l'hydrogramme de crue.
Dans le cadre d'une description des processus selon le principe établi par Horton, la
transformation de la pluie en hydrogramme de crue se traduit par l'application successive de
deux fonctions, nommées respectivement fonction de production – ou fonction d'infiltration
- et fonction de transfert. La fonction de production permet de déterminer le hyétogramme
de pluie nette à partir de la pluie brute. La fonction de transfert permet quant à elle de
déterminer l'hydrogramme de crue résultant de la pluie nette (la pluie nette est la fraction de
pluie brute participant totalement à l'écoulement). (Fig. 11.2).
Une averse, définie dans le temps et dans l'espace, tombant sur un bassin versant de
caractéristiques connues, et dans des conditions initiales données, provoque à l'exutoire du
bassin considéré un hydrogramme défini. La figure 11.3 définit quelques éléments essentiels
relatifs à l'hydrogramme résultant d'un hyétogramme spécifique.
L'hydrogramme de crue présente la forme générale d'une courbe en cloche dissymétrique que
l'on divise en quatre parties : tarissement (avant la pluie nette), crue, décrue et tarissement
(après la recension hydropluviométrique étudiée). On définit alors des temps
caractéristiques :
1Il s'agit là d'un constat identique à celui que l'on a fait lorsque l'on a présenté les notions de
coefficient d'écoulement de surface et coefficient de ruissellement au chapitre consacré à l'étude de
l'écoulement et de l'infiltration.
La nature et l'origine des crues ou hautes eaux sont liées aux régimes hydrologiques et à la
taille du bassin versant. Les bassins versants du Plateau suisse, par exemple, appartiennent
au régime pluvial ou au régime nivo-pluvial. Les crues auront ainsi pour origine les averses
(liquides et /ou solides) et/ou la fonte de neige. Les crues peuvent être groupées, selon les
causes qui les engendrent selon :
♣ les crues d'embâcle2 de glace (lorsque des blocs de glace d'un cours d'eau gelé
emportés lors du dégel s'accumulent et forment des barrages à l'amont desquels
les plaines s'inondent). La débâcle résulte de la brusque rupture de ces barrages,
provoquant ainsi des crues violentes mais brèves.
2Les termes d'embâcle et de débâcle sont aussi parfois utilisés pour désigner une accumulation
puis une rupture de matériaux autre que la glace (exemple : tronc d'arbres...).
La réponse hydrologique d'un bassin versant est influencée par une multitude de facteurs
tels que ceux liés :
Les facteurs liés aux précipitations ainsi qu'aux conditions climatiques sont des facteurs
externes au milieu tandis que la morphologie, les propriétés physiques du versant, la
structuration du réseau et les conditions antécédentes d'humidité sont des facteurs internes.
L'ensemble de ces éléments ayant été étudié tout au long des chapitres précédent, nous n'y
reviendrons pas ici à l'exception du rôle de la pluviosité qui sera développé dans les
prochains paragraphes ainsi que l'importance de l'état antécédent d'humidité.
Dans ce cadre, l'Institut d'Aménagement des Terres et des Eaux (IATE / HYDRAM) a
développé un didacticiel permettant en outre d'étudier la réponse hydrologique d'un bassin
versant lorsque celui-ci est soumis à une sollicitation pluvieuse dans quelques situations
simples. Cet outil, nommé DHYDRO est une application informatique pouvant être acquise
ici. La figure 11.4 ci-dessous donne un aperçu sommaire des fonctionnalités de cette
application. On y a représenté quatre fenêtres nommées A, B, C et D. La fenêtre A illustre le
contrôle essentiel du bassin versant, à savoir sa couverture du sol et la pente du versant
ainsi que le type de réseau hydrographique qu'il comporte, la fenêtre B donne la description
topographique du versant ainsi que celle du réseau hydrographique. La fenêtre C permet de
déterminer la pluie nette à partir d'un fichier de pluie brute selon différentes fonctions
d'infiltration que l'utilisateur peut choisir. Enfin, la dernière fenêtre (D) représente
graphiquement l'hydrogramme résultant d'une précipitation donnée sur un bassin versant
défini.
Pour illustrer l'influence de ces facteurs, on se base ci-dessous sur une schématisation très
ancienne (Linsley et Crowford, 1966, reprise par Réménérias, 1976) qui conserve malgré son
âge toute son acuité pédagogique.
Pour un événement pluvieux donné, la réponse hydrologique d'un bassin dépend du volume
précipité, mais également des variations d'intensité et de la durée de l'averse. Afin d'illustrer
tout d'abord l'influence de la durée de l'averse, nous emprunterons le raisonnement suivant :
Soit un bassin versant divisé en quatre zones A-B-C-D concentriques d'égales surfaces,
d'égales coefficients de ruissellement et délimitées par des lignes isochrones. Rappelons que
les lignes isochrones sont des lignes d'isovaleurs de temps de concentration entre le point
considéré et l'exutoire du versant. On assigne à chaque secteur un temps d'écoulement
variant de 1 heure pour la zone la plus proche de l'exutoire (A) à 4 heures pour la zone la
plus éloignée (D) (Fig. 11.5).
On considère un événement pluvieux de volume égal à 10 mm constant pour tous les cas de
figure examinés. La figure 11.6 présente trois cas distincts de durée :
Fig. 11.6 - Influence de la durée de l'averse sur la réponse hydrologique d'un bassin versant (d'après
Réménérias. Les hyétogrammes figurent sur la partie gauche de la figure tandis que les
hydrogrammes résultants figurent à droite.
On peut alors observer qu'il existe une durée critique de la précipitation pour laquelle le débit
de pointe est maximal. Cette durée critique est égale à la durée de concentration du bassin
versant.
Fig. 11.7 - Influence de la distribution spatiale de l'averse sur la réponse hydrologique d'un bassin
versant. Les hyétogrammes figurent sur la partie gauche de la figure tandis que les hydrogrammes
résultants figurent à droite. (d'après Réménérias ).
l'averse sur l'allure de l'hydrogramme résultant est illustrée dans la figure 11.8 dans laquelle
les hyétogrammes figurent sur la partie gauche de la figure tandis que les hydrogrammes
résultants figurent à droite.
Fig. 11.8 - Influence des variations d'intensité de l'averse sur la réponse hydrologique d'un bassin
versant). (D’après Réménérias
Avant de présenter l'analyse des crues ainsi que des événements pluies-débits, on insistera
une fois encore sur le fait que le rôle de l'ingénieur est d'identifier les facteurs dominant du
comportement hydrologique d'un bassin versant ainsi que de quantifier leur importance sur
la dite réponse. Dans cet esprit, l'étude du comportement hydrologique fait souvent appel à
des modèles permettant d'explorer plus ou moins systématiquement l'effet d'un ou de
plusieurs facteurs sur le comportement hydrologique. La modélisation est donc ici exploitée
comme outil de compréhension et d'analyse.
La réponse hydrologique dépend également fortement de l'état hydrique initial du bassin, lui-
même relié aux séquences de divers types de périodes pluvieuses et sèches qui ont précédé
l'événement pluie-débit étudié. Ainsi, une pluie tombant sur un sol sec servira d'abord à
combler le déficit d'humidité du sol. Cette même pluie sera entièrement disponible pour le
ruissellement si le sol est déjà saturé au début de l'événement.
Dans un état initial sec (nappe profonde, faible extension des surfaces saturées), la recharge
de la nappe sera beaucoup plus lente et le ruissellement par saturation beaucoup plus limité.
La figure 11.9 donne un exemple de débit simulé pour deux conditions antécédentes
particulières à savoir une condition humide et une condition sèche. La condition initiale
sèche se traduit par une première crue bien plus faible que dans la situation d'une condition
initiale humide puisque l'eau de pluie va commencer par remplir le réservoir sol avant de
générer de l'écoulement.
Fig. 11.9 - Influence des conditions antécédentes d'humidité sur le comportement hydrologique d'un
bassin versant (sol initialement humide en noir et sol initialement sec en gris).
11.5 Conclusion
Ce chapitre nous a permis d'effectuer un passage important allant des aspects descriptifs de
l'hydrologie par le biais de la compréhension des processus sous-jacents à la notion de
réponse hydrologique. Arrivant au terme de ce cours, nous avons successivement vu les
principaux facteurs qui conditionnent le comportement hydrologique d'un bassin versant. Du
rôle de la topographie à celui des conditions initiales d'humidité du sol en passant par
l'importance des facteurs météorologiques, nous avons ainsi souligné la diversité des
processus intervenant dans le cycle de l'eau à l'échelle du bassin versant mais aussi leur
complexité et leurs interrelations.
L'analyse fréquentielle est une méthode statistique de prédiction consistant à étudier les
événements passés, caractéristiques d'un processus donné (hydrologique ou autre), afin d'en
définir les probabilités d'apparition future.
Cette prédiction repose sur la définition et la mise en oeuvre d'un modèle fréquentiel, qui est
une équation décrivant le comportement statistique d'un processus. Ces modèles décrivent la
probabilité d'apparition d'un événement de valeur donnée.
L'analyse fréquentielle fait appel à diverses techniques statistiques et constitue une filière
complexe qu'il convient de traiter avec beaucoup de rigueur. Ses diverses étapes peuvent être
schématisées très simplement selon le diagramme suivant :
Avant de commencer tout travail, il est primordial de formuler clairement les buts de
l'analyse et d'adapter la démarche en conséquence. A cet égard, en hydrologie, l'un des
critères essentiels est certainement l'échelle spatio-temporelle : étudier le comportement des
crues dans un microbassin urbain (à très faible temps de concentration) avec des données de
pluie au pas de temps mensuel n'aurait pas de sens ! L'inverse est tout aussi vrai : il est
probablement inutile de disposer de pluies au pas de temps de la minute pour l'étude du
bassin versant de l'Amazone !
La validité des résultats d'une analyse fréquentielle dépend du choix du modèle fréquentiel et
plus particulièrement de son type. Diverses pistes peuvent contribuer à faciliter ce choix,
mais il n'existe malheureusement pas de méthode universelle et infaillible.
La loi normale se justifie, théoriquement par le théorème central-limite, comme la loi d'une
variable aléatoire formée de la somme d'un grand nombre de variables aléatoires. En
hydrologie fréquentielle des valeurs extrêmes, les distributions ne sont cependant pas
symétriques, ce qui constitue un obstacle à son utilisation. Cette loi s'applique toutefois
généralement bien à l'étude des modules annuels des variables hydro-météorologiques en
climat tempéré.
La loi log-normale est préconisée par certains hydrologues dont V.-T. Chow qui la justifient
en argumentant que l'apparition d'un événement hydrologique résulte de l'action combinée
d'un grand nombre de facteurs qui se multiplient. Dès lors la variable
E.-J. Gumbel postule que la loi double exponentielle, ou loi de Gumbel, est la forme limite de
la distribution de la valeur maximale d'un échantillon de valeurs. Le maximum annuel
d'une variable étant considéré comme le maximum de 365 valeurs journalières, cette loi doit
ainsi être capable de décrire les séries de maxima annuels.
Il est à remarquer que plus le nombre de paramètres d'une loi est grand, plus l'incertitude
dans l'estimation est importante. Pratiquement il est par conséquent préférable d'éviter
l'utilisation de lois à trois paramètres ou plus.
4. avec n > 0 (lois de type logarithme): loi log-normale (Galton), Pearson V, Fréchet, log-
gamma.
Cette approche suggère la plus grande prudence avec des lois de type logarithmique qui
peuvent largement surestimer les valeurs correspondant à des fréquences rares.
Le choix d'un type de modèle probabiliste est souvent basé sur des habitudes locales, qui, à
la longue, deviennent de véritables coutumes. Un tel choix résulte fréquemment de
l'expérience d'un ou de plusieurs auteurs, constituant peu à peu des écoles. L'attitude du
praticien, suivant telle ou telle école, est légitime dès lors qu'elle le fait bénéficier d'une
grande somme d'expérience, mais elle a aussi ses dangers. En effet elle peut conduire à
perpétuer un choix, parfois mal fondé, dans des situations, où il peut se révéler inadéquat.
Dans certains pays, ou dans certaines administrations, il existe en effet des règles ou normes
qui fixent la méthodologie d'une analyse fréquentielle. Pour l'étude des débits maximums, par
exemple, la loi log-Pearson III est recommandée aux Etats-Unis.
Beaucoup d'auteurs utilisent les tests d'adéquation (voir paragraphe contrôle de l'ajustement)
comme technique permettant de choisir le modèle fréquentiel approprié. Cependant il est à
remarquer qu'un test statistique ne permet que de conclure au rejet, ou à l'acceptation, de
rapport qui est une mesure d'aplatissement, peuvent être utilisés dans un
diagramme analogue à celui de K. Pearson. La figure ci-dessous illustre la position des lois
fréquemment utilisées en hydrologie.
Diagramme des L-moments : symétrie - aplatissement . Les abréviations, désignant les lois, utilisées
dans ce diagramme sont les suivantes : UNI=uniforme, NOR=normale, GUM=Gumbel, EXP=exponentielle,
GEV=extrêmes généralisées, PIII=Pearson III, LN=log-normale, GPA=Pareto généralisée.
Il est à remarquer que dans ce graphique les lois à deux paramètres sont représentées
comme un point, tandis que les lois à 3 paramètres sont figurées par une courbe.
La distribution des valeurs extrêmes provenant de n'importe quelle distribution converge vers
la loi des extrêmes généralisées (GEV). La distribution de cette loi s'exprime de la manière
suivante :
III Weibull
I Gumbel +
II Fréchet
l'expression d'un quantile est alors linéaire. En effet pour trouver la valeur d'un quantile,
Dans le cas d'un ajustement selon la loi de Gumbel, la méthode graphique repose sur le fait
que l'expression d'un quantile correspond à l'équation d'une droite. En conséquence, les
points de la série à ajuster peuvent être reportés dans un système d'axes ; il est alors
possible de tracer la droite qui passe le mieux par ces points et d'en déduire les deux
Dans la mesure où les points sont connus (ils font partie de la donnée du problème), il
suffit de définir les coordonnées correspondant à chaque point pour pouvoir le positionner
dans le graphique. Ces coordonnées se déterminent à partir de la relation inverse de la
chaque valeur .
Nom Formule
Weibull
Cunnane
Gringorten
Hazen
Des simulations ont montré que pour la loi de Gumbel, il est judicieux utiliser la distribution
empirique de Hazen.
L'ajustement graphique, bien qu'étant une méthode approximative, a le très grand avantage
de fournir une représentation visuelle des données et de l'ajustement. Celle-ci constitue un
aspect essentiel du jugement porté sur l'adéquation entre la loi choisie et les données
traitées, quelle que soit la méthode d'ajustement utilisée.
La méthode des moments consiste à égaler les moments échantillonnaux et les moments
On obtient donc les formules suivantes pour l'estimation par la méthode des moments :
Le but de cette méthode est de réaliser un ajustement lorsque les moments classiques ne
conviennent pas. Les deux paramètres et sont obtenus très simplement à partir des
valeurs des deux premiers L-moments de la loi de Gumbel et des estimations calculées sur
l'échantillon :
La solution des moindres rectangles conduit à trouver la droite bissectrice des solutions
la variable réduite de Gumbel et l'axe par celui de la variable hydrologique étudiée que
nous notons ici .
Nous obtenons par cette méthode les estimateurs suivants :
La première équation doit être résolue de façon itérative. Dans ce cas la solution de la
méthode des moments peut par exemple être utilisée comme première approximation.
L'examen visuel du graphique représentatif de l'ajustement réalisé, même s'il peut paraître
rudimentaire, reste un des bons moyens pour juger de la qualité d'un ajustement et devrait
toujours constituer un préambule à tout test statistique. La figure ci-dessous en présente un
exemple.
Ajustement de la série tronquée des débits de pointe [m3/s] du Nozon à Orny (1923-1931) à une loi
exponentielle.
A.4.2 Le test chi-carré de K. Pearson
Ce test est appliqué dans une situation où l'on observe la répartition de objets dans
classes. Il est utilisé pour tester l'hypothèse que la répartition des données s'effectue selon
une distribution théorique. On se pose donc la question de l'adéquation d'une distribution
théorique à des données.
Pour tester l'adéquation d'une répartition théorique, on dispose de deux éléments. D'une
Le score du test de Pearson est une mesure de la distance entre la répartition empirique et la
loi théorique. Elle se base sur la répartition des objets selon la loi théorique :
Ensuite, on calcule :
L'hypothèse nulle que l'on teste avec le test de Pearson est : " La distribution théorique
est la vraie distribution sous-jacente aux données ". On peut démontrer que : la distribution
de la statistique de Pearson sous l'hypothèse est bien approchée par une loi (chi-carré
Dans le cas où la variable aléatoire considérée est continue, il faut discrétiser, ce qui
introduit un élément d'ambiguïté. Pour le cas continu il existe un autre test qui utilise la
distribution empirique et qui, en règle générale, est plus puissant que le test de Pearson.
La statistique est tabulée dans plusieurs ouvrages. Le principe de ce test est illustré dans
la figure suivante :
qui permet de donner plus d'influence aux faibles et fortes fréquences. Cela conduit à la
statistique notée .
on obtient un test sensible au comportement pour des fréquences rares. Cette procédure de
test peut son se révéler particulièrement utile lorsqu'on s'intéresse, comme c'est
généralement le cas en hydrologie, aux valeurs extrêmes.
A ce stade de l'analyse nous disposons d'un modèle fréquentiel , obtenu après plusieurs
étapes. On est donc en droit de se poser la question de sa fiabilité ou degré de confiance que
l'on peut y accorder.
estimés par la méthode des moments l'expression d'un quantile peut s'écrire .
où en introduisant :
Lorsque les paramètres ont été estimés par la méthode du maximum de vraisemblance la
procédure de calcul de l'erreur-type d'un quantile se base sur la méthode delta (méthode de
linéarisation se basant sur le développement de Taylor). Pour la loi de Gumbel, on obtient :
Exercice : Estimation des débits de crue pour différents temps de retour par la méthode
statistique – Application au bassin versant de la Mentue à Yvonand (VD, Suisse)
Pour le bassin versant de la Mentue (station à Yvonand), localisé dans la région de Plateau,
en Suisse Romande, et se jetant dans le lac de Neuchâtel, on vous demande :
1) Ajuster la série des débits maximums annuels selon une distribution de Gumbel. Ajuster
les données graphiquement.
2) Estimer les débits de pointe de temps de retour, 5, 20, 50, 100 ans.
3) Ajuster les données par la méthode des moments. Estimer les débits de pointe de temps de
retour, 5, 20, 50, 100 ans.
Données :
L’exercice porte sur le bassin versant de la Mentue (station à Yvonand). Les données
nécessaires à la réalisation de cet exercice se trouvent dans le tableau ci-dessous. Il s’agit
d’une série de débits maximums annuels en [m3/s].
débit de débit de débit de
Année Année Année
pointe pointe pointe
1971 23.00 1980 33.25 1989 16.83
1972 15.40 1981 27.59 1990 27.99
1973 13.20 1982 52.66 1991 37.27
1974 19.08 1983 32.55 1992 37.99
1975 18.09 1984 30.47 1993 25.41
1976 20.81 1985 37.43 1994 21.83
1977 41.50 1986 35.47 1995 45.40
1978 30.82 1987 29.65
1979 43.67 1988 33.60
Réponse
Analyse fréquentielle des débits observés - Application au Bassin versant de la Mentue
Ajustement par la méthode graphique
nombre d'observations = 25
débit de fréquence variable
Année Tri Année Qpmax rang r
pointe empirique F réduite u
1971 23.00 1973 13.2 1 0.02 -1.364
1972 15.40 1972 15.4 2 0.06 -1.034
1973 13.20 1989 16.8 3 0.1 -0.834
1974 19.08 1975 18.1 4 0.14 -0.676
1975 18.09 1974 19.1 5 0.18 -0.539
1976 20.81 1976 20.8 6 0.22 -0.415
1977 41.50 1994 21.8 7 0.26 -0.298
1978 30.82 1971 23.0 8 0.3 -0.186
1979 43.67 1993 25.4 9 0.34 -0.076
1980 33.25 1981 27.6 10 0.38 0.033
1981 27.59 1990 28.0 11 0.42 0.142
1982 52.66 1987 29.7 12 0.46 0.253
1983 32.55 1984 30.5 13 0.5 0.367
1984 30.47 1978 30.8 14 0.54 0.484
1985 37.43 1983 32.6 15 0.58 0.607
1986 35.47 1980 33.3 16 0.62 0.738
1987 29.65 1988 33.6 17 0.66 0.878
50
40
y = 7.9767x + 25.524
it
2
éb
R = 0.9679
30
d
20
10
0
-2 -1 0 1 2 3 4 5
u : Variable centrée et réduite
Chow V. T., Applied Hydrology, Civil Engineering Series, McGraw-Hill International Edition,
1988.
Cosandey C., Robinson M., Hydrologie continentale, Armand Collin, Paris, 2000.
Gleick, P.H. (Editor),. Water in Crisis. A guide to the World's Fresh Water Resources. Oxford
Hillel D., L'eau et le sol -Principes et processus physiques, Vander-éditeur, Louvain, 1974.
Meylan P., Musy A., Hydrologie Fréquentielle, Edition HGA Bucarest, 1999
Romandes, 1997
Musy A., Higy, C., Hydrologie Appliquée, Edition HGA Bucarest, 1998
Musy A., Soutter M., Physique du sol, Presses Polytechniques Universitaires Romandes,
Lausanne, 1991.
Singh V.P., Computer models of watershed hydrology. Water Resources Publications, 1129
pp., 1995.
Viesman, W., Knapp, J., Lewis, G., Harbaugh, T., Introduction to Hydrology, Harper and