Rainfall may be considered as the most important source of drinking water and watering land in different areas all over the world. Therefore, simulation and estimation of the hydrological phenomenon is of paramount importance. In this study, for the first time, the long-term rainfall in Rasht city was simulated using an optimum hybrid artificial intelligence (AI) model over a 62 year period from 1956 to 2017. The gene expression programming (GEP) and wavelet transform (WT) were combined to develop the hybrid AI model (WGEP). Firstly, the most effective lags of time series data were identified by means of the autocorrelation function (ACF); then eight various GEP and WGEP models were defined. Next, the GEP models were analyzed and the superior GEP model as well as the most influenced lags was detected. For instance, the variance accounting for (VAF), correlation coefficient (R) and scatter index (SI) for the superior GEP model was calculated to be 0.765, 0.508 and 0.709, respectively. Additionally, lags (t-1), (t-2), (t-3) and (t-12) were the most influenced. Then, the different mother wavelets were examined, indicating that the demy mother wavelet was the most optimal one. Moreover, analyzing the numerical simulations showed that the mother wavelet enhanced the performance of the GEP model significantly. For example, the VAF index for the superior WGEP model was increased almost three times after using the mother wavelet. Furthermore, the R and MARE statistical indices for the WGEP model were computed to be 0.935 and 0.862, respectively.