- 1 http://www.amice-project.eu
1En 2009, l’EPAMA (Établissement Public d’Aménagement de la Meuse et de ses Affluents) a initié un projet concernant des solutions d’adaptation de la Meuse aux impacts des extrêmes hydrologiques liés au changement climatique. Ce projet a pris le nom d’AMICE pour : Adaptation of the Meuse to the Impacts of Climate Evolution1et a été retenu par le programme européen INTERREG IVB ENO. Il regroupe 17 partenaires provenant des 4 pays contrôlant une partie du bassin versant de la Meuse (France, Belgique, Allemagne et Pays-Bas).
2Nous présentons ici un des aspects d’AMICE relatif à l’estimation de l’évolution des débits d’étiage de la Meuse au cours du XXIème siècle (Drogue & al., 2010). La méthode choisie consiste à alimenter un modèle hydrologique empirique par des données de scénarios climatiques (un scénario sec et un scénario sec modéré). Elle tient compte également de l’incertitude liée à l'estimation des débits par le modèle. La valeur du débit mensuel d’étiage (QMNA) de fréquence 1/5 calculée sur 3 périodes trentenaires (1971-2000, 2021-2050 et 2071-2100) sert d’indicateur de cette évolution.
3La démarche mise en œuvre pour estimer les QMNA de la Meuse comporte 3 étapes (fig. 1) :
Figure 1 : Méthodes mises en oeuvre
- 2 Nous avons choisi cette période car nous disposions déjà d’un référentiel statistique sur les débi (...)
4Dans un premier temps, nous nous sommes intéressés aux données issues de scénarios climatiques sur la base de simulations réalisées par le modèle ARPEGE-Climat v.4 de Météo- France. Le modèle numérique ARPEGE est un modèle de circulation générale global et spectral développé par Météo-France en collaboration avec le CEPMMT pour la prévision numérique du temps. La version climat d'ARPEGE, nommée ARPEGE-Climat, a été développée dans les années 1990. Elle est dérivée du modèle de prévision opérationnelle à courte échéance. La grille d'ARPEGE peut être basculée pour changer la position du pôle et étirée pour augmenter la résolution dans une zone d'intérêt. Cette capacité de zoom a permis de développer des études régionales du climat avec ARPEGE-Climat notamment sur l’Europe (Déqué et Piedelievre, 1995). Dans la version climat d’ARPEGE, la résolution varie de 50 km au centre de la Méditerranée à 450 km dans le Pacifique Sud. Ceci assure une résolution d’au moins 60 km sur la France. La discrétisation verticale se fait sur 31 couches qui suivent le relief avec des tranches d’altitude variables. Le pas de temps de calcul est de 30 mn. Nous avons caractérisé le climat actuel (1971-20002) à partir des données provenant de E-OBS 2.0. Cette base de données hébergée sur le site du "European Climate Assessment & Dataset" regroupe des températures de l’air et des précipitations journalières spatialisées selon une grille de résolution de 0.25° (environ 17 km en longitude et 27 km en latitude dans la zone du bassin de la Meuse) sur la période 1950-2008 (Haylock & al., 2008). Les valeurs moyennes de températures et de précipitations aux stations hydrométriques résultent d'une spatialisation par des polygones de Thiessen. Les scénarios climatiques utilisés pour cette étude procèdent d'une estimation par Météo-France des tendances saisonnières pour les périodes 2021-2050 et 2071-2100 par rapport au climat actuel (méthode des anomalies). Ces tendances ont été appliquées globalement au bassin de la Meuse sans tenir compte de nuances locales qui apparaissent dans les données d’échelle fine (Commeaux et al., 2010). Elles concernent 2 scénarios d'émission de gaz à effet de serre siglés A2 et A1B appliqués au modèle ARPEGE Climat v.4. L’approche adoptée dans le cadre de cette analyse n’est donc pas ensembliste. Elle se fonde sur les simulations d’un seul modèle de circulation générale et ne prend en compte que les deux scénarios d’émission A1B et A2 produits dans le cadre des expériences numériques d’intercomparaison du 3e rapport d’avancement du GIEC paru en 2001. L’incertitude réflexive (Pagé et Terray, 2010) est donc prise en compte, mais de façon limitée. L’objectif est de fournir une première indication du signe et de l’intervalle de variation de la tendance affectant les débits d’étiage pour deux futurs scénarios climatiques répondant à deux trajectoires sociétales équiprobables: i) une trajectoire « médiane » (scénario A1B) tablant sur une forte croissance économique qui suppose des échanges mondiaux importants, ii) une trajectoire « tendancielle » (scénario A2). Ce dernier table sur une croissance économique qui met l’accent sur les échanges régionaux. Le scénario A2 est le plus pessimiste des scénarios utilisés et simule la situation qu'il adviendrait si aucune mesure n'était prise pour restreindre les émissions de gaz à effet de serre. Dans ce scénario la concentration en CO2 en 2100 est de 840 ppm contre 770 ppm pour le scénario A1B. Par souci de simplification nous désignons respectivement A2 et A1B, scénario "sec modéré" et scénario "sec". Pour obtenir des séries de données représentatives des périodes 2021-2050 et 2071-2100, les tendances saisonnières de chaque scénario exprimées en pourcentage pour la pluie et en écart absolu pour les températures de l’air ont été imputées aux données du climat actuel issues de la climatologie E-OBS. Les valeurs d'évapotranspiration potentielle proviennent de l'application de la formule à base thermique de Oudin (Oudin & al., 2005).
- 3 http://www.cemagref.fr/webgr/
5La deuxième phase a consisté à caler un modèle hydrologique pluie-débit journalier avec les données observées de la période 1971-2000. Nous avons utilisé le modèle GR4j (Perrin et al., 2003) car il est abondamment documenté3, il nécessite peu de paramètres et a déjà montré de bonnes performances (en comparaison avec d'autres modèles) sur la région d’étude. Son architecture est présentée sur la figure n°2. Pour son calage, nous avons utilisé comme fonction- objectif le critère de Nash-Sutcliffe calculé sur le logarithme des débits. Pour les 7 stations étudiées, les valeurs du critère de Nash-Sutcliffe s’échelonnent de 77 à 90 % et correspondent à des niveaux d’efficience satisfaisants pour simuler des QMNA. Les sorties du modèle pluie-débit font ensuite l'objet d'une correction et d'une estimation de leur incertitude selon une méthode empirique développée par Berthier (2005). Pour chaque classe de débits journaliers (d'amplitude exponentielle), on calcule le rapport médian débit journalier observé / débit journalier simulé (quantile 50 % période de référence). Celui-ci sert de "coefficient correcteur" (par classe de débit) et permet d'améliorer la qualité des simulations (par exemple la correction du biais sur les faibles valeurs à Chalaines, fig. 3). Afin de s’assurer de la représentativité des résultats, les quantiles d’incertitude figurant sur le graphique concernent la période de calage (15 premières années, en bleu) et la période de référence (30 ans en rouge). Les paramètres de calage du modèle et les intervalles d'incertitudes déterminés ici pour la période de référence 1971-2000 serviront par la suite pour les simulations résultant des scénarios climatiques 2021-2050 et 2071-2100.
6La troisième et dernière phase reprend les résultats et méthodes des deux premières. Les données des scénarios climatiques alimentent GR4j, les débits ainsi simulés passent l'opération de correction et au final, on extrait les valeurs de QMNA. Celles-ci sont enfin ajustées à une loi log-normale pour définir le quantile quinquennal sec (seuil fréquemment utilisé en France pour caractériser un étiage sévère).
Figure 2 : Le modèle GR4j
Figure 3 : Correction et incertitudes des débits simulés – La Meuse à Chalaines
7Les bassins versants étudiés se situent principalement en France (fig. 4). Ils sont drainés par la Meuse et ses affluents. Celle-ci prend sa source dans la région du Bassigny. Jusqu’à la station de Goncourt, elle draine essentiellement des terrains imperméables marneux ou argileux. En aval jusqu’à la station de Stenay, le fleuve entaille successivement les plateaux calcaires du Dogger puis le plateau de l’Oxfordien. Les premiers concentrent une densité relativement importante de phénomènes karstiques qui peuvent être à l’origine de pertes (en amont de Neufchateau) et de circulations souterraines vers des bassins contigus (la Moselle ou la Marne). Ces échanges s’effectuent dans la plupart des cas aux dépends du bassin de la Meuse. En aval de Stenay, la Chiers, principal affluent de la Meuse (qui apparait particulièrement démembrée), draine en rive gauche les terrains du Dogger et en rive droite les grès d’Ardenne-Luxembourg qui assurent un soutien important en période d’étiage. En aval de Charleville-Mézières, la Meuse entaille le massif ardennais. Celui-ci ne dispose de ressource aquifère que très localement. Par contre il constitue le pôle pluviométrique du bassin avec ponctuellement plus de 1300 mm /an. Dans la partie amont du bassin la pluviométrie demeure très homogène avec des hauteurs annuelles comprises entre 800 et 900 mm/an.
8Nous avons sélectionné 7 stations hydrométriques présentant des chroniques de débits relativement longues. La plus récente a débuté en 1973 (Chalaines) et toutes ces stations fonctionnent encore actuellement (Tableau 1). Elles sont toutes gérées par la DREAL Lorraine et leurs données proviennent de la Banque Hydro. Elles se répartissent assez régulièrement le long du cours de la Meuse avec des bassins versants relativement grands compris entre 364 et 10120 km² (fig. 4).
Figure 4 : Domaine d’étude
Tableau 1 : Stations hydrométriques étudiées
9Les sorties d'ARPEGE Climat pour les 2 scénarios indiquent une augmentation des températures pour toutes les saisons. Comparativement à 1971-2000 cette hausse se situe entre +1.2°C et +2.1°C pour 2021-2050 et +2.6°C à +5.6°C pour 2071-2100 (Tableau 2). Pour la modélisation hydrologique, l'ETP répercute cette évolution car son estimation s'appuie uniquement sur les données de température de l’air.
10A cette élévation de la température de l’air s'ajoute une diminution globale des précipitations saisonnières pouvant aller jusqu'à -39% (été 2071-2100, Tableau 2). Le scénario sec modéré envisage juste une petite augmentation (+4%) des précipitations printanières sur la période 2021-2050. L'écart entre les 2 scénarios s'avère peu important hormis pour les précipitations hivernales 2071-2100.
Tableau 2 : Scénarios climatiques, évolution des cumuls de précipitations et des températures moyennes
11Les principaux résultats sont présentés sous la forme de profils hydrologiques qui mettent en relation l’évolution des débits d'amont en aval avec la surface du bassin versant drainé.
12Le profil de la Meuse permet de juger la capacité du modèle hydrologique à simuler les débits mensuels d'étiage. Pour Goncourt (station la plus en amont) le modèle peine à produire des simulations d'étiage comparables aux valeurs observées (fig 5). Avec les valeurs simulées et corrigées, on atteint à peine la borne supérieure de l'intervalle de confiance statistique (90%) de l'ajustement des QMNA. Il s'agit du plus petit bassin versant (364 km²) pour lequel les précipitations moyennes issues de la climatologie E-OBS sont peut-être sur-estimées.
Figure 5 : Profil hydrologique de la Meuse française – Débits mensuels d’étiage quinquennal sec - Période 1971-2000
13En revanche, pour les stations aval avec des bassins versants supérieurs à 1000 km², le modèle présente des QMNA simulés et corrigés proches des QMNA observés (profils confondus à partir de Stenay). La procédure de correction exposée plus haut (Berthier, 2005) renforce parfois très avantageusement la qualité des simulations comme par exemple à Saint-Mihiel.
14Si les scénarios climatiques sec et sec modéré présentent, somme toute, assez peu d'écarts (cf. § 3.1), les QMNA produits par GR4j à partir de ces scénarios se différencient encore moins (Tableau 3). Quelle que soit l’échéance envisagée (2021-2050 ou 2071-2100), nous pouvons considérer qu'il n'existe qu'une seule hypothèse de travail, les deux scénarios aboutissant à des résultats similaires. L'écart significatif constaté pour les précipitations hivernales ne se répercute pas sur les débits d'étiage simulés.
15Nous avons choisi d’adopter comme référence les QMNA 1/5 simulés et corrigés pour calculer les écarts avec les QMNA des périodes 2021-2050 et 2071-2100 car ces derniers résultent des mêmes procédures (modélisation et correction). S’il existe un biais comme à Goncourt, on peut concevoir que les valeurs de QMNA projetées pour nos périodes-objectifs égalent le QMNA observé plus les écarts présentés dans le tableau 3.
Tableau 3 : Débits mensuels d’étiage quinquennaux sec (QMNA 1/5)
16Dans la partie amont du bassin, en raison des incertitudes du modèle, on ne peut pas clairement estimer une décroissance significative du QMNA 1/5 2021-2050 par rapport à la période 1971-2000 (fig. 6). Par contre de Chalaines à Montcy-Notre-Dame, cette baisse apparaît plus marquée avec des valeurs relatives comprises entre -28 et -44% (Tableau 3). Le profil hydrologique 1971-2000 ne chevauche pas le rail d'incertitude à 70 % des QMNA 1/5 simulés de la période 2021-2050. Par contre à l'extrémité aval du bassin, l'écart entre les QMNA des 2 périodes s'atténue (-26% à -25%) et le QMNA 1/5 2021-2050 dépasse la borne inférieure de l'intervalle de confiance à 90% du QMNA 1/5 1971-2000.
Figure 6 : Profil hydrologique de la Meuse française – Débits mensuels d’étiage quinquennal sec - Période 1971-2000 vs 2021-2050
17L'indication la plus spectaculaire concerne les stations de Stenay et Montcy où le QMNA 1/5 2021-2050 équivaudrait au débit mensuel d'étiage minimum connu sur la période de référence 1971-2000. Cet étiage correspond à l'année 1976 habituellement perçue comme particulièrement sévère voire "catastrophique". En fait, en suivant notre hypothèse de travail et pour ces 2 stations, d'ici le milieu du XXIème siècle, les niveaux d'étiage de 1976 seraient sous-passés en moyenne tous les 5 ans.
18Pour la fin du XXIème siècle, l'écart des QMNA 1/5 avec ceux de la période 1971-2000 s'amplifie (fig. 7). Il n'y a qu'à Goncourt où le rail d'incertitude dépasse encore la valeur du climat actuel. Partout ailleurs le profil de la période 2071-2100 apparaît très proche du débit minimum journalier connu sur la période de référence 1971-2000. En d'autres termes pour situer une valeur qui statistiquement n'est pas rarissime (le QMNA 1/5), nous avons dû nous référer au minimum absolu de la période 1971-2000. D'ici à la fin du siècle, les valeurs de QMNA 1/5 devraient globalement diminuer de moitié.
Figure 7 : Profil hydrologique de la Meuse française – Débits mensuels d’étiage quinquennal sec - Période 1971-2000 vs 2071-2100
19Actuellement on se situe en moyenne 5% (env. 18 jours/an) du temps sous le seuil du QMNA 1/5 (seuil de crise, fig. 8). On sous-passerait ce même seuil entre 6 (env. 22 jours/an) et 22 % (env. 80 jours/an) du temps en 2021-2050 et entre 15 (env. 55 jours/an) et 52% (env. 200 jours/an) en 2071-2100. Notre définition des seuils de crise relatifs aux étiages se voient ainsi remis en cause et corrélativement certains usages de l'eau.
Figure 8 : Proportion de débits journaliers inférieurs au QMNA 1/5 (intervalle d’incertitude de 70 % inclus)
20A l'évidence, ce travail n'a pas comme ambition de représenter une prévision des débits d'étiage de la Meuse pour le XXIème siècle. Il s'agit davantage d'un "exercice de style" où la contrainte est de produire des débits d'étiage à partir de scénarios climatiques hypothétiques et équiprobables. Nous ne disposons d'aucun moyen pour valider nos simulations, les seules garanties que nous pouvons proposer concernent la qualité de la simulation climatique de référence et des débits simulés sur une période connue. La crédibilité des scénarios climatiques demeure prépondérante, pour que les débits d'étiage "possibles" proposés ici soient considérés comme des valeurs "probables". Ainsi un certain nombre d'hypothèses fortes ont jalonné notre démarche.
21Pour la première phase, l'élaboration des scénarios climatiques résultant d'une modification globale des abats de précipitations journaliers par saison et de la température de l’air selon la même méthode peut paraître simpliste. Nous n'avons pas envisagé la possibilité d'apparition plus fréquente de longs épisodes secs aggravant les situations d'étiage. Ces épisodes peuvent également devenir plus durables.
22Les simulations pour le XXIème siècle proviennent d'un modèle hydrologique calé sur la période 1971-2000. Ce calage nous semble robuste pour notre période de référence avec l'emploi d'une procédure de calage/contrôle. D’autre part nous avons testé différentes périodes de calage (années sèches, années humides) et constaté une faible sensibilité des paramètres. Mais dans le cadre non-stationnaire de l'évolution du climat que nous étudions, les paramètres imposés au modèle hydrologique pour 2021-2050 et 2071-2100 ne correspondent peut-être pas aux paramètres optimaux pour estimer les écoulements de ces périodes. Enfin les rails d'incertitude et la correction des débits simulés par le modèle hydrologique résultent également de calculs effectués sur notre période de référence. Leur transposition au XXIème siècle constitue un autre postulat d'autant plus téméraire que la distribution des débits simulés en fonction de leurs erreurs relatives (fig. 3) sera immanquablement différente. Une classe de débit peu représentée dans le climat actuel (comme les étiages les plus sévères) verra son effectif augmenter durant les périodes-objectifs (2021-2050 et 2071-2100).
23Ces réserves s'appliquent à la plupart des travaux prospectifs. Notre démarche est surtout pragmatique et a pour souci de montrer l'incertitude parfois importante liée à des débits simulés malgré des critères d'ajustement du modèle pluie-débit satisfaisants. Elle montre aussi la nécessité de recourir à une démarche ensembliste permettant de probabiliser les futurs climatiques utilisés pour forcer les modèles d’impact en hydrologie.
24A l'échelle du bassin de la Meuse, l'évolution des QMNA 1/5 à partir d'un modèle hydrologique et en fonction des changements climatiques liés aux émissions additionnelles de gaz à effet de serre d’origine anthropique, n’apparaît pas toujours décelable. Pour le petit bassin versant de la Meuse à Goncourt l'incertitude des estimations issues du modèle pluie-débit ne permet pas de démontrer une diminution significative même pour la période 2071-2100.
25En revanche aux autres stations, si pour notre première période-objectif (2021-2050) la baisse des QMNA 1/5 peut sembler modérée, pour la seconde (2071-2100) la chute serait spectaculaire : le QMNA 1/5 se rapprocherait du débit journalier minimum connu sur la période 1971-2000.
26Les partenaires du projet AMICE vont proposer une stratégie d'adaptation sous-tendue par cette évolution "possible" des ressources en eau superficielle d'ici la fin du XXIe siècle.
27Du côté des géographes hydrologues et des hydrologues « globaux », ce type de travail prospectif lié aux évolutions du climat fait l'objet de nombreuses publications et ouvre des pistes de recherche particulièrement intéressantes (Brigode, 2010).