Rainfall is one of the most important components of the water cycle and plays a very important role in the measurement of climate characteristic in any area. Limitations such as lack of sufficient information about the amount of rainfall in time and space scale and complexity of the relationship between meteorological elements related to rainfall, causes the calculation of these parameters using the conventional method not to be implemented. One method of evaluating and forecasting of rainfall in each region is time series models. In this research, to predict the average annual rainfall synoptic station at Mahabad, Uromiya and Mako in West Azarbayejan provience during 1984-2013, linear time series ARIMA was used. To investigate model static, Auto Correlation Function (ACF) and Partial Auto Correlation Function (PACF) was applied and with differencing method, the non-static data transformed to static data. In next step, stochastic models to estimate the annual rainfall average were used. With regard to the evaluation criterion such as T, P-VALUE < 0.05 and Bayesian Information Creterion (BIC), ARIMA (1,0,0), ARIMA (0,1,1) and ARIMA (0,1,1) models was determined as a suitable model for predicting annual rainfall in the three selected stations at Uromiya, Makoo and Mahabad. In the following, the annual rainfall for 3 (2013-2016) years is forecasted which based on rainfall data in that time, the adjusted model was acceptable.