A hybrid model is presented for stochastic simulation of multiseason streamflows. This involves partial prewhitening of the streamflows using a parsimonious linear periodic parametric model, followed by resampling the resulting residuals using moving block bootstrap to obtain innovations and subsequently postblackening these innovations to generate synthetic replicates. This model is simple and is efficient in reproducing both linear and nonlinear dependence inherent in the observed streamflows. The first part of this paper demonstrates the hybrid character of the model through stochastic simulations performed using monthly streamflows of Weber River (Utah) that exhibit a complex dependence structure. In the latter part of the paper the hybrid model is shown to be efficient in simulating multiseason streamflows, through an example of the San Juan River (New Mexico). This model ensures annual-to-monthly consistency without the need for any adjustment procedures. Furthermore, the hybrid model is able to preserve both within-year and cross-year monthly serial correlations for multiple lags. Also, it is seen to be consistent in predicting the reservoir storage (validation) statistic at low as well as high demand levels. ¿ 2001 American Geophysical Union