The characteristics of water quality variation were predicted by stochastic model in Chungju dam, north Chungcheong province of south Korea, Monthly time series data of water quality from 1989 to 2001;temperature, BOD, COD and SS, were obtained from environmental yearbook and internet homepage of ministry of environment. Development of model was carried out with Box-Jenkins method, which includes model identification, estimation and diagnostic checking. ACF and PACF were used to model identification. AIC and BIC were used to model estimation. Seosonal multiplicative ARIMA(1, 0, 1)(1, 1, 0)$_{12}$ model was appropriate to explain stochastic characteristics of temperature. BOD model was ARMa(2, 2, 1), COD was seasonal multiplicative ARIMA(2. 0. 1)(1. 0, 1)$_{12}$, and SS was ARIMA(1, 0, 2) respectively. The simulated water quality data showed a good fitness to the observed data, as a result of model verification.ion.