Etude Des Relations Nappe-Rivière En Haute-Normandie, Bassin Versant De L’Austreberthe Modélisations Hydrologiques GARDENIA

Etude des relations nappe-rivière en Haute-Normandie, bassin versant de l’Austreberthe Modélisations hydrologiques GARDENIA

BRGM/RP-55813-FR Novembre 2007

Etude des relations nappe-rivière en Haute-Normandie, bassin versant de l’Austreberthe Modélisations hydrologiques GARDENIA

BRGM/RP-55813-FR Novembre 2007

Étude réalisée dans le cadre des opérations

de recherche EAU-R18

L. Arnaud

Vérificateur : Approbateur :

Nom : Stollsteiner P. Nom : Pasquet J.-F.

Date : 23 novembre 2007 Date : 23 novembre 2007

Signature : Signature :

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Mots clés : Aquifère crayeux, Bilan hydrologique, Modélisation Gardénia, Pluie efficace, Recharge, Relation nappe-rivière, Ruissellement, Austreberthe, Blacqueville, Seine-Maritime, Haute-Normandie.

En bibliographie, ce rapport sera cité de la façon suivante : Arnaud L. – Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA. Rapport BRGM/RP-55813-FR, 55 p., 26 ill., 2 ann.

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Synthèse

Dans le but d’approfondir la connaissance relative aux relations nappe-rivière en Haute-Normandie, l’Agence de l’Eau Seine-Normandie (AESN) et le BRGM se sont associés dans le cadre d’un projet de recherche. La méthodologie proposée par le brgm se décompose en deux approches complémentaires : • une approche par modélisation hydrologique globale du bassin versant, objet du présent rapport ;

• une approche par traçages isotopiques qui devrait permettre d’évaluer la nature des relations entre les nappes, les sources et la rivière (Rapport BRGM/RP- 55983-FR).

Le Comité de Pilotage du projet, constitué de l’AESN, de la DIREN Haute-Normandie et du BRGM a, dans un premier temps, décidé de tester la méthodologie sur le bassin versant de l’Austreberthe (Seine-Maritime).

Des modélisations satisfaisantes (pluie-débit et pluie-niveau) ont pu être obtenues pour le bassin versant de l’Austreberthe. Les différents termes du bilan hydrologique moyen calculé sur la station de jaugeage de Saint-Paër, permettent de tirer les enseignements suivants :

• les précipitations efficaces annuelles moyennes sont de l’ordre de 375 mm. Cette valeur est cohérente avec la cartographie nationale des pluies efficaces élaborée par le MEDAD sur la période 1946-2001 ;

• avec un débit de fuite négligeable (prélèvements non restitués au bassin versant + écoulement souterrain au droit de la station de jaugeage), le bilan hydrologique ne fait apparaître aucun déséquilibre quantitatif du bassin versant ;

• la contribution moyenne des écoulements souterrains lents au débit de l’Austreberthe est de l’ordre de 90 %. Cette contribution moyenne est en accord avec les travaux actuellement menés par l’Université de Rouen. En période d’étiage (juillet à octobre), le cours d’eau est quasi-exclusivement alimenté par les eaux souterraines. En période hivernale (décembre à mars), même si la contribution du ruissellement augmente, l’alimentation de l’Austreberthe par les eaux souterraines reste prépondérante (en moyenne entre 86 et 89 %). Des écoulements rapides significatifs ne sont observés sur le bassin versant que lors des hivers exceptionnellement pluvieux où leur contribution peut alors atteindre entre 20 et 30 %.

En ce qui concerne les deux campagnes de prélèvements réalisées dans le cadre des traçages isotopiques, en novembre 2006 et en mai 2007, le modèle semble indiquer une contribution des écoulements lents à l’alimentation de l’Austreberthe sensiblement identique : respectivement 98 et 96 %.

BRGM/RP-55813-FR – Rapport final 3

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Sommaire

1. Introduction...... 9

2. Contextes géologique et hydrogéologique du bassin versant de l’Austreberthe ...... 11

2.1. SITUATION GEOGRAPHIQUE ...... 11

2.2. CONTEXTE GEOLOGIQUE...... 12

2.3. RESSOURCES EN EAU SOUTERRAINE...... 15

2.3.1. L’aquifère crayeux...... 15 2.3.2. L’exploitation de la ressource...... 17

3. Modélisations GARDENIA...... 19

3.1. LE MODELE HYDROLOGIQUE GLOBAL GARDENIA DU BRGM...... 19

3.1.1. Principe du modèle ...... 19 3.1.2. Fonctionnement du modèle...... 20 3.1.3. Equations du bilan hydrologique ...... 22

3.2. APPLICATION AU BASSIN VERSANT DE L’AUSTREBERTHE...... 24

3.2.1. Les données d’entrée...... 25 3.2.2. Les données de sortie...... 27 3.2.3. Le calage du modèle...... 32 3.2.4. Bilan hydrologique moyen...... 37 3.2.5. Contribution des écoulements lents ...... 38

4. Conclusions ...... 43

5. Bibliographie...... 45

BRGM/RP-55813-FR – Rapport final 5 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Liste des illustrations

Illustration 1 – Carte de localisation du bassin versant de l’Austreberthe (en hachuré)...... 11 Illustration 2 – Modèle Numérique de Terrain pour le bassin versant de l’Austreberthe (source : IGN)...... 12 Illustration 3 – Carte géologique au 1 / 50 000 au droit du bassin versant de l’Austreberthe (carte n°76 – Yvetot, 1974, BRGM ©) ...... 14 Illustration 4 – Extension des terrains crétacés (extrait de la notice de la carte n°76 – Yvetot, 1974, BRGM©) ...... 15 Illustration 5 – Carte piézométrique sur le bassin versant de l’Austreberthe à l’échelle du 1 / 100 000 (BRGM, 1990) ...... 17 Illustration 6 – Evolution des prélèvements en eau sur le bassin versant de l’Austreberthe entre 1994 et 2005 (source AESN) ...... 18 Illustration 7 – Représentation schématique des éléments du bilan hydrologique d’un bassin versant ...... 19 Illustration 8 – Fonctionnement du modèle hydrologique GARDENIA pour la simulation du débit d'un cours d'eau ou d'un niveau piézométrique (Schéma classique)...... 21 Illustration 9 – Localisation des différentes stations de mesures disponibles sur le bassin versant ...... 24 Illustration 10 – Répartition annuelle des précipitations et de l’ETP sur la station climatique de Rouen-Boos (source : Météo France) ...... 25 Illustration 11 – Répartition de la pluie efficace annuelle moyenne (1946-2001) sur le bassin versant de l’Austreberthe (source : MEDAD) ...... 26 Illustration 12 – Graphe des double cumuls pour les stations climatiques d’Auzebosc et de Boos (source : Météo France) ...... 27 Illustration 13 – Caractéritiques des 2 stations de jaugeage suivies par la DIREN Haute- Normandie...... 28 Illustration 14 – Graphique de corrélation entre les lames d’eau journalières ruisselées à Saint-Paër et à Duclair ...... 29 Illustration 15 – Graphique de double cumul pour les lames d’eau ruisselées à Saint- Paër et à Duclair...... 29 Illustration 16 – Débits caractéristiques de l’Austreberthe à la station hydrologique de Saint-Paër (source : DIREN Haute-Normandie)...... 30 Illustration 17 – Hauteur d’eau précipitée à Rouen et débits de l’Austreberthe en 2001 (au pas de temps journalier)...... 30 Illustration 18 – Chronique piézométrique sur Blacqueville de 1968 à 2006 ...... 31 Illustration 19 – Visualisation sur la période 1997-2006 du calage pluie-débit de l’Austreberthe à Saint-Paër sur la période 1997-2006 (en bleu : la part des écoulements lents) ...... 33

6 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 20 – Reconstitution des débits de l’Austreberthe à Saint-Paër depuis 1979 (en bleu : la part des écoulements lents) ...... 34 Illustration 21 – Visualisation sur la période 1980-2006 du calage pluie-niveau à Blacqueville sur la période 1980-2006 ...... 36 Illustration 22 – Bilans hydrologiques annuels moyens calculés pour chacune des modélisations sur la période 1997-2006 ...... 37 Illustration 23 – Représentation schématique du bilan hydrologique à Saint-Paër ...... 37 Illustration 24 – Contributions des écoulements souterrains lents au débit de l’Austreberthe au pas de temps mensuel (1969-2006) ...... 38 Illustration 25 – Evolution des paramètres hydro-climatiques pour le cycle hydrologique 2006-2007...... 40 Illustration 26 – Synthèse des données piézométriques et hydrologiques pour les deux campagnes de prélèvements...... 41

Liste des annexes

Annexe 1 – Compte-rendu du calage de la modélisation pluie-débit sous GARDENIA ...... 47 Annexe 2 – Compte-rendu du calage de la modélisation pluie-niveau sous GARDENIA ...... 51

BRGM/RP-55813-FR – Rapport final 7

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

1. Introduction

Dans le but d’approfondir la connaissance relative aux relations nappe-rivière en Haute-Normandie, l’Agence de l’Eau Seine-Normandie (AESN) et le BRGM se sont associés dans le cadre d’un projet de recherche. La méthodologie proposée par le BRGM se décompose en deux approches complémentaires : • une approche par modélisation hydrologique globale du bassin versant, objet du présent rapport ;

