Hierarchy of climate and hydrological uncertainties in transient low flow projections

International audienceThis paper proposes a methodology for estimating the transient probability distribution of yearly hydrological variables conditional to an ensemble of projections built from multiple general circulation models (GCMs), multiple statistical downscaling methods (SDMs) and multiple hydrological models (HMs). The methodology is based on the quasi-ergodic analysis of variance (QE-ANOVA) framework that allows quantifying the contributions of the different sources of total uncertainty, by critically taking account of large-scale internal variability stemming from the transient evolution of multiple GCM runs, and of small-scale internal variability derived from multiple realizations of stochastic SDMs. The QE-ANOVA framework was initially developed for long-term climate averages and is here extended jointly to (1) yearly anomalies and (2) low flow variables. It is applied to better understand possible transient futures of both winter and summer low flows for two snow-influenced catchments in the southern French Alps. The analysis takes advantage of a very large dataset of transient hydrological projections that combines in a comprehensive way 11 runs from 4 different GCMs, 3 SDMs with 10 stochastic realizations each, as well as 6 diverse HMs. The change signal is a decrease in yearly low flows of around -20 % in 2065, except for the most elevated catchment in winter where low flows barely decrease. This signal is largely masked by both large- and small-scale internal variability, even in 2065. The time of emergence of the change signal on 30 year low-flow averages is however around 2035, i.e. for time slices starting in 2020. The most striking result is that a large part of the total uncertainty - and a higher one than that due to the GCMs - stems from the difference in HM responses. An analysis of the origin of this substantial divergence in HM responses for both catchments and in both seasons suggests that both evapotranspiration and snowpack components of HMs should be carefully checked for their robustness in a changed climate in order to provide reliable outputs for informing water resource adaptation strategies