Next Article in Journal
Transfer Learning with Convolutional Neural Networks for Rainfall Detection in Single Images
Previous Article in Journal
Optimization Methodology for Estimating Pump Curves Using SCADA Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Environmental Capacity Calculation Based on Control of Contamination Zone for Water Environment Functional Zones in Jiangsu Section of Yangtze River, China

1
Institute of Environment Planning and Evaluation, Hohai University, Nanjing 210098, China
2
Nanjing Water Facilities Management Center, Nanjing 210098, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(5), 587; https://doi.org/10.3390/w13050587
Submission received: 7 January 2021 / Revised: 11 February 2021 / Accepted: 18 February 2021 / Published: 24 February 2021
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
In recent years, due to unsustainable production methods and the demands of daily life, the water quality of the Yangtze River has deteriorated. In response to Yangtze River protection policy, and to protect and restore the ecological environment of the river, a two-dimensional model of the Jiangsu section was established to study the water environmental capacity (WEC) of 90 water environment functional zones. The WEC of the river in each city was calculated based on the results of the water environment functional zones. The results indicated that the total WECs of the study area for chemical oxygen demand (COD), ammonia nitrogen (NH3-N), and total phosphorus (TP) were 251,198 t/year, 24,751 t/year, and 3251 t/year, respectively. Among the eight cities studied, Nanjing accounted for the largest proportion (25%) of pollutants discharged into the Yangtze River; Suzhou (11%) and Zhenjiang (12%) followed, and Wuxi contributed the least (0.4%). The results may help the government to control the discharge of pollutants by enterprises and sewage treatment plants, which would improve the water environment and effectively maintain the water ecological function. This research on the WEC of the Yangtze River may serve as a basis for pollution control and water quality management, and exemplifies WEC calculations of the world’s largest rivers.

1. Introduction

Jiangsu Province, a Chinese province currently undergoing rapid economic development, is faced with serious problems with regard to its water environment [1,2,3]. In particular, drinking water sources, nature reserves, and control sections on the Yangtze River are at risk due to pollution. In the 21st century, water security issues have become of utmost importance [4], and mathematical models have been used to assist in water quality management [5].
In order to further protect the water ecological environment, it is necessary to study the water environmental capacity (WEC), which refers to the maximum amount of pollutants that can be contained in a water body while maintaining normal function [6,7], based on hydraulic modeling that can calculate the flow direction, flow volume, and water quality transport [8]. Widely used models include the Streeter-Phelos (S-P) model [9,10,11,12], the MIKE model [13], the Environmental Fluid Dynamics Code (EFDC) model [14], the Water Quality Analysis Simulation Program (WASP) model [15,16,17], the Quality Simulation Along River Systems (QUASAR) model [18], the Soil and Water Assessment Tool (SWAT) model, the Modular Three-dimensional Finite-difference Ground-water Flow Model (MODFLOW) [19], and the Stream Water Quality Model (QUAL2K) [20,21].
Some scholars have established mathematical models to calculate the WEC of the study areas. For instance, the WEC of the urban lakeside area of Lake Taihu was estimated by the total standard method, and the water quality was calculated with a zero-dimensional model [22]. Considering the synergetic influence of point and surface sources, a nonlinear optimization allocation model was used to calculate the WEC with the genetic algorithm based on controlled section water quality standards [23]. The subsection summation model was used to simulate the WEC of the river based on the demands and changes in the hydrological conditions [24]. An adequate characterization of the hydraulic state of a river is fundamental to the success of any water quality model, since it strongly affects various kinetic processes [25]. However, the above three WEC calculation methods do not fully consider the influence of topography on hydrodynamic simulations. As a professional modeling software, MIKE takes into account both topography and hydrological conditions to improve the hydrodynamic simulation accuracy and further improve the water quality simulation accuracy. The calculation of WEC is beneficial to resolve contradictions among administrative regions [26]. A MIKE 11 model was established to study the WEC of Qinhuai River, the main tributary of the Yangtze River [27]. Also, a WEC calculation framework based on MIKE 11 for the regional water environment functional zones was built in a plain river network area [28].
Some scholars have already studied the WEC of the Yangtze River. The WEC of the Chongqing section of the river was calculated using the coupled 1-D and 2-D models [29]. In addition, Chen et al. [30] calculated the allowable discharge loads in the Jiujiang section of the river under different hydrological conditions and in different water quality restricted zones. Li et al. [31] measured the WEC of the Nantong section of the river. However, there has been little research on the recent WEC of the whole Jiangsu section of the Yangtze River.
The Jiangsu section of Yangtze River is located at an estuary and has a wide river surface, so a 0-D or 1-D model cannot be easily applied. In contrast, the MIKE 21 model can generate a 2-D simulation and is able to simulate the water quantity and quality of sewage outlet-type problems. Jia et al. [32] simulated the increase in the concentration of pollutants using the MIKE 21 model in the Caofeidian sea area, which is significantly affected by the surrounding terrain. Therefore, a MIKE 21 model was built to calculate the pollution zone under different discharge amounts in each water environmental function zone in the Jiangsu section of the Yangtze River, considering the pollution sources, topography, hydrological conditions, and water quality objectives. When the pollution zone did not affect the drinking water sources, nature reserves, and control sections, the maximum discharge quantity was taken as the WEC. The Yangtze River is the third-longest river in the world and the longest river in Asia. Research on the WEC of the Yangtze River can serve as a basis for pollution control and water quality management of sewage outlet problems, and may play an exemplary role in WEC calculation of the world’s largest rivers.

