Due to limited direct measurements, regional or global terrestrial evapotranspiration (ET) is generally derived from a combination of meteorological and satellite observations. Although the inhomogeneity of the observed climate data has been widely reported, its impact on the calculated ET has not been adequately quantified. This study aimed to calculate ET using the modified Penman‐Monteith (MPM) model with raw and homogenized meteorological data. Additionally, we compared the calculated ET with those estimates from variable methods (water balance, satellite‐based, and reanalysis) in China and its six major river basins from 1982 to 2020. During the overlapping period of 1997–2018, ET calculated from raw input data decreased slightly at −0.39 mm yr−2 (p = 0.64) in China, whereas homogenized ET showed a significant increasing trend of 0.93 mm yr−2 (p = 0.02), with a better agreement with water balance ET (1.93 mm yr−2, p = 0). Global Land Evaporation Amsterdam Model (GLEAM) and Modern‐Era Retrospective Analysis for Research and Applications, version 2 (MERRA2) could reproduce the increasing trends with 2.08 mm yr−2 (p = 0) and 2.59 mm yr−2 (p = 0). The intercomparison of input variables (solar radiation, relative humidity, wind speed, precipitation, and air temperature) among ET products revealed substantial differences, which can account for the discrepancies in ET estimates. Homogenized ET, GLEAM and MERRA2 exhibited significant increasing trends in China and most river basins from 1982 to 2020. Our findings underscore the importance of utilizing homogenized input data for more accurate ET estimation.