Application of 1-D lake models coupled within earth-system prediction models will improve accuracy but requires accurate initialization of lake temperatures. Here, we describe a lake initialization method by cycling within a weather prediction model to constrain lake temperature evolution. We compared these lake temperature values with other estimates and found much reduced errors (down to 1-2 K). The lake cycling initialization is now applied to two operational US NOAA weather models.