Using Exceedance Probability to Determine Total Maximum Daily Loads for Reservoir Water Quality Management

Total maximum daily loads (TMDLs) are used to protect water bodies based on their assimilative ability and are transferred as a maximum allowable load, which is the sum of all pollution emissions in a watershed that cannot be exceeded. This allowable load is calculated from a target water quality concentration and a flow state. The target water quality concentration is typically consistent with water quality standards; however, it is difficult to determine which flow state to use, especially for lakes and reservoirs. In this study, an exceedance probability method is established for determining the TMDL for reservoirs. The SWMM (Storm Water Management Model) was used to understand the pollution loads from the watershed, and the Vollenweider model was used to simulate the total phosphorous (TP) concentrations in the reservoir. Using the validated Vollenweider model, the relationship between pollution loads and the target TP concentration is illustrated. This relationship is associated with real changes in the reservoir water volume and is presented as the exceedance probability. In the study area (i.e., the Shiman Reservoir in Taiwan), an exceedance probability of 50% is suggested, and the allowable TP load is 22,209 kg/year when considering a target TP concentration of 20 μg/L. When considering effective management, the pollution sources in three hot spots are priorities, and controlling their point and nonpoint pollution sources can decrease TP from 25 to 22 μg/L.


Introduction
Total maximum daily loads (TMDLs) were originally developed for the Clean Water Act in the United States, which was passed to protect waterbodies based on their assimilative abilities or carrying capacities.Wastewater effluents are required to comply with specified discharge pollution concentrations and total emission loadings.The maximum load of a receiving waterbody is allocated to pollution sources as their allowable emission loading quota [1,2].The TMDL provides complete protection at the watershed scale and includes point and nonpoint pollution sources.In addition, the TMDL considers water quality management from effluent-based control to ambient-based management [3], accounting for relationships between pollution sources and water quality conditions in the waterbody.The TMDL has been widely accepted and applied in the U.S. [4].It has also been applied in other countries; for example, Kang et al. [5] considered the TMDL in a small watershed in Korea, Boyacioglu and Alpaslan [3] demonstrated a TMDL program in the Tahtali Basin, Turkey, and Wang et al. [6] established TMDLs for Taihu Lake, China.
In Taiwan, the Water Pollution Control Act suggests that the TMDL be used when effluent standards fail to meet water-use standards and when special protection is needed for the water.A TMDL was established in Taiwan in 2016 that was aimed at limiting the maximum loads of heavy metal for protecting agricultural land.However, few other realistic TMDL cases in Taiwan have been addressed, and those that have are difficult to expand due to their high complexity, including the complexities involved in determining assimilative capacity and a feasible allocation policy.
The development and application of TMDLs are more complex and require more analysis tools than do traditional pollution control policies.At least three quantification process are required: the calculation of point and nonpoint pollution source emissions, the determination of the allowable maximum loads of the receiving waterbody, and the associations between the pollution source emissions and water quality in the receiving waterbody.Mathematical models often play an important role in TMDL processes [7,8].In general, dynamic models, load duration curves, general watershed models, and export coefficients are commonly used to develop TMDLs [9].Point and nonpoint source emissions and their effects on water quality can be analyzed by modeling with several different water quality models.Mishra and Deng [10] demonstrated a new 1-D sediment transport model that consisted of GIS, the universal soil loss equation (USLE) model, and Hydrologic Engineering Center's River Analysis System (HEC-RAS) to establish a sediment TMDL.Diffuse pollution is difficult to identify and quantify, and the transportation of diffuse pollutants across an entire watershed is difficult to determine due to high uncertainties [11].Gulati et al. [12] indicated that the flow and concentration of diffuse pollution do not follow a normal distribution.This result implies that although allowable pollutant loads can be calculated from waterbody requirements, accurately identifying the pollution sources and the water quality in receiving waterbodies is very challenging.However, before the load allocations of pollution sources, the allowable maximum loads of the receiving waterbody remain uncertain.
Generally, the objective of waterbody protection is to meet the designated use and associated water quality concentrations of the waterbody.The water quality standard is a fixed value, but the flow of the waterbody and the carrying capacity fluctuate.For a river or stream, the TMDL is the target water quality concentration multiplied by the low flow discharge.In the U.S., the critical flow considered for the TMDL of a river or stream is 7Q10 (the average flow of the 7 consecutive days with the lowest flow every 10 years).In Taiwan, the Q 75 (Q 75 refers to the daily flow value that is met or exceeded 75% of the days) is used, which implies a 75% protection level, and the acceptable probability of violating the water quality standard is 25%.Using lower flow to determine the TMDL results in greater protection.When using low flow, the allowable maximum load is low and the limitation to pollution sources is large.Unlike river or stream TMDLs, reservoir TMDLs are more complex.In a reservoir, the input flow is not equal to the outflow.Thus, the methods used to determine the TMDL for a reservoir are different from those used for a river, in which the flow and water quality concentration are identical at the input and output.The allowable maximum loads for a reservoir are decided by changes in the reservoir water volume.When the water volume decreases, the assimilative ability decreases and cannot handle the usual pollution inputs.All water entering the reservoir is then typically detained in the reservoir for a period of time, and the water quality changes.Thus, more details must be considered when determining the TMDL of a reservoir than when determining that of a river.In addition to the hydraulic effects, the pollution loads are complex.For lake P TMDLs, both internal and external P loads must be considered.Steinman and Ogdahl [13] reviewed previously established TMDLs for internal P loads and concluded that the internal loading target had already been met and that no management action was needed.Therefore, in the present study, reservoir TMDLs are determined by considering changes in the reservoir water volume and the external pollution loads of the watershed.The effects of internal pollution loads are excluded.
In this study, an exceedance probability method is used to establish reservoir TMDLs.The use of probability to illustrate water quality is not a new; for example, Havens and Walker [14] related total phosphorus (TP) concentrations to the occurrence probability of algae bloom and applied an empirical model to evaluate the level of protection of the waterbody from bloom occurrence under maximum allowable external TP loads.However, the use of exceedance probability for determining the TMDL of a reservoir has not been conducted previously.In addition to the novel evaluation of this method, a complete flowchart, including pollution source analysis, water quality modeling at the watershed and reservoir scales, the association of pollution sources with the resulting water quality, and pollution reduction strategies for hot spots are demonstrated in this study.This procedure is applied in a Taiwanese reservoir as a case study and is expected to be applicable in other cases.

