This paper presents a new method for estimating global land-surface evaporation using satellite observations. The approach uses a variety of satellite-derived data products to estimate daily evaporation at a global scale with a spatial resolution of 0.25 degrees. The method relies on the Priestley and Taylor (PT) evaporation model, which requires a small number of inputs, most of which can be obtained from satellite data. Key components of the method include the use of microwave-derived soil moisture, land surface temperature, vegetation density, and the estimation of rainfall interception loss. The model is validated against one year of eddy covariance measurements from 43 stations, showing good agreement with observed evaporation values (R = 0.80, N = 43) and a low average bias (-5%). The model performs well across different vegetation types and climate conditions, with an average correlation coefficient of 0.83 for daily time series, slightly lower than the 0.90 found for monthly time series. The first global map of annual evaporation derived from this methodology is also presented. The method uses satellite data to create a spatially coherent estimate of evaporative flux over land, with a focus on global validity and the use of satellite-based data sets. The model includes a detailed estimation of rainfall interception loss and uses a Kalman filter to assimilate satellite soil moisture data. The model is validated against FLUXNET observations, showing good agreement with measured evaporation values (R = 0.83 for daily and R = 0.90 for monthly time series). The model also shows good performance in dry regions, with no significant negative impact from the modeling of evaporation stress. The results indicate that the model can effectively capture the global distribution of annual evaporation.This paper presents a new method for estimating global land-surface evaporation using satellite observations. The approach uses a variety of satellite-derived data products to estimate daily evaporation at a global scale with a spatial resolution of 0.25 degrees. The method relies on the Priestley and Taylor (PT) evaporation model, which requires a small number of inputs, most of which can be obtained from satellite data. Key components of the method include the use of microwave-derived soil moisture, land surface temperature, vegetation density, and the estimation of rainfall interception loss. The model is validated against one year of eddy covariance measurements from 43 stations, showing good agreement with observed evaporation values (R = 0.80, N = 43) and a low average bias (-5%). The model performs well across different vegetation types and climate conditions, with an average correlation coefficient of 0.83 for daily time series, slightly lower than the 0.90 found for monthly time series. The first global map of annual evaporation derived from this methodology is also presented. The method uses satellite data to create a spatially coherent estimate of evaporative flux over land, with a focus on global validity and the use of satellite-based data sets. The model includes a detailed estimation of rainfall interception loss and uses a Kalman filter to assimilate satellite soil moisture data. The model is validated against FLUXNET observations, showing good agreement with measured evaporation values (R = 0.83 for daily and R = 0.90 for monthly time series). The model also shows good performance in dry regions, with no significant negative impact from the modeling of evaporation stress. The results indicate that the model can effectively capture the global distribution of annual evaporation.