Abstract. Evapotranspiration (ET) links the hydrological, energy, and carbon cycle on the land surface. Quantifying ET and its spatiotemporal changes is also key to understanding climate extremes such as droughts, heatwaves and flooding. Regional ET estimates require reliable observationally-based gridded ET datasets, and while many have been developed using physically-based, empirically-based and hybrid techniques, their efficacy, and particularly the efficacy of their uncertainty estimates, is difficult to verify. In this work, we extend the methodology used in Hobeichi et al. (2018) to derive a new version of the Derived Optimal Linear Combination Evapotranspiration (DOLCE) product, with observationally constrained spatiotemporally varying uncertainty estimates, higher spatial resolution, more constituent products and extended temporal reach (1980–2018). After successful evaluation of the efficacy of these uncertainty estimates out-of-sample, we derive novel ET climatology clusters for the land surface, based on the magnitude and variability of ET at each location. The verified uncertainty estimates and extended time period then allow us to examine the robustness of historical trends spatially and in each of these six ET climatology clusters. We find that despite robust decreasing ET trends in some regions, these do not correlate with behavioural ET clusters. Each cluster, and the vast majority of the Earth's surface, show clear robust increases in ET over the recent historical period.