Simulating the Impacts of an Upstream Dam on Pollutant Transport : A Case Study on the Xiangjiang River , China

Water pollution incidents can cause rapid deterioration of water quality, potentially leading to the illness and death of surrounding residents. Therefore, it is imperative to assess the impact of water pollution incidents, and to understand the mechanisms of pollutant transport and transformation. To address this problem, this paper constructed a hydraulic water quality model for the lower reaches of the Xiangjiang River, China, using the hydrodynamic module and convective diffusion module of MIKE21. Six pollution incident scenarios were simulated to investigate the transport process of pollutants, as affected by an upstream dam structure, the Changsha Comprehensive Control Project dam (CCCP). Analysis of the results suggests that the CCCP plays an essential role in controlling the transport and transformation of pollutants. With the CCCP, the process of transport is weakened, and the dispersion effect is strengthened. In particular, after the construction of the CCCP, the same amount of upstream discharge leads to lower peak pollutant concentrations and longer pollutant arrival times to each waterworks’ intake, thereby alleviating the impact of water pollution incidents. Further, comparative analysis suggests that the role of CCCP is much more significant with lower discharges (e.g., during the dry season), largely due to the higher amount of water quantity within the reach.


Introduction
Water pollution incidents are becoming major environmental events in the world [1].Water pollution incidents have caused rapid deterioration of water quality, subsequently leading to the illness and death of surrounding residents, bringing severe economic losses and negatively impacting water ecological environments in both the short and long terms [2].In China, approximately 1700 water pollution events have occurred annually, and about 140 million people have been impacted by polluted drinking water supplies [3].As an example, a sudden water pollution accident occurred in Xiangjiang River, in January 2006, with the cadmium concentration in the nearby water exceeding 25 times that of the national standard, thereby causing a serious threat to the safety of drinking water.Unlike other water pollution issues, the prediction and prevention of water pollution incidents are extremely difficult due to the unpredictable nature of human activities [4].The consequences of water pollution incidents, however, are always significant, largely due to their unpredictable nature and the rapid spread of pollutants in a very short time frame.Therefore, it is imperative to simulate the impact of water pollution incidents, and to understand the mechanisms of pollutant transport and transformation.
Due to the significant impact of water pollution incidents, scientists and professionals have attempted to understand the mechanisms of pollution migration and transformation, as well as estimate the geographic extent of the affected area with various pollutant concentrations.Earlier water quality models in the world can be traced back to 1925, when Streeter and Phelps established the Streeter-Phelps model, based on a study of Ohio River pollution [5].Since then, a large number of water quality models have been developed [6,7].For water pollution incidents, in particular, He et al. (2011) developed a water pollution emergency response system for monitoring water pollution incidents in the Three Gorges Reservoir, Yangtze River, China [4].Zhang et al. (2011) explored water pollution incidents using a spatio-temporal simulation system, based on system dynamic and geographic information systems (GIS) models applied to the Songhua River, China [2].In addition to the aforementioned studies, water quality models also have the potential to be applied in the study of water pollution incidents [8,9].In particular, Cox (2003) reviewed a number of water quality models, especially QUAL2E (The Enhanced Stream Water Quality Model), MIKE-11, SIMCAT (SIMulation of CATchments), TOMCAT, ISIS, and QUASAR (Quality simulation along river systems), for estimating dissolved oxygen (DO) transportation in lowland rivers, and concluded that a model of medium complexity might be applicable for examining water quality issues [6].For water pollution incidents, however, little research has been conducted to examine pollutant transportation within a river system, although a hydrodynamic model that can provide a spatio-temporal dynamics of pollutant transport is essential [10].Furthermore, although a number of studies have examined the impact of water projects (e.g., dams) on water quality and natural habitats [11][12][13][14], to our knowledge, no research has been conducted on examining the role of dams on pollutant transportation and diffusion, especially those resulting from water pollution incidents.
To address the aforementioned problem, this research, therefore, developed a hydrodynamic module and convective diffusion module to simulate the spatio-temporal dynamics of pollutants after the occurrence of a water pollution accident, with an emphasis on the impact of the Changsha Comprehensive Control Project (CCCP) dam on downstream parts of the river.To reach such an objective, we employed MIKE zero, a commercial software package developed by the Danish Hydraulic Institute (DHI), in order to model the water quality.MIKE zero is a hydraulic model that simulates water movements dynamically in a stream or river [6,15].The model can be used to simulate sediment transport, eutrophication, water quality, and rainfall-runoff [6].The hydrodynamic module and convective diffusion module of MIKE21, part of the MIKE zero products, are suitable for simulating pollutant dispersion as a full-hydrodynamic model.Three research questions will be answered: (1) What are the geographic extents and concentrations of pollutants within a given time frame, and how do they impact the water quality at the intakes of waterworks?(2) What are the impacts of the CCCP dam, downstream, on pollutant transport and transformation?(3) If the CCCP plays a role in the process of pollutant transport and transformation, does such a role vary with different upstream discharges?

