Climate change is causing rapid warming and altered precipitation patterns in mountain watersheds, both of which influence the timing of ice breakup in mountain lakes. To enable predictions of ice breakup in the future, we analyzed a dataset of mountain lake ice breakup dates derived from remote sensing and historical downscaled climate data. We evaluated drivers of ice breakup, constructed a predictive statistical model, and developed projections of mountain lake ice breakup date with global climate models. Using Random Forest analysis, we determined that winter and spring cumulative snow fraction (portion of precipitation falling as snow) and air temperature are the strongest predictors of ice breakup on mountain lakes. Interactions between precipitation, cumulative winter air temperature and lake surface area indicate that shifts in air temperature and precipitation affect smaller lakes (< 2 km2) more than larger lakes (> 2–10 km2). A linear mixed effects model (RMSE of 18 d), applied with an ensemble of 15 global climate models, projected that end-of-century ice breakup in mountain lakes will be earlier by 25 ± 4 and 61 ± 5 (mean ± SE) days for representative concentration pathways 4.5 and 8.5, respectively.
