Reservoir Operation Rules for Controlling Algal Blooms in a Tributary to the Impoundment of Three Gorges Dam

Since the first impoundment of Three Gorges Dam in 2003, algal blooms occur frequently in the near-dam tributaries. It is widely recognized that the impoundment-induced change in hydrodynamic condition with the lower current velocity will make the eutrophication problem even more severe when an excessive amount of nutrients is already loaded into a reservoir and/or its tributaries. Operation tests carried out by Three Gorges Corporation in 2010 point to some feasible reservoir operation schemes that may have positive impacts on reducing the algal bloom level. In our study, an attempt is made to obtain, through a numerical hydrodynamic and water quality modeling and analysis, the reservoir operation rules that would reduce the level of algal blooms in the Xiangxi River (XXR), a near-dam tributary. Water movements and algal blooms in XXR are simulated and analyzed under different scenarios of one-day water discharge fluctuation or two-week water level variation. The model results demonstrate that the reservoir operations can further increase the water exchange between the mainstream of the Three Gorges Reservoir (TGR) and the XXR tributary and thus move a larger amount of algae into the deep water where it will die. Analysis of the model results indicate that the water discharge fluctuation constituted of a lower valley-load flow and a larger flow difference for the short-term operation (within a day), the rise in water level for the medium-term operation (e.g., over 3201 weeks), and the combination of the above two for the long-term operation (e.g., over months) can be the feasible reservoir operation rules in the non-flood season for TGR.