Study Flowchart
This study focused on reservoir TMDL and required a reservoir water quality model and a watershed model.A flowchart of this study is shown in Figure 1.The watershed model was used to model the nonpoint source pollution and the transportation of pollutants across the watershed.The reservoir water quality model was used to predict water quality changes in reservoirs.The pollutants generated from the watershed became the inputs of the reservoir water quality model.Before combining the two models, model calibration and the validation of the watershed model were conducted.When the quality of the pollutant input was acceptable, model calibration and validation of the reservoir model were implemented sequentially.The validated reservoir model illustrates the relationship between pollutant inputs from the watershed and the water quality changes in the reservoir.It was used to determine the allowable pollutant loadings under a required water quality standard.The main contribution of this study is that an exceedance probability curve of the allowable pollution loads under fixed and ideal water quality concentrations in a reservoir is suggested as a TMDL value.In addition, the validated watershed model can demonstrate the spatial distributions of pollutants and identify pollution hot spots, which can inform the design of pollution control measures based on the determined TMDL value.The reservoir and watershed models in the study flowchart are not the only models possible when using this approach.Users can select models according to their preferences and the available data to ensure that the results of the model are reliable.
Water 2016, 8, 541 3 of 14 the TMDL of a reservoir has not been conducted previously.In addition to the novel evaluation of this method, a complete flowchart, including pollution source analysis, water quality modeling at the watershed and reservoir scales, the association of pollution sources with the resulting water quality, and pollution reduction strategies for hot spots are demonstrated in this study.This procedure is applied in a Taiwanese reservoir as a case study and is expected to be applicable in other cases.

