Precise determination of the reference evapotranspiration (ET0) is vital to studying the hydrological cycle. In addition, It plays a significant role in properly managing and allocating water resources in agriculture. The objective of this research was to examine the effectiveness of five different data-driven techniques, including artificial neural networks "multilayer perceptron" (ANN), gene expression programming (GEP), random forest (RF), support vector machine "radial basis function" (SVM), and multiple linear regression (MLR) to model the daily ET0. These methods were also compared with Hargreaves-Samani (HS), Oudin, Ritchie, Makkink (MAK), and Jensen Haise (JH) empirical models and their calibrated versions. The empirical models JH and MAK performed better than the models HS and Oudin after being calibrated by linear regression. All data-driven methods with four inputs were superior to the original and calibrated empirical models. The RF and ANN methods generally demonstrated better estimation accuracy than other data-driven methods. The performance of the RF and ANN models that utilized Tmax, Tmin, and Rs inputs, as well as those that incorporated Tmax, Tmin, Rs, and U2 inputs, proved to be superior to their corresponding MLR-based and GEP-based models for predicting ET0 in the Adana plain, which is characterized by a Mediterranean climate. Nevertheless, the GEP and MLR methods have the advantage of utilizing explicit algebraic equations, making them more convenient to apply, especially in the context of agricultural irrigation practices.