Emission Control in River Network System of the Taihu Basin for Water Quality Assurance of Water Environmentally Sensitive Areas

As pollution incidents frequently occurred in the functional water areas of the Taihu Basin, Yangtze Delta, effective emission control to guarantee water quality in the Taihu Basin became the priority for environmental management. In this study, a new total emission control (TEC) method was proposed with an emphasis on the concept of water environmentally sensitive areas (WESAs). This method was verified in Wujiang District and the techniques can be concluded in three steps: (1) a 1-D mathematical model for the study area was established and the model was calibrated using field measurement data; (2) based on an analysis of administrative planning and regulations, WESAs were identified as the main controlling objectives for emission control calculations. The weighting coefficient of local pollution sources was investigated to discuss the effectiveness of TEC on water quality improvement at WESAs; and (3) applying the river network mathematical model, water quality along the river segments was simulated under different pollution control plans. The results proved the effectiveness of TEC in the study area and indicated that a 14.6% reduction in the total amount of ammonia-nitrogen (NH3-N), as well as a 31.1% reduction in the total amount of chemical oxygen demand (CODcr), was essential in order to meet the water quality standard in the WESAs.


Introduction
The Taihu Basin, located in the estuary of the Yangtze River Delta, is the core economic area in East China.As it is experiencing rapid economic development and industrialization, rivers are overloaded with pollutants, resulting in deterioration of water quality and potential threats to the safety of functional water zones, especially for drinking water sources [1][2][3].As such, efficient pollution control in this complex river network area is the key to water environment management of the Taihu Basin.
Total emission control (TEC), firstly proposed in the 1970s, has been applied in basins all over the world as an effective method for water quality management [4][5][6].In Article 18 of the 2008 Law on the Prevention and Control of Water Pollution promulgated by the National People's Congress of China, the TEC system was established as an important mechanism to shift the incentives of subnational authorities.Total emissions of regulated pollutants are allocated by the national government in order to achieve desired water environment quality, and the TEC system authorizes the imposition of significant penalties to discipline authorities who discharge higher levels than the allocated amount [7].While the overall goal of TEC is to maintain a healthy water environment, the controlling objectives and techniques of TEC vary differently according to the target management areas.A large number of research studies have examined TEC in the Taihu Basin.Generally, water quality in either all river segments or several cross-sections of the river networks were identified as the control factors for TEC [8][9][10].However, the control based on all river segments was fairly strict and the implementation was quite difficult.Oppositely, it was comparably difficult to ensure a sufficiently high level of water quality in important functional or environmentally sensitive river segments if only very limited or unrepresentative sections were identified.In addition, the effectiveness of local TEC has rarely been evaluated in previous studies.The effect from the reduction of local emissions on water quality improvement may be limited in many circumstances when non-local sources contribute most of the pollutants in study area.If so, the estimated local emission reduction rate may fail to achieve the relevant water quality standards.
The concept of environmentally sensitive areas (ESAs) was first proposed in 2008.According to the "Environmental Impact Assessment of Construction Project and Classification Management List" [11], ESAs include many different types of natural and cultural protection areas, as well as some important functional areas, which are particularly sensitive to some types of contaminants or ecological impacts.In recent research of the Taihu Basin river networks, water environmentally sensitive areas (WESAs) refer to special river areas, including all of the river protection zones proposed by officials and key river segments with important ecological and water quality functions.
In this study, a new method of TEC based on the WESAs in Wujiang District was proposed.The overall objective is to estimate the local pollution reduction rate in order to reach desired water quality standards in the study area.Specifically, the research focused on the following objectives: (i) validate the 1-D river network model for the study area; (ii) explore the effectiveness of the TEC for Wujiang; and (iii) estimate the feasibility of the potential pollution control plans in the research area.It is expected that this method can provide much information and guidance for the government in environmental management.

Study Area
Wujiang District is located in the south of Jiangsu Province.It also neighbors two other highly developed areas, Zhejiang province and Shanghai Municipality (Figure 1).The total population for this district is over 1.2 million.Wujiang currently ranks as one of the most economically successful cities in China and the major industries include electronics, printing and dyeing, and textile manufacturing, which produce large amounts of sewage.The surface water area of Wujiang District is approximately 267.1 km 2 , accounting for 22.7% of the total area.The river networks in this district are mainly composed by three rivers: Taipu River, Grand River, and Lanxi River (Figure 1).The Taipu river flows easterly into Shanghai and the other two rivers flow southward into Zhejiang province [12,13].

