AbstractThe lakewide seasonal carbon cycle of Lake Michigan is poorly quantified and lacks the mechanistic links necessary to determine impacts upon it from eutrophication, invasive species, and climate change. A first step toward a full appreciation of Lake Michigan's carbon cycle is to quantify the dominant mechanisms of its internal carbon cycle. To achieve this, we use the MIT general circulation model configured to the bathymetry of Lake Michigan and coupled to an ecosystem model to simulate the seasonal cycle of productivity, temperature, circulation, and the partial pressure of CO2 in water (pCO2). This biogeochemistry is designed to be appropriate for the prequagga mussel state of the lake. The primary mechanism behind the seasonal cycle of primary productivity is lake physics. The offshore spring phytoplankton bloom begins following a reduction in deep vertical mixing and ends with the depletion of nutrients via thermal stratification. The exception is the western shoreline, where summer winds drive coastal upwelling, providing hypolimnetic nutrients and generating significant productivity. Surface pCO2 is controlled by the net effect from temperature on solubility, and is modulated by biological uptake of dissolved inorganic carbon (DIC) and isothermal mixing of DIC‐rich water in winter. Temperature tends to have the greatest seasonal impact in nearshore regions, while local DIC has the greatest impact in offshore regions. Lakewide, the model suggests that carbon is absorbed from the atmosphere during the spring bloom and released to the atmosphere during winter mixing and when summer surface temperatures are at their maximum.