Inaccuracies in air–sea heat fluxes severely degrade the accuracy of ocean numerical simulations. Here, we use artificial neural networks to correct air–sea heat fluxes as a function of oceanic and atmospheric state predictors. The correction successfully improves surface and subsurface ocean temperatures beyond the training period and in prediction experiments.

This study introduces a new 3D lake–ice–atmosphere coupled model that significantly improves winter climate simulations for the Great Lakes compared to traditional 1D lake model coupling. The key contribution is the identification of critical hydrodynamic processes – ice transport, heat advection, and shear-driven turbulence production – that influence lake thermal structure and ice cover and explain the superior performance of 3D lake models to their 1D counterparts.