A number of models have been suggested •or hydrologic time series in general and streamflow series in particular. Most of them are normal autoregressive (AR) of order 1 with either constant or periodic parameters. Since generally hydrologic time series are nonnormal (skewed), transformations have been suggested to make the series approximately normal. A new class of univariate models is proposed herein which incorporates skewed and correlation properties within the model structure without the necessity of transformations. Such models assume a gamma marginal distribution and a constant or periodic aut0regressive structure. The models may be additive gamma, multiplicative. gamma, or a mixed model which incorporates properties of both additive and multiplicative models. The gamma models were tested and compared in relation to (transformed) normal AR models by computer simulation studies based on five weekly streamflow series with samples varying from 35 to 40 years of record. The results show that the new class of gamma models compares favorably with respect to the normal models in reproducing the basic statistics usually analyzed for streambow simulation. It is expected that the proposed gamma models will be of interest to other researchers for further developments and applications to hydrologic and geophysical time series. INTRODUCTION Synthetic generation of streamflow sequence_s is commonly utilized in planning and operational studieg of water resources systems. Such synthetic sequences are generated so as to rep•roduce a set of statistical properties which are found in the historical streamflow data. For such purpose, a number of time series models have been suggested since Hennan [1955] and Thomas and Fiering [1962] applied G•tussian auto•egressire (AR) models to seasonal hydrologic series. In modeling these series it is usually assumed that they are seasonally stationary and that seasonality is reflected in the mean, variance, covariance, and the skewness. However, in general, one may assume seasonality in the marginal distribution. Annual streamflow series may be considered approximately normal' however, for shorter time interval series such as monthly, weekly, and daily series, departures from normality become important. The problem of skewed streamflows has been handled by a number of methods. A Widely used technique is to use transformations to render a series close to normal [Box and Cox, 1965]. Examples of this approach are given by Matalas [1967], Mejia and Rodri•iuez-lturbe [1974], Lettenrnaier and Burqes [1977], Hirsch [1979], and Stedinqer [1981]. Another approach is to find the statistical properties of the noise in such a way as to reproduce the skew of the series. See, for instance, the Wilson-Hilferty (W-H) transformations used by Thomas and Fierinq [1962], Beard [1965], Payne et al. [1969], McGinnis and Sammens [1970], McMahon and Miller [1971], and O'Connell [1974]. Although the coefficient of skewness can be reproduced by this approach, the underlying variable is not gamma. A related problem ...