Improving the simulation of large lakes in regional climate modeling: Two-way lake-atmosphere coupling with a 3D hydrodynamic model of the great lakes

Accurate representations of lake–ice–atmosphere interactions in regional climate modeling remain one of the most critical and unresolved issues for understanding large-lake ecosystems and their watersheds. To date, the representation of the Great Lakes two-way interactions in regional climate models is achieved with one-dimensional (1D) lake models applied at the atmospheric model lake grid points distributed spatially across a 2D domain. While some progress has been made in refining 1D lake model processes, such models are fundamentally incapable of realistically resolving a number of physical processes in the Great Lakes. In this study, a two-way coupled 3D lake-ice–climate modeling system [Great Lakes–Atmosphere Regional Model (GLARM)] is developed to improve the simulation of large lakes in regional climate models and accurately resolve the hydroclimatic interactions. Model results are compared to a wide variety of observational data and demonstrate the unique skill of the coupled 3D modeling system in reproducing trends and variability in the Great Lakes regional climate, as well as in capturing the physical characteristics of the Great Lakes by fully resolving the lake hydrodynamics. Simulations of the climatology and spatiotemporal variability of lake thermal structure and ice are significantly improved over previous coupled, 1D simulations. At seasonal and annual time scales, differences inmodel results are primarily observed for variables that are directly affected by lake surface temperature (e.g., evaporation, precipitation, sensible heat flux) while no significant differences are found in other atmospheric variables (e.g., solar radiation, cloud cover). Underlying physical mechanisms for the simulation improvements using GLARM are also discussed