Precise estimation of reference evapotranspiration (ET0) is of significant importance in hydrologic processes. In this study, a genetic algorithm (GA) optimized back propagation (BP) neural network model was developed to estimate ET0 using different combinations of meteorological data across various climatic zones and seasons in China. Fourteen climatic locations were selected to represent five major climates. Meteorological datasets in 2018–2020, including maximum, minimum and mean air temperature (Tmax, Tmin, Tmean, °C) and diurnal temperature range (∆T, °C), solar radiation (Ra, MJ m−2 d−1), sunshine duration (S, h), relative humidity (RH, %) and wind speed (U2, m s−1), were first subjected to correlation analysis to determine which variables were suitable as input parameters. Datasets in 2018 and 2019 were utilized for training the models, while datasets in 2020 were for testing. Coefficients of determination (r2) of 0.50 and 0.70 were adopted as threshold values for selection of correlated variables to run the models. Results showed that U2 had the least r2 with ET0, followed by ∆T. Tmax had the greatest r2 with ET0, followed by Tmean, Ra and Tmin. GA significantly improved the performance of BP models across different climatic zones, with the accuracy of GABP models significantly higher than that of BP models. GABP0.5 model (input variables based on r2 > 0.50) had the best ET0 estimation performance for different seasons and significantly reduced estimation errors, especially for autumn and winter seasons whose errors were larger with other BP and GABP models. GABP0.5 model using radiation/temperature data is highly recommended as a promising tool for modelling and predicting ET0 in various climatic locations.