2. Materials and Methods

2.1. Study Area

The study area is the section of the Yangtze River, the longest river in China, that is located in central Jiangsu Province (Figure 1). The area extends from the junction of Anhui Province and Jiangsu Province to the estuary of the river. The area is south of Yangzhou, Taizhou, and Nantong, north of Zhenjiang, Changzhou, Wuxi, and Suzhou, and crosses Nanjing. Jiangsu Province has a subtropical monsoon climate, which is characterized by significant temperature fluctuations and distinct seasons. Due to the influence of the monsoon climate, the area has abundant precipitation, i.e., 724–1210 mm annually. However, there are obvious regional differences in precipitation, with more in the east than in the west, and more in the south than in the north.
There are 90 water environment functional zones in the Jiangsu section of the Yangtze River (as per the functional zoning of the surface water environment in Jiangsu Province, People’s Government of Jiangsu Province, 2003), including 72 class II functional zones, 17 class III functional zones, and 1 functional zone that maintains current water quality. The chemical oxygen demand (COD), ammonia nitrogen (NH3-N), and total phosphorus (TP) concentrations of class II zones are less than 15 mg/L, 0.5 mg/L and 0.1 mg/L, respectively, and those of class III zones are less than 20 mg/L, 1.0 mg/L and 0.5 mg/L, respectively. In addition, there are 40 sensitive targets (Figure 1).

2.2. Data Collection

Water quality data were collected from 26 sampling sites along the banks of the Yangtze River in 2015. The water samples were quickly transferred in a refrigerated state to the laboratory for analysis. All water quality indices were measured according to the China National Standard approach or the Environment Standard approach. Herein, the water quality indices consisted of COD, NH3-N, and TP.
Data on the water level in the study area were obtained from the hydrological statistics of the Yangtze River Basin in the Annual Hydrological Report, People’s Republic of China, in 2015 [33]. The underwater elevation of the model was extracted from the actual topographic computer-assisted design (CAD) graphic documents of the Yangtze River from Nanjing to Gaoqiao. Information on 199 industrial enterprises (IE) and 64 sewage treatment plants (STP) that discharge sewage directly into the Yangtze River was obtained from environmental statistics from 2017 for Jiangsu Province. Data on the water quality and quantity of the 111 main tributaries into the Yangtze River were provided by the Jiangsu Environmental Protection Bureau.

2.3. Generalization of Sewage Outlet

A statistical investigation revealed that the main sources of pollutants were industrial enterprises, sewage treatment plants, and the main tributaries along the Yangtze River. A total of 134 sewage outlets were generalized based on the environmental statistics of Jiangsu Province, including the geographical location and emission amount in 2017 (Figure 1). The equation for the discharge amount of each generalized outlet is as follows:
E m i = e E m i e + f E m i f + g ( C m i g × Q m i g )
where E m i is the total emission amount, E m i e is the industrial enterprise emission amount, and E m i f is the sewage treatment plant emission amount of the ith generalized sewage outlet in the mth water functional zone (t/year). C m i g is the water quality concentration in the tributary of the ith generalized sewage outlet in the mth water functional zone (mg/L). Q m i g is the tributary flow of the ith generalized sewage outlet in the mth water functional zone (m3/s).

2.4. Water Environmental Capacity Calculation Methods

2.4.1. Water Environmental Mathematical Model

