Abstract. The regional climate model REMO was coupled with the FLake lake model to include an interactive treatment of lakes. Using this new version, the Fenno-Scandinavian climate and lake characteristics were studied in a set of 35-year hindcast simulations. Additionally, sensitivity tests related to the parameterization of snow albedo were conducted. Our results show that overall the new model version improves the representation of the Fenno-Scandinavian climate in terms of 2 m temperature and precipitation, but the downside is that an existing wintertime cold bias in the model is enhanced. The lake surface water temperature, ice depth and ice season length were analyzed in detail for 10 Finnish, 4 Swedish and 2 Russian lakes and 1 Estonian lake. The results show that the model can reproduce these characteristics with reasonably high accuracy. The cold bias during winter causes overestimation of ice layer thickness, for example, at several of the studied lakes, but overall the values from the model are realistic and represent the lake physics well in a long-term simulation. We also analyzed the snow depth on ice from 10 Finnish lakes and vertical temperature profiles from 5 Finnish lakes and the model results are realistic.