The lake-rich Qinghai-Tibet Plateau (QTP) has significant impacts on regional and global 15 water cycles and monsoon systems through heat and water vapor exchange. The lake-atmosphere 16 interactions have been quantified over open-water periods, yet little is known about the lake ice 17 thermodynamics and heat and mass balance during ice-covered season due to a lack of field data. 18Modeling experiments on ice evolution and energy balance were performed in a shallow lake with a 19 high-resolution snow and ice thermodynamic model. The bottom ice growth and decay dominated the 20 seasonal evolution of the thickness of lake ice. Strong surface sublimation was a crucial pattern of ice 21 loss, which was up to 40% of the maximum ice thickness. Simulation results matched well with the 22 observations with respect to ice mass balance components, net ice thickness, and ice temperature. 23Strong solar radiation, consistent freezing air temperature, and low air moisture were the major driving 24 forces controlling the seasonal ice mass balance. Energy balance was estimated at the ice surface and 25 bottom, and within the ice interior and under-ice water. Particularly, almost all heat fluxes showed 26 significant diurnal variations including short-and long-wave radiation, turbulent heat fluxes, water heat 27 fluxes at ice bottom, and absorbed and penetrated solar radiation. The calculated ice surface 28 temperature indicated that the atmospheric boundary layer was consistently stable and neutral over the 29 ice-covered period. The turbulent heat fluxes between the lake ice and air and the net heat gain by the 30 lake were much lower than those during open-water period. Ice surface sublimation (vapor flux) was 31 demonstrated to be a vital seasonal water balance component, accounting for 41% of lake water loss 32 during the ice seasons. 33