This study examined the evaporation from Lake Erie during the record lake effect snowfall event, November 17th-20th, 2014, by reconstructing heat fluxes and evaporation rates over Lake Erie using the unstructured grid, Finite-Volume Community Ocean Model (FVCOM). Nine different model runs were conducted using combinations of three different flux algorithms: the Met Flux Algorithm (COARE), a method routinely used at NOAA’s Great Lakes Environmental Research Laboratory (SOLAR), and the Los Alamos Sea Ice Model (CICE); and three different meteorological forcings: the Climate Forecast System version 2 Operational Analysis (CFSv2), Interpolated observations (Interp), and the High Resolution Rapid Refresh (HRRR). A few non-FVCOM model outputs were also included in the evaporation analysis from an atmospheric reanalysis (CFSv2) and the large lake thermodynamic model (LLTM). Model-simulated water temperature and meteorological forcing data (wind direction and air temperature) were validated with buoy data at three locations in Lake Erie. The simulated sensible and latent heat fluxes were validated with the eddy covariance measurements at two offshore sites; Long Point Lighthouse in north central Lake Erie and Toledo water crib intake in western Lake Erie. The evaluation showed a significant increase in heat fluxes over three days, with the peak on the 18th of November. Snow water equivalent data from the National Snow Analyses at the National Operational Hydrologic Remote Sensing Center showed a spike in water content on the 20th of November, two days after the peak heat fluxes. The ensemble runs presented a variation in spatial pattern of evaporation, lake-wide average evaporation, and resulting cooling of the lake. Overall, the evaporation tended to be larger in deep water than shallow water near the shore. The lake-wide average evaporations from CFSv2 and LLTM are significantly smaller than those from FVCOM. The variation among the nine FVCOM runs resulted in the 3D mean water temperature cooling in a range from 3 degrees C to 5 degrees C (6-10 EJ loss in heat content), implication for impacts on preconditioning for the upcoming ice season.