An Automated and Improved Methodology to Retrieve Long-time Series of Evapotranspiration Based on Remote Sensing and Reanalysis Data

The large-scale quantification of accurate evapotranspiration (ET) time series has substantially been developed in recent decades using automated approaches based on remote sensing data. However, there are still several model-related uncertainties that require precise assessment. In this study, the Surface Energy Balance Algorithm for Land (SEBAL) and meteorological data from the Global Land Data Assimilation System (GLDAS) were used to estimate long-term daily actual ET based on three endmember selection procedures: two land cover-based models, one with (WF) and the other without (WOF) morphological functions, and the Allen method (with the default percentiles) for 2270 Landsat images. Models were evaluated for 23 flux tower sites with four main vegetation cover types as well as different climate types. Results showed that endmember selection with morphological functions (WF_ET) generally performed better than the other endmember approaches. Climate-based classification assessment provided the clearest discrimination between the performance of the different endmember selection approaches for the humid category. For humid zones, the land cover-based methods, especially WF, appropriately outperformed Allen. However, the performance of the three approaches was similar for sub-humid, semi-arid and arid climates together; the Allen approach was therefore recommended to avoid the need for dependency on land cover maps. Tower-by-tower validation also showed that the WF approach performed best at 12 flux tower sites, the WOF approach best at 5 and the Allen approach best at 6, suggesting that the use of land cover maps alone does not explain the differences between the performance of the land cover-based models and the Allen approach. Additionally, the satisfactory error metrics results when comparing the EC estimations with EC measurements, with root mean square error (RMSE) ≈ 0.91 and 1.59 mm·day−1, coefficient of determination (R2) ≈ 0.71 and 0.41, and bias percentage (PBias) ≈ 2% and 60% for crop and non-crop flux tower sites, respectively, supports the use of GLDAS meteorological forcing datasets with the different automated ET estimation approaches. Overall, given that the thorough evaluation of different endmember selection approaches at large scale confirmed the validity of the WF approach for different climate and land cover types, this study can be considered an important contribution to the global retrieval of long time series of ET.