Study Area and Data
The lower reaches of the Xiangjiang River (see Figure 1), located in the province of Hunan, China, was selected as the study area.The Xiangjiang River runs through a number of large cities of Hunan (e.g., Changsha, Zhuzhou, and Xiangtan), and a number of industrial and residential wastes are major pollution sources of the Xiangjiang River.The reach starts from the Heishipu Bridge to the CCCP, which was constructed in October 2012, with a length of approximately 36.8 km, a width of 800-1600 m, and an average gradient of 0.0824‰ (see Figure 2).Within this reach, several waterworks, including the Second Waterworks, the Eighth Waterworks, the Third Waterworks, the First Waterworks, and the Fourth Waterworks, exist, from the upstream to the downstream.It is essential to simulate pollutant concentrations at the intakes of these waterworks.In order to implement the MIKE21 model, we employed river terrain information, hydrological data, and water quality information as basic data.In particular, the river terrain map of 2008 (see Figure 2) was obtained through in situ measurements, and hydrological data of 2007 and 2008 were also obtained for model calibration and validation.
Water2016, 8, 516 3 of 14 essential to simulate pollutant concentrations at the intakes of these waterworks.In order to implement the MIKE21 model, we employed river terrain information, hydrological data, and water quality information as basic data.In particular, the river terrain map of 2008 (see Figure 2) was obtained through in situ measurements, and hydrological data of 2007 and 2008 were also obtained for model calibration and validation.Water2016, 8, 516 3 of 14 essential to simulate pollutant concentrations at the intakes of these waterworks.In order to implement the MIKE21 model, we employed river terrain information, hydrological data, and water quality information as basic data.In particular, the river terrain map of 2008 (see Figure 2) was obtained through in situ measurements, and hydrological data of 2007 and 2008 were also obtained for model calibration and validation.

Model Construction
For model construction, a two dimensional hydrodynamic model was implemented using MIKE21.In particular, with the parameters and boundary conditions, the contaminant advection diffusion process associated with water movement was described by the Navier-Stokes equations and the advection diffusion equation, listed as follows: where h is the water depth, u is the average flow velocity at the x direction, v is the average flow velocity at the y direction, f is the coefficient of the Coriolis force, η is the water level, τ sx and τ sy are the wind stress components at the x and y directions, T xx , T xy , T yx , and T yy are the side stresses, p a is the air pressure, S is the point source emissions (due to the water pollution accident), D x and D y are the pollutant diffusion coefficient at the x and y directions, c is the concentration of pollutants, and finally F is the pollutants degradation coefficient.

Model Parameters and Boundary Conditions
The model entry point is located at the Heishipu Bridge, and the exit point is at the CCCP dam, with a total length of 36.8 km (see Figure 2).For the reach, we constructed the riverbed using triangular unstructured grids based on the measured topographic information.The entire reach has 41,998 grids, 21,771 grid nodes, with a minimum grid size of 50 m, and a maximum grid size of 100 m.The Manning's roughness of the sections within the reach were obtained through calibrating and validating hydraulic models with water depth, riverbed morphology, vegetation conditions, etc., as input parameters, and the resultant values were between 0.025 and 0.045.For different scenarios, the boundary condition of the upstream was set as the discharge at the Heishipu Bridge, and the downstream boundary condition was set as the water level at the dam site of the CCCP.

