The WATEQ4F-based "multiple step method", or the method of running WATEQ4F repeatedly, was used to evaluate the error propagation in computations of water-mineral interactions under the circumstances of constant and varied temperatures due to the uncertainties of input variables. The results show the following: the errors of water chemistry analysis can strongly affect the modelling results of water-mineral reactions; different input variables (errors) have different effects on the saturation indices (S.I.) of different minerals; in many cases, the S.I. errors of minerals change with temperatures.