Modeling of pan evaporation (Ep) is of paramount importance in the evaluation of drinking water supplies, planning of regional water resources and reservoir management. The main aim of this study is to investigate the accuracy of linear and non-linear ensemble approaches for monthly Ep modeling in Erbil and Salahaddin meteorological stations of Iraq. For this purpose, sensitivity analysis was performed to determine the dominant input parameters. The results showed that Tmean, Tmax and Tmin are the most effective parameters. Thereafter, two scenarios were involved for the Ep modeling. In scenario 1, the ability of artificial neural network, least-squares support-vector machine and multiple linear regression models was examined for the estimation of Ep. The results demonstrated that different input combinations led to different performance, model 3 (which has Tmean, Tmax, Tmin, RH) for Erbil station and model 2 (which has Tmean, Tmax, Tmin) for Salahaddin station provided the best performance among several models developed. In scenario 2, linear and non-linear ensemble approaches were employed as simple linear average, weighted linear average and non-linear ensemble (NLE) models to improve predictions of the single models. The results reported that ensemble modeling could improve performance of single models and NLE model provided the best results due to its non-linear nature. The general results demonstrated that the proposed ensemble models could improve predictions of single models up to 5% and 16% for Erbil and Salahaddin stations, respectively.