Precise estimation of pan evaporation is necessary to manage available water resources. In this study, the capability of three hybridized models for modeling monthly pan evaporation (Epan) at three stations in the Dongting lake basin, China, were investigated. Each model consisted of an adaptive neuro-fuzzy inference system (ANFIS) integrated with a metaheuristic optimization algorithm; i.e., particle swarm optimization (PSO), whale optimization algorithm (WOA), and Harris hawks optimization (HHO). The modeling data were acquired for the period between 1962 and 2001 (480 months) and were grouped into several combinations and incorporated into the hybridized models. The performance of the models was assessed using the root mean square error (RMSE), mean absolute error (MAE), Nash–Sutcliffe Efficiency (NSE), coefficient of determination (R2), Taylor diagram, and Violin plot. The results showed that maximum temperature was the most influential variable for evaporation estimation compared to the other input variables. The effect of periodicity input was investigated, demonstrating the efficacy of this variable in improving the models’ predictive accuracy. Among the models developed, the ANFIS-HHO and ANFIS-WOA models outperformed the other models, predicting Epan in the study stations with different combinations of input variables. Between these two models, ANFIS-WOA performed better than ANFIS-HHO. The results also proved the capability of the models when they were used for the prediction of Epan when given a study station using the data obtained for another station. Our study can provide insights into the development of predictive hybrid models when the analysis is conducted in data-scare regions.