A water environmental mathematical model of the Yangtze River was developed based on the MIKE 21 program by DHI. The governing equations of the 2-D hydrodynamics in a Cartesian coordinate system are continuity equations and momentum equations of integral Navier–Stokes mean equations of incompressible fluid in three dimensions along the direction of water depth. The continuity equation was expressed as:
h t + h u ¯ x + h v ¯ y = h Q
The momentum equations were expressed as:
h u ¯ t + h u ¯ 2 x + h ν u ¯ y = f ν ¯ h g h η t h ρ p a x g h 2 2 ρ ρ x + τ s x ρ τ b x ρ 1 ρ ( s x x x + s x y y ) + x ( h T x x ) + y ( h T x y ) + h u s Q
h v ¯ t + h u v ¯ x + h v ¯ 2 y = f u ¯ h g h η y h ρ P a y g h 2 2 ρ ρ y + τ s y ρ τ b y ρ 1 ρ ( s y x x + s y y y ) + x ( h T x y ) + y ( h T y y ) + h ν s Q
where t is time (s), x , y are Cartesian coordinates (m), h is the total depth of the water (m), η is water density (kg/m3), u and v are the velocity components in the x and y directions, respectively (m/s). f = 2 Ω sin φ is Coriolis factor. τ s x , τ s y are surface wind stress in the x and y directions, respectively (kg·m/s2). τ b x , τ b y are base drag in the x and y directions, respectively (kg·m/s2). S x x , S x y , S y y are radiation stress tensors. P a is atmospheric pressure (pa). Q is point source emission flow (m3/s). g is the acceleration of gravity (m/s2). ρ is the density of water (kg/m3). u s , v s are the rates of external discharge into an environmental body of water (m/s).
Transverse stress T i j includes viscous resistance, turbulent friction resistance, and differential advection friction resistance, which can be calculated by the vorticity equation of the average vertical velocity:
T x x = 2 A u ¯ x , T x y = A ( u ¯ y + v ¯ x ) , T y y = 2 A v ¯ x
where A is the area of contact (m2).
The basic equation of convective diffusion of pollutants in two-dimensional, nonuniform flow can be expressed as:
h C ¯ t + h u ¯ C ¯ x + h v ¯ C ¯ y = h [ x ( D x x ) + y ( D y y ) ] C ¯ + S + S k
where C ¯ is the average pollutant concentration (mg/L), u ¯ and v ¯ are the velocity components in the x and y directions, respectively (m/s), and D x and D y are the diffusion coefficients in the x and y directions, respectively (m2/s). S is the source of pollution (g/m2/s). S k is the dynamic conversion (g/m2/s).

2.4.2. Designed Hydrological Conditions

According to the principle of selecting the most unfavorable hydrological conditions for the water environment, this calculation selected the water flow at Datong Hydrographic Station of the Yangtze River with a 90% hydrological guarantee rate after the construction of the Three Gorges Dam, comprehensively considering the influence of the Middle Route, the East Route of the South-to-North Water Diversion Project, and the water diversion project from Yangtze River to Huaihe River.

2.4.3. Model Parameter Sensitivity Analysis

A sensitivity analysis of the Manning number (N), COD degradation coefficient (KC), ammonia nitrogen degradation coefficient (KN), and total phosphorus degradation coefficient (KP) in the model was carried out using the Morris screening method. The model was perturbed several times with a 5% step size, and the average value of the Morris coefficient was calculated using the following equation [34,35,36,37]:
S = i = 0 n 1 ( Y i + 1 Y i ) / Y 0 P i + 1 P i n 1 ,
where S is the sensitivity discriminant factor, Yi+1 is the output value of the i+1-th run of the model, Yi is the output value of the i-th run of the model, is the initial value of the calculated result after parameter calibration, Pi+1 is the percentage of the parameter value after the i+1-th model operation relative to the calibration parameter, Pi is the percentage of the parameter value after the i-th model operation relative to the calibration parameter, and n is the total number of times the model ran.

2.4.4. Percentage Bias and the Nash–Sutcliffe Efficiency Coefficient

The calibration and validation of the model were carried out using the technique of graphic evaluation and error evaluation.
As a graphic evaluation technique, the time series diagram compared the measured value with the simulated value to judge the accuracy of the latter.
The model error was evaluated by calculating the percentage bias (PBIAS) and Nash–Sutcliffe efficiency coefficient (NSE) between the calculated value and the measured value. The PBIAS and NSE were calculated as follows [9,38,39,40]:
P B I A S = i = 1 n ( Y i s Y i o ) i = 1 n Y i o × 100 % ,
N S E = 1 i = 1 n ( Y i o Y i s ) 2 i = 1 n ( Y i o Y i o ¯ ) 2 ,
where P B I A S is the percentage bias (%), N S E is the Nash-Sutcliffe efficiency coefficient, Y i s is the simulated value of the model with time i (the units of water flow is m3/s, water quality is mg/L), Y i o is the observed value of the model with time i (the units of water flow is m3/s, water quality is mg/L), and Y i o ¯ is average the observed value of the model with time i (the units of water flow is m3/s).

2.4.5. Water Environmental Capacity Calculation

Taking the emission amount of each generalized outlet as the background value, the contamination zone was calculated by importing the pollution source at the generalized outlet in the established mathematical model of the water environment. Then, the WEC of each water function was calculated as follows in Figure 2.
Where P m i is the pollution source of the ith generalized sewage outlet in the mth water functional zone (t/year), C m s , C m o are the water quality target and current water quality, respectively, in the mth water functional zone (mg/L), L ( P m i ) is the maximum longitudinal length of the contamination zone in the P m i sewage outlet (m), L 1 is 10% of the length of the water functional zone (m), L 2 is the distance between the P m i sewage outlet and the nearest sensitive target, which is a drinking water source, nature reserve, or control section (m), and W E C m is the WEC in the mth water functional zone (t/year).

3. Results

3.1. Generalization of Direct Discharge Pollution Sources

