The interannual variability of the annual maximum ice cover (AMIC) of the Great Lakes is strongly influenced by large-scale atmospheric circulations that drive regional weather patterns. Based on statistical analyses from 1980 to 2020, we identify a reduced number of accumulated freezing degree days across the winter months in recent decades, a step-change decrease of AMIC after the winter of 1997/98, and an increased interannual variability of AMIC since 1993. Our analysis shows that AMIC is significantly correlated with El Niño-Southern Oscillation (ENSO), the North Atlantic Oscillation (NAO), and the Pacific/North American (PNA) before the winter of 1997/98. After that, the AMIC is significantly correlated with the Tropical/Northern Hemisphere (TNH) and Eastern Pacific Oscillation (EPO). Singular value decomposition of the 500 hPa geopotential height and surface air temperature shows a dipole pattern over the northeast Pacific and North America, demonstrating the ridge-trough system. This dipole pattern shifts northward to the northern Rocky Mountains, placing the Great Lakes region in the trough after 1997/98. This shift coincides with the increased interannual variability of the EPO index, as well as the change in the sea surface temperature (SST) over the northeast Pacific, where the second mode of the empirical orthogonal function (EOF) on SST shows a warm blob-like feature manifested over the Gulf of Alaska. The regression of wave activity flux onto the SST EOF shows that the source of upward and eastward propagation of a stationary Rossby wave shifts to the west coast of North America, likely moving the ridge-trough system eastward after the winter of 1997/98.