Introduction
Development of reservoirs with large regulation capacity is an effective way to utilize water resources, produce hydropower energy and reduce flood risk.However, the hydrodynamic condition after impoundment significantly changes from the natural state with lower water velocity and longer residence time, resulting in deteriorated water quality and eutrophication problems [1,2].For the largest man-made reservoir system in the world, the Three Gorges Reservoir (TGR), according to the annual report of eco-environmental monitoring of TGR released by Executive Office of State Council TGP Construction Committee [3], more than 30% of the monitoring sections in the major tributaries are in the state of high eutrophication, and algal blooms occur in over 20 tributaries from March to October every year, such as Xiangxi River (XXR), Tongzhuang River and Daning River.Algal blooms occurred in XXR seriously affect local residents' life with turbid and stinky water [4,5].Harmful algal blooms with toxins as a prime agent of water quality deterioration also result in a loss of water for recreation and drinking [6][7][8].In order to alleviate the increasing algal bloom problems in TGR, more and more attention is being paid to the feasible solutions.
In the last decade, long-term field monitoring and researches of algal blooms in the near-dam tributaries of TGR have been carried out by government, corporations, colleges and research institutions [9][10][11][12][13][14]. Agreement is reached that the significant change in hydrodynamic condition is a compounding factor for the more frequent occurrence of algal blooms after the impoundment of Three Gorges Dam (TGD) with its effect on the transport of nutrient loads, water temperature structure and environment for algae growth [1,9,[15][16][17].An operation test is carried out by Three Gorges Corporation in 2010, and the result points to some feasible reservoir operation schemes that may have positive impacts on reducing the algal bloom level in the near-dam tributaries.With the help of measured data, influencing factors on water movement of the tributaries and their influence on algal blooms are found.Ji [18] and Yang [2] point out that the reverse density flow produced by water level fluctuation and water temperature difference both have positive impacts on water exchange and interaction between the mainstream and the tributaries.Jiang [19] verifies the result and emphasizes that thermal stratified flow by water temperature difference should be considered for the algal bloom control.Then, analysis is conducted to validate that a wave-type water release operation is beneficial in reducing algal bloom in the tributary [20].Furthermore, a feasible water level process of TGR is also demonstrated to be helpful for controlling phytoplankton blooms in the tributary bays by the analysis of mixing depth [21].Besides above, effects of environment factors and operation of TGR on algal blooms in the near-dam tributaries are also discussed.Several numerical models have been developed to predict the algal blooms in tributaries of TGR with analysis of hydrodynamic condition, nutrient concentration, temperature, and light illumination [22][23][24].Zhou [25,26] proposes that water exchange and interaction are enhanced by daily peak-load regulation between the mainstream and the XXR.Ma [27,28] describes water movement in the tributaries under multi-timescale operations of TGR and develops multi-objective operation rules based on one-dimensional numerical model in non-flood season for TGR with water quality considered.With above pieces of work, some puzzles have been clarified for a clearer understanding of the mechanism of algal blooms in TGR.However, the relationship among operation of TGR, hydrodynamics condition and algal bloom still should be further studied and quantified that will form the basis for proposing a feasible operation rule for TGR.Therefore, an attempt is made in this study to obtain, through a numerical hydrodynamic and water quality modeling and analysis, the reservoir operation rules that would reduce level of algal blooms in the Xiangxi River (XXR), a near-dam tributary.
The paper is arranged as follows.In Section 2, the study area is introduced in detail.From consideration of remarkable vertical flow structure of TGR as a river-type reservoir, a hydrodynamic and water quality model for the mainstream of TGR and tributary XXR (TGR-XXR model) is developed and calibrated based on the CE-QUAL-W2 model that has been applied successfully to many lakes, estuaries and reservoirs [29][30][31][32][33][34].In Section 3, water movement and algal bloom process in XXR are simulated and analyzed under 11 different scenarios of short-term daily water discharge fluctuations or medium-term water level variations using chlorophyll a as the indicator for algal blooms.Then, operation rules for reducing level of algal blooms are proposed in the concluding Section 4.

Study Area
The TGD is located in the middle reach of the Yangtze River in China, creating TGR with a total storage capacity of 393 × 10 8 m 3 and a regulation storage capacity of 221.5 × 10 8 m 3 .Backwater zone length of TGR is 660 km from the dam, and the reservoir has a total drainage basin area of 55,742 km 2 .Average annual runoff is 14,000 m 3 /s, and inflow ranges from 3000 m 3 /s to 30,000 m 3 /s normally within a year.The normal reservoir level is 175.0 m and the flood control level is 145.0 m.XXR is one of the largest near-dam tributaries of TGR in Hubei Province, located only 31.3 km upward from the TGD (110°25'-111°06' E, 31°04'-31°34' N), which flows southwardly into the Yangtze River.The mainstream of XXR is 33 km long and the total length of XXR is 97.4 km with a drainage basin area of 3099 km 2 .Average annual runoff is 65.5 m 3 /s, and inflow ranges from 20 m 3 /s to 200 m 3 /s normally within a year.Figure 1 shows location of TGD and sketch map of the study area (river reaches in red).

Model Introduction
CE-QUAL-W2 (W2) model, developed by the U.S. Army Corps of Engineers' Waterways Experiment Station, is a two-dimensional (longitudinal-vertical) hydrodynamic and water quality model [35].Because the model assumes lateral homogeneity, it is best suited for relatively long and narrow water bodies exhibiting longitudinal and vertical water quality gradients, similar to the TGR.Water surface elevations, velocities and temperatures can be predicted well by the hydrodynamic module of W2 model.Furthermore, 28 constituents and over 60 derived variables can be included in the water quality module, by which the whole process of algal blooms and the transport of algae can be simulated with integrity and clarity.

Model Construction
The TGR-XXR model is set up with two water bodies (the mainstream is wb_1 and XXR is wb_2) and eight reaches (seven reaches for the mainstream and one for XXR) according to the variation of bottom slope.The mainstream is 598.5 km long upwards from the TGD, and the XXR is 32 km long upwards from the confluence.Based upon the reservoir geometry and the computational efficiency and desirable accuracy, the model is simulated with 213 longitudinal segments, each 4.5 km long for wb_1 and 0.5 km long for wb_2.Depending on the local water depth, up to 107 vertical layers are used at each segment with a layer thickness of 2 m. Figure 2 shows the model segmentations, and the typical segments that could be used in the model calibration and application are also labeled in it.A maximum time step of 48 s is adopted for the model simulation.Input data for the model simulations include reservoir or river bathymetry, meteorological data (air temperature, dew point temperature, wind speed, wind direction, and percentage of cloud cover), inflow and outflow characterized by flow rate or water level, water temperature and nutrient concentration.The important calibration coefficients for the hydrodynamic module and water quality module are listed in Table 1 and their values are based upon those recommended in the past studies [35][36][37][38].Among these calibration coefficients, Manning's friction factor, maximum algal growth rate, ammonium decay rate, and nitrate decay rate are most sensitive for the model results.

Field Data Collection
The TGR-XXR model is calibrated with field data collected from 10 October 2008 to 15 June 2009 (Julian Day: 284-532).The meteorological data is the daily average of monitoring data at Zigui station obtained from the Meteorological Administration of Hubei Province.The inflow, outflow and water temperature are the daily data, which are provided by China Three Gorges Corporation obtained from the long-term field monitoring.China Three Gorges University carries out the field observation in XXR monthly and provides the constituent concentrations over the model calibration period.

Hydrodynamic Model Calibration
Simulation results of water surface elevation and water temperature are used to calibrate the hydrodynamic module of model, as shown in Figures 3 and 4. The statistic errors, mean absolute error (MAE) and Nash-Sutcliffe efficiency (NSE) [39] are listed in Table 2. Figure 3 and the statistical errors show that the modeled water surface elevation matches the measured data well, which rises from the initial 155.59-172.76m and then falls to 145.93 m.The comparisons of vertical profiles of water temperature at the three typical stations from November to April depicted in Figure 4 as well as the associated statistical errors show good model skill in reproducing the thermal structure of the water column.Thermal stratification is weak in the mainstream of TGR but relatively strong in XXR, which matches well with the actual situation of TGR.

Water Quality Model Calibration
The procedure used for the water quality module calibration searches for the best fit between the observed and predicted constituents including nitrogen (ammonia and nitrate), phosphorus, dissolved silica and chlorophyll a (Chla).The dominant species of algal blooms in XXR varies seasonally, and diatom is typically the dominant species of the spring algal bloom [3].Because algal blooms do not occur in the mainstream, the calibration results shown here are only for the XXR.Comparisons of temporal variations between modeled nutrient and observed data at a typical station (Segment 176) in the upper middle reach are shown in Figure 5, and the calibration results of Chla at different stations are shown in Figure 6, with the statistic errors listed in Table 2. Visual observations and the statistic errors indicate that the model performed well in simulating the water quality of the mainstream of TGR and XXR.

Generalization of Water Discharge Process of TGHS
According to the actual regulation rules, the Three Gorges Hydropower Station (TGHS) plays roles in both electricity base generation and peak-load regulation.For the flood-control, environment-maintenance, and other purposes, reservoir level is set to fall from the normal level of 175 m in November to the flood control level of 145 m in May or June of the following year, but it should be higher than 155.0 m before mid-May.After the impoundment of TGD, sedimentation rate at from the upstream end increases.Clear water released from TGHS and Gezhouba Hydropower Station aggravates the river bed erosion, especially in the downstream of Gezhouba Reservoir [40].In addition, low water level has a negative effect on navigation and water supply.In order to solve the problem, average daily water discharge of TGHS is constrained to be more than 5000 m 3 /s in the on-flood season.During the non-flood season as well as the algal bloom season from March to June, inflow to TGR is larger than 5000 m 3 /s, thus it is feasible to utilize the excess flow/volume to implement a reservoir regulation with a water level fluctuation for weeks.
Water discharge processes of TGHS are generalized into short-term operation (daily operation), medium-term operation (approximate two weeks) and long-term operation (for the non-flood season).For the short-term operation, due to an important role of power generation and water supply to downstream river, daily water discharge process of TGHS in non-flood season is planned according to daily inflow, release requirement and load demand with peak-load regulation.In this study, 24 h of a day are divided into four typical periods and four transition periods with a total of 48 time intervals (each lasts 0.5 h), including valley-load period (when the discharge is the lowest) (0:00-6:00), morning peak-load period (8:00-11:00), medium-load period (12:00-17:00), and evening peak-load period (19:00-22:00).Water discharge does not change over a typical period but changes linearly over a transition period.Water discharge in the evening peak-load period (peak-load flow) is equal to that in the valley-load period (valley-load flow) plus flow difference (Figure 7).Water discharges in the morning peak-load period or the medium-load period are obtained from its ratio over that in evening peak-load period, set to be 0.85 and 0.65, respectively.For the medium-term operation, the water level variation over a period of time is set to be three different scenarios: constant, rise or fall (Figure 8).Additionally, the long-term operation is planned step by step with periodical medium-term operation and short-term operation for a month approximately in the non-flood season.Here, the long-term operation rules are obtained from the inference based on the above two operation rules without individual simulation.

Scenario Construction
The calibrated model is applied to predict the influence of different timescale operation on algal blooms in XXR using Chla concentration as an indicator that is common for algae biomass evaluation [41].Forty days (Julian Day: 0-40) are set as the simulation period in order to model an entire algal bloom process probably occurred in spring (from March to May) in XXR based on the actual situation.To separate the effects of reservoir regulation from the effects of other variables, all other variables are set to have constant values.The air temperature and dew point temperature are set to be constant at 18 °C and 13 °C that are the average values used during the model calibration when the algal blooms occurred.The effects of wind and cloud cover are not considered in the model.The initial reservoir water level is 160 m, the inflow of XXR is 65.5 m 3 /s, and the inflow of TGR and discharge of TGHS are set to be corresponding to the model demand for flow or water level in scenarios (9041.875m 3 /s is used as the average daily discharge in most cases).Long-term observation shows that the density current due to water temperature difference between mainstream and XXR has seasonal variation over a year [42].The corresponding plunging position in which water flows into XXR from mainstream is different.When algal blooms occur in spring, the density current maintains that water from the mainstream plunges into the XXR from upper layers weakly and middle layers strongly.In order to maintain the same flow plunging/entrainment type as that describe above, initial and boundary water temperatures for the mainstream and XXR are set to be that as depicted in Table 3.For XXR, Zhang [43] proposes that nitrogen is mainly provided by the TGR mainstream and phosphorous is mainly from the XXR inflow.Nutrients are not the limiting factor for algal growth in TGR [44].Thus, the boundary conditions for nutrients should be conformed to the above situation and the concentrations should be set to be sufficient for algal growth.The average Chla concentration in surface water of XXR is 10 μg/L, and reduces gradually with the depth.The threshold value of high eutrophication of TGR is 32.59 ~ 62.81 μg/L [45], and the average value 47.7 μg/L is assumed as the threshold value of algal bloom occurrence in XXR in this study.The initial conditions and boundary conditions are depicted in Table 3.
Influence of two operation factors on algal bloom process are considered, including peak-load regulation type (short-term operation) and water level variation (medium-term operation) (Table 4).S1 ~ S5 scenarios are operation solutions with different daily water discharge processes, in which S2 to S5, when compared with S1, aim at detecting the influence of continuous daily peak-load regulation with the same valley-load flow but of different "flow difference" (peak-load flow) (Figure 7).S4 and S6 to S9 are operation solutions of the same daily average discharge flow but different combinations of the valley-load flow and the flow difference (Figure 7).S10 and S11 are operation solutions without daily peak-load regulation in which a rise or a fall in water level lasts for 2 weeks (Figure 8).These two scenarios can be used to show the influence of water level variation on algal blooms when compared with S1.Notes: a Distances from the confluence of the two boundaries where the algal blooms occurred; b The algal bloom in the scenario of water level fall is still in the occurrence at the end of the simulation, thus the duration is more than 24 days.

Results and Analysis
Because the 11 scenarios have the same nutrient, illumination and meteorology conditions, the differences of simulation results on algal blooms are mainly caused by different hydrodynamic conditions due to the change in reservoir operation.The simulated basic flow field of XXR indicates that water from the mainstream plunges from the upper layers into XXR and outflow of XXR discharges from the lower layers stably as a result of the water temperature difference between the two river branches.The water flowing from the upstream towards the downstream is defined as positive flow and water in the opposite direction is defined as negative flow.The negative flow can only reach the distance of 29 km from the confluence approximately within the backwater zone, because the reach of further upstream represents the river-type characteristic and the positive velocity is rather large that plunging flow cannot go forward any more.
The baseline conditions: Take S1 as an example to describe a simulated algal bloom process for this study as it does not have any intervening reservoir operation.The contour plots of Chla concentration on four different Julian days, from the beginning to the end day of the 40-day simulation period, are shown in Figure 9.The vertical distribution of velocity vectors at various locations along the longitudinal direction are also shown in the figure, which reflect the above-described flow field and also indicate the transport process of algae.(Note the distance shown on the horizontal axis is from upstream to downstream ends.)Algae primarily grow in the upper layers and emerge in the water surface which is influenced mainly by the negative flow.The negative flow pushes the algae towards the upstream of XXR, and the concentration of algae at the confluence is diluted at first since Chla concentration of the mainstream of TGR is always lower than that in XXR [2].When the algae reaches the upper boundary of the backwater zone, the combined effects of continuous pushing of negative flow and large flow of river-type reach at the farthest upstream force the algae move downwards into the deep water where they die gradually without light penetrating into.Thus, the whole transport process of algae is an anticlockwise motion.
The temporal variations of Chla concentration on the water surface at typical stations under S1 are shown in Figure 10.The longitudinal distributions of Chla on the water surface at different times during the simulation period are shown in Figure 11.Chla at confluence is diluted from the beginning, but Chla in the downstream, midstream, and upstream all go through the increase and decrease phases.At the upper reach, the increase time is longer and the peak value is larger.From the longitudinal distribution, Chla in the downstream decreases throughout the simulation but increases in the upper stream with a coverage length of high concentration shrinking before it reaches the peak value.The peak value of Chla 117.76 μg/L occurs at Julian Day 25.75 around the reach 26 km from the confluence.Then, the reach of 17.5-28.5km from the confluence is entirely covered by the algal bloom.The algal bloom in S1 lasts for 20 days all together.So, the Chla concentration along the upstream reach reflects the highest level of algal bloom in XXR, where it occurs more easily, is more severe and lasts longer.The reservoir operation: Based on the flow field and algal bloom process in S1, reservoir operation is used to accelerate the transport process of algae by enhancing the water exchange between the two river branches, in order to reduce the algal bloom level.As described above, the strength of flow field in S1, S10 and S11 is constant but fluctuates in S4 with peak-load regulation throughout a day.The contour plots of Chla concentration on Julian Day 14 (considered as the end day of operation under the scenarios of water level variation) under the four primary scenarios S1, S4, S10 and S11 are shown in Figure 12.The vertical distribution of velocity vectors at various locations along the longitudinal direction are also shown in Figure 12.These plots reflect the same flow field and algae transport process of algae, but at different magnitudes to be specifically described below.
Specific and quantitative effects of the short-term peak-load regulation: The plunging flow and discharge flow, defined as the flow rate through the segment of XXR confluence (Segment 211) can be calculated by the negative and positive velocity multiplied by river section area, which are obvious indicators of water exchange between the two river branches.Figure 13 shows the flow value of each hour during a day in S4 respectively in comparison to S1.The plunging flow and discharge flow in S1 (the baseline condition) are constant and their difference is also a constant that is equal to the XXR inflow (65.5 m 3 /s), but they fluctuate strongly in S4 with the reservoir's daily peak-load (flow) regulation.For the most time, the plunging flow in S4 is much larger and the discharging flow is at the same time much smaller than S1.For the entire day, water from the mainstream plunges into XXR and discharges out of XXR in S4 are both much larger than S1.Thus, flow fluctuation through the peak-load regulation is an effective way of strengthening the water exchange between the two branches (between the mainstream and the tributary).The velocity vectors shown in Figure 12 are at the moment (at 0 o'clock) when the plunging flow is the largest during the day.It can be obviously seen that all the velocities of negative (plunging) flow in S4 are larger than S1.Then, the more intensive plunging flow will push the algae towards the upstream and moved into deep water faster.Just as the contour plot shows, the coverage area of algae in S4 is smaller than S1 and the Chla concentration is also lower than S1.Therefore, the daily peak-load regulation is an effective way of reducing the level of algal bloom.Specific and quantitative effects of the medium-term water level variation: Compared with S1, water level variations in S10 and S11 primarily change the vertical velocity distribution of XXR in the upper layers, as shown in Figure 12.The thickness of plunging flow layer is much larger with larger velocity in S10 than S1, but the thickness is much smaller in S11 with lower velocity than S1.The discharge flows under all the three scenarios (S1, S10 and S11) have minor differences, however, a much larger amount of water from the mainstream plunges into XXR in S10 resulting in a stronger dilution effect as well as faster upstream transport of algae and leading to lower Chla concentration.On the contrary, a much smaller amount of water from the mainstream plunges into XXR in S11, resulting in the slower transport process, decreased water volume and weaker dilution effect, higher Chla concentration and larger coverage area of algae.At the same time, as indicated before, the velocity of the river-type reach is rather high that algae cannot be transported all the way upstream.In conclusion, the water level rise, similar to the peak-load regulation demonstrated above, is also effective for reducing the algal bloom level by strengthening the flow plunging into XXR to accelerate the upstream transport process of algae.
Different levels of the effects resulted from various values of reservoir operation factors: Different values of operation factors have different levels of effect on the reduction of algal blooms.The temporal variations of Chla concentration at Segment 163 of various operational scenarios are compared in Figure 14.The major indicators of algal bloom levels are listed in Table 4, and all the Chla concentrations refer to that on the water surface.The time when the algal bloom reaches the highest level is listed in the fourth column.The maximum Chla concentration and the average Chla concentration along the outbreak reach at this time are shown in Figure 15 with values listed in the fifth and ninth columns of the table.The respective reduced percentages of Chla concentration are calculated using Chla concentration of S1 as a benchmark and they are listed in the sixth and 10th columns.The coverage length and outbreak duration of algal blooms are also shown in Figure 16 with their values listed in the table.
From the whole growth process of algae shown in Figure 14, it can be observed that differences in Chla concentration under different operation scenarios reveal themselves mostly from the time when the algal bloom begins to occur.The Chla concentration decreases with the larger flow difference under the same valley-load flow, with the lower valley-load flow plus the larger flow difference under the same daily average discharge flow, and as the water level rises.The four quantitative indicators of the algal bloom level are shown in Figure 15 and Figure 16 for various operation scenarios.They demonstrate the rather similar impacts.Combined with the understanding of the hydrodynamic process with the specific values of indicators in the table, it can be observed that the scenario with a larger flow difference under the same valley-load flow has a larger plunging flow, higher velocities and a shorter transport process.Thus, the algal bloom reaches the peak value earlier but the values are lower, and the algal bloom level is also lower with a shorter coverage length and shorter outbreak duration.Owing to the same reason, the scenario with a lower valley-load flow and larger flow difference under the same daily discharge flow works better on the reduction of algae concentration.
Among all the operation factors, the effect of water level rise is the most significant.Compared with the scenario without any regulation (S1), the coverage length is shrunk by 1.5 km, the outbreak duration is shortened by 9 days (an almost one half reduction, reduced from 20 to 11 days), and the average and maximum Chla concentrations are reduced by 22.38% and 31.95%(almost one third reduction), respectively.However, the water level fall makes the algal bloom problem worse as the maximum Chla concentration on the contrary increases by 33.24% (about one third increase) and the algal bloom lasts beyond the end of the simulation (with the algal bloom duration longer than 24 days).Reservoir operation rules for reducing algal bloom level: As demonstrated through the hydrodynamic and water quality modeling analysis in this study, the daily peak-load regulation will help to reduce the level of algal bloom.It should be conducted for short-term operation with a larger flow difference under the same valley-load flow or a combination of the lower valley-load flow and the higher flow difference under the same daily average discharge flow.Also, the water level rise can be conducted for the medium-term operation (e.g., over two weeks).With both of these types of reservoir operation being favorable to the algal bloom level reduction, they can be combined together periodically as a long-term operation (e.g., monthly) in the non-flood season for TGR.Multi-objective reservoir operation: It has been demonstrated that the peak-load regulation can improve the utilization efficiency of electricity generated and the raising reservoir level is beneficial for energy storage and hydropower generation [40].Therefore, an integration of the daily operation and medium-term operation is not only a feasible way for improving the hydrodynamic conditions and water quality, but also a probable way to achieve multi-objective regulation for both hydropower generation and water environment improvement.
For the future work, with the runoff, electricity generation demand and the regulations of cascade reservoirs considered, influence of multi-time scale operation on the eutrophication and algal blooms in the tributaries can be further investigated to propose more feasible operation rules.Then, an ecology-friendly, multi-objective, and multi-scale optimal operation for TGR can be achieved.