• une approche par traçages isotopiques qui devrait permettre d’évaluer la nature (qualitativement et éventuellement quantitativement) des relations entre les nappes (aquifère crayeux, nappe d’accompagnement), les sources et la rivière. Dans le but d’étudier l’impact des variations saisonnières sur ces relations (contribution de la nappe à l’écoulement de la rivière à l’étiage), deux campagnes de mesures ont être menées (Rapport BRGM/RP-55983-FR).

Les modélisations hydrologiques (pluie-débit et pluie-niveau) ont été réalisées à l’aide du logiciel GARDENIA, développé par le BRGM. L’objectif principal de ces modélisations est ici de différencier les deux composantes de l’écoulement : • celle rapide que l’on peut assimiler au ruissellement superficiel et aux écoulements karstiques ; • celle lente que l’on peut assimiler à la composante lente des écoulements souterrains (milieux poreux et fissurés), c’est-à-dire à la contribution des aquifères au débit des rivières.

Après une présentation des contextes géologique et hydrogéologique du bassin versant étudié, il sera procédé à la description et à l’analyse critique des données utilisées. La présente étude repose sur les données issues des différentes stations de mesures disponibles (données climatiques, piézométriques, hydrologiques et de prélèvements). Enfin, les résultats obtenus seront commentés et plus particulièrement détaillés pour les périodes correspondant aux deux campagnes de prélèvements relatives à l’approche isotopique.

BRGM/RP-55813-FR – Rapport final 9

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

2. Contextes géologique et hydrogéologique du bassin versant de l’Austreberthe

2.1. SITUATION GEOGRAPHIQUE

En rive droite de la Seine, le point de confluence de l’Austreberthe et du fleuve est situé sur la commune de Duclair, au sommet de la troisième boucle de la Seine à partir de l’estuaire (Illustration 1).

Illustration 1 – Carte de localisation du bassin versant de l’Austreberthe (en hachuré)

L’Austreberthe prend sa source sur la commune de Sainte-Austreberthe. Son seul affluent pérenne est le Saffimbec (rive droite), ce dernier prend naissance au niveau de résurgences karstiques, sur la commune de Limésy. Le point de confluence des deux cours d’eau est situé sur la commune de Pavilly.

La superficie du bassin versant hydrographique de l’Austreberthe est de 215 km².

L’altitude moyenne du bassin versant est de 126 m NGF avec une altitude minimale de 3 m NGF et une altitude maximale de 179 m NGF (Illustration 2).

BRGM/RP-55813-FR – Rapport final 11 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 2 – Modèle Numérique de Terrain pour le bassin versant de l’Austreberthe (source : IGN)

2.2. CONTEXTE GEOLOGIQUE

D’après la carte géologique à 1/50 000 n°76 – Yvetot (1974, BRGM ©) et les données BSS (Banque de données du sous-sol gérée par le BRGM et accessible au public ), le sous-sol du bassin versant de l’Austreberthe est composé d’une superposition de plusieurs formations, avec du sommet vers la base (Illustration 3) : • des limons non différenciés (notation LP ; Illustration 3) couvrant le sommet des plateaux ; • un deuxième ensemble de formations superficielles, constitué d’argiles résiduelles à silex (notation RS ; Illustration 3) provenant de l’altération de la craie sous-jacente. Ces dépôts semblent présenter une épaisseur variant de

12 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

quelques mètres à 10 mètres sur le bassin versant étudié. Les Argiles à silex affleurent en tête des vallons secs du bassin versant ainsi que sur les différents versants ; • un substratum composé de craies : - la craie blanche du Santonien (notation C5) affleurant sur les versants de l’Austreberthe dans sa partie aval (au Sud de Pavilly). Sa puissance est variable, elle est de 30 à 40 m près de Barentin et Pavilly ; - la craie blanche grisâtre du Coniacien (notation C4) affleurant également dans la basse vallée de l’Austreberthe ; - la craie grise blanchâtre à rares silex du Turonien (notation C3) affleurant dans les vallées du Saffimbec et de l’Austreberthe en amont de Pavilly. Son épaisseur varie entre 40 et 80 m sur la partie amont du bassin versant ; - craie grise glauconieuse à silex du Cénomanien moyen et supérieur (notation C2) qui n’affleure que dans la vallée du Saffimbec. L’épaisseur de cet ensemble est de 30 à 35 mètres.

D’un point de vue structural, la partie amont du bassin versant de l’Austreberthe est traversé par un système complexe de failles et de flexures prolongeant la faille de la Seine (Illustrations 3 et 4) : une faille reliant Sainte-Austreberthe et Pavilly (de direction varisque) recoupe la haute vallée de l’Austreberthe et présente un rejet de 70 m.

Trois failles, de direction armoricaine, ont été reconnues au nord-ouest de Pavilly, la faille nord recoupe la vallée du Saffimbec et présente un panneau nord-est affaissé et un rejet de 20 à 30 m près de Pavilly et probablement de 70 à 80 m près de Mesnil- Panneville. Les deux failles méridionales présentent au contraire des panneaux sud- ouest affaissés avec un rejet total évalué à environ 200 mètres. Ce dernier secteur est finalement intercalé entre deux régions séparées par la faille de la Seine et son prolongement : une région nord relevée et une région sud abaissée (Cf. Illustration 4).

BRGM/RP-55813-FR – Rapport final 13 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 3 – Carte géologique au 1 / 50 000 au droit du bassin versant de l’Austreberthe (carte n°76 – Yvetot, 1974, BRGM ©)

14 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 4 – Extension des terrains crétacés (extrait de la notice de la carte n°76 – Yvetot, 1974, BRGM©)

2.3. RESSOURCES EN EAU SOUTERRAINE

2.3.1. L’aquifère crayeux

La principale ressource en eau souterraine du bassin versant de l’Austreberthe est l’aquifère crayeux. Quelques rares ouvrages captent tout de même la nappe alluviale (§ 2.3.2).

La craie est une roche carbonatée résultant du dépôt marin de boues et de coquilles de micro-organismes. Pour différentes raisons géologiques (tectoniques, variations climatiques ...), la craie a été affectée par de nombreux accidents (failles, fissures ....) et soumise à des altérations dont l'un des résultats est le développement de réseaux karstiques. La craie est ainsi un milieu mixte : poreux, fissuré, et « karstique ». Lorsque la craie est compacte (sous les plateaux et à grande profondeur), les eaux s’écoulent dans les pores intergranulaires qui ont des dimensions de l’ordre du micron. Les valeurs de perméabilité sont égales ou inférieures à 1.10 -6 m/s. Les vitesses de circulation de la nappe sont alors de l’ordre de 1 à 10 mètres par an. Dans les vallées sèches et humides, des fissures d’ouverture millimétrique voire centimétrique affectent ce milieu, mais il est encore assimilable à un milieu poreux car les écoulements ont un

BRGM/RP-55813-FR – Rapport final 15 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

régime laminaire et respectent la loi de Darcy. Les perméabilités ont alors des valeurs comprises entre 1.10 -5 et 1.10 -3 m/s. Les vitesses d’écoulement peuvent atteindre plusieurs dizaines et même plusieurs centaines de mètres par an. Lorsque l’ouverture des fissures s’agrandit et que les conditions de charge le permettent, le régime d’écoulement devient turbulent, avec des vitesses atteignant plusieurs centimètres à la seconde. Les écoulements obéissent à d’autres lois que celles de Darcy et parfois à d’autres conditions (flux) que les conditions de pression. Ces réseaux sont connus dans toute la région recouverte par la craie, mais ils n’intéressent que certains talwegs, certaines vallées sèches et humides, et n’occupent qu’une petite partie de l’aquifère.

Des opérations de traçages menées sur le bassin versant de l’Austreberthe et du Saffimbec ont mis en évidence de telles circulations karstiques. Sur l’Illustration 5 (carte piézométrique), les flèches vertes représentent les traçages positifs et les points rouges représentent les bétoires (pertes karstiques). Ces données sont issues de l’inventaire pilote bétoires-traçages réalisé par le BRGM en 2004 (Rapport BRGM/RP- 52423-FR). Il est rappelé que cet inventaire pilote, préalable à un futur inventaire régional, n’est pas exhaustif.

Sur le bassin versant, la surface piézométrique de la nappe épouse fortement la morphologie du sol (Illustration 5) qui dépend en partie de la répartition de la fissuration de la craie. Elle forme des dômes d’alimentation sous les plateaux où l’aquifère peu fissuré a une fonction capacitive, et des dépressions dans les zones fissurées drainantes à fonction transmissive (vallées humides et sèches, réseaux « karstiques »).

La nappe s’écoule par gravité depuis le sommet des plateaux vers les sources et les cours d’eau du Saffimbec et de l’Austreberthe (Cf. figures de drainance sur la carte piézométrique ; Illustration 5). Ces deux cours d’eau sont donc alimentés en grande partie par l’apport des eaux souterraines, notamment en période d’étiage (règle générale en Haute-Normandie).

L’Illustration 5 présente également la délimitation du bassin versant souterrain de l’Austreberthe, tracé à partir de la carte piézométrique départementale (BRGM, 1990). Aucun échange souterrain (via le réseau karstique notamment) n’étant reconnu avec les bassins versants souterrains, le bassin versant hydrogéologique reste très proche du bassin versant hydrographique (délimitation plus lissée pour une superficie de 203 km²).

