Abstract. This paper presents the analysis results of temperatures collected at three monitoring stations which were used to study ice freezing and thawing processes on a reservoir along Irtysh River. The measured temperatures were comprehensively analyzed and correlated with air temperature measured at a meteorological station. The results showed that air temperature was closely related to temperature at the ice surface, e.g., T40, T20 and T0, and temperatures in ice almost increased linearly with depth. In addition, ice thicknesses were calculated based on measured temperatures along arrays of temperature and compared with that calculated using a simplified Stefan’s model. The results indicated that ice thickness varied spatially and temporally, and Stefan’s model overestimated the ice thicknesses with a maximum discrepancy of 12 cm. Moreover, the calculated ice thickness was correlated with temperatures, variations of temperature and accumulated freezing degree days (AFDD) based on Pearson correlation analysis, showing that ice thickness was proportional to AFDD with a coefficient of 0.89, and negatively correlated with T0, T−3, and Tice(c) with coefficients of −0.69, −0.73 and −0.62, respectively. Therefore, linear and non-linear models were proposed, which were validated using datasets from three stations in Russia and Finland, demonstrating that the linear model incorporating AFDD and T0 can capture local ice freezing and thawing processes with a relatively minor discrepancy, and the results were consistent at different stations. The paper provides an approach to comprehensively study the ice formation process and a practical model to calculate local ice thickness.