A B S T R A C TThe seasonal evolution of snow and ice on Lake Oraja¨rvi, northern Finland, was investigated for three consecutive winter seasons. Material consisting of numerical weather prediction model (HIRLAM) output, weather station observations, manual snow and ice observations, high spatial resolution snow and ice temperatures from ice mass balance buoys (SIMB), and Moderate Resolution Imaging Spectroradiometer (MODIS) lake ice surface temperature observations was gathered. A snow/ice model (HIGHTSI) was applied to simulate the evolution of the snow and ice surface energy balance, temperature profiles and thickness. The weather conditions in early winter were found critical in determining the seasonal evolution of the thickness of lake ice and snow. During the winter season (Nov.ÁApr.), precipitation, longwave radiative flux and air temperature showed large inter-annual variations. The uncertainty in snow/ice model simulations originating from precipitation was investigated. The contribution of snow to ice transformation was vital for the total lake ice thickness. At the seasonal time scale, the ice bottom growth was 50Á70% of the total ice growth. The SIMB is suitable for monitoring snow and ice temperatures and thicknesses. The Mean Bias Error (MBE) between the SIMB and borehole measurements was (0.7 cm for snow thicknesses and 1.7 cm for ice thickness. The temporal evolution of MODIS surface temperature (three seasons) agrees well with SIMB and HIGHTSI results (correlation coefficient, R 00.81). The HIGHTSI surface temperatures were, however, higher (2.88C5MBE53.98C) than the MODIS observations. The development of HIRLAM by increasing its horizontal and vertical resolution and including a lake parameterisation scheme improved the atmospheric forcing for HIGHTSI, especially the relative humidity and solar radiation. Challenges remain in accurate simulation of snowfall events and total precipitation.