<p>Unlike global climate models, hydrological models cannot simulate the feedbacks among atmospheric processes, vegetation, water, and energy exchange at the land surface. This severely limits their ability to quantify the impact of climate change and the concurrent increase of atmospheric CO<sub>2</sub> concentrations on evapotranspiration and thus runoff. Hydrological models generally calculate actual evapotranspiration as a fraction of potential evapotranspiration (PET), which is computed as a function of temperature and net radiation and sometimes of humidity and wind speed. Almost no hydrological model takes into account that PET changes because the vegetation responds to changing CO<sub>2</sub> and climate. This active vegetation response consists of three components. With higher CO<sub>2</sub> concentrations, 1) plant stomata close, reducing transpiration (physiological effect) and 2) plants may grow better, with more leaves, increasing transpiration (structural effect), while 3) climatic changes lead to changes in plants growth and even biome shifts, changing evapotranspiration. Global climate models, which include dynamic vegetation models, simulate all these processes, albeit with a high uncertainty, and take into account the feedbacks to the atmosphere.</p><p>Milly and Dunne (2016) (MD) found that in the case of RCP8.5 the change of PET (computed using the Penman-Monteith equation) between 1981- 2000 and 2081-2100 is much higher than the change of non-water-stressed evapotranspiration (NWSET) computed by an ensemble of global climate models. This overestimation is partially due to the neglect of active vegetation response and partially due to the neglected feedbacks between the atmosphere and the land surface.</p><p>The objective of this paper is to present a simple approach for hydrological models that enables them to mimic the effect of active vegetation on potential evapotranspiration under climate change, thus improving computation of freshwater-related climate change hazards by hydrological models. MD proposed an alternative approach to estimate changes in PET for impact studies that is only a function of the changes in energy and not of temperature and achieves a good fit to the ensemble mean change of evapotranspiration computed by the ensemble of global climate models in months and grid cells without water stress. We developed an implementation of the MD idea for hydrological models using the Priestley-Taylor equation (PET-PT) to estimate PET as a function of net radiation and temperature. With PET-PT, an increasing temperature trend leads to strong increases in PET. Our proposed methodology (PET-MD) helps to remove this effect, retaining the impact of temperature on PET but not on long-term PET change.</p><p>We implemented the PET-MD approach in the global hydrological model WaterGAP2.2d. and computed daily time series of PET between 1981 and 2099 using bias-adjusted climate data of four global climate models for RCP 8.5. We evaluated, computed PET-PT and PET-MD at the grid cell level and globally, comparing also to the results of the Milly-Dunne study. The global analysis suggests that the application of PET-MD reduces the PET change until the end of this century from 3.341 mm/day according to PET-PT to 3.087 mm/day (ensemble mean over the four global climate models).</p><p>Milly, P.C.D., Dunne K.A. (2016). DOI:10.1038/nclimate3046.</p>