River water temperature is a key control of many physical and bio-chemical processes in river systems, which theoretically depends on multiple factors. Here, four different machine learning models, including multilayer perceptron neural network models (MLPNN), adaptive neuro-fuzzy inference systems (ANFIS) with fuzzy c-mean clustering algorithm (ANFIS_FC), ANFIS with grid partition method (ANFIS_GP), and ANFIS with subtractive clustering method (ANFIS_SC), were implemented to simulate daily river water temperature, using air temperature (T), river flow discharge (Q), and the components of the Gregorian calendar (CGC) as predictors. The proposed models were tested in various river systems characterized by different hydrological conditions. Results showed that including the three inputs as predictors (T, Q, and the CGC) yielded the best accuracy among all the developed models. In particular, model performance improved considerably compared to the case where only T is used as predictor, which is the typical approach of most of previous machine learning applications. Additionally, it was found that Q played a relevant role mainly in snow-fed and regulated rivers with higher-altitude hydropower reservoirs, while it improved to a lower extent model performance in lowland rivers. In the validation phase, the MLPNN model was generally the one providing the highest performances, although in some river stations ANFIS_FC and ANFIS_GP were slightly more accurate. Overall, the results indicated that the machine learning models developed in this study can be effectively used for river water temperature simulation.