16 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 5 – Carte piézométrique sur le bassin versant de l’Austreberthe à l’échelle du 1 / 100 000 (BRGM, 1990)

2.3.2. L’exploitation de la ressource

La ressource en eau souterraine concernait 99 % des prélèvements sur le bassin versant de l’Austreberthe en 2005 (Illustration 6). Aujourd’hui, il ne subsiste qu’un seul captage en rivière : celui de la société Gailliard à Barentin.

BRGM/RP-55813-FR – Rapport final 17 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Après une décroissance des prélèvements entre 1994 et 1998, on constate une nette augmentation en 1999/2000, à attribuer aux prélèvements industriels. En 2005, les prélèvements destinés à l’alimentation en eau potable représentaient de l’ordre de 60 % des prélèvements du bassin. Ils se font par l’intermédiaire de 3 forages :

• le forage de Limesy (indice BSS n° 00767X0021) ;

• le forage de Duclair (indice BSS n° 00992X0197) ;

• le forage de Sierville (indice BSS n° 00767X0020).

Certains des ouvrages industriels exploitent la nappe alluviale (ex. de la société BADIN).

Enfin, aucun prélèvement agricole n’est connu par l’AESN dans l’emprise du bassin versant de l’Austreberthe.

3500000

3000000

2500000

2000000

1500000

Volume annuel (m3) 1000000 Prélèvements totaux Prélèvements en nappe 500000 Prélèvements AEP Prélèvements industriels 0 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 Année

Illustration 6 – Evolution des prélèvements en eau sur le bassin versant de l’Austreberthe entre 1994 et 2005 (source AESN)

18 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3. Modélisations GARDENIA

3.1. LE MODELE HYDROLOGIQUE GLOBAL GARDENIA DU BRGM

3.1.1. Principe du modèle

Le modèle hydrologique global GARDENIA (Thiery, 2003) est un modèle hydrologique global développé par le BRGM. Ce logiciel simule le cycle de l’eau (Illustration 7), depuis les précipitations sur le bassin versant jusqu’au débit à l’exutoire d’une rivière (ou d’une source) ou au niveau ponctuel d’un aquifère (niveau piézométrique).

Illustration 7 – Représentation schématique des éléments du bilan hydrologique d’un bassin versant

Ce modèle est global, car il considère des « données d’entrée » globales en pluie et évapotranspiration potentielle et une sortie unique qui est, suivant le cas, le débit à l’exutoire ou le niveau piézométrique en un point de la nappe sous-jacente. GARDENIA simule le cycle de l’eau par plusieurs réservoirs en cascade.

BRGM/RP-55813-FR – Rapport final 19 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.1.2. Fonctionnement du modèle

Plusieurs schémas hydrauliques peuvent être utilisés : • le schéma type GARDENIA au sens strict avec son module de production propre. Plusieurs configurations sont possibles avec notamment 1 (Illustration 8) ou 2 réservoirs souterrains. Dans la présente étude, nous avons utilisé la configuration avec un seul réservoir. Dans la configuration à 1 réservoir souterrain, les 3 réservoirs représentent schématiquement : - la zone superficielle du sol sujette à la reprise par évaporation ; - la zone non saturée siège des écoulements rapides (Ruissellement et composante rapide des écoulements karstiques) ; - la zone saturée siège des écoulements lents de l’aquifère ; • le schéma mixte type GARDENIA avec le module de production type GR4 ; • le schéma type GR4 développé par Edijatno et Michel (1989).

Les échanges entre les réservoirs sont définis par des fonctions de transfert faisant intervenir une dizaine de paramètres globaux (réserve utile, temps de tarissement, etc.) définis pour un bassin versant ou une entité homogène au sein d’un bassin versant.

Pour initialiser l’état de remplissage des réservoirs, il est nécessaire, en modélisation hydrologique globale, de prendre en compte une période d’initialisation. Celle-ci doit être d’autant plus longue que l’inertie ou « mémoire » de l’aquifère est grande. Cela suppose de disposer de chroniques de pluie et d'ETP sur une période aussi longue que possible antérieurement à la période de calage.

Le calage du modèle consiste à ajuster au mieux les données d’entrée (pluie, ETP) et les données de sortie (débit ou niveau piézométrique) avec un jeu de paramètres aussi réaliste que possible, calé sur l'exploitation de l'ensemble des données disponibles. Ce calage s’effectue par une méthode semi-automatique. L’utilisateur fournit un jeu de paramètres initiaux, donne des bornes de variations plausibles pour chacun d'eux et indique ceux qui doivent être optimisés. L'optimisation concomitante des différents paramètres est effectuée selon un algorithme non-linéaire adapté de la méthode de Rosenbrock. Le logiciel teste des valeurs de ces paramètres jusqu’à trouver un jeu fournissant la meilleure adéquation entre séries calculées et observées. Ce calage s’effectue sous le contrôle de l’utilisateur qui peut le contraindre, notamment par la fixation des valeurs de certains paramètres ou de bornes plausibles.

Une fois calé, le modèle peut être utilisé pour : • Etablir un bilan hydrologique, et/ou vérifier la cohérence entre les différentes données (pluie, ETP, débit et/ou niveau), en permettant notamment, dans une certaine mesure, de différencier les deux composantes de l’écoulement : - celle rapide que l’on peut assimiler au ruissellement superficiel et à la composante rapide des écoulements karstiques ; - celle lente que l’on peut assimiler à la composante lente des écoulements souterrains ;

20 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

• Simuler des débits ou des niveaux résultant d'épisodes climatiques exceptionnels. Ces simulations permettent le dimensionnement d'ouvrages tels que retenues pour l'irrigation, aménagements de sources, fondations de bâtiments, niveaux de sous-sol, dimensionnement d'ouvrages hydrauliques pour l'évacuation des crues ou pour la régulation des débits (barrages), etc ; • Faire une extension de débits ou de niveaux piézométriques sur une période pendant laquelle on possède les données d’entrée mais pas de débit et/ou de niveau piézométrique ; • Réaliser une prévision de débits ou de niveaux piézométriques. Il est possible de prolonger une série observée pour les semaines, mois ou années à venir en considérant différents scénarii possibles (étiage, niveau piézométrique de hautes eaux, changement climatique, ……) ; • Reconstituer des valeurs manquantes dans une chronique de débits (rivière ou source), ou de niveau piézométrique (nappe). Les séries ainsi complétées peuvent être utilisées dans d'autres calculs ou modélisations.

Pluie - Niveau piézométrique BASSIN DE PROVINS 150

Charge 142

134

126

118

110 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002 Date GARDENIA R= 0.968 ------=Observé - - - =Simulé Pluie - Débit

…. Ecoulement souterrain lent

Illustration 8 – Fonctionnement du modèle hydrologique GARDENIA pour la simulation du débit d'un cours d'eau ou d'un niveau piézométrique (Schéma classique)

BRGM/RP-55813-FR – Rapport final 21 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.1.3. Equations du bilan hydrologique

L’exploitation des modèles hydrologiques globaux GARDENIA pluie – débit et pluie – niveau piézométrique permet de quantifier les différents termes du bilan hydrologique.

Les équations du bilan hydrologique global pluie – débit sont les suivantes :

PLUIE = E TR + PEF

PEF = QRAP + QSOUT + dSTO

QRivière = ()QRAP + QSOUT + QEchange avec :

PLUIE Lame d’eau précipitée sur le bassin

ETR Lame d’eau évapotranspirée sur le bassin

PEF Hauteur de pluie efficace ou écoulement global

QRAP Lame d’eau écoulée rapide (Ruissellement superficiel + écoulement karstique rapide)

QSOUT Lame d’eau écoulée correspondant à l’écoulement souterrain lent

QEchange Lame d’eau écoulée correspondant au débit d’échange avec d’autres systèmes hydrologiques

Q > 0 si apport extérieur au bassin Echange

QEchange < 0 si sortie d'eau bassin

Dans les modélisations hydrologiques globales pluie – niveau piézométrique le niveau piézométrique est déduit de la hauteur d’eau dans le réservoir souterrain (G) par la relation :

G NP = + NB SG avec :

NP Niveau piézométrique

G Hauteur d’eau dans le réservoir souterrain

SG Coefficient d’emmagasinement global

NB Niveau de base (niveau qui serait atteint en l’absence totale d’alimentation au bout d’un temps infini)

22 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

On peut ainsi différencier deux types d’écoulements :

o les écoulements rapides, c'est-à-dire le ruissellement superficiel et la composante rapide des écoulements karstiques ;

o les écoulements lents constitués principalement des écoulements souterrains. C’est ce que nous assimilerons à la contribution des aquifères au débit des rivières (aux écoulements karstiques près) .

Ces termes « lent » et « rapide » doivent être relativisés en fonction notamment de la taille du bassin versant et du pas de temps de calcul.

Les termes du bilan hydrologique peuvent être calculés à des pas de temps différents.

La signification et l’utilisation des différents paramètres de GARDENIA sont explicités dans le tableau ci-dessous :

Paramètres Unité Dénomination BV km² Surface du bassin versant Débit d'échange m3/s Débit d’échange avec d’autres systèmes hydrologiques Q > 0 si apport extérieur au bassin Echange