In 2017, 33,547 tons of COD, 1627 tons of NH3-N, and 178 tons of TP were discharged directly into the Yangtze River water environment functional zones in eight cities along the river by industrial enterprises and sewage treatment plants (Figure 1). Nanjing directly discharged the most pollution into the Yangtze River; the amounts of COD, NH3-N, and TP were 8924 t/year, 464 t/year, and 56 t/year, accounting for 26.6%, 28.5%, and 31.5% of the total, respectively (Table 1). Nantong discharged the second-highest amount of pollution, while Wuxi discharged no pollutants, directly.

3.2. Model Parameter Sensitivity Analysis Results

According to a model parameter sensitivity analysis, N was an insensitive parameter, while KC, KN, and KP were high sensitivity parameters [34,35,37] (Table 2).

3.3. Model Calibration and Validation Results

3.3.1. Model Calibration Results

The measured water level data of Nanjing Station (NJ) and Zhenjiang Station (ZJ) from 15 October to 22 October 2015, and the measured water quality data of the drinking water source of Zhangjiagang (ZJG) and Honggang (HG) from 17 October to 19 October 2015, were used to calibrate the model.
According to the model calibration results, N ranged from 29.3 to 50.6, and the degradation coefficients of COD, NH3-N, and TP were 0.2 day−1, 0.15 day−1, and 0.06 day−1, respectively.
The PBIAS and NSE of the simulated water level were calculated according to the measured water level data of NJ and ZJ. The PBIAS and NSE of NJ were 8.3% and 0.910 (Figure 3a). The PBIAS and NSE of ZJ were 12.4% and 0.905 (Figure 3b).
The PBIAS of the simulated water quality was calculated according to the measured water quality data of the drinking water source, e.g., ZJG and HG. The average PBIAS of ZJG and HG are shown in Table 3. Using ZJG as an example, the time series diagrams are shown in Figure 4.

3.3.2. Model Validation Results

The water level data of Jiangyin Station (JY) and Yingchuangang Station (YCG) from 15 October 2015 to 22 October 2015, and the measured water quality data of the drinking water source of Langgang (LG) and Liuhe (LH) from 17 October 2015 to 19 October 2015, were used to validate the model.
The PBIAS and NSE of the simulated water level were calculated according to the measured water level data of JY and YCG. The PBIAS and NSE of JY were 7.8% and 0.904 (Figure 3c). The PBIAS and NSE of YCG were 13.8% and 0.894 (Figure 3d).
The PBIAS of the simulated water quality was calculated according to the measured water quality data of the drinking water source, e.g., LG and LH. The average PBIAS of LG and LH are shown in Table 3. Using LH as an example, the time series diagrams are shown in Figure 4.

3.4. Calculation Results of Water Environmental Capacity

In accordance with the WEC calculation method, the capacity of 90 water environment functional zones was calculated using the established mathematical model (Figure 5). There were 251,198 tons, 24,751 tons, and 3251 tons of total capacities for COD, NH3-N, and TP per year, respectively. Nanjing ranked first, with total capacities for COD, NH3-N, and TP of 59,537 t/year, 8099 t/year, and 1008 t/year, accounting for 23.7%, 32.7%, and 31.0% of the total, respectively. Wuxi had the smallest capacities, with total capacities for COD, NH3-N, and TP of 831 t/year, 121 t/year, and 15 t/year, accounting for 0.3%, 0.5%, and 0.5% of the total, respectively (Table 4).

4. Discussion

4.1. Sewage Outlet Generalization

The environmental statistics of Jiangsu Province showed that Nanjing had 32 industrial enterprises and 9 sewage treatment plants along the Yangtze River. In contrast, Wuxi had no industrial enterprises or sewage treatment plants along the river. Moreover, the industrial enterprises and sewage treatment plants in Nanjing produced a large amount of sewage discharge, which explained why Nanjing directly discharged the most pollutants while Wuxi discharged no pollutants (Table 1).

4.2. Model Performance Rating

Model performance can be judged on the grounds of general performance ratings [38,39] (Table 5). Based on the PBIAS of the water level in NJ and JY, the model performance can be rated “very good” and “good” for ZJ and YCG, respectively. Based on the NSE of the water level, the model performance can be rated “very good” for NJ, ZJ, JY, and YCG.
In order to be consistent with the hydrodynamic calibration and validation time (15 October 2015 to 22 October 2015), the water quality data of the same period were selected in this study. The PBIAS for water quality in ZJG, LG, and LH received a “very good” rating, and the model performance can be rated “good” on the basis of the PBIAS in HG. Therefore, the established model was shown to be accurate and reliable, and can be used as the basis for WEC calculations.

4.3. Water Environmental Capacity

Based on the calculation of direct discharge pollutants earlier in the paper, Nanjing was found to be the city with the largest pollution discharge (Table 1). Nanjing was also the city with the largest WEC among the eight cities (Table 4), despite the fact that it had stricter water quality standards within the Yangtze River water environment functional area. According to the WEC calculation method, WEC was related to the water quality target, the maximum longitudinal length of the functional zone (L), sensitive targets, and the current sewage discharge. In order to prove that the calculation of WEC was reasonable, the length of the water functional areas was superimposed with the WEC and the discharged contaminants of each city. The results showed that cities with longer water functional zones generally had a larger WEC (Figure 6).
Furthermore, under the condition of calculating the main pollution sources, the WEC calculation result was consistent with that of the Yangtze River Water Conservancy Commission, with a gap of less than 5%. The results of this calculation were scientific and reliable, and can provide technical support for pollution load reduction and water environment treatment and improvement in each water environment functional zone.

