<h4>The global annual lake ice phenological dataset includes the freeze-up date, break-up date for 74,245 lakes in the northern hemisphere. The dataset is divided into three parts: 1) current data, obtained from MODIS productions through a DLRM model (with parameters provided), covering the period of 2001 to 2020; 2) historical and 3) future simulation data, obtained from the temperature-based lake-specific models, for the periods of 1861-2005 and 2006-2099, respectively. The historical and future simulations were only performed for 30,063 lakes that meet the model conditions and are presented in the dataset.</h4>