The topic of satellite remote sensing of lake ice has gained considerable attention in recent years. Optical satellite data from the Moderate Resolution Imaging Spectroradiometer (MODIS) allow for the monitoring of lake ice cover (an Essential Climate Variable or ECV), and dates associated with ice phenology (freeze-up, break-up, and ice cover duration) over large areas in an era where ground-based observational networks have nearly vanished in many northern countries. Ice phenology dates as well as dates of maximum and minimum ice cover extent (for lakes that do not form a complete ice cover in winter or do not totally lose their ice cover in summer) are useful for assessing long-term trends and variability in climate, particularly due to their sensitivity to changes in near-surface air temperature. Existing knowledge-driven (threshold-based) retrieval algorithms for lake ice cover mapping that use top-of-atmosphere (TOA) reflectance products do not perform well under lower solar illumination conditions (i.e. large solar zenith angles), resulting in low TOA reflectance. This research assessed the capability of four machine learning classifiers (i.e. multinomial logistic regression, MLR; support vector machine, SVM; random forest, RF; gradient boosting trees, GBT) for mapping lake ice cover, water and cloud cover during both break-up and freeze-up periods using the MODIS/Terra L1B TOA (MOD02) product. The classifiers were trained and validated using samples collected from 17 large lakes across the Northern Hemisphere (Europe and North America); lakes that represent different characteristics with regards to area, latitude, freezing frequency, and ice duration. Following an accuracy assessment using random k-fold cross-validation (k = 100), all machine learning classifiers using a 7-band combination (visible, near-infrared and shortwave-infrared) were found to be able to produce overall classification accuracies above 94%. Both RF and GBT provided overall and class-specific accuracies above 98% and a more visually accurate depiction of lake ice, water and cloud cover. The two tree-based classifiers offered the most robust spatial transferability over the 17 lakes and performed consistently well across ice seasons. However, only RF was relatively insensitive to the choice of the hyperparameters compared to the other three classifiers. The results demonstrate the potential of RF for mapping lake ice cover globally from MODIS TOA reflectance data. • This study assessed the capability of ML classifiers for lake ice mapping from MOIDS. • RF and GBT produced the best performance in terms of classification accuracies. • RF and GBT offered the most robust spatial and temporal transferability. • RF was insensitive to the choice of the hyperparameters compared to other classifiers. • The results show the potential of RF for mapping lake ice cover globally from MODIS.