This dataset contains modeled daily lake area, volume, constituent mass, and biogeochemical rates for 3,692 lakes in the Northern Highlands Lake District (NHLD) for one retrospective model run (1986-2010) and 12 model runs under future climate scenarios. This dataset was created using published tools developed to simulate detailed hydrological and biogeochemical fluxes for thousands of lakes and reservoirs over large spatiotemporal scales. The lake hydrology model utilized a computationally-efficient integrated surface water and groundwater modeling framework that informed a lake water budget model incorporating daily hydrologic inputs and exports from individual lakes within the modeling domain. The lake biogeochemical model was informed by the hydrologic information and was built upon a simple lake energy budget, constituent loading, and lake biogeochemical model to track carbon storage and processing for all lakes within the NHLD modeling domain. Our one retrospective model run was driven by historic meteorological data and the projected model runs were driven by projected future climate scenario periods that are representative through the year 2100. For more details on the historic and projected driver data and model set up, please see Zwart et al. (year and DOI to be entered once MS is published).