Data Collection
To investigate the status of the river networks, 15 cross-sections at the major rivers were selected for field measurement (Figure 1).Hydrology and water chemistry measurement were conducted twice in both March 2013 and April 2014, with three days for each field measurement.Hydrological parameters included cross-section discharge, water depth, and velocity.Water chemistry parameters included chemical oxygen demand (CODCr) and ammonia-nitrogen (NH3-N), which are both typical parameters to reflect water quality in the Taihu Basin, as guided in the Jiangsu provincial "Twelve-Five" comprehensive program of energy conservation and emission reduction [14].Sampling and measurement techniques for hydrology/water chemistry were based on the "Technical Specification Requirements for Monitoring of Surface Water" (HJ/T91-2002) [15].Specifically, water depth was evenly recorded at five points of each cross-section and the horizontal velocity was measured with a flow meter from surface to bottom at each point.Discharge at the cross-section was derived as a product of velocity and cross-section area.Water samples were collected with a 2 L plastic water sampler at the surface and bottom for chemical analysis.Samples were immediately delivered to Wujiang Environment Monitoring Center for Measurement.
According to Environmental Water Quality Standards for Surface Water, promulgated in 2002 [16], water quality conditions can be classified as six grades, from best to worst: grade I, II, III, IV, V and V-.The corresponding criteria for each grade are presented in Table 1.

Data Collection
To investigate the status of the river networks, 15 cross-sections at the major rivers were selected for field measurement (Figure 1).Hydrology and water chemistry measurement were conducted twice in both March 2013 and April 2014, with three days for each field measurement.Hydrological parameters included cross-section discharge, water depth, and velocity.Water chemistry parameters included chemical oxygen demand (COD Cr ) and ammonia-nitrogen (NH 3 -N), which are both typical parameters to reflect water quality in the Taihu Basin, as guided in the Jiangsu provincial "Twelve-Five" comprehensive program of energy conservation and emission reduction [14].Sampling and measurement techniques for hydrology/water chemistry were based on the "Technical Specification Requirements for Monitoring of Surface Water" (HJ/T91-2002) [15].Specifically, water depth was evenly recorded at five points of each cross-section and the horizontal velocity was measured with a flow meter from surface to bottom at each point.Discharge at the cross-section was derived as a product of velocity and cross-section area.Water samples were collected with a 2 L plastic water sampler at the surface and bottom for chemical analysis.Samples were immediately delivered to Wujiang Environment Monitoring Center for Measurement.
According to Environmental Water Quality Standards for Surface Water, promulgated in 2002 [16], water quality conditions can be classified as six grades, from best to worst: grade I, II, III, IV, V and V-.The corresponding criteria for each grade are presented in Table 1.Pollutants to the river networks of the Taihu Basin are mainly from the discharge of industrial, residential, livestock, and farming sewage.According to the statistics of emission investigation from thousands of sites from the Jiangsu Provincial Environmental Monitoring Center in 2013, the total pollutants discharge from these sources, mainly including COD Cr and NH 3 -N, can be roughly estimated.The average annual emissions of COD Cr and NH 3 -N from all the pollution source types were up to 26,791.4 t/a and 1684.2 t/a, respectively.The contribution proportion from different sources are displayed in Figure 2. The results illustrated that industry and urban residential sewage were dominant in the contributions of COD Cr emissions, and industry and rural residential sewage were the main sources for NH 3 -N.
Sustainability 2017, 9, 301 4 of 13 pollutants discharge from these sources, mainly including CODCr and NH3-N, can be roughly estimated.The average annual emissions of CODCr and NH3-N from all the pollution source types were up to 26,791.4 t/a and 1684.2 t/a, respectively.The contribution proportion from different sources are displayed in Figure 2. The results illustrated that industry and urban residential sewage were dominant in the contributions of CODCr emissions, and industry and rural residential sewage were the main sources for NH3-N.

