A regional water quality control model is developed by linking a steady-state water quality simulation model with an optimization model. The water quality simulation model can be applied to complex river systems with both point and nonpoint loads using multiple interdependent pollution parameters described by either linear or nonlinear equations. Twelve water quality parameters can be modeled simultaneously: four non conservative constituents (or conservative constituents if the decay rate is set equal to zero); coliform bacteria (MPN); phosphorus; biochemical oxygen demand (BOD); ammonia (NH3); nitrate (NO3); dissolved oxygen (DO); temperature (°C); and algae. The water quality model is used to generate constraint equation for the optimization model. The optimization model is formulated as an integer programming problem in which the integer decision variables are wastewater treatment levels or diffuse source management practices to be determined for each load. The model considers the addition or upgrading of wastewater treatment with structural and nonstructural schemes for both point and diffuse pollution sources. A least cost solution is found subject to water quality standards at surveillance points. Additional constraints can include uniform and zoned uniform treatment. Low flow augmentation and bypass piping can be considered with slight water quality simulation model modification. A simulation model-optimization model iteration procedure is used to find an optimum solution. The regional water quality model is applied to two problems, a hypothetical problem and the Jordan River in Utah. The hypothetical consists of four pollution discharge points, at which seven possible treatment levels are available for six quality constituents: phosphorus; biochemical oxygen demands; ammonia; nitrate; dissolved oxygen; and algae. Water quality standards for three constituents are imposed at five surveillance points along the river. The portion of the Jordan River examined consists of seven pollution discharge points, at which seven treatment levels are available for the same six quality constituents as in the hypothetical problem. Water quality standards for two constituents are imposed at three surveillance points. The cost minimization problem for the Jordan River (1975 flows) required tertiary sand filters at all point loads at an increase from current costs of 1,75,881peryeartomeetstreamwaterqualitystandards.Toassistingainingfamiliaritywiththelinkedsimulationoptimizationmodel,severalsensitivitystudiesareperformed.Thesensitivityoftheoptimalsolutiontotwomodelinputparametersisinvestigated.Theseparametersarethewaterqualityequationcoefficientsandthewaterqualitystreamstandards.Substantialreductionsintreatmentcostswerepossiblebymakingminorchangesinsomeoftheinputparameters.Inthehypotheticalproblem,10percentincreaseintheammoniadecayrateora15percentincreaseinthemaximumspecificalgaegrowthratewouldresultina16percentreductionintheminimumsystemtreatmentcost.A10percentrelaxationofthestreamstandardsatallsurveillancepointswouldresultina54percentreductionintheminimumsystemtreatmentcosts.Theoptimalsolution’ssensitivitytochangesinheadwaterandpointdischargeflowisalsoinvestigated.Theoptimaltreatmentschemefortheprojected1995flowsintheJordanRiverwasthesameasforthe1975flows.Theincreasefromcurrentcostsforthe1995flowswas2,407,092 per year