Abstract. The increasing air temperature in a changing climate will impact actual evaporation and have consequences for water resources management in energy-limited regions. In many hydrological models, evaporation is assessed by a preliminary computation of potential evaporation (PE) representing the evaporative demand of the atmosphere. Therefore, in impact studies the quantification of uncertainties related to PE estimation, which can arise from different sources, is crucial. Indeed, a myriad of PE formulations exist and the uncertainties related to climate variables cascade into PE computation. So far, no consensus has emerged on the main source of uncertainty in the PE modelling chain for hydrological studies. In this study, we address this issue by setting up a multi-model and multi-scenario approach. We used seven different PE formulations and a set of 30 climate projections to calculate changes in PE. To estimate the uncertainties related to each step of the PE calculation process (namely Representative Concentration Pathways, General Circulation Models, Regional Climate Models and PE formulations), an analysis of variance decomposition (ANOVA) was used. Results show that PE would increase across France by the end of the century, from +40 to +130 mm/year. In ascending order, uncertainty contributions by the end of the century are explained by: PE formulations (below 10 %), then RCPs (above 20 %), RCMs (30–40 %) and GCMs (30–40 %). Finally, all PE formulations show similar future trends since climatic variables are co-dependent to temperature. While no PE formulation stands out from the others, in hydrological impact studies the Penman-Monteith formulation may be preferred as it is representative of the PE formulations ensemble mean and allows accounting for climate and environmental drivers co-evolution.