1-D Mathematical River Network Model
Mathematical models have become important tools in water environmental research.Currently, a sizeable number of models, such as MIKE developed by DHI [17], the QUAL series model developed by US EPA [18], and EFDC developed by VIMS [19], have been widely used for water environment simulation.Among them, MIKE has shown favorable functionality and strong applicability in a number of river network studies.For example, Zeng et al. used MIKE11 to simulate water quality trends in the Wenrui river networks [20].Lian et al. also applied this model to predict the hydrology and water quality of the river networks in Fuzhou City [21].In our study, the MIKE11 model was applied to simulate the hydrodynamic status and transport of contaminants in Wujiang District.
In MIKE11, the Saint-Venant equation was used to simulate the hydrology of the river networks: where, q is the discharge of a tributary; Q is the discharge of the main stream; A, B, and Z represent cross-sectional area, width, and depth, respectively; VX is the flow velocity along the river in x direction; K is the volume correction coefficient; and α is the momentum correction coefficient.The equations above were discretized with the four-point linear implicit scheme.The linear water elevation and velocity were calculated based on the matrix mark method.
Water quality was predicted under the mass balance equation as below: where, C is the concentration of the substance, u is the average stream velocity, Ex is convective diffusion coefficient, and K is the attenuation coefficient.
The convective diffusion coefficient can be estimated through the empirical formula:

1-D Mathematical River Network Model
Mathematical models have become important tools in water environmental research.Currently, a sizeable number of models, such as MIKE developed by DHI [17], the QUAL series model developed by US EPA [18], and EFDC developed by VIMS [19], have been widely used for water environment simulation.Among them, MIKE has shown favorable functionality and strong applicability in a number of river network studies.For example, Zeng et al. used MIKE11 to simulate water quality trends in the Wenrui river networks [20].Lian et al. also applied this model to predict the hydrology and water quality of the river networks in Fuzhou City [21].In our study, the MIKE11 model was applied to simulate the hydrodynamic status and transport of contaminants in Wujiang District.
In MIKE11, the Saint-Venant equation was used to simulate the hydrology of the river networks: where, q is the discharge of a tributary; Q is the discharge of the main stream; A, B, and Z represent cross-sectional area, width, and depth, respectively; V X is the flow velocity along the river in x direction; K is the volume correction coefficient; and α is the momentum correction coefficient.The equations above were discretized with the four-point linear implicit scheme.The linear water elevation and velocity were calculated based on the matrix mark method.
Water quality was predicted under the mass balance equation as below: where, C is the concentration of the substance, u is the average stream velocity, E x is convective diffusion coefficient, and K is the attenuation coefficient.
The convective diffusion coefficient can be estimated through the empirical formula: where, V is the velocity predicted from the hydrodynamic component of the model, and a and b are empirical parameters.The 1-D plain river network model for Wujiang District included 19 river segments in total, with a length from 8.2 km to 39.6 km.The dimensions for each segment was obtained from Wujiang Environmental Monitoring Center and were correctly implemented into the model as geometry data.The model consisted of 104 h-points and 65 Q-points (Figure 3).Surface elevation and water quality (e.g., COD Cr and NH 3 -N) parameters were predicted on h-points while river discharge was simulated on Q-points.

Weighting Coefficient
All of the emission sources that influence water quality in WESAs can be divided into two sections based on the scope of TEC: local and non-local pollution sources.Non-local pollution sources contribute major pollutants from outside of the study area, and local pollution sources directly emit contaminants into rivers within the study area.In general, only local pollution sources are chosen as control targets in environmental management.
To guarantee the effectiveness of TEC on the water quality improvement for WESAs, it is crucial to evaluate the effect of local pollution sources on target water bodies.The effectiveness of local emission control in the study area can be quantified by calculating the weighting coefficient for local pollution sources.The weighting coefficient of local pollution sources was first proposed by Zhu et al. [10] in the study of the Hangjiahu River networks, and it is mainly to evaluate whether the TEC in the scope of the control area had effective influence on water quality improvement of WESAs.The equation can be expressed as follows: where, is the weighting coefficient of local pollution sources; is the average water

Weighting Coefficient
All of the emission sources that influence water quality in WESAs can be divided into two sections based on the scope of TEC: local and non-local pollution sources.Non-local pollution sources contribute major pollutants from outside of the study area, and local pollution sources directly emit contaminants into rivers within the study area.In general, only local pollution sources are chosen as control targets in environmental management.
To guarantee the effectiveness of TEC on the water quality improvement for WESAs, it is crucial to evaluate the effect of local pollution sources on target water bodies.The effectiveness of local emission control in the study area can be quantified by calculating the weighting coefficient for local pollution sources.The weighting coefficient of local pollution sources was first proposed by Zhu et al. [10] in the study of the Hangjiahu River networks, and it is mainly to evaluate whether the TEC in the scope of the control area had effective influence on water quality improvement of WESAs.The equation can be expressed as follows: where, α i is the weighting coefficient of local pollution sources; c i is the average water concentration of the target water body, under the condition that local pollutants were considered; c 0 is the average water concentration of the target water body, under the condition that local pollutants were ignored.c i and c 0 were derived from the validated 1-D river network model.