Conclusions
The occurrence of algal blooms in near-dam tributaries of TGR in the studied non-flood season is caused by the slow water motion that provides suitable environment for algae growth, together with the unfavorable meteorological conditions (e.g., warm air, low wind) and the high nutrient load.To regulate the hydrodynamic conditions to reduce algal blooms, a two-dimensional, laterally averaged hydrodynamic and water quality model based on the CE-QUAL-W2 model for TGR and XXR is developed and applied to study the influences of reservoir operations on the algal blooms.
The formulation of reservoir operation rules is based on the current scheduling rules of TGHS as well as the scenario simulations.The reservoir operation will be helpful for water quality improvement as it will enhance water exchange between the mainstream of TGR and the XXR tributary resulting in more water from the mainstream plunging into XXR leading to the higher dilution as well as the acceleration of the upstream transport process of algae.Through the comparisons of algal bloom levels under different values of operation factors, the operation rules oriented to water environment improvement are formulated.For the short-term operation, daily peak-load regulation should be conducted with a larger flow difference under the same valley-load flow, or a combination of the lower valley-load flow and a higher flow difference under the same daily average discharge flow.For the medium-term operation (e.g., over two weeks), water level rise can be conducted.For the long-term operation (e.g., over months), the short-term flow regulation and medium-term water level rise can be combined periodically in the non-flood season for TGR.

