Reference evapotranspiration (ET0) from FAO-Penman-Monteith equation is highly sensitive to the surface incoming solar radiation (SISR) and therefore accurate estimate of this parameter would result in more accurate estimation of ET0. In this study, the accuracy of three main approaches for SISR estimation including empirical models (Angstrom and Hargreaves-Samani), physically-based data assimilation models (Global Land Data Assimilation System-Noah, GLDAS/Noah, and National Centers of Environmental Predictions/National Center for Atmospheric Research, NCEP/NCAR), and a satellite observation model (Satellite Application Facility on Climate Monitoring, CM-SAF) were evaluated using ground-based measurements from 2012 to 2015. Then SISR outputs from introduced approaches were implemented in FAO-Penman-Monteith equation for ET0 estimation on daily and monthly basis. The Angstrom calibrated model was the most accurate model with a coefficient of determination (R2) of 0.9 and standard error of estimate (SEE) of 2.58 MJ. m−2. d−1, and GLDAS/Noah, Hargreaves-Samani, NCEP/NCAR, and CM-SAF, had lower accuracy, respectively. However, the lack of the meteorological data and required empirical coefficients are the main limitations of applying the empirical models, however, satellite-based approaches are more practical for operational purposes. The results indicated that, in spite of slight overestimation in warm months, GLDAS/Noah model had better performance with R2=0.87 and SEE = 3.5 MJ. m−2. d−1 in case of lack of meteorological data. The accuracy of ET0 derived from FAO-Penman-Monteith equation was directly depended on the accuracy of SISR estimation. The ET0 estimation error was related to SISR estimation error with a fourth-degree function and had a linear relationship with SISR error at daily and monthly scales, respectively.