Results
Field water quality measurement at selected sites in March 2013 displayed little temporal variation.Averaged monitoring data is presented in Table 2. Generally, according to Environmental Water Quality Standards for Surface Water, promulgated in 2002 [14], water quality conditions for the study river segments were between grade III and IV.Poor water quality sites (grade IV or worse) occupied a portion up to 47% and most of these sites are located in south of the study area.Water quality deterioration from upstream (e.g., site 1, site 7) to downstream (e.g., site 6, site 10) was clearly noted due to the transport of pollutants.As presented in the table for the field measurement of the Grand Canal, COD Cr concentrations along the Grand Canal have increased by 55%, and NH 3 -N has increased by 94%.Water quality for Taipu River was comparably good since its upstream directly connects to the west Lake Taihu where the water quality condition was fairly good.Particularly, NH 3 -N concentration of upstream sites 7 and 8 were of very low values.However, a five-fold increase of NH 3 -N concentration was noted from sites 8 to 9 after being mixed with the water from the Grand Canal.

Model Calibration and Validation
To build the mathematical hydrodynamic model for Wujiang District, model parameters were calibrated using a trial and error method through which the parameters were continually adjusted until appropriate numerical results were derived.This can efficiently ensure a stable calculation process and good agreement between measured and simulated results.In our study for the plain river networks, the bottom roughness and degradation coefficient were determined through calibrating the hydrodynamic and water quality models, respectively.
Field hydrology data were obtained from 17 to 19 March 2013 in Wujiang District.The canal roughness of the river networks in Wujiang District was calibrated in a range of 0.023-0.032.According to the study of the river networks in the lower reach of the Yangtze River of China by Han et al. the bottom roughness for this area was estimated between 0.019 and 0.031 [22] and, hence, our calibrated roughness was completely reasonable.As is shown in Figure 4, the simulated values were close to the measured ones at river segment 2, 4, and 9. Relative errors for discharge at these three sites were between 3.9% and 18.2%.Average errors for the entire river networks was 12%. ) derived in the study of north Lake Taihu by Wang et al. [23].As is illustrated in Figures 5 and 6, the relative errors for CODCr were between 3.8% and 14.6%, while the results for NH3-N were between 5.2% and 25%.Average errors for the study area was 10.3% for CODCr and 18.4% for NH3-N.The parameters of the model were further validated with field measurement data during April 2014 in order to guarantee the stability and precision of the model.The validation results for the river discharge are shown in Figure 7.The results from the model was in good agreement with the measured values, with relative errors between 3.0% and 14.8% at the three major sites.An average error of 8% for the whole river networks suggested a confident validity of the model performance in hydrology predictions.) derived in the study of north Lake Taihu by Wang et al. [23].As is illustrated in Figures 5 and 6, the relative errors for COD Cr were between 3.8% and 14.6%, while the results for NH 3 -N were between 5.2% and 25%.Average errors for the study area was 10.3% for COD Cr and 18.4% for NH 3 -N.Based on the field water quality measurement (28-30 March 2013), the degradation parameters (CODCr and NH3-N) in the model were calibrated in a range of 0.09-0.15d −1 and 0.07-0.13d −1 , respectively.These values were close to the results (0.09-0.11 d −1 and 0.05-0.07d −1 ) derived in the study of north Lake Taihu by Wang et al. [23].As is illustrated in Figures 5 and 6, the relative errors for CODCr were between 3.8% and 14.6%, while the results for NH3-N were between 5.2% and 25%.Average errors for the study area was 10.3% for CODCr and 18.4% for NH3-N.The parameters of the model were further validated with field measurement data during April 2014 in order to guarantee the stability and precision of the model.The validation results for the river discharge are shown in Figure 7.The results from the model was in good agreement with the measured values, with relative errors between 3.0% and 14.8% at the three major sites.An average error of 8% for the whole river networks suggested a confident validity of the model performance in hydrology predictions.Based on the field water quality measurement (28-30 March 2013), the degradation parameters (CODCr and NH3-N) in the model were calibrated in a range of 0.09-0.15d −1 and 0.07-0.13d −1 , respectively.These values were close to the results (0.09-0.11 d −1 and 0.05-0.07d −1 ) derived in the study of north Lake Taihu by Wang et al. [23].As is illustrated in Figures 5 and 6, the relative errors for CODCr were between 3.8% and 14.6%, while the results for NH3-N were between 5.2% and 25%.Average errors for the study area was 10.3% for CODCr and 18.4% for NH3-N.The parameters of the model were further validated with field measurement data during April 2014 in order to guarantee the stability and precision of the model.The validation results for the river discharge are shown in Figure 7.The results from the model was in good agreement with the measured values, with relative errors between 3.0% and 14.8% at the three major sites.An average error of 8% for the whole river networks suggested a confident validity of the model performance in hydrology predictions.The parameters of the model were further validated with field measurement data during April 2014 in order to guarantee the stability and precision of the model.The validation results for the river discharge are shown in Figure 7.The results from the model was in good agreement with the measured values, with relative errors between 3.0% and 14.8% at the three major sites.An average error of 8% for the whole river networks suggested a confident validity of the model performance in hydrology predictions.The water quality results between the model and field observations are shown in Figures 8 and  9.It was noted that the model-predicted results matched well with the measured ones, with relative errors between 3.2% and 16%.The average relative error for the entire river networks was 9.2%, which indicated that the transport process of pollution could be well represented by our model as well.