QEchange < 0 si sortie d'eau bassin Coef. PLUIE % Coefficient correcteur de la pluie Coef. ETP % Coefficient correcteur de l'ETP RUMAX mm Capacité de rétention maximale du sol RUMAX_PROG mm Capacité de la réserve superficielle Progressive, cette alternative à la RUMAX permet de rendre compte des variations de l’ETP en fonction du taux d’humidité du sol. Dans notre cas, de meilleurs calages ont été obtenus avec la prise en compte de ce paramètre. RUIPER mm Hauteur d'égal ruissellement percolation THG Pas de temps Temps de demi-percolation TG1 Pas de temps Temps de demi tarissement rapide du réservoir souterrain 1 TG1_2 Pas de temps Temps de demi-transfert entre les réservoirs souterrains 1 et 2 TG2 Pas de temps Temps de demi-tarissement lent du réservoir souterrain 2 Coef. Ajustement Sans unité Coefficient d'ajustement Pondération Sans unité Pondération des débits (0=non, 1=oui)

BRGM/RP-55813-FR – Rapport final 23 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.2. APPLICATION AU BASSIN VERSANT DE L’AUSTREBERTHE

L’Illustration 9 permet de localiser les différentes stations de mesures utilisées dans le cadre de la présente étude.

Illustration 9 – Localisation des différentes stations de mesures disponibles sur le bassin versant

24 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.2.1. Les données d’entrée

Il est rappelé que les données d’entrée du modèle hydrologique global GARDENIA sont, pour cette étude, les données de précipitations et d’ETP. Les données climatologiques utilisées ont été récupérées auprès de Météo France pour la station de Rouen-Boos sur la plus longue période possible, à savoir du 01/03/1968 au 30/06/2007. Cette station enregistre la mesure des précipitations et des paramètres climatologiques permettant le calcul des valeurs de l’ETP. Le calcul de l’ETP, réalisé par Météo France, repose sur l’utilisation de la méthode de Penman-Monteih, au pas de temps journalier. Il s’agit du poste climatique avec mesure de l’ETP le plus proche du bassin étudié.

L’Illustration 10 représente les cumuls annuels (année civile) des précipitations et de l’ETP obtenus à la station de Rouen Boos sur la période 1969/2006. Pour ce qui est des précipitations annuelles, on enregistre une valeur minimale de 482 mm en 1976, une valeur maximale de 1 053 mm en 2000 pour une moyenne inter-annuelle de 821 mm. Sur les dernières décennies, ce graphique montre :

• des années particulièrement sèches en 1971, 1976 et 1989 avec des précipitations annuelles inférieures à 600 mm ;

• et des années très humides en 1981, 1993, 1999, 2000 et 2001 avec des précipitations annuelles supérieures à 1 000 mm.

1200 Précipitations ETP 1000

800

600 Hauteur(mm) 400

200

0 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 Année

Illustration 10 – Répartition annuelle des précipitations et de l’ETP sur la station climatique de Rouen-Boos (source : Météo France)

BRGM/RP-55813-FR – Rapport final 25 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Quelques rares lacunes d’information ont été constatées au niveau des mesures quotidiennes de l’ETP, la série a été reconstituée à partir de la moyenne inter-annuelle du pas de temps (aucune lacune pour les précipitations).

Dans le cadre de la gestion d’une ressource en eau souterraine, il est nécessaire de différencier les précipitations d’hiver (entre octobre et avril) des précipitations d’été (entre mai et septembre) : les précipitations d’hiver donnent lieu à une précipitation efficace qui alimente les nappes. Les précipitations d’été ne produiront pas de précipitations efficaces (sauf évènement pluvieux exceptionnels) mais seront utiles pour l’agriculture dont les besoins sont alors maximums.

Les modélisations pluie-débit et pluie-niveau permettront ultérieurement d’évaluer mois par mois la pluie efficace sur le bassin versant.

Néanmoins, en première approximation, on peut se référer à la carte de France des précipitations efficaces annuelles moyennes élaborée par le MEDAD sur la période 1946-2001 (Illustration 11). Sur le bassin versant de l’Austreberthe, la lame d’eau moyenne annuelle des précipitations efficaces se répartit en 2 grands domaines : • de 300 à 350 mm sur la partie aval du bassin ; • de 350 à 400 mm sur la partie amont.

Illustration 11 – Répartition de la pluie efficace annuelle moyenne (1946-2001) sur le bassin versant de l’Austreberthe (source : MEDAD)

26 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Il apparaît que les précipitations sont significativement plus importantes sur l’ensemble du bassin versant de l’Austreberthe qu’au droit de la commune de Boos (station climatique utilisée).

Les précipitations mensuelles moyennes des postes d’Auzebosc (implanté à l’ouest du bassin et plus proche de la zone d’étude) et de Boos ont été comparées par la méthode des doubles cumuls (Illustration 12). Cette dernière méthode, propre à l’hydrologie, compare plusieurs échantillons acquis à différentes stations afin de déceler une éventuelle hétérogénéité. Le principe de la méthode consiste à vérifier la proportionnalité de la somme des valeurs mesurées aux deux stations.

950 y = 0.8914x + 3.8944 R2 = 0.9963 850

750

650

550

(mm) 450

350

250

150

Cumul des précipitations mensuellesmoyennesBoos à desprécipitations Cumul 50 50.0 250.0 450.0 650.0 850.0 1050.0 Cumul des précipitations mensuelles moyennes à Auzebosc (mm)

Illustration 12 – Graphe des double cumuls pour les stations climatiques d’Auzebosc et de Boos (source : Météo France)

L’utilisation de la station de Boos reste justifiée compte tenu des bonnes corrélations obtenues avec la station d’Auzebosc (Illustration 12). Comme annoncé précédemment, un coefficient de correction global compris entre 10 et 15% devra, cependant, être affecté aux valeurs de pluie (coefficient optimisé par le modèle).

3.2.2. Les données de sortie

Données de débits :