Model Calibration and Validation
We employed the average daily flow and water level measures in 2007 to calibrate the hydrodynamic model, and relevant in situ measurements in 2008 to validate the model.Through analysis of flows in the Xiangjiang River, and those from its tributaries, we found a simultaneous pattern, and the contributing flows from the tributaries are considered negligible.For model calibration and validation, we used the flow and water level data of Heishipu Bridge (upstream) and CCCP (downstream).On the basis of this hydrodynamic mathematical model, the water quality module can be superimposed to form a two-dimensional hydrodynamic water quality mathematical model.With this hydrodynamic model, the pollutant diffusion process can be simulated by applying the convection-diffusion equation and the degradation coefficients of individual pollutants.

Water Pollution Accident Simulation
Six different scenarios with different discharges at the upstream (Heishipu Bridge) and water levels at the downstream (CCCP) were evaluated.For discharges at the Heishipu Bridge, three scenarios were examined according to the monthly average discharges at the Xiangtan Hydrological Station (on the upstream), with values of 2400, 1300, and 500 m 3 /s (see Table 1).Water levels at the CCCP were set as the normal water level of the reservoir (29.70 m) after the construction of the CCCP, and three natural water levels (25.12 m, 23.93 m, and 22.30 m, respectively) without the project.This study assumes that, at 12:00 a.m., 1 December 2012, a spill of pollutants in the Changsha Heishipu Bridge section occurred, and that the discharge of pollutants went into the center of the Heishipu Bridge.Contamination concentration was set to 10,000 mg/L, with the maximum leakage discharge as 5 m 3 /s for a period of 30 min.We also assumed that the concentrations of other pollutants were zero; the pollutant degradation coefficient was 0, and the diffusion coefficient was 10 m 2 /s by following the guidance provided in [10].

MIKE Model
With the model parameters and boundary conditions, we constructed the hydrodynamic model using MIKE 21 module.The model was calibrated using the average daily flow and water level information at Heishipu Bridge and CCCP in 2007, and the in situ measurements using data from 2008 were for validation.Results (Figure 3) suggest that the developed model performed well in simulating water flow at the CCCP, with a relative error below 3.0%.The model performs better with a high water level and is relatively poor when the water level is low.Two factors may contribute to the discrepancy with low water levels.First, dredging activities in this reach lead to irregular shapes of riverbed with dispersed currents during the dry season, resulting in the modeled water levels being higher than the measured water levels.On the other hand, with the significant topographic variations of the riverbed, the terrain difference between the defined riverbed and the actual riverbed increases, thereby leading to modeling errors.In general, though, the estimation error is satisfactory, and the model can be employed for further simulations of the impact of water pollution incidents.

Water Pollution Accident Simulation
With the calibrated hydrodynamic model, the transport and transformation of pollutants were simulated under six different scenarios, listed in Table 1.The first three scenarios (A, B, C) were generated after the construction of the CCCP (with water level H = 29.70 m at CCCP) with different discharges at the Heishipu Bridge (Q = 2400, 1300, and 500 m 3 /s).The last three scenarios (D, E, F) were to simulate the water pollutant transport and transformation using natural waters (without the CCCP).With the pollutant quantity of 2600 m 3 at the Heishipu Bridge, the pollutant transport and transformation temporal dynamics were obtained.As an example, Figures 4 and 5 illustrate the pollutant concentration dynamics at 0.5 h, 1 h, 2 h, 4 h, 8 h, and 12 h under the scenarios A and D. Results suggest that the developed two-dimensional hydrodynamic model can clearly demonstrate the spatial-temporal trend of pollutants dynamically.This two-dimensional model is significantly better than the one-dimension model, as it can illustrate detailed information of pollutant dispersion within the river channel.As an example, under Scenario A, the simulation results show that the peak pollutant concentration is very high (approximately 10 mg/L) after 0.5 h, and this value decreases to about 7 mg/L after 1 h, 5 mg/L after 2 h, 3 mg/L after 4 h, and 1 mg/L after 8 h.Moreover, the pollutant dispersion pattern of is also very clear, and it is feasible to calculate the pollutant arrival time for a specific location, especially for waterworks intake spots.This is essential, as necessary emergency response should be taken before pollutants reach waterwork intakes.