Effectiveness Analysis
"Taihu Lake Basin Water Environmental Function Zoning" approved by the State Council of China and "Regional planning of protection ecological red line in Jiangsu Province" approved by the Peoples' Government of Jiangsu Province [24,25] have both made regulations about functional and sensitive river segments in study area.According to the regulations by government, four river segments with important functions or sensitive properties (Table 3) in Wujiang District were identified as WESAs.The control factors for TEC were established based on achieving water quality standards in these four WESAs.The locations of the WESAs are shown in Figure 10.
WESA1 and WESA4 belong to drinking water protection areas on the Taipu River, and the water quality at these two areas is required to reach grade II-III standard according to "Taihu Lake Basin Water Environmental Function Zoning".WESA 2 belongs to an interregional buffer area of the Former Grand Canal and WESA3 is part of the interregional buffer area of the Lanxi River.Both The water quality results between the model and field observations are shown in Figures 8 and 9.It was noted that the model-predicted results matched well with the measured ones, with relative errors between 3.2% and 16%.The average relative error for the entire river networks was 9.2%, which indicated that the transport process of pollution could be well represented by our model as well.The water quality results between the model and field observations are shown in Figures 8 and  9.It was noted that the model-predicted results matched well with the measured ones, with relative errors between 3.2% and 16%.The average relative error for the entire river networks was 9.2%, which indicated that the transport process of pollution could be well represented by our model as well.

Effectiveness Analysis
"Taihu Lake Basin Water Environmental Function Zoning" approved by the State Council of China and "Regional planning of protection ecological red line in Jiangsu Province" approved by the Peoples' Government of Jiangsu Province [24,25] have both made regulations about functional and sensitive river segments in study area.According to the regulations by government, four river segments with important functions or sensitive properties (Table 3) in Wujiang District were identified as WESAs.The control factors for TEC were established based on achieving water quality standards in these four WESAs.The locations of the WESAs are shown in Figure 10.
WESA1 and WESA4 belong to drinking water protection areas on the Taipu River, and the water quality at these two areas is required to reach grade II-III standard according to "Taihu Lake Basin Water Environmental Function Zoning".WESA 2 belongs to an interregional buffer area of the Former Grand Canal and WESA3 is part of the interregional buffer area of the Lanxi River.Both  The water quality results between the model and field observations are shown in Figures 8 and  9.It was noted that the model-predicted results matched well with the measured ones, with relative errors between 3.2% and 16%.The average relative error for the entire river networks was 9.2%, which indicated that the transport process of pollution could be well represented by our model as well.

Effectiveness Analysis
"Taihu Lake Basin Water Environmental Function Zoning" approved by the State Council of China and "Regional planning of protection ecological red line in Jiangsu Province" approved by the Peoples' Government of Jiangsu Province [24,25] have both made regulations about functional and sensitive river segments in study area.According to the regulations by government, four river segments with important functions or sensitive properties (Table 3) in Wujiang District were identified as WESAs.The control factors for TEC were established based on achieving water quality standards in these four WESAs.The locations of the WESAs are shown in Figure 10.
WESA1 and WESA4 belong to drinking water protection areas on the Taipu River, and the water quality at these two areas is required to reach grade II-III standard according to "Taihu Lake Basin Water Environmental Function Zoning".WESA 2 belongs to an interregional buffer area of the Former Grand Canal and WESA3 is part of the interregional buffer area of the Lanxi River.Both

