Abstract. There is growing evidence that climate change will alter water availability in Europe. Here, we investigate how hydrological low flows are affected under different levels of future global warming (i.e. 1.5, 2, and 3 K with respect to the pre-industrial period) in rivers with a contributing area of more than 1000 km 2 . The analysis is based on a multi-model ensemble of 45 hydrological simulations based on three representative concentration pathways (RCP2.6, RCP6.0, RCP8.5), five Coupled Model Intercomparison Project Phase 5 (CMIP5) general circulation models (GCMs: GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, NorESM1-M) and three state-of-theart hydrological models (HMs: mHM, Noah-MP, and PCR-GLOBWB). High-resolution model results are available at a spatial resolution of 5 km across the pan-European domain at a daily temporal resolution. Low river flow is described as the percentile of daily streamflow that is exceeded 90 % of the time. It is determined separately for each GCM/HM combination and warming scenario. The results show that the low-flow change signal amplifies with increasing warming levels. Low flows decrease in the Mediterranean region, while they increase in the Alpine and Northern regions. In the Mediterranean, the level of warming amplifies the signal from −12 % under 1.5 K, compared to the baseline period 1971-2000, to −35 % under global warming of 3 K, largely due to the projected decreases in annual precipitation. In contrast, the signal is amplified from +22 (1.5 K) to +45 % (3 K) in the Alpine region due to changes in snow accumulation. The changes in low flows are significant for regions with relatively large change signals and under higher levels of warming. However, it is not possible to distinguish climate-induced differences in low flows between 1.5 and 2 K warming because of (1) the large inter-annual variability which prevents distinguishing statistical estimates of periodaveraged changes for a given GCM/HM combination, and (2) the uncertainty in the multi-model ensemble expressed by the signal-to-noise ratio. The contribution by the GCMs to the uncertainty in the model results is generally higher than the one by the HMs. However, the uncertainty due to HMs cannot be neglected. In the Alpine, Northern, and Mediterranean regions, the uncertainty contribution by the HMs is partly higher than those by the GCMs due to different representations of processes such as snow, soil moisture and evapotranspiration. Based on the analysis results, it is recommended (1) to use multiple HMs in climate impact studies and (2) to embrace uncertainty information on the multimodel ensemble as well as its single members in the adaptation process.