Abstract Continuous simulation of monthly evapotranspiration rates for 1979–2015 was performed by the latest, calibration-free version of the complementary relationship of evaporation over the conterminous United States. The results were compared to similar estimates of the WREVAP program and the North American Regional Reanalysis (NARR) project. Validation of the three methods was performed by the Parameter-Elevation Regressions on Independent Slopes Model precipitation and Hydrologic Unit Code level-6 runoff data. The present method outperforms the WREVAP and NARR estimates with a root-mean-square error (RMSE) of 89 mm yr−1, an R2 value of 0.87, an absolute bias (σ) of −5 mm yr−1, and slope (m) and intercept (c) values of 0.97 and 22 mm yr−1, respectively, for the best-fit line, in comparison to similar values (RMSE = 161 mm yr−1, R2 = 0.8, σ = 124 yr–1 mm yr−1, m = 0.88, c = 191 mm yr−1; and RMSE = 195 mm yr−1, R2 = 0.81, σ = 146 mm yr−1, m = 1.05, c = 120 mm yr−1) of the latter two methods. The value of the Priestley–Taylor (PT) coefficient was determined by inversion of the PT-equation via a model-independent identification of wet cells and their estimated surface temperatures.