Les valeurs de débits moyens journaliers de l’Austreberthe sont disponibles au droit de deux stations hydrométriques de la DIREN Haute-Normandie : une station à Duclair et une à Saint-Paër. Les fiches signalétiques de ces stations sont consultables sur internet ( http://www.hydro.eaufrance.fr ).

BRGM/RP-55813-FR – Rapport final 27 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Superficie BV Station Début suivi Fin de suivi (km²)

Duclair 01/06/1964 04/02/1999 209

Saint-Paër 27/02/1997 En cours 157 Illustration 13 – Caractéritiques des 2 stations de jaugeage suivies par la DIREN Haute- Normandie

Une analyse critique de ces données a été effectuée par l’intermédiaire d’études de corrélation et des doubles cumuls sur la période commune d’observation (du 27/02/1997 au 04/02/1999).

Le graphique de corrélation des lames d’eau ruisselées montre un nuage élargi par suite d’une série de valeurs. Le graphique des doubles cumuls montre, lui, une cassure entre le 17/06/1997 et le 06/01/1998 (Illustrations 13 et 14).

La DIREN a pu nous fournir l’origine de cette discordance entre les 2 séries de données : la crue du 16 juin 1997 a provoqué la rupture du pont situé à l’amont immédiat de la station de Duclair et donc une modification importante de l’écoulement local de l’Austreberthe (détarage de plus de 10 cm). Au-delà de cet incident, la DIREN accorde plus de fiabilité à la station de Saint-Paër, en particulier en période de crue où le pont était fréquemment submergé à Duclair.

En définitive, compte tenu des observations faites précédemment et des précisions apportées par la DIREN, il a été décidé de réaliser la modélisation pluie-débit à partir des données issues de la station de jaugeage de Saint-Paër.

28 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.5

3.0

2.5

2.0

(mm) 1.5

1.0

0.5

Lame d'eau journalières ruisselées Lamejournalières ruisselées d'eau à Saint-Paër 0.0 0.0 0.5 1.0 1.5 2.0 2.5 Lame d'eau journalières ruisselées à Duclair (mm)

Illustration 14 – Graphique de corrélation entre les lames d’eau journalières ruisselées à Saint- Paër et à Duclair

600

500

400

06/01/1998 300

Paër (mm) Paër 17/06/1997 200

100

Cumul à Saint- lames d'eau des ruisselées 0 0 100 200 300 400 500 600 Cumul des lames d'eau ruisselées à Duclair (mm)

Illustration 15 – Graphique de double cumul pour les lames d’eau ruisselées à Saint-Paër et à Duclair

BRGM/RP-55813-FR – Rapport final 29 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Les débits caractéristiques de l’Austreberthe à Saint-Paër sont reportés dans le tableau suivant :

Définition Valeur (m3/s)

Débit journalier moyen 1,88

Débit journalier minimal 0,95 (03/10/1997)

Débit journalier maximal 11.2 (26/12/99)

Débit instantané maximal 18.7 (11/05/00) Illustration 16 – Débits caractéristiques de l’Austreberthe à la station hydrologique de Saint-Paër (source : DIREN Haute-Normandie)

janv.-01 avr.-01 juil.-01 oct.-01 févr.-02 mai-02 août-02 déc.-02 mars-03 juin-03 sept.-03 janv.-04 avr.-04 juil.-04 nov.-04 févr.-05 mai-05 août-05 déc.-05 0 17000

10 15000

13000 20

11000 30 9000 40 7000

50 5000 Hauteur d'eau précipitée d'eau (mm) précipitée Hauteur 60 3000 Débit de l'Austreberthe à Saint-Paër (l/s) Débit de l'Austreberthe 70 1000

Illustration 17 – Hauteur d’eau précipitée à Rouen et débits de l’Austreberthe en 2001 (au pas de temps journalier)

Comme conclu par la DIREN Haute-Normandie lors de l’étude de la vallée sèche du Paulu (bassin versant de l’Austreberthe), l’Illustration 16 montre que la précipitation journalière n’est pas un critère suffisant pour expliquer le ruissellement sur ce bassin versant. Le taux de saturation du sol semble, en effet, jouer un rôle déterminant dans

30 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

le déclenchement de crues importantes. Ce graphique permet également de constater que le débit observé est composé d’un débit de base (écoulement souterrain) très prononcé et de débits de crue de durée très courte. Ces derniers correspondent au ruissellement superficiel se produisant lors d’épisodes pluvieux importants, notamment sur les zones imperméabilisées. Les précipitations ponctuellement importantes en 2004 et 2005 n’ont engendré aucune crue significative, ceci peut s’expliquer par une non atteinte de la saturation des sols du bassin versant lors de cette succession d’hivers secs. Dans ces périodes d’étiage sévère, l’Austreberthe ne semble être alimentée que par la vidange des aquifères, hormis quelques petits pics engendrés par du ruissellement au niveau des zones urbanisées.

Données de niveaux piézométriques :

Le seul point de suivi du réseau piézométrique du bassin Seine-Normandie implanté dans le bassin versant souterrain de l’Ausreberthe est l’ouvrage de Blacqueville (en limite ouest, Illustration 9). Ce puits de 85 m de profondeur capte l’aquifère crayeux sous plateau.

Les fluctuations piézométriques (Illustration 17) enregistrées au droit de Blacqueville sont de type pluri-annuel, c’est-à-dire caractérisées par des fluctuations saisonnières très atténuées (voire inexistantes) et des fluctuations pluriannuelles de grande amplitude. Un tel régime de fluctuation est, localement, caractéristique d’un aquifère à fonction capacitive (peu transmissif), d’un recouvrement argileux et limoneux important et d’exutoires éloignés. Il est, en effet, rappelé que le puits de Blacqueville est localisé en sommet de plateau.

Illustration 18 – Chronique piézométrique sur Blacqueville de 1968 à 2006

BRGM/RP-55813-FR – Rapport final 31 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.2.3. Le calage du modèle

Modélisation pluie-débit :

La série de débits de la station de Saint-Paër a été calée sur la période 1997-2006 (station hydrologique mise en service en 1997). A partir des données climatiques disponibles, une période d’initialisation de 28 ans a pu être retenue, ce qui est intéressant compte tenu de l’inertie importante du système. Cette longue période d’initialisation permet de créer des conditions initiales de simulation très voisines des conditions réelles.

Différentes hypothèses ont été retenues pour réaliser le calage :

• la superficie du bassin versant a été prise égale à celle du bassin versant topographique contrôlée par la station de jaugeage de Saint- Paër. Celle-ci, déterminée à partir de la carte topographique au 1 / 25 000 de l’IGN, est de 157 km². Cette valeur a bien entendu été imposée sans possibilité d’optimisation ;

• comme vu précédemment, les précipitations de Boos ont été laissés libres à l’optimisation pour permettre l’application d’un coefficient multiplicateur (+ 12%). Les valeurs d’ETP n’ont en revanche pas été affectées d’un coefficient de correction ;

• les paramètres laissés libres à l’optimisation sont les différents paramètres de calage standard du modèle GARDENIA c'est-à-dire : la capacité de la hauteur d’égal ruissellement percolation (mm), la réserve superficielle progressive (mm), le temps de percolation (mois), le temps de demi-tarissement (mois) et le débit de fuite. L’optimisation est cependant contrainte pour obtenir un jeu de données cohérent et réaliste.

Le meilleur calage a été obtenu pour un pas de temps décadaire. Plusieurs essais de calage ont été menés au pas de temps journalier, mais l’ajustement s’est avéré moins bon, notamment en raison de la difficulté du modèle à reproduire les épisodes de forte crue. Vu la très courte durée des écoulements rapides, il aurait été nécessaire de fournir des données d’entrée avec un pas de temps encore inférieur à la journée pour espérer les reproduire de manière satisfaisante. L’objectif principal de la présente modélisation étant d’estimer la contribution des écoulements lents au débit de l’Austreberthe, un calcul au pas de temps décadaire s’est avéré suffisant.

Le coefficient d’ajustement ainsi obtenu est de 0.958 (Illustration 19). Les paramètres de calage sont reportés en Annexe 1.

32 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Illustration 19 – Visualisation sur la période 1997-2006 du calage pluie-débit de l’Austreberthe à Saint-Paër sur la période 1997-2006 (en bleu : la part des écoulements lents)

Le graphique résultat de l’Illustration 19 montre que les grandes tendances des variations de débits de l’Austreberthe à Saint-Paër sont bien reconstituées par le modèle. Les augmentations progressives des débits sont bien reproduites pendant l’hiver et les corrélations des vitesses de tarissement sont bonnes.

En revanche, d’une manière générale, la simulation rencontre des difficultés pour représenter :

• les évènements de crue soudains liés à des précipitations hivernales importantes (octobre-novembre 1998, décembre 1999, mai et novembre 2000, janvier et mars 2001). Les débits calculés sont systématiquement inférieurs à ceux observés. Par suite du pas de temps des données d’entrée et de la grande différence de temps de réponse entre les deux composantes de l’écoulement (rapide et lente), le modèle ne parvient pas à reproduire tout à fait les ruissellements importants qui surviennent suite à de fortes précipitations en période de saturation des sols ;

• l’étiage observé depuis 2004 : le modèle réagit, en effet, aux précipitations hivernales alors qu’en réalité le cours d’eau semble avoir atteint son débit de base (globalement stationnaire sur 2004 et 2005). Ceci traduit la quasi-absence de précipitations efficaces sur ces deux années.

BRGM/RP-55813-FR – Rapport final 33 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

A partir des paramètres retenus à l’issue du calage, la chronique a été reconstituée sur la période 1979-1996 (Illustration 20). Cette reconstitution n’est qu’indicative, car elle est basée sur le niveau d’anthropisation actuel du bassin versant.

Illustration 20 – Reconstitution des débits de l’Austreberthe à Saint-Paër depuis 1979 (en bleu : la part des écoulements lents)

Finalement, le calage pluie-débit sur l’Austreberthe à Saint-Paër a permis d’identifier les données non mesurées telles que la pluie efficace et le débit de fuite. On obtient sur la période de simulation (1997-2006) un débit de fuite de 26 L/s (négligeable) et une pluie efficace annuelle moyenne de 377 mm .

Même si les périodes d’observation ne sont pas les mêmes, les ordres de grandeur de la pluie efficace annuelle moyenne calculée par GARDENIA (377 mm) et par le MEDAD (environ 360 mm, cf. Illustration 11) sont identiques.

Le débit de fuite, ajusté automatiquement par le modèle, est supposé constant tout au long de l’année et correspond au flux non contrôlé à l’exutoire (débit souterrain + prélèvements en nappe). Dans le cas présent, le débit de fuite calculé peut être considéré comme négligeable, car il correspond à une lame d’eau annuelle de 5.2 mm ce qui est, en effet, très faible comparé aux incertitudes liées à la pluie moyenne et à l’ETP sur le bassin.

34 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Le volume annuel des prélèvements pour l’irrigation, l’alimentation en eau potable et l’industrie sur le bassin versant de l’Austreberthe en amont de Saint-Paër qui s’élevait à 1.78 Mm3 en 2005, correspond à 57 L/s soit 3 % du débit moyen. Une importante partie de ces prélèvements retournant au milieu, il n’est pas surprenant que le débit de fuite calculé soit inférieur aux prélèvements, en particulier pour de tels ordres de grandeur.

Modélisation pluie-niveau piézométrique :

La chronique piézométrique du puits de Blacqueville a été calée sur la période 1980- 2006 (station piézométrique mise en service en 1968). A partir des données climatiques disponibles, une période d’initialisation de 11 ans a ainsi pu être retenue.

Comme dans le cas de la modélisation pluie-débit, les précipitations à Boos ont été affectées d’un coefficient correcteur (+ 12%). Dans le même temps, les paramètres laissés libres à l’optimisation sont la capacité de la hauteur d’égal ruissellement- percolation (mm), le temps de percolation (mois), le temps de demi-tarissement (mois) et le coefficient d’emmagasinement (optimisation contrainte).

Le calage a également été réalisé au pas de temps décadaire.

Le coefficient d’ajustement finalement obtenu est de 0.931 (Illustration 21). Le compte-rendu complet de la modélisation est fourni en Annexe 2.

Sur la base de ce calage, la pluie efficace annuelle moyenne calculée par le modèle sur la période 1980-2006 est de 400 mm .

Le coefficient d’emmagasinement a été laissé libre à l’optimisation. Sa valeur calculée de 0.0097 semble supérieure aux valeurs connues pour l’aquifère crayeux en sommet de plateau (0.0001 à 0.001). Attention, le coefficient d’emmagasinement GARDENIA, correspondant plus à un coefficient d’influence globale des fluctuations de stock sur le niveau piézométrique à Blacqueville, ne peut être relié à une interprétation de pompage d’essai dont la validité reste locale et très souvent réalisé sur une courte période. Le coefficient d’emmagasinement GARDENIA peut éventuellement être égal au coefficient d’emmagasinement moyen de la nappe si le point d’observation est situé loin de tout cours d’eau comme dans le cas présent. En tout état de cause, le coefficient d’emmagasinement global calculé par GARDENIA apparaît réaliste.

En ce qui concerne la chronique piézométrique simulée, on constate, une nouvelle fois, que les grandes tendances sont bien reconstituées. Les recharges et les vidanges pluri-annuelles sont notamment très bien retranscrites.

En revanche, d’une manière générale, la simulation rencontre des difficultés pour représenter :

• les années de très hautes eaux où la piézométrie calculée reste toujours inférieure à celle observée (1995 et 2001 par exemple). Ceci peut s’expliquer par la forte hétérogénéité de l’aquifère crayeux dans sa profondeur (coefficient d’emmagasinement plus important dans sa partie supérieure). Pour tenter de

BRGM/RP-55813-FR – Rapport final 35 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

tenir compte de cette hétérogénéité de l’aquifère, une modélisation avec deux réservoirs souterrains a été menée mais sans amélioration significative du calage ;

• certaines variations saisonnières (fluctuations annuelles de 1999 et 2000 par exemple). La modélisation semble donner plus de poids à la fonction capacitive de l’aquifère, déjà très importante en réalité.

Illustration 21 – Visualisation sur la période 1980-2006 du calage pluie-niveau à Blacqueville sur la période 1980-2006

36 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.2.4. Bilan hydrologique moyen

Avant toute analyse des bilans hydrologiques calculés par la modélisation globale, il est rappelé que les résultats sont à prendre avec précaution : il s’agit d’une approche tendancielle et non d’une approche quantitative.

Les bilans hydrologiques annuels moyens, obtenus pour chacune des modélisations, sont sensiblement équivalents si on les compare sur la période 1997-2006 (Illustration 22). Type deType modélisation Station hydrométrique / Piézomètre Pluie(mm) (mm)ETR Pluieefficace (mm) (mm) Qrap. Qsout.(mm) Stock Diff. (mm) Qsout. %

Pluie-Débit Saint-Paër 972 598 377 36 339 -2 90

Pluie-Niveau Blacqueville 972 589 386 45 345 +4 88.5

Illustration 22 – Bilans hydrologiques annuels moyens calculés pour chacune des modélisations sur la période 1997-2006

Sur cette période, la contribution moyenne des écoulements lents au débit de l’Austreberthe est de l’ordre de 90 % pour les deux modélisations. Ce résultat est en accord avec les travaux menés par l’Université de Rouen (Dupuis).

L’Illustration 23 présente les différents termes du bilan hydrologique sur le bassin versant. Il s’agit d’un bilan annuel calculé à partir des données moyennes résultant du calage pluie-débit à Saint-Paër (période de simulation 1997-2006).

Illustration 23 – Représentation schématique du bilan hydrologique à Saint-Paër

BRGM/RP-55813-FR – Rapport final 37 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

3.2.5. Contribution des écoulements lents

L’étude des variations saisonnières du bilan hydrologique a été réalisée à partir des résultats de la modélisation pluie-débit, celle-ci est, en effet, plus représentative du fonctionnement hydrologique du bassin versant (calcul plus intégrateur).

Le modèle permet de décomposer l’écoulement à la rivière en écoulement rapide (ruissellement rapide + écoulement karstique rapide) et en écoulement lent souterrain (contribution des aquifères au débit de la rivière). Ces résultats doivent être pris avec précautions, car il n’y a pas de certitude dans le fait d’assimiler les écoulements identifiés comme lents aux écoulements souterrains. Cette décomposition permet, cependant, de dégager des informations tendancielles sur l’alimentation du cours d’eau par la nappe.

L’Illustration 24 représente les variations saisonnières de la contribution des écoulements lents au débit de l’Austreberthe sur la période 1969-2006.

100

80 Minimum Moyenne Maximum 70

60 lents au débit de l'Austreberthe (%) lents au débit de l'Austreberthe Contribution mensuelle écoulements des 50 Jan. Fév. Mar. Avr. Mai Juin Juil. Août Sept. Oct. Nov. Déc. Mois

Jan. Fév. Mar. Avr. Mai Juin Juil. Août Sept. Oct. Nov. Déc. Min 73.7 65.6 69.9 72.0 78.4 82.2 78.0 88.3 86.2 66.3 59.5 62.4 Max 98.2 96.3 98.5 99.5 100.0 100.0 100.0 100.0 100.0 100.0 100.0 99.3 Moyenne 86.4 87.1 88.9 91.1 92.2 94.5 96.2 97.4 97.7 95.4 92.0 87.6 Ecart-Type 7.9 7.9 7.2 6.1 5.0 4.4 4.0 2.6 2.6 6.7 8.5 9.8

Illustration 24 – Contributions des écoulements souterrains lents au débit de l’Austreberthe au pas de temps mensuel (1969-2006)

38 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

A titre indicatif, plusieurs constats peuvent être faits :

• sur la période 1969-2006, la modélisation GARDENIA calcule une contribution mensuelle minimale des écoulements souterrains lents au débit de l’Austreberthe de 59 %, une contribution maximale de 100 % et une contribution moyenne de 90 % ;

• en période d’étiage (juillet à octobre), le débit de l’Austreberthe est essentiellement assuré par les écoulements souterrains : de 95 à 98 % en moyenne d’après la modélisation globale. L’analyse statistique montre que les écarts interannuels sont les plus faibles pour ces périodes : quelque soient les précipitations, la part du ruissellement reste, en effet, très faible (absence de pluies efficaces) ;

• en période de recharge hivernale, la contribution des écoulements souterrains au débit de l’Austreberthe est logiquement plus faible mais reste largement prépondérante : de 86 à 89 % en moyenne. Les écarts interannuels sont, en revanche, plus importants : la part du ruissellement est, en effet, fortement dépendante des précipitations hivernales. La part du ruissellement peut ainsi, par exemple, varier entre 4 % en février 1976 (hiver sec) et 34 % en février 1995 (hiver exceptionnellement pluvieux) ;

• des ruissellements significatifs ne sont observés que lors des hivers exceptionnellement pluvieux avec en général plusieurs mois très humides consécutifs. A titre indicatif, on peut noter l’hiver 1974-1975 avec une contribution mensuelle minimale des écoulements souterrains lents de 59 %, l’hiver 1979- 1980 (62 %), 1992-1993 (76 %), 1995 (66 %) ou encore 2000-2001 (70 %).

Mise en situation des campagnes de prélèvements

Dans le cadre de l’approche isotopique (Rapport BRGM/RP-55983-FR), deux campagnes de prélèvements ont été menées : du 9 au 10 novembre 2006 pour représenter une situation de basses eaux et du 9 au 11 mai 2007 pour une situation de hautes eaux (en fin de période de recharge).

Le graphique de l’Illustration 25 présente les évolutions des précipitations à Boos, des niveaux piézométriques de l’aquifère crayeux à Blacqueville et des débits journaliers moyens de l’Austreberthe à Saint-Paër pour le cycle hydrologique en cours (de septembre 2006 à juin 2007).

En novembre 2006, la situation hydrogéologique de l’aquifère crayeux correspondait à un état de basses eaux marqué, et la recharge hivernale n’avait toujours pas débuté sous les plateaux. Le piézomètre de Blacqueville indiquait alors un taux de remplissage de l’aquifère de 28 % pour un niveau piézométrique de 51 m NGF. Les précipitations cumulées sur le mois d’octobre 2006 à Boos ont été de 55 mm (moyenne mensuelle interannuelle de 75 mm) et de 6 mm sur la première décade de novembre.

BRGM/RP-55813-FR – Rapport final 39 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Ce déficit pluviométrique n’a pas permis de générer de pluies efficaces (ruissellement + infiltration) significatives en fin d’automne. Une légère augmentation du débit de l’Austreberthe à Saint-Paër n’est, en effet, observée qu’à partir de la deuxième décade de novembre. Cette augmentation de débit s’amplifie ensuite sur la troisième décade de novembre et surtout sur la première décade de décembre.

3500 59 Pluie à Boos

Débit à Saint-Paër 3000 57

Niveau piézométrique à Blacqueville

2500 55

2000 53

1500 51

1000 49 Débit(10² mm) (L/s) Pluie / Niveau piézométrique Niveau (mNGF)

500 47

0 45 sept.-06 oct.-06 nov.-06 déc.-06 janv.-07 févr.-07 mars-07 avr.-07 mai-07 juin-07 Date

Illustration 25 – Evolution des paramètres hydro-climatiques pour le cycle hydrologique 2006- 2007

Lors de la première campagne de prélèvement, l’étiage était ainsi marqué pour la nappe comme pour le cours d’eau. Le modèle calcule une contribution des écoulements souterrains lents de 98 % pour le mois de novembre (part des ruissellements négligeable).

Début mai 2007, le piézomètre de Blacqueville indiquait un taux de remplissage de l’aquifère de 36 % avec un niveau de l’ordre de 56 m NGF (+ 5 m par rapport à novembre 2006). Sur le cycle hydrologique 2006-2007, la recharge de l’aquifère crayeux s’est étalée de décembre 2006 à juillet 2007 (Illustration 25). En revanche, en l’absence de précipitations sur le mois d’avril (excepté le 29 avril), l‘Austreberthe avait déjà entamé sa décrue estivale. En définitive, le printemps 2007, exceptionnellement sec, a empêché tout ruissellement significatif sur le bassin de l’Austreberthe les semaines précédant la deuxième campagne de prélèvement. Le modèle calcule une contribution des écoulements souterrains lents de 96 % . L’ordre de grandeur reste donc identique à celui calculé pour novembre

40 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

2006 même si les conditions hydrologiques sont différentes (recharge de l’aquifère en cours et débit de l’Austreberthe supérieur).

Des prélèvements réalisés à la fin mars 2007 auraient sans doute permis la comparaison de situations plus contrastées (le modèle calcule une contribution des écoulements lents de 87 % en mars).

Pour chacune des campagnes de prélèvements, les conditions hydrologiques sont reportées dans le tableau de l’Illustration 26.

Débit décadaire de Piézométrie de l’aquifère l’Austreberthe à Saint-Paër crayeux à Blacqueville (L/s) Contribution Campagne de des Date prélèvements écoulements Taux de lents (%) NP (m NGF) remplissage Observé Simulé (%)

Basses eaux 9 au 10/11/06 51 28 1 155 1 052 98

Hautes eaux 9 au 11/05/07 56 36 1 779 1 286 96

Illustration 26 – Synthèse des données piézométriques et hydrologiques pour les deux campagnes de prélèvements

BRGM/RP-55813-FR – Rapport final 41

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

4. Conclusions

Des modélisations GARDENIA satisfaisantes ont pu être obtenues pour le bassin versant de l’Austreberthe. Il a, par conséquent, été possible de calculer un bilan hydrologique fiable, notamment à partir de la modélisation pluie-débit. Le calage pluie- débit est, en effet, représentatif de l’ensemble du bassin versant contrôlé par la station de jaugeage, tandis que le calage pluie-niveau reste une représentation plus ponctuelle.

Après analyse des différents termes du bilan hydrologique moyen calculé sur le bassin versant de l’Austreberthe en amont de Saint-Paër, les principaux enseignements sont les suivants :

En ce qui concerne les deux campagnes de prélèvements réalisées en novembre 2006 et en mai 2007, le modèle semble indiquer une contribution des écoulements lents à l’alimentation de l’Austreberthe sensiblement identique : respectivement 98 et 96 %.

BRGM/RP-55813-FR – Rapport final 43

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

5. Bibliographie

Bassompierre P., Martin P. et Roux J.-C. – 1967 – BRGM – Données géologiques et hydrogéologiques acquises à la date du 30/01/1967 sur le territoire de la feuille topographique au 1 / 50 000 d’Yvetot. Rapport BRGM/DSGR.67.A24.

Conrad G. – Septembre 1980 – Avis d’Hydrogéologue Agréé – Définition des périmètres de protection, données techniques sur l’ouvrage et les installations annexes (indice BRGM 00767X0021).

Eberentz P. – Novembre 1984 – BRGM – Recherche des causes de pollution de captages d’eau potable en Seine-Maritime. Identification des points de pollution par colorations. Propositions d’actions. Rapport BRGM/84 AGI 310 HNO.

De la Querière P. – Octobre 1989 – Avis d’Hydrogéologue Agréé – Révision des périmètres de protection du forage de Limesy (indice BRGM 00767X0021). Additif au rapport de G. Conrad.

Holé J.-P., De la Querière P., Pernel F. et Peckre M. – 1990 –Carte hydrogéologique du département de la Seine-Maritime.

Rico G., Peckre M. – Septembre 1990 – BRGM – Impact des rejets des stations d’épuration sur les eaux souterraines captées pour l’AEP. Département de la Seine-Maritime. Rapport BRGM/R31417 HNO 4S 90.

De la Querière P. – Juillet 1994 – Avis d’Hydrogéologue Agréé – Définition des périmètres de protection du forage de Limesy (indice BRGM 00767X0021). Actualisation.

De la Querière P. – Octobre 1997 – Avis d’Hydrogéologue Agréé – Définition des périmètres de protection du forage de Limesy (indice BRGM 00767X0021). Actualisation.

Clermonte J. – Juillet 2001 – Avis d’Hydrogéologue Agréé – Définition des périmètres de protection du forage de Duclair (indice BRGM 00992X0197).

Equilbey E. – Juin 2003 – BRGM – Pilote de l’inventaire historique régional Haute-Normandie des bétoires, itinéraires souterrains des eaux (traçages) et des exutoires. Rapport BRGM/RP- 52423-FR.

BURGEAP – Février 2004 – Bassin versant de l’Austreberthe et du Saffimbec. Définition d’une méthode de diagnostic de bassin versant pour les produits phytosanitaires. Rapport Burgeap/R4263/A.11474/C.3R3581.

Clermonte J. – Avril 2005 – Avis d’Hydrogéologue Agréé – Définition des périmètres de protection du forage de Duclair (indice BRGM 00992X0197).

BRGM/RP-55813-FR – Rapport final 45

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Annexe 1 -

Compte-rendu du calage de la modélisation pluie-débit sous GARDENIA

BRGM/RP-55813-FR – Rapport final 47

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Modélisation de l'Austreberthe à Saint-Paër au pas de temps décadaire (1997-2006)

*** Pré-Options Générales *** N=Prise en compte de Neige (Déf=Non ; 1=Oui ; G=Oui + générer les param.) 0=Débit Souterrain le plus Lent (Déf=Normal ; 1=Perte ; -1=Uniquement) 0=Reprise complém. par Évapotr. (Déf=Non ; 1=Res H ; 2=Res G1 ; 3=Res G2) N=Analyse de Sensibilité (Déf=Non ; 1=Oui uniquement analyse de Sensibilité) N=Mode d'exécution (C = Contrôle sur écran ; Déf=Rapide ; D=Direct ; M=Muet) *** Options Générales *** 1=Nombre de Bassins à traiter successivement 0=Type de calcul (0=Débits , 1=Niveaux) 1=Existence d'une série de Débits ou Niveaux Observés (0=Non ; 1=Oui) 1=Édition Pluie Efficace ou Recharge (0=Non ; 1=Pluies Efficaces 2=Recharge) 1=Édition des Débit/Niveaux simulés : (0=Non ; 1=Oui) 2=Édition des termes du Bilan : (0=Non 1=Bilans annuels 2=Bilans Mensuels) 0=Allègement du Listing (0=Complet ; 1=Allégé ; 2=Supprimé) 0=Schéma de calcul (0=Gardénia ; 1=GR4 ; 2=Product GR4 ; 3=Transfert GR4) 2=Dessin de la série simulée (0=non ; 1=Oui ; 2=Oui avec décomposition) 0=Pondération des écarts pour le calage (0=non ; 2=Favorise les étiages) 0=Prise en compte de Pompage (0=non ; 1=Oui ; 2=Oui en rivière) 0=Prévision (0=Non ; 1=Préparation ; 2=émission de Prévision) 0=Méthode de Prévision (0=Méthode précise [réadaption] ; 1=Méthode simplifiée) *** Pas de temps du Fichier Pluie Neige Pompage ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier Température ************** 0= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier ETP ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier Observations ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Durée des pas de temps si non standard ************** Standard=Unité de durée des Pas si non standard (sec,min,heu,jou,moi,ann) 1.000=Durée d'un pas de temps (dans l'unité) *** Options du Bassin *** 38=Nombre d'Années des séries de données (Pluie, ETP, Observ.) [0 => Toutes] 28=Nombre d'Années démarrage (-n pour générer n année moy fictives de démarrage) 1969=Date de la Première Année des données (par ex. 2012) 0=Décalage dans la série des Pluies [+5 => Retarde de 5 pas ; -4 Avance de 4 pas] 0=Décalage de la série des Débits/Niveaux observés [ex: -2 => Avance 2 pas ] 0=État initial :0=Pl Effic. ou équil ;1=Niv ou Deb. -1=Vide ; -2=RuMax Vide aussi 300=Nombre maxi. d'itérations pour le calage (-1 = aucune itération; pas de calage)