5. Conclusions

1.
In order to respond to the national protection policy for the Yangtze River, a two-dimensional mathematical model of the water environment in the Jiangsu section was established to study the WEC of the Yangtze River. Based on the PBIAS of water level and quality, the model performance can be rated “good” or “very good,” and the model can be well applied to the calculation of WEC.
2.
The WEC was calculated on the condition that 10% of the length of the water environment functional zone was used as the sewage mixing zone, and the sensitive targets were not affected. There were 251,198 tons, 24,751 tons, and 3251 tons of total WEC for COD, NH3-N, and TP per year, respectively. The calculation of the WEC in the 90 water environment functional zones of the river can provide a scientific basis for the protection and management of water resources, and lay a foundation for water resources planning.
3.
The main tributaries, industrial enterprises, and sewage treatment plants that discharge sewage directly into the Yangtze River were considered in the process of sewage outlet generalization. Among them, the emissions of the main tributaries into the river changed with the water quality and quantity, and were not easily controlled. In contrast, the emissions from industrial enterprises and sewage treatment plants were less than the allowable emissions determined dialectically by national government departments, and were generally fixed. Controlling the discharge of pollutants in the main tributaries will be the most important and difficult task in protecting the Yangtze River. Linking the allowable discharge of the main tributaries to the WEC of the Yangtze River has become an urgent challenge.

Author Contributions

Conceptualization, Q.M. and Y.P.; methodology, Q.M. and Y.P.; software, Q.M. and R.M.; resources, Y.P.; writing—original draft preparation, Q.M.; writing—review and editing, Q.M. and R.M.; visualization, Q.M.; project administration, Y.P.; funding acquisition, Y.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Science and Technology Major Project of China, grant number 2018ZX07208-005, and National Natural Science Foundation of China, grant number 51879070.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data was obtained from Jiangsu Environmental Protection Bureau and are available from the authors with the permission of Jiangsu Environmental Protection Bureau.

Acknowledgments