Water Pollution Accident Simulation
With the calibrated hydrodynamic model, the transport and transformation of pollutants were simulated under six different scenarios, listed in Table 1.The first three scenarios (A, B, C) were generated after the construction of the CCCP (with water level H = 29.70 m at CCCP) with different discharges at the Heishipu Bridge (Q = 2400, 1300, and 500 m 3 /s).The last three scenarios (D, E, F) were to simulate the water pollutant transport and transformation using natural waters (without the CCCP).With the pollutant quantity of 2600 m 3 at the Heishipu Bridge, the pollutant transport and transformation temporal dynamics were obtained.As an example, Figures 4 and 5 illustrate the pollutant concentration dynamics at 0.5 h, 1 h, 2 h, 4 h, 8 h, and 12 h under the scenarios A and D. Results suggest that the developed two-dimensional hydrodynamic model can clearly demonstrate the spatial-temporal trend of pollutants dynamically.This two-dimensional model is significantly better than the one-dimension model, as it can illustrate detailed information of pollutant dispersion within the river channel.As an example, under Scenario A, the simulation results show that the peak pollutant concentration is very high (approximately 10 mg/L) after 0.5 h, and this value decreases to about 7 mg/L after 1 h, 5 mg/L after 2 h, 3 mg/L after 4 h, and 1 mg/L after 8 h.Moreover, the pollutant dispersion pattern of is also very clear, and it is feasible to calculate the pollutant arrival time for a specific location, especially for waterworks intake spots.This is essential, as necessary emergency response should be taken before pollutants reach waterwork intakes.Based on the simulation results, we calculated pollutant arrival time (see Table 2), pollutant peak concentration, and arrival time (see Table 3) for the five waterworks along the reach, including the Second Waterworks, the Eighth Waterworks, the Third Waterworks, the First Waterworks, and the Fourth Waterworks, sequentially, following the river direction (Figure 6).Table 2 indicates that the arrival time is positively associated with the distance from the Heishipu Bridge for each scenario.As an example, under Scenario A, the pollutant reaches the Second Waterworks after 0.33 h, the Eighth Waterworks after 3 h, the Third Waterworks after 3.17 h, the First Waterworks after 5.5 h, and the Fourth Waterworks after 9.67 h.Similarly, Table 3 indicates a decreasing trend of Based on the simulation results, we calculated pollutant arrival time (see Table 2), pollutant peak concentration, and arrival time (see Table 3) for the five waterworks along the reach, including the Second Waterworks, the Eighth Waterworks, the Third Waterworks, the First Waterworks, and the Fourth Waterworks, sequentially, following the river direction (Figure 6).Table 2 indicates that the arrival time is positively associated with the distance from the Heishipu Bridge for each scenario.
As an example, under Scenario A, the pollutant reaches the Second Waterworks after 0.33 h, the Eighth Waterworks after 3 h, the Third Waterworks after 3.17 h, the First Waterworks after 5.5 h, and the Fourth Waterworks after 9.67 h.Similarly, Table 3 indicates a decreasing trend of pollutant peak concentrations of waterworks from upstream to downstream.Taking Scenario A as an example, the pollutant arrives at the Second Waterworks with a peak concentration of 7.28 mg/L, at the Eighth Waterworks with a peak concentration of 1.91 mg/L, at the Third Waterworks with a peak concentration of 1.84 mg/L, at the first Waterworks with a peak concentration of 1.54 mg/L, and finally at the Fourth Waterworks with a peak concentration of 0.92 mg/L.Water2016, 8, 516 9 of 14 pollutant peak concentrations of waterworks from upstream to downstream.Taking Scenario A as an example, the pollutant arrives at the Second Waterworks with a peak concentration of 7.28 mg/L, at the Eighth Waterworks with a peak concentration of 1.91 mg/L, at the Third Waterworks with a peak concentration of 1.84 mg/L, at the first Waterworks with a peak concentration of 1.54 mg/L, and finally at the Fourth Waterworks with a peak concentration of 0.92 mg/L.In addition to the general trend of pollutant dispersion, we also performed comparative analyses of different scenarios to examine (1) the impact of the CCCP on the pollutant arrival times and concentrations; (2) the impact of discharges under CCCP scenarios; and (3) the impact of discharges with the natural water scenarios.