Figure 1 .
Figure 1.Location of the Three Gorges Dam (TGD) and sketch map of the study area.

Figure 2 .
Figure 2. Model segmentations of TGR-XXR.(a) Sketch map of grid cross-cut of XXR; (b) Grid straight-cut of XXR; (c) Model grid plan and location of typical grids.(Number in the inset of Figure 2c represents the model segment number, e.g., 137 represents Segment 137).

Figure 4 .
Figure 4. Comparison of modeled and observed vertical temperature profiles.

Figure 7 .
Figure 7. Water discharge scenarios for short-term operation.

Figure 8 .
Figure 8. Four types of water level variation scenarios for medium-term operation.

Figure 9 .
Figure 9. Transport process for algae in XXR under the baseline condition S1.

Figure 10 .
Figure 10.Temporal variation of Chla concentration in water surface of XXR at typical stations on S1.

Figure 11 .
Figure 11.Longitudinal distribution of Chla concentration in water surface of XXR on S1.

Figure 12 .
Figure 12.Contour plot of Chla concentration and vector plot of longitudinal velocity at Julian Day 14.00.

Figure 13 .
Figure 13.Daily plunging flow and discharge flow at XXR confluence.

Figure 14 .
Figure 14.Temporal variations of Chla concentrations of upstream of XXR (Segment 163) under different operation factors.(a) the same valley-load flow but different flow differences; (b) the same daily water discharge but different combinations of valley-load flows and peak-load flows; (c) different water level variations.

Table 2 .
Results of statistic errors.
Figure 3.Comparison of modeled and observed water surface elevation for TGR (location before TGD site, Segment 146 in model).

Table 3 .
Initial and boundary conditions of water temperature and constituent concentrations.

Table 4 .
Comparisons of Chla concentration and algal bloom levels of the outbreak reach.