Effectiveness Analysis
"Taihu Lake Basin Water Environmental Function Zoning" approved by the State Council of China and "Regional planning of protection ecological red line in Jiangsu Province" approved by the Peoples' Government of Jiangsu Province [24,25] have both made regulations about functional and sensitive river segments in study area.According to the regulations by government, four river segments with important functions or sensitive properties (Table 3) in Wujiang District were identified as WESAs.The control factors for TEC were established based on achieving water quality standards in these four WESAs.The locations of the WESAs are shown in Figure 10.
WESA1 and WESA4 belong to drinking water protection areas on the Taipu River, and the water quality at these two areas is required to reach grade II-III standard according to "Taihu Lake Basin Water Environmental Function Zoning".WESA 2 belongs to an interregional buffer area of the Former Grand Canal and WESA3 is part of the interregional buffer area of the Lanxi River.Both WSEA 2 and 3 are required to reach grade III water quality standards.Details of these WESAs are shown in Table 3. WSEA 2 and 3 are required to reach grade III water quality standards.Details of these WESAs are shown in Table 3.It is important to explore the feasibility of TEC in the scope of Wujiang District by comparing the pollution contributions from local and non-local sources before the control rate is estimated.A low local weighting coefficient suggested local TEC might not work to achieve better water quality and an alternate method in improving water environmental condition was desired.The model was simulated in two scenarios: with and without local emission sources.The weighting coefficient was then calculated following Equation (4).The coefficient results for 4 WESAs are shown in Table 4.The average CODCr weighting coefficients of local pollution sources on the four WESAs was over 60.7%, which implies the local pollution source for CODCr has a relatively stronger influence compared to non-local sources for Wujiang District.It is important to explore the feasibility of TEC in the scope of Wujiang District by comparing the pollution contributions from local and non-local sources before the control rate is estimated.A low local weighting coefficient suggested local TEC might not work to achieve better water quality and an alternate method in improving water environmental condition was desired.The model was simulated in two scenarios: with and without local emission sources.The weighting coefficient was then calculated following Equation (4).The coefficient results for 4 WESAs are shown in Table 4. would affect the existing water quality even in the unfavorable hydraulic condition.As for that, field water quality measurement data in March, 2013 was directly applied as the initial condition for model simulation.The model was run separately under the three proposed pollution control plans.Since the length of river reaches in the WESAs varied from 4.2 km to 17.5 km, water quality changes along the reaches under the effect of emissions from both riversides.Based on the water environmental capacity calculation method [29], the average water quality at the computational nodes in the WESAs was used to evaluate the status of the WESAs.
Different water quality results for the WESAs under relevant emission control conditions are presented in Figure 11a The bar chart indicates the changes in water quality for different WESAs.After analysis, the results were as follows: (1) Scheme I: under the fixed hydrological conditions of a P = 90% guarantee rate, maintaining the current emission status without any reduction, WESA1, WESA2, and WESA3 all exceeded the corresponding standards of water quality; therefore, the current total emission amounts should be immediately reduced; (2) Scheme II: according to regional water management planning, it has been recommended that local CODCr and NH3-N pollution load should be reduced by 5.4% and 12.7%, respectively.After the implementation of this reduction planning, water quality for these four WESAs was improved to some extent.However, water quality of WESA2 and WESA3, especially in terms of NH3-N, still cannot meet the requirements.A further reduction would be necessary on the basis of integrated water resource planning in Wujiang District; (3) Scheme III: the total emissions of CODCr and NH3-N in Wujiang District were required to be reduced by 14.6% and 31.1%.The reduction result showed all the water quality in environmentally sensitive areas reached the corresponding requirements, except for NH3-N in WESA3.However, a 9.4% exceedance in WESA3 appears to be in an acceptable range.This scheme can basically reach the requirements of TEC in Wujiang District.