BRGM/RP-55813-FR – Rapport final 49 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

20=Durée des pluies en moyenne par pas (%) (utilisations avancées)[défaut = 100 %] 1=Nombre de Réservoirs Souterrains (1 ou 2 ou 3=Double + seuil) [def = 1] 0=Utilisation de la même ETP (moyenne pour toutes les années [0=non ; 1=Oui] 0=Numéro du Jour initial [Def=1] (si durée non standard) ; ex. 31 pour 31 Déc. 0=Numéro du Mois initial [Def=1] (si durée non standard) ; ex. 12 pour 31 Déc. 0=Heure initiale [Def=0] (si durée non standard) ; par ex. 15 pour 15h30 0=Minute initiale [Def=0] (si durée non standard) ; Par ex. 30 pour 15h30 *** Paramètres de Prévision *** 0=Nombre d'Années de données du fichier de pluies pour la Prévision 0=Jour d'émission de la prévision (1-31) (Si pas de temps journalier sinon : 0) 0=Numéro du Mois [si journal ou mensuel] (ou du pas) d'émission de la Prévision) 0=Portée de la Prévision (Nombre de pas de la prévision) *** Paramètres Physiques *** 12.00000=Coefficient de correction globale des Pluies (%) Opti= 0 1.000E-04=Capacité de la réserve superficielle (RuMax) (mm) Opti= 0 0.00000=Coefficient de correction globale de l'ETP (%) Opti= 0 339.36166=Hauteur de répartition Ruissellement-Percolation (mm) Opti= 1 0.56317=Temps de 1/2 montée (temps de 1/2 Percolation) (mois) Opti= 1 4.84014=Temps de 1/2 tarissement Souterrain 1 (mois) Opti= 1 0.15010=Temps de 1/2 transfert Souterr. 1--> Souterr. 2 (mois) Opti= 0 0.15010=Temps de 1/2 tarissement Souterrain 2 (mois) Opti= 0 1.100E-03=Hauteur de réduction de reprise complém. par ETP (mm) Opti= 0 428.35513=Capacité de la réserve superficielle Progressive (mm) Opti= 1 -0.02610=Débit extérieur éventuel (m3/s) Opti= 1 157.00000=Surface du bassin versant (km2) Opti= 0 250.00000=Valeur démarrage : Pluie Effic annuelle/Niveau/Débit *** Bornes des paramètres physiques *** -10.00000=Min : Coefficient de correction globale des Pluies (%) Max = 12.00000 0.00000=Min : Capacité de la réserve superficielle (RuMax) (mm) Max = 500.00000 -15.00000=Min : Coefficient de correction globale de l'ETP (%) Max = 15.00000 1.00000=Min : Hauteur de répartition Ruissellement-Percolation (mm) Max =9999.00000 0.15000=Min : Temps de 1/2 montée (temps de 1/2 Percolation) (mois) Max = 10.00000 0.15000=Min : Temps de 1/2 tarissement Souterrain 1 (mois) Max = 80.00000 0.15000=Min : Temps de 1/2 transfert Souterr. 1--> Souterr. 2 (mois) Max = 50.00000 0.15000=Min : Temps de 1/2 tarissement Souterrain 2 (mois) Max = 50.00000 1.000E-03=Min : Hauteur de réduction de reprise complém. par ETP (mm) Max =9999.00000 0.00000=Min : Capacité de la réserve superficielle Progressive (mm) Max = 650.00000

50 BRGM/RP-55813-FR – Rapport final Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Annexe 2 -

Compte-rendu du calage de la modélisation pluie-niveau sous GARDENIA

BRGM/RP-55813-FR – Rapport final 51

Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

Modélisation de l'Austreberthe à Blacqueville au pas de temps décadaire (1980-2006)

*** Pré-Options Générales *** 0=Prise en compte de Neige (Déf=Non ; 1=Oui ; G=Oui + générer les param.) 0=Débit Souterrain le plus Lent (Déf=Normal ; 1=Perte ; -1=Uniquement) 0=Reprise complém. par Évapotr. (Déf=Non ; 1=Res H ; 2=Res G1 ; 3=Res G2) 0=Analyse de Sensibilité (Déf=Non ; 1=Oui uniquement analyse de Sensibilité) D=Mode d'exécution (C = Contrôle sur écran ; Déf=Rapide ; D=Direct ; M=Muet) *** Options Générales *** 1=Nombre de Bassins à traiter successivement 1=Type de calcul (0=Débits , 1=Niveaux) 1=Existence d'une série de Débits ou Niveaux Observés (0=Non ; 1=Oui) 1=Édition Pluie Efficace ou Recharge (0=Non ; 1=Pluies Efficaces 2=Recharge) 1=Édition des Débit/Niveaux simulés : (0=Non ; 1=Oui) 2=Édition des termes du Bilan : (0=Non 1=Bilans annuels 2=Bilans Mensuels) 0=Allègement du Listing (0=Complet ; 1=Allégé ; 2=Supprimé) 0=Schéma de calcul (0=Gardénia ; 1=GR4 ; 2=Product GR4 ; 3=Transfert GR4) 2=Dessin de la série simulée (0=non ; 1=Oui ; 2=Oui avec décomposition) 0=Pondération des écarts pour le calage (0=non ; 2=Favorise les étiages) 0=Prise en compte de Pompage (0=non ; 1=Oui ; 2=Oui en rivière) 0=Prévision (0=Non ; 1=Préparation ; 2=émission de Prévision) 0=Méthode de Prévision (0=Méthode précise [réadaption] ; 1=Méthode simplifiée) *** Pas de temps du Fichier Pluie Neige Pompage ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier Température ************** 0= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier ETP ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Pas de temps du Fichier Observations ************** 2= Pas de temps : 0=Journalier 1=Pentadaire 2=Décadaire 3=Mensuel 4=Autre 1= Format : 0=Gardénia_Sequentiel 1=Gardénia_Annuaire 2=Libre 3=Excel *** Durée des pas de temps si non standard ************** Standard=Unité de durée des Pas si non standard (sec,min,heu,jou,moi,ann) 1.000=Durée d'un pas de temps (dans l'unité) *** Options du Bassin *** 38=Nombre d'Années des séries de données (Pluie, ETP, Observ.) [0 => Toutes] 11=Nombre d'Années démarrage (-n pour générer n année moy fictives de démarrage) 1969=Date de la Première Année des données (par ex. 2012) 0=Décalage dans la série des Pluies [+5 => Retarde de 5 pas ; -4 Avance de 4 pas] 0=Décalage de la série des Débits/Niveaux observés [ex: -2 => Avance 2 pas ] 0=État initial :0=Pl Effic. ou équil ;1=Niv ou Deb. -1=Vide ; -2=RuMax Vide aussi 150=Nombre maxi. d'itérations pour le calage (-1 = aucune itération; pas de calage) 20=Durée des pluies en moyenne par pas (%) (utilisations avancées)[défaut = 100 %] 1=Nombre de Réservoirs Souterrains (1 ou 2 ou 3=Double + seuil) [def = 1] 0=Utilisation de la même ETP (moyenne pour toutes les années [0=non ; 1=Oui] 0=Numéro du Jour initial [Def=1] (si durée non standard) ; ex. 31 pour 31 Déc. 0=Numéro du Mois initial [Def=1] (si durée non standard) ; ex. 12 pour 31 Déc. 0=Heure initiale [Def=0] (si durée non standard) ; par ex. 15 pour 15h30