•
Impact of CCCP on the pollutant arrival time and concentration.
As shown in Figure 6, it can be discerned that the CCCP has a significant impact on the dispersion of the pollutant.Taking Scenarios A (with CCCP) and D (without CCCP) as examples, with the same discharge (e.g., 2400 m 3 /s), it takes a much longer time for the pollutant to reach the intakes of each waterworks, and the peak concentration of the pollutant is much lower.As an example, under Scenario A, it takes approximately 4.33 h (with a peak concentration of 1.84 mg/L) for the pollutants to reach the Third Waterworks, while it only takes 3.17 h (with a peak concentration of 2.98 mg/L) under Scenario D. Similar patterns can be discerned with different discharges (e.g., when comparing B and E, and C and F).This result is unsurprising, as the water level in the downstream is significantly higher due to the construction of the CCCP, which subsequently leads to decreased water velocity.Therefore, the transport velocity of pollutant is decreased, and the pollutant is diluted and dispersed in all directions.
• Impact of discharge under CCCP scenarios.
Figure 7 shows the pollutant concentration characteristics with different discharges (e.g., 2400 m 3 /s, 1300 m 3 /s, and 500 m 3 /s) under CCCP scenarios (e.g., water level is 29.7 m).Comparative analysis indicates that, when discharge decreases, the arrival time to waterworks intakes increases with decreased pollution concentration.Taking the Eighth Waterworks as an example, the arrival time of pollutant is 3 h with a peak concentration of 1.91 mg/L with the discharge of 2400 m 3 /s; the arrival time increases to 5 h with a peak concentration of 1.56 mg/L when the discharge decreases to 1300 m 3 /s; and, finally, the arrival time increases to 10.5 h with a peak concentration of 1.08 mg when the discharge further decreases to 500 m 3 /s.As a summary, under CCCP scenarios, with the same downstream water level (29.70 m), discharge at the upstream is negatively associated with the arrival time and positively associated with the peak concentrations of each waterwork intake.This result is reasonable, as the higher value of discharges leads to a higher velocity of pollutant transport, and a higher pollutant concentration.
• Impact of discharge with the natural water scenarios.
Under the natural water scenarios, the impacts of different discharges on pollutant transport are shown in Figure 8. Comparative analyses suggest that, when discharge decreases, the arrival time to waterwork intakes increases, although with increased pollution concentrations.Taking the Eighth Waterworks as an example, the arrival time of the pollutant is 2.17 h with a peak concentration of 3.09 mg/L when the discharge is 2400 m 3 /s; the arrival time increases to 3.17 h with a peak concentration of 3.13 mg/L when the discharge decreases to 1300 m 3 /s; and, finally, the arrival time increases to 5.33 h with a peak concentration of 5.33 mg when the discharge further decreases to 500 m 3 /s.As a summary, under the natural water scenarios, discharge is negatively associated with the arrival time, as well as the peak concentrations of each waterworks' intake.This result is likely be explained by the water quantity within the reach.That is, under natural water conditions, lower discharge in the dry season is associated with lower water quantity, thereby leading to higher peak pollutant concentrations.
Water2016, 8, 516 11 of 14  Impact of discharge with the natural water scenarios.
Under the natural water scenarios, the impacts of different discharges on pollutant transport are shown in Figure 8. Comparative analyses suggest that, when discharge decreases, the arrival time to waterwork intakes increases, although with increased pollution concentrations.Taking the Eighth Waterworks as an example, the arrival time of the pollutant is 2.17 h with a peak concentration of 3.09 mg/L when the discharge is 2400 m 3 /s; the arrival time increases to 3.17 h with a peak concentration of 3.13 mg/L when the discharge decreases to 1300 m 3 /s; and, finally, the arrival time increases to 5.33 h with a peak concentration of 5.33 mg when the discharge further decreases to 500 m 3 /s.As a summary, under the natural water scenarios, discharge is negatively associated with the arrival time, as well as the peak concentrations of each waterworks' intake.This result is likely be explained by the water quantity within the reach.That is, under natural water conditions, lower discharge in the dry season is associated with lower water quantity, thereby leading to higher peak pollutant concentrations.