We are very grateful to the anonymous reviewers and editors for their constructive comments.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Howarth, R.B.; Kennedy, K. Economic growth, inequality, and well-being. Ecol. Econ. 2016, 121, 231–236. [Google Scholar] [CrossRef]
  2. Zuo, Y.; Shi, Y.-L.; Zhang, Y.-Z. Research on the Sustainable Development of an Economic-Energy-Environment (3E) System Based on System Dynamics (SD): A Case Study of the Beijing-Tianjin-Hebei Region in China. Sustainability 2017, 9, 1727. [Google Scholar] [CrossRef] [Green Version]
  3. Xu, M.; Hu, W.-Q. A research on coordination between economy, society and environment in China: A case study of Jiangsu. J. Clean. Prod. 2020, 258, 120641. [Google Scholar] [CrossRef]
  4. McBean, E.A.; Huang, J. Water Security Implications in the 21st Century for Coastal Cities: The Imperative Need for Action. J. Water Resour. Plan. Manag. 2020, 146, 02520003. [Google Scholar] [CrossRef]
  5. De Pinto, J.V.; Freedman, P.L.; Dilks, D.M.; Larson, W.M. Models Quantify the Total Maximum Daily Load Process. J. Environ. Eng. 2004, 130, 703–713. [Google Scholar] [CrossRef]
  6. Wagner, R.C.; Dillaha, T.A.; Yagow, G. An assessment of the reference watershed approach for TMDLs with biological impairments. Water Air Soil Pollut. 2007, 181, 341–354. [Google Scholar] [CrossRef]
  7. Hu, K.; Wang, Y.; Feng, B.; Wu, D.; Tong, Y.; Zhang, X. Calculation of water environmental capacity of large shallow lakes—A case study of Taihu Lake. Hydrol. Res. 2020, 22, 223–236. [Google Scholar] [CrossRef] [Green Version]
  8. Berglund, E.Z.; Pesantez, J.; Rasekh, A.; Shafiee, M.E.; Sela, L.; Haxton, T. Review of Modeling Methodologies for Managing Water Distribution Security. J. Water Resour. Plan. Manag. 2020, 146, 03120001. [Google Scholar] [CrossRef]
  9. Gotovtsev, A.V. Determination of biochemical oxygen demand and the oxidation rate based on a modified Streeter-Phelps system. Dokl. Earth Sci. 2015, 460, 195–197. [Google Scholar] [CrossRef]
  10. Hiramatsu, K.; Harada, M.; Cuu, N.T.; Lan, T.T. Estimation of Water Environment Capacity in the Cau River Basin, Vietnam us-ing the Streeter-Phelps Model. J. Fac. Agric. Kyushu Univ. 2017, 62, 163–169. [Google Scholar] [CrossRef]
  11. Arifin, A.; Mohamed, R.M.S.R.; Al-Gheethi, A.; Kassim, A.H.M.; Yaakob, M.A. Assessment of household greywater discharge from village houses using Streeter-Phelps model in stream. Desalination Water Treat. 2020, 179, 8–18. [Google Scholar] [CrossRef]
  12. Long, B.T. Inverse algorithm for Streeter–Phelps equation in water pollution control problem. Math. Comput. Simul. 2020, 171, 119–126. [Google Scholar] [CrossRef]
  13. Nigussie, T.A.; Altunkaynak, A. Modeling the effect of urbanization on flood risk in Ayamama Watershed, Istanbul, Turkey, using the MIKE 21 FM model. Nat. Hazards 2019, 99, 1031–1047. [Google Scholar] [CrossRef]
  14. Shi, Y.; Li, Y.; Cheng, Y.; Tang, C.; Craig, P.M.; Chung, D.H. A novel sigma coordinate system to simulate abrupt changes of underwater terrain in a hydrodynamic model: Application to Lake Mead, USA. Environ. Sci. Pollut. Res. 2020, 27, 9261–9273. [Google Scholar] [CrossRef]
  15. Gordillo, G.; Morales-Hernández, M.; García-Navarro, P. Finite volume model for the simulation of 1D unsteady river flow and water quality based on the WASP. J. Hydroinform. 2019, 22, 327–345. [Google Scholar] [CrossRef]
  16. Ramirez, L.A.R.; Luna, M.A.G.; Chiriví-Salomón, J.S.; Nieto, D.M.M. Simulation of Cr-III dispersion in the High Bogotá River Basin using the WASP model. Rev. Fac. Ing. Univ. Antioq. 2019, 97, 30–40. [Google Scholar] [CrossRef]
  17. Wool, T.; Ambrose, J.R.B.; Martin, J.L.; Comer, A. WASP 8: The Next Generation in the 50-Year Evolution of USEPA’s Water Quality Model. Water 2020, 12, 1398. [Google Scholar] [CrossRef]
  18. Wang, Q.; Li, S.; Jia, P.; Qi, C.; Ding, F. A Review of Surface Water Quality Models. Sci. World J. 2013, 2013, 1–7. [Google Scholar] [CrossRef] [Green Version]
  19. Sith, R.; Watanabe, A.; Nakamura, T.; Yamamoto, T.; Nadaoka, K. Assessment of water quality and evaluation of best management practices in a small agricultural watershed adjacent to Coral Reef area in Japan. Agric. Water Manag. 2019, 213, 659–673. [Google Scholar] [CrossRef]
  20. Zhang, R.; Qian, X.; Li, H.; Yuan, X.; Ye, R. Selection of optimal river water quality improvement programs using QUAL2K: A case study of Taihu Lake Basin, China. Sci. Total Environ. 2012, 431, 278–285. [Google Scholar] [CrossRef]
  21. Gurjar, S.K.; Tare, V. Spatial-temporal assessment of water quality and assimilative capacity of river Ramganga, a tributary of Ganga using multivariate analysis and QUEL2K. J. Clean. Prod. 2019, 222, 550–564. [Google Scholar] [CrossRef]
  22. Yan, R.; Gao, Y.; Li, L.; Gao, J. Estimation of water environmental capacity and pollution load reduction for urban lakeside of Lake Taihu, eastern China. Ecol. Eng. 2019, 139, 105587. [Google Scholar] [CrossRef]
  23. Chen, L.; Han, L.; Ling, H.; Wu, J.; Tan, J.; Chen, B.; Zhang, F.; Liu, Z.; Fan, Y.; Zhou, M.; et al. Allocating Water Environmental Capacity to Meet Water Quality Control by Considering Both Point and Non-Point Source Pollution Using a Mathematical Model: Tidal River Network Case Study. Water 2019, 11, 900. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, X.; Luo, J.; Zhang, X.; Xie, J. Dynamic simulation of river water environmental capacity based on the subsection sum-mation model. Water Environ. Res. 2020, 92, 278–290. [Google Scholar] [CrossRef]
  25. Keupers, I.; Willems, P. Development and testing of a fast conceptual river water quality model. Water Res. 2017, 113, 62–71. [Google Scholar] [CrossRef]
  26. Avisse, N.; Tilmant, A.; Rosenberg, D.; Talozi, S. Quantitative Assessment of Contested Water Uses and Management in the Conflict-Torn Yarmouk River Basin. J. Water Resour. Plan. Manag. 2020, 146, 05020010. [Google Scholar] [CrossRef]
  27. Song, W.; Pang, Y.; Li, Y.; Zhang, P.; Xu, Q.; Fu, X. Calculation of water environmental capacity and pollutant sharing rate with water diversion: A case study of Qinhuai River. Water Supply 2018, 19, 1026–1035. [Google Scholar] [CrossRef]
  28. Huang, J.; Zhang, X.; Tong, Y.; Pang, Y. Water Environmental Capacity Calculation and Allocation of the Taihu Lake Basin in Jiangsu Province Based on Control Unit. Int. J. Environ. Res. Public Health 2019, 16, 3774. [Google Scholar] [CrossRef] [Green Version]
  29. Zhong, C.H.; Xing, Z.G.; Jiang, L.W.; Wang, P. Study on water environmental capacity of Chongqing Section of Yangtze River and Jialing River. Chongqing Environ. Sci. 1994, 3, 43–50. (In Chinese) [Google Scholar]
  30. Chen, Y.; Yanwu, L.; Mou, H.; Yu, C.; Peng, M. Study on water environmental capacity of Jiujiang Section of Yangtze River. Res. Environ. Sci. 1994, 1, 24–29. [Google Scholar] [CrossRef]
  31. Li, X.; Chen, X.M.; Cui, P. Research on the confirmation of water environmental capacity and total control of the Yagntze River in Nantong city. Environ. Monit. Forewarn 2012, 4, 10–14. [Google Scholar] [CrossRef]
  32. Jia, P.; Wang, Q.; Lu, X.; Zhang, B.; Li, C.; Li, S.; Li, S.; Wang, Y. Simulation of the effect of an oil refining project on the water environment using the MIKE 21 model. Phys. Chem. Earth Parts A/B/C 2018, 103, 91–100. [Google Scholar] [CrossRef]
  33. Hydrology Bureau of Ministry of Water Resources of China. The lower reaches of the Yangtze River. In Annual Hydrological Report, People’s Republic of China, 1st ed.; Hydrology Bureau of Ministry of Water Resources of China: Beijing, China, 2015; Volume 6, pp. 1–75. [Google Scholar]
  34. Lenhart, T.; Eckhardt, K.; Fohrer, N.; Frede, H.-G. Comparison of two different approaches of sensitivity analysis. Phys. Chem. Earth Parts A/B/C 2002, 27, 645–654. [Google Scholar] [CrossRef]
  35. Feng, K.; Lu, Z.; Yang, C. Enhanced Morris method for global sensitivity analysis: Good proxy of Sobol’ index. Struct. Multidiscip. Optim. 2018, 59, 373–387. [Google Scholar] [CrossRef]
  36. Ho, L.V.; Khatir, S.; Roeck, G.D.; Bui-Tien, T.; Wahab, M.A. Finite element model updating of a cable-stayed bridge using metaheuristic algorithms combined with Morris method for sensitivity analysis. Smart Struct. Syst. 2020, 26, 451–468. [Google Scholar] [CrossRef]
  37. Xu, H.; Badea, A.F.; Cheng, X. Sensitivity analysis of thermal-hydraulic models based on FFTBM-MSM two-layer method for PKL IBLOCA experiment. Ann. Nucl. Energy 2020, 147, 107732. [Google Scholar] [CrossRef]
  38. Moriasi, D.N.; Arnold, J.G.; van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  39. Hydrologic and Water Quality Models: Performance Measures and Evaluation Criteria. Trans. ASABE 2015, 58, 1763–1785. [CrossRef] [Green Version]
  40. Almagro, A.; Oliveira, P.T.S.; Rosolem, R.; Hagemann, S.; Nobre, C.A. Performance evaluation of Eta/HadGEM2-ES and Eta/MIROC5 precipitation simulations over Brazil. Atmos. Res. 2020, 244, 105053. [Google Scholar] [CrossRef]
