In the face of rapidly compounding climate change impacts, including ocean acidification (OA), it is critical to understand present-day stress exposure and to anticipate the biogeochemical conditions experienced by vulnerable ecosystems like coral reefs. To meaningfully predict nearshore carbonate chemistry, we must account for the complexity of the local benthic community, as well as connectivity between habitats and relevant endmember carbonate chemistry. Here, we adopt a system-scale approach to predict site-scale effects of benthic metabolism on the carbonate system of the Florida Reef Tract (FRT). We utilize bimonthly carbonate chemistry data from ten cross-shelf transects spanning 250 km of the FRT to model changes in dissolved inorganic carbon (DIC) and total alkalinity (TA). Benthic habitat maps were used to broadly classify communities known to impact carbonate chemistry. A SLIM 2D hydrodynamic model with mesh resolution reaching 100 m over reefs and along the coastline was used to determine the relevant water mass histories and identify the upstream benthic communities shaping local carbonate chemistry. These historical metabolic footprints, or “flowsheds”, were used to build predictive models of the change in DIC and TA at each station. The best predictive models included the chemical impacts of benthic ecosystem metabolism, as defined by water mass trajectories, weighted endmember chemistry, volume, time, and other environmental parameters (light, temperature, salinity, chlorophyll-a, and nitrate). Considering water mass for 5 days prior to sample collection yielded the highest model skill.