Discussion
To understand the mechanism of the transport and transformation of pollutants in a river system, this study developed a hydrodynamic model based on the convective diffusion module of MIKE21 and simulated the pollutant transport under different scenarios based on the options of discharge and the construction of the CCCP.Two fundamental questions have been answered through this simulation: (1) When will the pollutants arrive at a given point (e.g., waterwork intakes)?(2) What are the peak pollutant concentrations (and when) for a given point?
Comparative analyses suggest that the CCCP plays a central role in the process of pollutant transport.With the CCCP, it takes a longer time for the pollutants to reach a waterwork intake, and the pollutant concentration is also much lower for all scenarios.This is expected; as the water level at the CCCP is raised to a normal water level (29.70 m in this study), the velocity of the water (and pollutants) decreases significantly, which subsequently leads to a longer pollutant arrival time for a specific spot, and a lower peak pollutant concentration due to dispersion.

Discussion
To understand the mechanism of the transport and transformation of pollutants in a river system, this study developed a hydrodynamic model based on the convective diffusion module of MIKE21 and simulated the pollutant transport under different scenarios based on the options of discharge and the construction of the CCCP.Two fundamental questions have been answered through this simulation: (1) When will the pollutants arrive at a given point (e.g., waterwork intakes)?(2) What are the peak pollutant concentrations (and when) for a given point?
Comparative analyses suggest that the CCCP plays a central role in the process of pollutant transport.With the CCCP, it takes a longer time for the pollutants to reach a waterwork intake, and the pollutant concentration is also much lower for all scenarios.This is expected; as the water level at the CCCP is raised to a normal water level (29.70 m in this study), the velocity of the water (and pollutants) decreases significantly, which subsequently leads to a longer pollutant arrival time for a specific spot, and a lower peak pollutant concentration due to dispersion.
The analyses associated with upstream discharge also suggest the importance of the CCCP.Although a general pattern that higher discharges will lead to shorter arrival times can be discerned for all scenarios, a contradictory trend of peak pollutant concentration with CCCP and without CCCP has been found.That is, with the CCCP, lower discharge leads to lower peak pollutant concentration, while, without the CCCP, lower discharge leads to higher peak pollutant concentration.This difference can be explained by the water quantity within the reach.Therefore, with the CCCP, the water quantity does not vary significantly, as the water level of the downstream has been set to 29.7 m.Under the natural condition, however, the water quantity has been significantly reduced during dry season; moreover, with a lower discharge at the upstream, the peak pollutant concentration is higher, largely due to lower water quantity.As a result, the role of the CCCP in alleviating water pollution is more significant with lower upstream discharges (e.g., during the dry season).
While this simulation model has satisfactorily examined the role of the CCCP dam on the transport and transformation of pollutants resulted from a water pollution accident, a real-world pollutant transport analysis is essential for model validation.In this study, however, we could not perform such an analysis due to the lack of in situ measurement data, and the scarcity of water pollution incidents.Real-world measurements and analyses, however, are essential for a better understanding of the pollutant transport and transformation processes within a river system.For the generation of the simulation model, we did not consider a number of factors that may impact pollutant transport and transformation.For example, we did not incorporate water degradation values into the analyses.For a particular water pollution accident, it is necessary to consider the effect of pollutant degradation.In addition, we did not consider the impact of effluents.Within this reach, no major effluents could be found, and this omission should not significantly affect our results.In future research, it may be necessary to consider more factors in order to build a more realistic model.

Conclusions
This paper constructed a hydraulic water quality model to simulate the transport and transformation of pollutants caused by water pollution incidents using the hydrodynamic module and convective diffusion module of MIKE21.Six scenarios, three after the construction of the CCCP and three under natural conditions, were built in order to examine the transport and transformation processes of pollutants.Especially, the spatio-temporal dynamics of pollutant distribution, as well as the arrival times and peak concentrations at waterwork intakes, were simulated and analyzed.Analysis of the results suggests two major conclusions.
First, the CCCP plays a central role in controlling the transport and transformation of pollutants.After the construction of the CCCP, the water level in the downstream area is raised, therefore reducing the velocity of water and pollutants.As a result, the process of transport is weakened, and the dispersion effect is strengthened.Therefore, the arrival time of the pollutants at each waterwork intake becomes longer with a decreased peak pollution concentration.Therefore, the construction of the CCCP reduces the adverse effect of pollutants on human beings and natural environments.
Second, although, in general, the CCCP serves a role in reducing pollutant concentrations and increasing the pollutant arrival time for each waterwork intake, such a function is more significant with lower discharges (e.g., during the dry season), largely due to the higher amount of water quantity within the reach with the CCCP.