Study Flowchart
This study focused on reservoir TMDL and required a reservoir water quality model and a watershed model.A flowchart of this study is shown in Figure 1.The watershed model was used to model the nonpoint source pollution and the transportation of pollutants across the watershed.The reservoir water quality model was used to predict water quality changes in reservoirs.The pollutants generated from the watershed became the inputs of the reservoir water quality model.Before combining the two models, model calibration and the validation of the watershed model were conducted.When the quality of the pollutant input was acceptable, model calibration and validation of the reservoir model were implemented sequentially.The validated reservoir model illustrates the relationship between pollutant inputs from the watershed and the water quality changes in the reservoir.It was used to determine the allowable pollutant loadings under a required water quality standard.The main contribution of this study is that an exceedance probability curve of the allowable pollution loads under fixed and ideal water quality concentrations in a reservoir is suggested as a TMDL value.In addition, the validated watershed model can demonstrate the spatial distributions of pollutants and identify pollution hot spots, which can inform the design of pollution control measures based on the determined TMDL value.The reservoir and watershed models in the study flowchart are not the only models possible when using this approach.Users can select models according to their preferences and the available data to ensure that the results of the model are reliable.

Reservoir Model: Vollenweider Model
The reservoir water quality model used in this study is the Vollenweider model [15,16], which has been used widely to assess phosphorous (P) concentrations in lakes and reservoirs.The Vollenweider model is a one-dimensional mass-balance model that assumes that the change in the P concentration over time is equal to the input P concentration minus the output P concentration and the amount of P lost inside the waterbody.Based on this assumption, the P is completely mixed in the reservoir and the output P concentration is the same as the P concentration in the reservoir.In addition, P is lost when it settles in the sediment.The Vollenweider equation is as follows.
where PV = total mass of P in reservoir (g); P = P concentration (g/m 3 ); V = reservoir water volume (m 3 ); t = time; M = annual mass of P input (g/year); Q = annual volume of water outflow (m 3 /year); The main problem encountered when using the Vollenweider model is that the settling coefficient is difficult to measure in different waterbodies.One of the most common solutions was developed by Kirchner and Dillon [17], who suggested that the proportion of P lost to sediment (R p ) can be estimated as R p = v/ (v + q s ), where v is the settling velocity (m/year) and q s is the areal hydraulic load (m/year).The parameter q s is calculated from the total flow entering the reservoir (m 3 /year) divided by the surface area of the reservoir (m 2 ).When incorporating this equation into the Vollenweider equation under the steady-state assumption, Equation ( 1) can be transformed to Equations ( 2) and (3).
In Equation ( 3), the P concentration is calculated based on the input P loads (M), outflow discharge (Q), and lost proportion (R p ). R p is obtained by considering the settling velocity (v) and the actual inflow and reservoir surface area (A) data.In this study, M was generated from a validated watershed model.Inflow and outflow data and water surface area were obtained from official monitoring stations.Under ideal conditions, as the water volume changes, the water surface area also changes.However, we did not have precise surface area data for calculating q s , which is used to calculate R p .Instead, a maximum water surface area was used.Using a fixed A affects the accuracy of q s , but its influence might be less than the influences of M and Q in the final TP concentration equation.This is a limitation of field-obtained data and is one factor affecting the precision of model simulation.The only unknown parameter was settling velocity (v).The parameter is adjusted and determined via model calibration and validation.Once we verified the settling velocity, the allowable M for a target P concentration was generated considering different inflow and outflow conditions.Finally, we obtained a curve of the allowable M for a target P concentration and determined the distribution of the curve using the exceedance probability.

