Simple Method to Extract Lake Ice Condition From Landsat Images

Ice plays key roles in regulating hydrological, ecological, biogeochemical, and socioeconomic functions of lakes. Long-term in situ lake ice phenological records indicate that lake ice is trending toward later freeze-up, earlier breakup, and a shorter ice duration. Parallel to study of lake ice using in situ records and process-based models, satellite remote sensing can expand our understanding of lake ice change over large spatial scales. However, most remote sensing studies have focused on large lakes or short periods of time, which may not robustly represent changes over multidecadal time periods or in the much more numerous small lakes. Here, we present a random forest model, Sensitive Lake Ice Detection (SLIDE), to accurately extract ice conditions from Landsat TM, ETM+, and OLI images. We trained the model using a manually labeled lake ice dataset (1089 labeled areas over 995 lakes globally). Our results show that our model achieves accurate classification between ice/snow and water (accuracy: 97.8%, kappa coefficient: 95.5%). Comparing Landsat-derived ice cover with in situ ice conditions, we show that our model produces less bias, lower RMSE, and higher kappa than does the Landsat snow/ice flag from the quality assessment band. This is especially true during the transitional period surrounding the ice on and off dates reported from in situ (mean bias −7.3% from our model, −17.3% from the Landsat quality band). Our results demonstrate the feasibility of mining the rich Landsat archive to study lake ice dynamics and of better flagging ice-affected lake observations.