BRGM/RP-55813-FR – Rapport final 53 Etude des relations nappe-rivière en Haute-Normandie – Modélisations GARDENIA

0=Minute initiale [Def=0] (si durée non standard) ; Par ex. 30 pour 15h30 *** Paramètres de Prévision *** 0=Nombre d'Années de données du fichier de pluies pour la Prévision 0=Jour d'émission de la prévision (1-31) (Si pas de temps journalier sinon : 0) 0=Numéro du Mois [si journal ou mensuel] (ou du pas) d'émission de la Prévision) 0=Portée de la Prévision (Nombre de pas de la prévision) *** Paramètres Physiques *** 12.00000=Coefficient de correction globale des Pluies (%) Opti= 1 0.00000=Capacité de la réserve superficielle (RuMax) (mm) Opti= 0 0.00000=Coefficient de correction globale de l'ETP (%) Opti= 0 2370.31372=Hauteur de répartition Ruissellement-Percolation (mm) Opti= 1 6.66714=Temps de 1/2 montée (temps de 1/2 Percolation) (mois) Opti= 1 8.10000=Temps de 1/2 tarissement Souterrain 1 (mois) Opti= 1 0.15010=Temps de 1/2 transfert Souterr. 1--> Souterr. 2 (mois) Opti= 0 0.15010=Temps de 1/2 tarissement Souterrain 2 (mois) Opti= 0 1.100E-03=Hauteur de réduction de reprise complém. par ETP (mm) Opti= 0 350.10001=Capacité de la réserve superficielle Progressive (mm) Opti= 1 29.01145=Niveau de base (m NGF) Opti= 1 0.01044=Coefficient d'emmagasinement équivalent (-) Opti= 1 0.00000=Valeur démarrage : Pluie Effic annuelle/Niveau/Débit *** Bornes des paramètres physiques *** -10.00000=Min : Coefficient de correction globale des Pluies (%) Max = 12.00000 0.00000=Min : Capacité de la réserve superficielle (RuMax) (mm) Max = 500.00000 -15.00000=Min : Coefficient de correction globale de l'ETP (%) Max = 15.00000 0.00000=Min : Hauteur de répartition Ruissellement-Percolation (mm) Max = 1.000E+05 0.00000=Min : Temps de 1/2 montée (temps de 1/2 Percolation) (mois) Max = 20.00000 8.10000=Min : Temps de 1/2 tarissement Souterrain 1 (mois) Max = 30.00000 0.15000=Min : Temps de 1/2 transfert Souterr. 1--> Souterr. 2 (mois) Max = 50.00000 0.15000=Min : Temps de 1/2 tarissement Souterrain 2 (mois) Max = 50.00000 1.000E-03=Min : Hauteur de réduction de reprise complém. par ETP (mm) Max =9999.00000 350.00000=Min : Capacité de la réserve superficielle Progressive (mm) Max = 500.00000

54 BRGM/RP-55813-FR – Rapport final

Centre scientifique et technique Service géologique régional Haute-Normandie 3, avenue Claude-Guillemin Parc de la Vatine BP 36009 10 rue A. Sakharov 45060 – Orléans Cedex 2 – France 76130 – Mont Saint Aignan - France Tél. : 02 38 64 34 34 Tél. : 02 35 60 12 00