Figure 1 .
Figure 1.Channel plan in the lower reaches of the Xiangjiang River (the left panel shows the map of Hunan Province, and the right panel illustrates the reach).

Figure 2 .
Figure 2. The riverbed terrain of the lower reaches of the Xiangjiang River (from Heishipu Bridge to the Changsha Comprehensive Control Project (CCCP)).

Figure 1 .
Figure 1.Channel plan in the lower reaches of the Xiangjiang River (the left panel shows the map of Hunan Province, and the right panel illustrates the reach).

Figure 1 .
Figure 1.Channel plan in the lower reaches of the Xiangjiang River (the left panel shows the map of Hunan Province, and the right panel illustrates the reach).

Figure 2 .
Figure 2. The riverbed terrain of the lower reaches of the Xiangjiang River (from Heishipu Bridge to the Changsha Comprehensive Control Project (CCCP)).

Figure 2 .
Figure 2. The riverbed terrain of the lower reaches of the Xiangjiang River (from Heishipu Bridge to the Changsha Comprehensive Control Project (CCCP)).

Figure 3 .
Figure 3. Simulated and measured daily runoff for model calibration and validation at CCCP ((a) Calibration; (b) Validation).

Figure 3 .
Figure 3. Simulated and measured daily runoff for model calibration and validation at CCCP ((a) Calibration; (b) Validation).

Figure 4 .
Figure 4. Pollution diffusion dynamics after the occurrence of the water pollution accident under Scenario A (Q = 2400 m 3 /s, H = 29.70 m).Figure 4. Pollution diffusion dynamics after the occurrence of the water pollution accident under Scenario A (Q = 2400 m 3 /s, H = 29.70 m).

Figure 4 .Figure 5 .
Figure 4. Pollution diffusion dynamics after the occurrence of the water pollution accident under Scenario A (Q = 2400 m 3 /s, H = 29.70 m).Figure 4. Pollution diffusion dynamics after the occurrence of the water pollution accident under Scenario A (Q = 2400 m 3 /s, H = 29.70 m).

Figure 5 .
Figure 5. Pollution diffusion dynamics after the occurrence of the water pollution accident under Scenario D (Q = 2400 m 3 /s, H = 25.12 m).

Figure 6 .
Figure 6.Pollutant concentration at different waterworks intakes along the reach between Heishipu Bridge and CCCP.Figure 6. Pollutant concentration at different waterworks intakes along the reach between Heishipu Bridge and CCCP.

Figure 6 .
Figure 6.Pollutant concentration at different waterworks intakes along the reach between Heishipu Bridge and CCCP.Figure 6. Pollutant concentration at different waterworks intakes along the reach between Heishipu Bridge and CCCP.

Figure 7 .
Figure 7. Pollutant concentrations with different discharges at different waterwork intakes along the Xiangjiang River under the scenarios with CCCP.

Figure 7 .
Figure 7. Pollutant concentrations with different discharges at different waterwork intakes along the Xiangjiang River under the scenarios with CCCP.
(a).The Second Water Intake (b).The Eighth Water Intake (c).The Third Water Intake (d).The First Water Intake (e).The Fourth Water Intake

Figure 8 .
Figure 8. Pollutant concentrations with different discharges at different waterwork intakes along the Xiangjiang River under natural water conditions.

Figure 8 .
Figure 8. Pollutant concentrations with different discharges at different waterwork intakes along the Xiangjiang River under natural water conditions.

Table 1 .
Different Scenarios in flow and water level conditions.
Note: * For water level calculation, 1985 National Height datum was employed.

Table 2 .
Pollutant arrival time at the intakes of waterworks.

Table 3 .
Pollutant peak concentration and arrival time at the intakes of waterworks.