Figure 1. Study area and pollution source generalization map: industrial enterprises (IE); sewage treatment plants (STP).
Figure 1. Study area and pollution source generalization map: industrial enterprises (IE); sewage treatment plants (STP).
Water 13 00587 g001
Figure 2. Calculation flow chart of the water environmental capacity (WEC).
Figure 2. Calculation flow chart of the water environmental capacity (WEC).
Water 13 00587 g002
Figure 3. Calibration and validation for water level: (a) Nanjing (NJ); (b) Zhenjiang (ZJ); (c) Jiangyin (JY); (d) Yingchuangang (YCG).
Figure 3. Calibration and validation for water level: (a) Nanjing (NJ); (b) Zhenjiang (ZJ); (c) Jiangyin (JY); (d) Yingchuangang (YCG).
Water 13 00587 g003
Figure 4. Calibration and validation for water quality: (a) Zhangjiagang (ZJG)-chemical oxygen demand (COD); (b) ZJG-ammonia nitrogen (NH3-N); (c) ZJG-total phosphorus (TP); (d) Liuhe (LH)-COD; (e) LH-NH3-N; (f) LH-TP.
Figure 4. Calibration and validation for water quality: (a) Zhangjiagang (ZJG)-chemical oxygen demand (COD); (b) ZJG-ammonia nitrogen (NH3-N); (c) ZJG-total phosphorus (TP); (d) Liuhe (LH)-COD; (e) LH-NH3-N; (f) LH-TP.
Water 13 00587 g004aWater 13 00587 g004b
Figure 5. Water environmental capacity (WEC) of 90 water environmental function zones: (a) Chemical oxygen demand (COD); (b) Ammonia nitrogen (NH3-N); (c) Total phosphorus (TP).
Figure 5. Water environmental capacity (WEC) of 90 water environmental function zones: (a) Chemical oxygen demand (COD); (b) Ammonia nitrogen (NH3-N); (c) Total phosphorus (TP).
Water 13 00587 g005aWater 13 00587 g005b
Figure 6. Relationship among the water environmental capacity (WEC) (left bars), the discharged amount of contaminants (right bars), and the maximum longitudinal length of the functional zone (L, solid line). COD Chemical oxygen demand; NH3-N Ammonia nitrogen; TP Total phosphorus.
Figure 6. Relationship among the water environmental capacity (WEC) (left bars), the discharged amount of contaminants (right bars), and the maximum longitudinal length of the functional zone (L, solid line). COD Chemical oxygen demand; NH3-N Ammonia nitrogen; TP Total phosphorus.
Water 13 00587 g006
Table 1. Emissions of pollutants from industrial enterprises and sewage treatment plants.
Table 1. Emissions of pollutants from industrial enterprises and sewage treatment plants.
CityCOD (t/year)NH3-N (t/year)TP (t/year)
Nanjing892446456
Yangzhou26664717
Zhenjiang345727615
Taizhou16105519
Changzhou60592589
Wuxi000
Suzhou364620319
Nantong718532343
SUM33,5471626178
COD: chemical oxygen demand; NH3-N: Ammonia nitrogen; TP: total phosphorus.
Table 2. Analysis of parameter sensitivity.
Table 2. Analysis of parameter sensitivity.
Parameter | S | Sensitivity LevelGrading Standard
N0.012–0.046insensitivity0 ≤ | S | < 0.05, insensitivity
0.05 ≤ | S | < 0.2, medium sensitivity
0.2 ≤ | S | < 1, sensitivity
| S | ≥ 1, high sensitivity
KC2.714–4.059high sensitivity
KN2.304–3.346high sensitivity
KP0.937–1.133high sensitivity
N: Manning number; KC: chemical oxygen demand degradation coefficient; KN: ammonia nitrogen degradation coefficient; KP: total phosphorus degradation coefficient; S: sensitivity.
Table 3. The percentage bias (PBIAS) of the simulated water quality.
Table 3. The percentage bias (PBIAS) of the simulated water quality.
Drinking Water SourceCOD (%)NH3-N (%)TP (%)
ZJG5.823.814.1
HG5.129.717.9
LG3.223.024.2
LH1.33.816.4
COD: chemical oxygen demand; NH3-N: Ammonia nitrogen; TP: total phosphorus. HG Honggang; LG Langgang; LH Liuhe; ZJG Zhangjiagang.
Table 4. Water environmental capacity (WEC) of eight cities.
Table 4. Water environmental capacity (WEC) of eight cities.
CityCOD (t/year)NH3-N (t/year)TP (t/year)
Nanjing59,53780991008
Yangzhou26,8302668356
Zhenjiang44,6833344480
Taizhou37,3512630388
Changzhou668953876
Wuxi83112115
Suzhou42,3844365546
Nantong32,8932986382
SUM251,19824,7513251
COD: chemical oxygen demand; NH3-N: Ammonia nitrogen; TP: total phosphorus.
Table 5. General performance ratings for recommended statistics.
Table 5. General performance ratings for recommended statistics.
Performance RatingPBIAS (%)NSE
Water LevelWater QualityWater Level
Very Good[−10, 10][−25, 25](0.75, 1]
Good[−15, −10) ∪ (10, 15][−40, −25) ∪ (25, 40](0.60, 0.75]
Satisfactory[−25, −15) ∪ (15, 25][−70, −40) ∪ (40, 70](0.40, 0.60]
Unsatisfactory[−100, −25) ∪ (25, +∞)[−100, −70) ∪ (70, +∞)(−∞, 0.40]
PBIAS: percentage bias; NSE: Nash–Sutcliffe efficiency coefficient.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, Q.; Pang, Y.; Mu, R. Water Environmental Capacity Calculation Based on Control of Contamination Zone for Water Environment Functional Zones in Jiangsu Section of Yangtze River, China. Water 2021, 13, 587. https://doi.org/10.3390/w13050587

AMA Style

Ma Q, Pang Y, Mu R. Water Environmental Capacity Calculation Based on Control of Contamination Zone for Water Environment Functional Zones in Jiangsu Section of Yangtze River, China. Water. 2021; 13(5):587. https://doi.org/10.3390/w13050587

Chicago/Turabian Style

Ma, Qiuxia, Yong Pang, and Ronghua Mu. 2021. "Water Environmental Capacity Calculation Based on Control of Contamination Zone for Water Environment Functional Zones in Jiangsu Section of Yangtze River, China" Water 13, no. 5: 587. https://doi.org/10.3390/w13050587

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop