Results/Conclusions: In comparison to multi-objective optimization method at a single site, the multi-site weighted average approach reduced the relative mean absolute error by 3.9-7.6% for monthly stream flow, 1-6.5% for monthly total nitrogen (TN), and 1.1-7.5% for monthly total phosphorus (TP), respectively. The GWLF model predicted a slight decrease in average annual stream flow and slight increases in average annual TN and TP discharges in the year 2030 compared with a baseline of 2000. By 2095, all stream flow and nutrient discharges were predicted to increase. The relative annual increase rate for TN was greater than for TP; and the relative annual rates of increase for both TN and TP were greater than the relative annual rate of stream flow increase. By 2095, annual discharges of water, TN, and TP were predicted to increase by 18%, 36%, 33%, respectively, compared to the 2000 baseline. In our simulations, the influence of climate change on discharges was much greater than that of land cover change; but only small changes in land cover were predicted for the particular watersheds we modeled (0.6-4.0% increases in urban development in 2030). Our results do indicate that the influence of climate change must be considered in management efforts to reducing nutrient loading to the Chesapeake Bay.