Discussion
In our study, WESAs were, for the first time, selected as control targets in the new TEC method which not only efficiently guarantee water quality conditions of important river segments, but also reasonably suggested an emission reduction rate and potentially maintained economic development.Moreover, effectiveness analysis was conducted as a pre-condition of the new method and a high weighting coefficient for Wujiang District was derived.Hence, our results indicated that, with a similar density of pollution emission sources, the new method could be sufficiently applied for an area larger than Wujiang (10,000 km 2 ) in the Taihu Basin.
The established 1-D model displayed good performance in reproducing hydrodynamics and water quality conditions of the study river networks.Relative errors between the model simulations and field measurements were comparably small for the model validation.However, it should be noted that river sediment in the study area was not directly considered in the model, given that little field data for sediment nutrient content was available.During the past several decades with high economic development, nutrients were significantly accumulated in the bottom sediment and it was The bar chart indicates the changes in water quality for different WESAs.After analysis, the results were as follows: (1) Scheme I: under the fixed hydrological conditions of a P = 90% guarantee rate, maintaining the current emission status without any reduction, WESA1, WESA2, and WESA3 all exceeded the corresponding standards of water quality; therefore, the current total emission amounts should be immediately reduced; (2) Scheme II: according to regional water management planning, it has been recommended that local COD Cr and NH 3 -N pollution load should be reduced by 5.4% and 12.7%, respectively.After the implementation of this reduction planning, water quality for these four WESAs was improved to some extent.However, water quality of WESA2 and WESA3, especially in terms of NH 3 -N, still cannot meet the requirements.A further reduction would be necessary on the basis of integrated water resource planning in Wujiang District; (3) Scheme III: the total emissions of COD Cr and NH 3 -N in Wujiang District were required to be reduced by 14.6% and 31.1%.The reduction result showed all the water quality in environmentally sensitive areas reached the corresponding requirements, except for NH3-N in WESA3.However, a 9.4% exceedance in WESA3 appears to be in an acceptable range.This scheme can basically reach the requirements of TEC in Wujiang District.

Discussion
In our study, WESAs were, for the first time, selected as control targets in the new TEC method which not only efficiently guarantee water quality conditions of important river segments, but also reasonably suggested an emission reduction rate and potentially maintained economic development.Moreover, effectiveness analysis was conducted as a pre-condition of the new method and a high weighting coefficient for Wujiang District was derived.Hence, our results indicated that, with a similar density of pollution emission sources, the new method could be sufficiently applied for an area larger than Wujiang (10,000 km 2 ) in the Taihu Basin.
The established 1-D model displayed good performance in reproducing hydrodynamics and water quality conditions of the study river networks.Relative errors between the model simulations and field measurements were comparably small for the model validation.However, it should be noted that river sediment in the study area was not directly considered in the model, given that little field data for sediment nutrient content was available.During the past several decades with high economic development, nutrients were significantly accumulated in the bottom sediment and it was found that internal sediment release can contribute as high as 25% of the nutrients in water bodies under certain hydrodynamic conditions [30].Even though calibration on the degradation coefficient for our model can reflect the role of sediment as sink or source to some extent, it is still desired to incorporate the sediment module into the model, especially for long-term simulations with highly varied hydrodynamic conditions.

Conclusions
As pollution of river networks is a serious environmental problem all over the world, total emission control has become the preferred choice for management of water quality by many governments, and achieved a great deal of success in environmental water management.To maintain the health status of important functional areas of rivers, strict restrictions and regulations should be applied in pollution emissions control.In addition, a reasonable reduction rate should be carefully determined.
In this paper, four WESAs in Wujiang District were identified to be the control factors of TEC.Based on weighting coefficient analysis methods for local pollution sources, the average weighting coefficient of local emissions (COD Cr ) was determined as 60.7%, and the corresponding value for NH 3 -N reached 72.5%, proving that TEC of local pollution sources in Wujiang is effective.Finally, simulations of the water environment in Wujiang District were conducted using a 1-D model under different pollution conditions, according to proposed pollution control plans.Analysis showed that to meet the water quality requirements of the WESAs, the total amount of NH 3 -N should be reduced by 14.6% and the total amount of COD Cr should be reduced by 31.1%.
The new method provided a practical way to plan a both effective and reasonable pollution reduction program in a river network of the Taihu Basin.This method, potentially applicable in other river networks, can be of great value in helping governments design water management policies.

Figure 1 .
Figure 1.River networks and distribution of monitoring sites in Wujiang District.

Figure 1 .
Figure 1.River networks and distribution of monitoring sites in Wujiang District.

