This study introduces the first use of Global Navigation Satellite System Reflectometry (GNSS-R) for monitoring lake ice phenology. This is demonstrated using Qinghai Lake, Tibetan Plateau, as a case study. Signal-to-Noise Ratio (SNR) values obtained from the Cyclone GNSS (CYGNSS) constellation over four ice seasons (2018 to 2022) were used to examine the impact of lake surface conditions on reflected GNSS signals during open water and ice cover seasons. A moving t-test algorithm was applied to time-varying SNR values allowing for the detection of lake ice at daily temporal resolution. Good agreement was achieved between ice phenology records derived from CYGNSS data and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. The CYGNSS timings for freeze-up, i.e., the period starting with the first appearance of ice on the lake (freeze-up start; FUS) until the lake becomes fully ice covered (freeze-up end; FUE), as well as those for breakup, i.e., the period beginning with the first pixel of open water (breakup start; BUS) and ending when the whole lake becomes ice-free (breakup end; BUE), were validated against the phenology dates derived from MODIS images. Mean absolute errors are 7, 5, 10, 4 and 5 days for FUS, FUE, BUS, BUE and ice cover duration, respectively. Observations revealed the sensitivity of GNSS reflected signals to surface melt prior to the appearance of open water conditions as determined from MODIS, which explains the larger difference of 10 days for BUS.

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.