Improving the quantification of climate change hazards by hydrological models: a simple ensemble approach for considering the uncertain effect of vegetation response to climate change on potential evapotranspiration

Abstract. Almost no hydrological model takes into account that changes in evapotranspiration are affected by how vegetation responds to changing CO2 and climate. This severely limits their ability to quantify the impact of climate change on evapotranspiration and, thus, water resources. As the simulation of vegetation responses is both complex and very uncertain, we recommend a simple approach to considering (in climate change impact studies with hydrological models) the uncertainty that the vegetation response causes with respect to the estimation of future potential evapotranspiration (PET). To quantify this uncertainty in a simple manner, we propose running the hydrological model in two variants: with its standard PET approach and with a modified approach to compute PET. In the case of PET equations containing stomatal conductance, the modified approach can be implemented by adjusting the conductance. We introduce a modified approach for hydrological models that computes PET as a function of net radiation and temperature only, i.e., with the Priestley–Taylor (PT) equation. The new PT-MA approach is based on the work of Milly and Dunne (2016) (MD), who compared the change in non-water-stressed actual evapotranspiration (NWSAET) as computed by an ensemble of global climate models (GCMs), which simulate vegetation response as well as interactions between the atmosphere and the land surface, with various methods to compute PET change. Based on this comparison, MD proposed estimating the impact of climate change on PET as a function of only the change in net energy input at the land surface. PT-MA retains the impact of temperature on daily to interannual as well as spatial PET variations but removes the impact of the long-term temperature trend on PET such that long-term changes in future PET are driven by changes in net radiation only. We implemented PT-MA in the global hydrological model WaterGAP 2.2d and computed daily time series of PET between 1901 and 2099 using the bias-adjusted output of four GCMs. Increases in GCM-derived NWSAET between the end of the 20th and the end of the 21st century for Representative Concentration Pathway 8.5 (RCP8.5) are simulated well by WaterGAP if PT-MA is applied but are severely overestimated with the standard PT method. Application of PT-MA in WaterGAP results in smaller future decreases or larger future increases in renewable water resources (expressed as the variable RWR) compared with the standard PT method, except in a small number of grid cells where increased inflow from upstream areas due to increased upstream runoff leads to enhanced evapotranspiration from surface water bodies or irrigated fields. On about 20 % of the global land area, PT-MA leads to an increase in RWR that is more than 20 % higher than in the case of standard PT, while on more than 10 % of the global land area, the projected RWR decrease is reduced by more than 20 %. While the modified approach to compute PET is likely to avoid the overestimation of future drying in many if not most regions, the vegetation response in other regions may be such that the application of the standard PET leads to more likely changes in PET. As these regions cannot be identified with certainty, the proposed ensemble approach with two hydrological model variants serves to represent the uncertainty in hydrological changes due to the vegetation response to climate change that is not represented in the model.