Figure 2 .
Figure 2. Emissions composition of different sources in the study area for 2013.
The model time step was set as 60 s, considering the model stability and calculation efficiency.Measured river discharge data at the boundary node was used to drive the model.An emission source in the study area was added at the nearest node point as external nutrient loading.Model calibration was conducted with field data in March 2013, while the validation of the model was done through comparing with the field measurement in April 2014.Sustainability 2017, 9, 301 5 of 13 where, V is the velocity predicted from the hydrodynamic component of the model, and a and b are empirical parameters.The 1-D plain river network model for Wujiang District included 19 river segments in total, with a length from 8.2 km to 39.6 km.The dimensions for each segment was obtained from Wujiang Environmental Monitoring Center and were correctly implemented into the model as geometry data.The model consisted of 104 h-points and 65 Q-points (Figure 3).Surface elevation and water quality (e.g., CODCr and NH3-N) parameters were predicted on h-points while river discharge was simulated on Q-points.The model time step was set as 60 s, considering the model stability and calculation efficiency.Measured river discharge data at the boundary node was used to drive the model.An emission source in the study area was added at the nearest node point as external nutrient loading.Model calibration was conducted with field data in March 2013, while the validation of the model was done through comparing with the field measurement in April 2014.

Figure 3 .
Figure 3. Diagram of the study area for the 1-D model.

Figure 3 .
Figure 3. Diagram of the study area for the 1-D model.

Sustainability 2017, 9 , 301 7 of 13 Figure 4 .
Figure 4. Comparison of discharge at three river segments between field measurements and the model.

Figure 5 .
Figure 5.Comparison of CODCr at three river segments between field measurements and the model.

Figure 6 .
Figure 6.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 4 .
Figure 4. Comparison of discharge at three river segments between field measurements and the model.

Figure 4 .
Figure 4. Comparison of discharge at three river segments between field measurements and the model.

Figure 5 .
Figure 5.Comparison of CODCr at three river segments between field measurements and the model.

Figure 6 .
Figure 6.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 5 .
Figure 5.Comparison of COD Cr at three river segments between field measurements and the model.

Figure 4 .
Figure 4. Comparison of discharge at three river segments between field measurements and the model.

Figure 5 .
Figure 5.Comparison of CODCr at three river segments between field measurements and the model.

Figure 6 .
Figure 6.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 6 .
Figure 6.Comparison of NH 3 -N at three river segments between field measurements and the model.

Figure 7 .
Figure 7.Comparison of discharge at three river segments between field measurements and the model.

Figure 8 .
Figure 8.Comparison of CODCr at 3 river segments between field measurement and model.

Figure 9 .
Figure 9.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 7 .
Figure 7.Comparison of discharge at three river segments between field measurements and the model.

Figure 7 .
Figure 7.Comparison of discharge at three river segments between field measurements and the model.

Figure 8 .
Figure 8.Comparison of CODCr at 3 river segments between field measurement and model.

Figure 9 .
Figure 9.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 8 .
Figure 8.Comparison of COD Cr at 3 river segments between field measurement and model.

Figure 7 .
Figure 7.Comparison of discharge at three river segments between field measurements and the model.

Figure 8 .
Figure 8.Comparison of CODCr at 3 river segments between field measurement and model.

Figure 9 .
Figure 9.Comparison of NH3-N at three river segments between field measurements and the model.

Figure 9 .
Figure 9.Comparison of NH 3 -N at three river segments between field measurements and the model.

10 .
Location of water environmentally sensitive areas in Wujiang.

Figure 11 .
Figure 11.Comparison of average water quality in sensitive areas under different conditions.(a) CODCr; (b) NH3-N.

Figure 11 .
Figure 11.Comparison of average water quality in sensitive areas under different conditions.(a) COD Cr ; (b) NH 3 -N.

Table 1 .
Surface water quality grade criteria.
Pollutants to the river networks of the Taihu Basin are mainly from the discharge of industrial, residential, livestock, and farming sewage.According to the statistics of emission investigation from thousands of sites from the Jiangsu Provincial Environmental Monitoring Center in 2013, the total

Table 1 .
Surface water quality grade criteria.

Table 2 .
Averaged water quality data at 15 cross-section.

Table 3 .
Details of water environmentally sensitive areas in Wujiang.

Table 3 .
Details of water environmentally sensitive areas in Wujiang.
Figure 10.Location of water environmentally sensitive areas in Wujiang.

Table 4 .
Weighting coefficient of local pollution for WESAs.

Table 4 .
Weighting coefficient of local pollution for WESAs.