The present paper describes an effort for developing the total maximum daily load (TMDL) for phosphorus and a load reduction strategy for the Feitsui Reservoir in Northern Taiwan. BASINS model was employed to estimate watershed pollutant loads from nonpoint sources (NPS) in the Feitsui Reservoir watershed. The BASINS model was calibrated using field data collected during a 2-year sampling period and then used to compute watershed pollutant loadings into the Feitsui Reservoir. The simulated results indicate that the average annual total phosphorus (TP) loading into the reservoir is 18,910 kg/year, which consists of non-point source loading of 16,003 kg/year, and point source loading of 2,907 kg/year. The Vollenweider mass balance model was used next to determine the degree of eutrophication under current pollutant loading and the load reduction needed to keep the reservoir from being eutrophic. It was estimated that Feitsui Reservoir can becoming of the oligotrophic state if the average annual TP loading is reduced by 37% or more. The results provide the basis on which an integrated control action plan for both point and nonpoint sources of pollution in the watershed can be developed.