It is of the most urgent necessity to get hydrological time series of long duration for the establishment of rational design and operation criterion for the Agricultural hydraulic structures. This study was conducted to select best fitted frequency distribution for the monthly runoff and to simulate long series of generated flows by multi-season first order Markov model with comparison of statistical parameters which are derivated from observed and sy- nthetic flows in the five watersheds along Geum river basin. The results summarized through this study are as follows. 1. Both two parameter gamma and two parameter lognormal distribution were judged to be as good fitted distributions for monthly discharge by Kolmogorov-Smirnov method for goodness of fit test in all watersheds. 2. Statistical parameters were obtained from synthetic flows simulated by two parameter gamma distribution were closer to the results from observed flows than those of two para- meter lognormal distribution in all watersheds. 3. In general, fluctuation for the coefficient of variation based on two parameter gamma distribution was shown as more good agreement with the observed flow than that of two parameter lognormal distribution. Especially, coefficient of variation based on two parameter lognormal distribution was quite closer to that of observed flow during June and August in all years. 4. Monthly synthetic flows based on two parameter gamma distribution are considered to give more reasonably good results than those of two parameter lognormal distribution in the multi-season first order Markov model in all watersheds. 5. Synthetic monthly flows with 100 years for eack watershed were sjmulated by multi- season first order Markov model based on two parameter gamma distribution which is ack- nowledged to fit the actual distribution of monthly discharges of watersheds. Simulated sy- nthetic monthly flows may be considered to be contributed to the long series of discharges as an input data for the development of water resources. 6. It is to be desired that generation technique of synthetic flow in this study would be compared with other simulation techniques for the objective time series.