Quantifying forcing uncertainties in the hydrodynamics of the Gironde estuary

Vanessya Laborie (1, 2) , Sophie Ricci (3) , Matthias de Lozzo (4) , Nicole Goutal (5) , Yoann Audouin (5) , Philippe Sergent (2)

1 LHSV - Laboratoire d'Hydraulique Saint-Venant / Saint-Venant laboratory for Hydraulics
2 Cerema Equipe-projet HA - Centre d'Etudes et d'Expertise sur les Risques, l'Environnement, la Mobilité et l'Aménagement - Equipe-projet HA
3 CERFACS - Centre Européen de Recherche et de Formation Avancée en Calcul Scientifique
4 DER - Département Etude des Réacteurs
5 EDF R&D LNHE - Laboratoire National d’Hydraulique et Environnement

High tide combined with high meteorological surge levels and discharges in the Garonne and Dor-dogne rivers in the Gironde estuary (south-west France), may lead to high water levels and flooding near the Blayais nuclear power plant and the city of Bordeaux, with significant economic and social impacts. A global sensitivity analysis (GSA) was performed with a Telemac2D numerical model currently used for operational water level forecasts. The major sources of uncertainties were identified by computing the Sobol' indices for uncertain inputs with an analysis of variance (ANOVA) approach for a 7-day storm event in 2003. The generation of the GSA ensemble of simulations consists of sampling scalar and field random variables: constant and uniform friction coefficients, as well as time-varying hydrological and maritime forcings. The temporal perturbation of time-dependent upstream hydrological and down-Vanessya LABORIE Saint-Venant Hydraulics laboratory, Cerema, 134, route de Beauvais, CS 60039, 60280 Margny Lès Philippe SERGENT Cerema, 134, route de Beauvais, CS 60039, 60280 Margny Lès Compiègne, France stream maritime forcings is assumed to be represented by a Gaussian Process characterized by a correlation time scale calculated from observations. A Karhunen-Loève decomposition was then applied to retain a limited number of eigenmodes. The GSA is performed for 20 random variables using GENCI HPC computational resources for task parallelism and domain decomposition. This requires the use of 250 000 runs for an elapsed simulation time of 101 days on 32768 cores. The performance of the ensemble was assessed with a rank diagram and a reliability curve in comparison to a set of measured water levels at 12 observing stations along the estuary. It was shown that, for this event, the maritime boundary conditions and the Strickler coefficients have a predominant role along the length of the estuary with an influence driven by the tidal cycle. In the upstream fluvial areas, the friction coefficient and hydrological inputs are predominant.

Domaines: Hydrologie; Modélisation et simulation; Ingénierie de l'environnement; Recherche opérationnelle [math.OC]

Vanessya Laborie, Sophie Ricci, Matthias de Lozzo, Nicole Goutal, Yoann Audouin, et al.. Quantifying forcing uncertainties in the hydrodynamics of the Gironde estuary. Computational Geosciences, 2019, 24 (1), pp.181-202. ⟨10.1007/s10596-019-09907-7⟩. ⟨hal-02535233⟩