In this study, five hydrological models, including the soil and water assessment tool (SWAT), identification of unit hydrograph and component flows from rainfall, evapotranspiration, and streamflow (IHACRES), Hydrologiska Byråns Vattenbalansavdelning (HBV), Australian water balance model (AWBM), and Soil Moisture Accounting (SMA), were used to simulate the flow of the Hablehroud River, north-central Iran. All the models were validated based on the root mean square error (RMSE), coefficient of determination (R2), Nash-Sutcliffe model efficiency coefficient (NS), and Kling-Gupta efficiency (KGE). It was found that SWAT, IHACRES, and HBV had satisfactory results in the calibration phase. However, only the SWAT model had good performance in the validation phase and outperformed the other models. It was also observed that peak flows were generally underestimated by the models. The sensitivity analysis results of the model parameters were also evaluated. A hybrid model was developed using gene expression programming (GEP). According to the error measures, the ensemble model had the best performance in both calibration (NS = 0.79) and validation (NS = 0.56). The GEP combination method can combine model outputs from less accurate individual models and produce a superior river flow estimate.