Multiple climate regimes in an idealized lake-ice-atmosphere model

In recent decades, the Laurentian Great Lakes have undergone rapid surface warming with the summertime trends substantially exceeding the warming rates of surrounding land. Warming of the deepest (Lake Superior) was the strongest, and that of the shallowest (Lake Erie)-the weakest of all lakes. To investigate the dynamics of accelerated lake warming, we considered single-column and multi-column thermodynamic lake-ice models coupled to an idealized two-layer atmosphere. The variable temperature of the upper atmospheric layer-a proxy for the large-scale atmospheric forcing-consisted, in the most general case, of a linear trend mimicking the global warming and atmospheric interannual variability, both on top of the prescribed seasonal cycle of the upper-air temperature. The atmospheric boundary layer of the coupled model exchanged heat with the lake and exhibited lateral diffusive heat transports between the adjacent atmospheric columns. In simpler single-column models, we find that, for a certain range of periodic atmospheric forcing, each lake possesses two stable equilibrium seasonal cycles, which we call "regimes"-with and without lake-ice occurrence in winter and with corresponding cold and warm temperatures in the following summer, respectively, all under an identical seasonally varying external forcing. Deeper lakes exhibit larger differences in their summertime surface water temperature between the warm and cold regimes, due to their larger thermal and dynamical inertia. The regime behavior of multi-column coupled models is similar but more complex, and in some cases, they admit more than two stable equilibrium seasonal cycles, with varying degrees of wintertime ice-cover. The simulated lake response to climate change in the presence of the atmospheric noise rationalizes the observed accelerated warming of the lakes, the correlation between wintertime ice cover and next summer's lake-surface temperature, as well as higher warming trends of the (occasionally wintertime ice-covered) deep-lake vs. shallow-lake regions, in terms of the corresponding characteristics of the forced transitions between colder and warmer lake regimes. Since the regime behavior in our models arises due to nonlinear dynamics rooted in the ice-albedo feedback, this feedback is also the root cause of the accelerated lake warming simulated by these models. In addition, our results imply that if Lake Superior eventually becomes largely ice-free (< 10% maximum ice cover every winter) under continuing global warming, the surface warming trends of the deeper regions of the lake will become modest, similar to those of the shallower regions of the lake.