Using a suite of numerical models, we investigated the influence of legacy P on water quality in the Yahara Watershed of southern Wisconsin, USA. The suite included Agro-IBIS, a terrestrial ecosystem model; THMB, a hydrologic and nutrient routing model; and the Yahara Water Quality Model which estimates water quality indicators in the Yahara chain of lakes. Using five alternative scenarios of antecedent P storage (legacy P) in soils and channels under historical climate conditions, we simulated outcomes of P yield from the landscape, lake P loading, and three lake water quality indicators. Data and code may also be found in https://github.com/SRCarpen/Yahara2070LakeModel_Fits2data.git