Evapotranspiration (ET) is one of the key components of the global hydrological cycle. Many models have been established to obtain an accurate estimation of ET, but the uncertainty of each model has not been satisfactorily addressed, and the weight determination in multi-model simulation methods remains unclear. In this study, the Bayesian model averaging (BMA) method was adopted to tackle this issue. We explored the combination of four surface energy balance (SEB) models (SEBAL, SSEB, S-SEBI and SEBS) with the BMA method by using Landsat 8 images over two study areas in China, the Huailai flux station (semiarid region) and the Sidaoqiao flux station (arid/semiarid region), and the data from two stations were used as validation for this method. The performances of SEB models and different BMA methods is revealed by three statistical parameters (i.e., the coefficient of determination (R2), root mean squared error (RMSE), and the Nash-Sutcliffe efficiency coefficient (NSE)). We found the best performing SEB model was SEBAL, with an R2 of 0.609 (0.672), RMSE of 1.345 (0.876) mm/day, and NSE of 0.407 (0.563) at Huailai (Sidaoqiao) station. Compared with the four individual SEB models, each of the BMA methods (fixed, posterior inclusion probability, or random) can provide a more accurate and reliable simulation result. Similarly, in Huailai (Sidaoqiao) station, the best performing BMA random model provided an R2 of 0.750 (0.796), RMSE of 0.902 (0.602) mm/day, and NSE of 0.746 (0.793). We conclude that the BMA method outperformed the four SEB models alone and obtained a more accurate prediction of ET in two cropland areas, which provides important guidance for water resource allocation and management in arid and semiarid regions.