This study was conducted to simulate long seres of synthetic monthly flows by multi-season first order Markov model with selection of best fitting frequency distribution, harmonic synthetic and harmonic regression models and to make a comparison of statistical parameters between observes and synthetic flows of five watersheds in Geum river system. The results obtained through this study can be summarized as follow. 1. Both gamma and two parameter lognormal distributions were found to be suitable ones for monthly flows in all watersheds by Kolmogorov-Smirnov test. 2. It was found that arithmetic mean values of synthetic monthly flows simulated by multi-season first order Markov model with gamma distribution are much closer to the results of the observed data in comparison with those of the other models in the applied watersheds. 3. The coefficients of variation, index of fluctuation for monthly flows simulated by multi-season first order Markov model with gamma distribution are appeared closer to those of the observed data in comparison with those of the other models in Geum river system. 4. Synthetic monthly flows were simulated over 100 years by multi-season first order Markov model with gamma distribution which is acknowledged as a suitable simulation modal in this study.