Abstract. Water temperatures in three morphometrically different lakes are simulated using a one-dimensional hydrodynamic lake model over the century (1911–2014) to elucidate the effects of increasing air temperature and decreasing wind speed on lake thermal variables (water temperature, stratification dates, strength of stratification, and surface heat fluxes). During the study period, epilimnetic temperatures increased, hypolimnetic temperatures decreased, and the length of the stratified season increased for the study lakes due to earlier stratification onset and later fall overturn. Additionally, there was an abrupt change in epilimnion temperature after 1930 in both Lake Mendota and Lake Wingra, and three changes, after 1934, 1995, and 2008 for Fish Lake. There was a significant change in the slope of trend of stratification duration after 1940 in Lake Mendota and a significant change in trend after 1981 for Fish Lake. Schmidt stability showed a statistically significant increasing trend for both deep lakes, with the larger trend and greater variability in the larger surface area lake. Sensible heat flux in all three lakes increases over the simulation period while longwave heat flux decreases. The shallow study lake had a greater change in latent heat flux and net heat flux, illustrating the role of lake depth to surface heat fluxes. Sensible heat flux in all three lakes had similar timing of abrupt changes, but the magnitude of the change increased with increasing depth. Abrupt changes in latent heat flux appear to be independent of lake morphometry, indicating that the timing of change may be primarily driven by climate. Perturbing drivers showed that increasing air temperature and decreasing wind speed caused earlier stratification onset and later fall overturn. For hypolimnetic water temperature, however, increasing air temperature warmed bottom waters while decreasing wind speed cooled bottom waters, indicating that the change of hypolimnetic temperatures globally may be influenced by local changes in wind speed. Overall, lake depth impacts the presence of stratification and magnitude of Schmidt stability, while lake surface area drives differences in hypolimnion temperature, hypolimnetic heating, variability of Schmidt stability, and stratification onset and fall overturn dates.