Modeling changes in ice dynamics and subsurface thermal structure in Lake Michigan-Huron between 1979 and 2021

The world’s largest lakes, including the Laurentian Great Lakes, have experienced significant surface warming and loss of ice cover over the last several decades. Although changing surface conditions have received substantial research interest, changes below the surface remain largely unexplored, despite their importance for turbulent mixing, nutrient cycling, and primary production. In this study, we investigate changes in subsurface thermal structure and timing in Lake Michigan-Huron related to ongoing climate warming. This work utilizes atmospheric reanalysis data to drive the Great Lakes Finite Volume Community Ocean Model (GL-FVCOM), providing three-dimensional hydrodynamic and ice simulations between 1979 and 2021. Results are used to analyze trends in ice and temperature dynamics, revealing significant changes in annually averaged ice cover (− 2.1– − 5.2