Watershed Model: SWMM Model
M can be obtained using several methods, such as the simple export coefficient method or via a complex watershed model.In this study, the SWMM (Storm Water Management Model) was used because the monitored climate, flow, and water quality data were sufficient and because this model has been successfully applied in several cases in Taiwan [18,19].In addition, the SWMM model includes an LID (low impact development) module, and it can assess the performance of different structural best management practice (BMP) facilities.The SWMM model was established by the U.S. Environmental Protection Agency (USEPA) and is freely available for downloading from their website.Before running the SWMM model, the GIS layers, including the watershed boundary, land use, digital elevation model (DEM), and river layers, are required.These GIS data are uploaded and combined using the BASINS (Better Assessment Science Integrating Point and Nonpoint Sources) platform to generate the necessary properties of the watershed for the SWMM model calculations.While the characteristics of the watershed are prepared, the rainfall data are inputted to simulate the hydrological performance.Model parameters or coefficients are determined by local references or through a calibration and validation process.
In the SWMM model, subwatersheds are hydrologic units of land in which the topography and drainage system elements direct surface runoff to a single discharge point.Each subwatershed is divided into pervious and impervious areas.The percentages of pervious and impervious area of the subwatershed are decided by the land use.Surface runoff can infiltrate the pervious subarea but not the impervious subarea.In the calculation of surface runoff, each subwatershed surface is treated as a nonlinear reservoir.Surface runoff occurs when the depth of water in the subwatershed exceeds the maximum depression storage, in which case the outflow is calculated by Manning's equation.Therefore, Manning's n for overland flow for both pervious and impervious areas is required.SWMM offers three infiltration methods: Horton infiltration, Green-Ampt infiltration, and Curve Number infiltration.In the present case study, we selected the Horton infiltration method.If users consider the effects of snow and groundwater, the snow-pack object and groundwater parameters need to be assigned.In this study, we did not include these two objects.
In the SWMM model, nonpoint source pollution is calculated from buildup and wash-off functions.The two functions calculate how many pollutants accumulate on dry days and wash away on wet days.Different functions are used for different land-use types, and the coefficients of the functions are confirmed through model validation.The detailed partitions in different phases and the transportation of pollutants are not included in this model.A more precise watershed model, GSSHA (Gridded Surface Subsurface Hydrologic Analysis), is suggested if more data are available [20].
The SWMM model should be calibrated and validated to ensure its reliability before it is used to obtain input P loads for the Vollenweider model.Realistic flow and water quality data obtained from monitoring stations were collected to compare with the simulated data at the same locations.Several statistical analyses can be used, such as r (Pearson's correlation coefficient), R 2 (coefficient of determination), MAPE (mean absolute percentage error), and NSE (Nash-Sutcliffe efficiency) [21].We used R 2 and MAPE, which are commonly used in Taiwan, to determine the fitness of the model.Various model parameters need to be adjusted, and the parameters for flow and water quality are different.It is helpful to begin by considering the most sensitive model parameters.The most sensitive model parameters of the SWMM model are described by Chen et al. [19].When the final model parameters are determined by model calibration and validation, the parameter values only fit the observed data.Although more observations might be helpful to obtain more realistic parameter values, in the present case, the parameter values introduce unavoidable uncertainty.Based on the results of model evaluation, we can adopt sufficiently accurate simulations, resulting in an acceptable level of model uncertainty.The unacceptable uncertainty should not be applied for decision-making.

Study Area
The Shiman Reservoir was constructed in 1964 and has a total water storage volume of 201.24 million m 3 , a water surface area of approximately 8 km 2 and a watershed area of 736.4 km 2 .The Shiman Reservoir is a multiple-function reservoir that is used for drinking water, irrigation, electricity generation, flood control, and recreation.The reservoir provides drinking water for nearly 2 million people in Taoyuan City, Taiwan.It is located upstream of Dahan Creek in northern Taiwan (Figure 2).The elevation of the watershed ranges from 100 m to 3500 m; accordingly, the flow rate is high, and serious soil erosion occurs.In this reservoir, several efforts have been made to mitigate soil erosion and reduce the amount of sediment entering the reservoir.Previous reservoir management The monitoring stations in the watershed of the reservoir include 6 weather stations, 5 flow monitoring stations, and 13 water quality monitoring stations.The locations of these stations are shown in Figure 2. Rainfall data are required to apply a watershed model.Thus, the rainfall data from the 6 stations were merged using the Thiessen polygon method, in which the contributing area is divided and used as a weighting factor.The land-use and river distribution of the watershed is shown in Figure 3.The land use of the Shiman Reservoir watershed is 78.3% forest and 14.8% agricultural land.Only 0.4% of the total land area is used for residential use.The potential pollution sources are nonpoint source pollution from agricultural lands and point-source pollution from untreated domestic sewers.The increasing tourist population is also a threat to water quality because the sewage system is not distributed throughout the entire watershed.
The average total nitrogen (TN)/TP ratio in the reservoir is 13.67, which indicates that TP is the most restricting factor in the reservoir.Thus, to control eutrophication, the TP input should be reduced.The standard TP concentration of the Shiman Reservoir is 20 μg/L.However, in the last 10 years, the maximum TP concentration reached as high as 48.4 μg/L, and the average TP concentration was 25.0 μg/L (Figure 4).The monitoring stations in the watershed of the reservoir include 6 weather stations, 5 flow monitoring stations, and 13 water quality monitoring stations.The locations of these stations are shown in Figure 2. Rainfall data are required to apply a watershed model.Thus, the rainfall data from the 6 stations were merged using the Thiessen polygon method, in which the contributing area is divided and used as a weighting factor.The land-use and river distribution of the watershed is shown in Figure 3.The land use of the Shiman Reservoir watershed is 78.3% forest and 14.8% agricultural land.Only 0.4% of the total land area is used for residential use.The potential pollution sources are nonpoint source pollution from agricultural lands and point-source pollution from untreated domestic sewers.The increasing tourist population is also a threat to water quality because the sewage system is not distributed throughout the entire watershed.
The average total nitrogen (TN)/TP ratio in the reservoir is 13.67, which indicates that TP is the most restricting factor in the reservoir.Thus, to control eutrophication, the TP input should be reduced.The standard TP concentration of the Shiman Reservoir is 20 µg/L.However, in the last 10 years, the maximum TP concentration reached as high as 48.4 µg/L, and the average TP concentration was 25.0 µg/L (Figure 4).

Watershed Model Calibration and Validation Results
The observed flow data consisted of continuous daily data, and the observed water quality consisted of monthly data.The water quality was sampled and analyzed once a month, yielding a total of 12 observed data sample sets in each year.The data from 2012 were used for model calibration to identify adequate model parameters, and the data from 2013 were used for model verification to verify the feasibility of the model parameters.Figure 5 presents the results of the flow simulation, and Figure 6 presents the results of the TP simulation.

Watershed Model Calibration and Validation Results
The observed flow data consisted of continuous daily data, and the observed water quality consisted of monthly data.The water quality was sampled and analyzed once a month, yielding a total of 12 observed data sample sets in each year.The data from 2012 were used for model calibration to identify adequate model parameters, and the data from 2013 were used for model verification to verify the feasibility of the model parameters.

Watershed Model Calibration and Validation Results
The observed flow data consisted of continuous daily data, and the observed water quality consisted of monthly data.The water quality was sampled and analyzed once a month, yielding a total of 12 observed data sample sets in each year.The data from 2012 were used for model calibration to identify adequate model parameters, and the data from 2013 were used for model verification to verify the feasibility of the model parameters.Figure 5 presents the results of the flow simulation, and Figure 6 presents the results of the TP simulation.In the model simulation process, flow is simulated before water quality.The R 2 values obtained from the calibration and verification tests were 0.83 and 0.64, respectively.After the simulation of flow was deemed acceptable, the water quality simulation was implemented, and the model parameters for simulating flow no longer changed during the water quality simulation.The MAPEs for the TP simulation were 30.0%for calibration and 31.0%for validation.Although the MAPE values were considered acceptable, the R 2 values for TP simulation were low: 0.38 for calibration and 0.18 for validation.However, the R 2 value is focused on the correlation.Because the water quality data are not continuous like the flow data, the R 2 value might not be appropriate for representing the  Loss percentage of P in the reservoir Calculated from v/(v + q s ) q s Areal hydraulic load (m/year) Calculated from Q in /A M Total TP input mass (g/year) Obtained from the verified SWMM P TP concentration in the reservoir (g/m 3 ) Calculation results from the Vollenweider model v Settling velocity of P (m/year) 10 (after model verification)

Reservoir Model Calibration and Validation Results
In the Vollenweider model, all equation parameters are denoted from realistic data except for the settling velocity (v), which differs among different reservoirs, and no experimental value for the Shiman Reservoir was obtained.Therefore, this parameter must be verified before it is applied.Observed TP concentrations were compared with the simulated TP concentrations to determine the appropriate v value.The input TP mass was obtained using verified SWMM.The input data used in the Vollenweider model are summarized in Table 2.The TP concentration in the reservoir was sampled each season for a total of four observations per year.Therefore, we used 5-year data to calibrate and verify the model (Figure 7).Consequently, the verified settling velocity for the Shiman Reservoir was 10 m/year.

Determination of Maximum Allowable Daily Loads of TP
The main objective of this study was to determine the maximum allowable daily load of TP for reservoirs.When the target TP concentration is determined, the total TP input mass is determined from the input and outflow of the reservoir using the previously verified reservoir model.Thus, when we input realistic discharge values for the reservoir into the model, the allowable TP input mass that meets the target TP concentration is produced.
The required P concentration is 20 μg/L for the Shiman Reservoir.The allowable TP input mass at the fixed TP concentration of 20 μg/L is obtained using the previous 5 years of flow data, and its accumulative probability distribution is shown in Figure 8.In addition to the TP of 20 μg/L, a TP concentration of 24 μg/L is shown in the figure because a TP concentration greater than 24 μg/L is recognized as a eutrophic state based on the CTSI (Carlson Trophic State Index), which is officially used to represent eutrophication states in Taiwan.The allowable TP input mass at the exceedance probability of 50% means that the input mass could protect the reservoir by avoiding the TP concentration target at the probability level of 50%; we call this input mass M50.Therefore, M75

Determination of Maximum Allowable Daily Loads of TP
The main objective of this study was to determine the maximum allowable daily load of TP for reservoirs.When the target TP concentration is determined, the total TP input mass is determined from the input and outflow of the reservoir using the previously verified reservoir model.Thus, when we input realistic discharge values for the reservoir into the model, the allowable TP input mass that meets the target TP concentration is produced.
The required P concentration is 20 µg/L for the Shiman Reservoir.The allowable TP input mass at the fixed TP concentration of 20 µg/L is obtained using the previous 5 years of flow data, and its accumulative probability distribution is shown in Figure 8.In addition to the TP of 20 µg/L, a TP concentration of 24 µg/L is shown in the figure because a TP concentration greater than 24 µg/L is recognized as a eutrophic state based on the CTSI (Carlson Trophic State Index), which is officially used to represent eutrophication states in Taiwan.The allowable TP input mass at the exceedance probability of 50% means that the input mass could protect the reservoir by avoiding the TP concentration target at the probability level of 50%; we call this input mass M 50 .Therefore, M 75 indicates that there is a 75% probability that this input mass could protect the reservoir and prevent TP concentrations greater than the target value.M 90 provides a higher level of protection because the allowable input mass is more limited.The allowable TP input masses for different exceedance probabilities are listed in Table 3.We suggest that M 50 can be used as the TMDL of the reservoir as a short-term strategy and that M 75 or M 90 can be used as a long-term strategy.The particular exceedance probability applied in any given case will depend on policy decision-making.For the Shiman Reservoir, the proposed allowable TP input mass is 22,209 kg/year (M 50 ).
indicates that there is a 75% probability that this input mass could protect the reservoir and prevent TP concentrations greater than the target value.M90 provides a higher level of protection because the allowable input mass is more limited.The allowable TP input masses for different exceedance probabilities are listed in Table 3.We suggest that M50 can be used as the TMDL of the reservoir as a short-term strategy and that M75 or M90 can be used as a long-term strategy.The particular exceedance probability applied in any given case will depend on policy decision-making.For the Shiman Reservoir, the proposed allowable TP input mass is 22,209 kg/year (M50).A simple average water volume should not be used to determine the reservoir TMDL because the water volume is controlled by manual operation and is not normally distributed (see Figure 8).For the Shiman Reservoir, considering an average water volume could indicate that the allowable pollution load is approximately 70,000 kg/year.This load would meet the TP standard at the 10% probability level, with the reservoir experiencing pollution or eutrophication >90% of the time.When considering the actual hydraulic conditions, the exceedance probability would provide feasible and objective information for determining the assimilative capacity of a reservoir.

Pollution Control Strategies and Their Performances for Improving Reservoir Water Quality
When determining the allowable TP input mass, the following procedure was used to allocate the input mass into the watershed and decrease the pollution content of the reservoir.Because it is a large watershed, the pollution hot spots in the reservoir are first defined as priority pollution control areas.The three following principles are used for defining these hot spots: the total contributed pollution loads, the percentages of point pollution sources, and the pollution load per land area.Of the 25 subwatersheds considered in this study, three subwatersheds were ranked as hot spots (Figure 9).The subwatersheds S1, S3, and S8 were characterized as areas of high pollution potential because  A simple average water volume should not be used to determine the reservoir TMDL because the water volume is controlled by manual operation and is not normally distributed (see Figure 8).For the Shiman Reservoir, considering an average water volume could indicate that the allowable pollution load is approximately 70,000 kg/year.This load would meet the TP standard at the 10% probability level, with the reservoir experiencing pollution or eutrophication >90% of the time.When considering the actual hydraulic conditions, the exceedance probability would provide feasible and objective information for determining the assimilative capacity of a reservoir.

Pollution Control Strategies and Their Performances for Improving Reservoir Water Quality
When determining the allowable TP input mass, the following procedure was used to allocate the input mass into the watershed and decrease the pollution content of the reservoir.Because it is a large watershed, the pollution hot spots in the reservoir are first defined as priority pollution control areas.The three following principles are used for defining these hot spots: the total contributed pollution loads, the percentages of point pollution sources, and the pollution load per land area.Of the 25 subwatersheds considered in this study, three subwatersheds were ranked as hot spots (Figure 9).The subwatersheds S1, S3, and S8 were characterized as areas of high pollution potential because they have high total pollution loads, a high percentage of point pollution sources, and high pollution loads per land area.At this stage, the subwatersheds with high potential pollution were identified.When detailing the sites and locations for best management practices (BMPs) and when determining load allocation, the advanced spatial optimization model might be needed, which includes cost-effectiveness and uncertainty considerations [22,23].
The current average TP concentration in the studied reservoir is 25 µg/L, and the current average TP input mass is 32,475 kg/year.If the target allowable TP input mass is 22,209 kg/year (M 50 ), pollution should be reduced by 31.6%.Because the locations of the pollution hot spots have been defined, the required pollution reduction can be implemented in these areas first.The potential reduction of TP from point sources in the three hot spots is 1018 kg/year, and the potential reduction of TP nonpoint sources is 2694 kg/year.The reduction of TP from point pollution sources could be achieved by improving domestic wastewater treatment.The reduction of TP from nonpoint sources could be achieved by using various BMPs, including structural and nonstructural measures.Therefore, controlling pollution at the three hot spots could decrease the TP from its current value of 25 µg/L to 22 µg/L.The TMDL and predicted TP concentration are shown in Figure 10.The control of hot spots can be regarded as a short-term management goal.If the goal of 20 µg/L is to be reached, pollution control measures must be implemented in the remaining subwatersheds.
Water 2016, 8, 541 12 of 14 they have high total pollution loads, a high percentage of point pollution sources, and high pollution loads per land area.At this stage, the subwatersheds with high potential pollution were identified.When detailing the sites and locations for best management practices (BMPs) and when determining load allocation, the advanced spatial optimization model might be needed, which includes costeffectiveness and uncertainty considerations [22,23].The current average TP concentration in the studied reservoir is 25 μg/L, and the current average TP input mass is 32,475 kg/year.If the target allowable TP input mass is 22,209 kg/year (M50), pollution should be reduced by 31.6%.Because the locations of the pollution hot spots have been defined, the required pollution reduction can be implemented in these areas first.The potential reduction of TP from point sources in the three hot spots is 1018 kg/year, and the potential reduction of TP from nonpoint sources is 2694 kg/year.The reduction of TP from point pollution sources could be achieved by improving domestic wastewater treatment.The reduction of TP from nonpoint sources could be achieved by using various BMPs, including structural and nonstructural measures.Therefore, controlling pollution at the three hot spots could decrease the TP from its current value of 25 μg/L to 22 μg/L.The TMDL and predicted TP concentration are shown in Figure 10.The control of hot spots can be regarded as a short-term management goal.If the goal of 20 μg/L is to be reached, pollution control measures must be implemented in the remaining subwatersheds.

Conclusions
Determining allowable load of a waterbody is important for implementing TMDLs and affects the subsequent pollution control practices.However, it is difficult to determine the flow status in a reservoir because the inflow and outflow are not equal.In this study, we used the exceedance probability to demonstrate the relationships between water quality and allowable loads.The results of this study imply that allowable loads can be met (without exceeding the target water quality concentration).The concept of this study is simple and can be applied to any lake or reservoir.In this study, we suggest using M 50 , which represents an exceedance probability of 50%.For stricter limitations, M 75 or M 90 can be considered.Notably, using the average reservoir volume significantly overestimates the assimilative ability and might fail to protect the reservoir.The manual operation of the reservoir water level affects the water volume and TP concentration.Thus, we suggest that the average reservoir water volume should not be used to determine the allowable loads.
To determine the allowable loads, it is necessary to combine the watershed and reservoir models.The choice of model is not limited to the two models used in this study.We chose the SWMM and Vollenweider models because they have been successfully applied in other cases in Taiwan.Users of this approach can follow the flowchart in this study and select the appropriate models for their case studies.In this study, we provide a total pollution reduction and do not show the water improvement performance of each pollution control practice, such as different BMPs.If no feasible watershed model exists, the exceedance probability method can be used to determine the maximum allowable load of a reservoir by replacing the watershed results with actual observations of inflow pollution loads.

Figure 1 .
Figure 1.Flowchart showing how the exceedance probability curve was produced for reservoir total maximum daily loads (TMDLs).

Figure 1 .
Figure 1.Flowchart showing how the exceedance probability curve was produced for reservoir total maximum daily loads (TMDLs).

14 Figure 2 .
Figure 2. The location of the Shiman Reservoir and the locations of the monitoring stations within the reservoir.

Figure 2 .
Figure 2. The location of the Shiman Reservoir and the locations of the monitoring stations within the reservoir.

Figure 3 .
Figure 3.The land use and rivers in the Shiman Reservoir watershed.The watershed was divided into 25 subwatersheds, which were used in the watershed model simulations.(a) The distribution of land use in the watershed; and (b) the results of watershed delineation.

Figure 3 .Figure 3 .
Figure 3.The land use and rivers in the Shiman Reservoir watershed.The watershed was divided into 25 subwatersheds, which were used in the watershed model simulations.(a) The distribution of land use in the watershed; and (b) the results of watershed delineation.
Figure 5  presents the results of the flow simulation, and Figure6presents the results of the TP simulation.

Figure 5 .Figure 5 .Figure 5 .
Figure 5.The flow simulation results.(a) The model calibration using the 2012 data; (b) the model validation using the 2013 data.

Figure 6 .
Figure 6.The TP simulation results.(a) The model calibration using the 2012 data, and (b) the model validation using the 2013 data.

Figure 6 .
Figure 6.The TP simulation results.(a) The model calibration using the 2012 data, and (b) the model validation using the 2013 data.

Figure 7 .
Figure 7. Validation of the Vollenweider model.

Figure 7 .
Figure 7. Validation of the Vollenweider model.

Figure 8 .
Figure 8.The exceedance probability of the input TP loads.This probability is regarded as the protection level to ensure that the reservoir TP concentration remains below the target concentration under certain input loadings.

Figure 8 .
Figure 8.The exceedance probability of the input TP loads.This probability is regarded as the protection level to ensure that the reservoir TP concentration remains below the target concentration under certain input loadings.

Figure 9 .
Figure 9.The hot spots (the areas with the highest pollution potential) of the Shiman Reservoir subwatersheds: S1, S3, and S8.

Figure 9 .
Figure 9.The hot spots (the areas with the highest pollution potential) of the Shiman Reservoir subwatersheds: S1, S3, and S8.

Figure 10 .
Figure 10.The predicted TMDL performance of the Shiman Reservoir and the relationships of the watershed TP loads and reservoir TP concentrations.For example, controlling the point source (PS) pollution in hot spots could reduce the TP load by 1018 kg/year and reduce the TP concentration from 25 to 24 µg/L.Subsequent control of nonpoint source (NPS) pollution in the hot spots could reduce the TP load by 2694 kg/year and result in a TP concentration of TP 22 µg/L.The target TP concentration is 20 µg/L, and the maximum allowable load is 22,209 kg/year.

Table 2 .
The data and their sources for the Vollenweider model applied to the Shiman Reservoir, Taiwan.

Table 2 .
The input data and their sources for the Vollenweider model applied to the Shiman Reservoir, Taiwan.