Next Article in Journal
Preliminary Study of Computational Time Steps in a Physically Based Distributed Rainfall–Runoff Model
Previous Article in Journal
Irrigation Management Based on Reservoir Operation with an Improved Weed Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment of Upper-Middle Reaches of Luanhe River Basin in Sudden Water Pollution Incidents Based on Control Units of Water Function Areas

1
Nankai University Key Laboratory of Pollution Processes and Environmental Criteria of the Ministry of Education, Remediation and Pollution Control for Urban Ecological Environment, College of Environmental Science and Engineering Numerical Simulation Group for Water Environment, Tianjin 300350, China
2
School of Computing, Mathematics and Computation, Manchester Metropolitan University, Manchester M1 5GD, UK
3
Water Resources Protection Bureau of Haihe Water Conservancy Commission of Ministry of Water Resources, Tianjin 300171, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(9), 1268; https://doi.org/10.3390/w10091268
Submission received: 8 August 2018 / Revised: 8 September 2018 / Accepted: 10 September 2018 / Published: 17 September 2018
(This article belongs to the Section Urban Water Management)

Abstract

:
Upper-middle reaches of Luanhe River Basin belongs to Haihe River Basin and it is a pretty important water source to North China, especially to Tianjin and Tangshan Cities, Hebei. Based on control units of the water function areas and the sub-basins the working units were produced. The index system for environmental risk source hazard was constructed by adopting the pressure state response (PSR) environmental analysis model. The environment risk sources are identified, and their hazard grade assessment is performed. In the environmental risk source hazard an assessment index system, namely “downstream characteristics of environmental risk sources” is added by taking the fact into account that environmental risk sources themselves are affected by different functional areas (working units) downstream of the rivers. Through collecting hazard data, determining their standards and weights for environment risk sources, the fuzzy comprehensive evaluation method is used to calculate the risk source hazard grades and the vulnerability grades of the working units. Using the one-dimensional exponential decay river model for pollutants in rivers, the hazard grade evaluation method of working units is established. This consists of two parts: (1) The risk source hazard grade of the working unit itself, and (2) the impact of the risk sources upstream on the working unit downstream of the rivers. Combining the hazard grade with the vulnerability grade of the working unit, the risk grade of the working unit is evaluated through the risk matrix. The risk zones of the watershed are realized by merging working units in the same control units of the water function areas with the same risk grades. The risk zoning of sudden water pollution incidents in the upper and middle reaches of Luanhe River Basin is obtained by applying the above risk zoning method. It is found that there are 55 risk zones in total, including three highest risk zones, 15 higher risk zones, 14 lower risk zones, 23 lowest risk zones. These results indicate that the upper and middle reaches of River Luanhe are overall at low risk. The corresponding management methods for the different risk zones are suggested.

1. Introduction

In recent years, due to the fast development of economy, more frequent incidents in water pollution occur in China, for instance, the 2005 Songhua River serious water pollution incident, the Taihu lake water pollution incident in 2007, Changzhi, Shanxi aniline leakage in 2012 [1]. Although these sudden incidents occur at low probability, it can cause irreparable consequences due to its sudden strong destructive characteristics at large scale. This together with its characteristics of multiple sources and receptors makes it very difficult to research.
There are two definitions of risk [2]. The first is defined as the product of the probability of occurrence of an event and the result or consequence of the event [3,4], and it is used for risk assessment in insurance, security fields in most cases [5,6], and also used for environmental risk sources, such as enterprises, construction projects, which is called “accident risk” [7,8,9]. The second is defined as the combination of hazard and vulnerability [10,11,12]. We construct environmental risk assessment model or system for environmental risk assessment of regional or watershed scale according to the second definition [13,14,15,16,17].
There are many types of basic units in risk zoning. The grid is used as the basic units [15], and all the risk evaluation indexes need to be obtained by interpolation. Administrative districts can also be used as basic units [10]. The social economic evaluation index can be obtained directly, but due to the existence of several scattered reaches in the basic unit the regional natural factors are not fully considered in the region.
There is a little study on risk assessment and regionalization of sudden water pollution events in the basin scale. Paul, et al. (2004) used 2-D surface water flow model CE-QUAL-W2 Version 3.1 to effectively study the risks associated with accidental or intentional releases of a hazardous material into an inland waterway [18]. Liu, et al. (2012) use a fuzzy fault tree model and the comprehensive influence index model to carry out sudden water pollution risk assessment in the region, which is lack of data, and establish multi-scale evaluation model to the river unit of the whole region [19]. In this study, the water function area is examined as an index of the receptor, but the interaction of the environmental risk between the upstream and downstream units of the river is not considered. Mehmet, et al. (2013) used modified pollution indices to evaluate the pollution status in the middle section of the Lower Seyhan River Basin and the water quality index (WQI) was utilized to evaluate water pollution levels [20]. Zhang, et al. (2014) introduced relative risk model (RRM) to quantify model evaluation factors, and undertook risk assessment to the old river basin. This study treats the sub basin as an evaluation unit, taking full account of the natural features of the watershed [21]. Most importantly, these studies have identified “Zoning units”, and the risk of the units are evaluated, but it is not made in the true sense of the “risk zoning”, because the zoning unit should satisfy the maximum risk similarity in the unit, while the risk difference is the biggest between the units [22].
Luanhe River Basin belongs to Haihe River basin, China, and Luanhe River is the second largest river into the Bohai Sea of north China, and belongs to typical semi-arid and seasonal wind climate region. It is an important natural barrier and drinking water resource for Beijing, Tianjin and Tangshan Region with more than 50 million populations. In the past, most researches were focused on the temporal and spatial distributions of water and sediment yield in the Luanhe River Basin [23] and landscape ecological risk response to land use change in the basin [24]. Little attention is paid to assess the risk of Luanhe River basin in sudden water pollution incidents. So this paper aims to improve the current risk assessment method and apply this new method to study the risk of the upper and middle reaches of Luanhe River Basin in sudden water pollution incidents in view of increasingly sudden water pollution incidents.

2. Research Area

The upper and middle reaches of Luanhe River Basin is the local basin in front of the dam of Daheiting Resivior, which is mainly across two provinces (autonomous regions): The Inner Mongolia Autonomous Region and Hebei Province, 18 counties/cities/districts/banners (CD County means Chengde County, FN County means Fengning County, KC County means Kuancheng County, LH County means Longhua County, LP County means Luanping County, PQ County means Pingquan County, Sl district means Shuangluan district, SQ district means Shuangqiao District, YSYZ Mining District means Yingshouyingzi Mining District, WC County means Weichang County, XL County means Xinglong County, QX County means Qianxi County, DL County means Duolun County, GY County means Guyuan County) (Figure 1).

3. Research Method

3.1. Data Source

Digital Elevation Model (DEM) data comes from the United States Geological Survey (USGS), the resolution is 90 m × 90 m, and the Albert equal projection, and GCS_Bejing_1954 geographic coordinate system is used. The data of the basin landform, water function areas, the environmental risk sources and regional natural economic and social information are obtained through the statistical yearbook, literature search, field investigation, visiting to relevant departments and expert consultation.

3.2. Water Function Area and Control Units

Water function area is the waters delimited according to its dominant function and carried out the corresponding water environment quality standard, which is based on natural conditions and development and utilization of water resources, comprehensive river basin planning, water resources and water ecosystem protection and economic and social development requirements to meet the demand of rational development, utilization, saving, and protection of water resources.
Control unit is obtained by superimposing the spatial data of the catchment area, water system distribution and its flow direction, boundary of water function area, control section distribution, administrative boundary and other indexes on a digital map, which is based on the principle of division of pollution control unit.
The relationship between the water functional area and control unit: (1) The water function area is the premise and basis for dividing the control unit, which is combined with three other elements (see Figure 2) to obtain a control unit; and (2) The water function area includes only the water area, and the control unit comprises the water area and the corresponding land area.
According to the control unit division method for water function areas, the upper-middle reaches of Luanhe River Basin includes 34 water function areas, and the corresponding control units, as shown in Figure 3 and Table 1.

3.3. The Division of Working Units and its Advantages

Each sub basin within the water function area control unit is defined as a working unit for risk assessment transition units based on the water function area control units, which are determined by the river basin management institution. All indexes are represented on the working units, and the indexes that are not in the working unit need to be interpolated to the working units. The advantages are that: (1) The risk assessment indexes associated with rivers within the working unit are complete, and there is no artificial cutting water bodies; (2) it is generally believed that once a sudden water pollution event occurs, the pollutant will diffuse in the sub basin where the environmental risk source lies and transmit through the river to the downstream; and (3) it is consistent with the current water management unit in China and easy to manage.
The upper, middle reaches of Luanhe River Basin are obtained by entering the DEM and using the ArcGIS hydrological analysis model and they are divided into 107 working units by combining the water function area control unit with the sub basin, as shown in Figure 4.

3.4. Environmental Risk Sources Hazard Grade Assessment

The environmental risk sources of sudden water pollution incidents are divided into four grades: Major risk sources, big risk sources, general risk sources and slight risk sources by referring to the National Emergency Plan for Environmental Emergencies [25]. Environmental risk source hazard assessment is identified through environmental risk source identification, hazard index system construction, quantification, index weight determination and grade calculation.

3.4.1. Establishing Environmental Risk Sources Hazard Index System

Environmental risk source data are obtained through the establishment of questionnaires, field research and visiting relevant departments. Referring to previous work [26], the pressure state response (PSR) environment analysis model is adopted, and the environmental risk sources hazard evaluation index system is constructed from three aspects: (1) Environmental risk sources hazard (pressure), (2) hazard receptor (state), and (3) risk control ability (response). As shown in Figure 5, it includes target, topic and index layers.

3.4.2. Quantifying Environmental Risk Sources Hazard Index System

As the target, it consists of three topics, as shown in Figure 5, and they are Environmental risk sources hazard S1, Hazard Receptor S2 and Hazard Control Ability S3. Following is an explanation about the system.

A Environmental Risk Sources Hazard S1

A1. Industry Categories I1
Risk sources in different industry categories have varying degrees of hazard. For example, “chemical raw materials and chemical products manufacturing” and petroleum processing, coking and nuclear fuel processing” have a higher score of 0.8; and “the food manufacturing industry”, “water production and supply industry” and “beverage manufacturing industry” have a low score of 0.2. Corresponding to the four grades of the risk sources, this paper divides the risk sources into four grades, and determines quantitative standards as listed in Table 2.
A2. Risk Source Size I2
Generally, the bigger the size of the risk sources, the higher the risk. Because most companies’ employees, sales and total assets that can measure risk sources size can’t be known, and the scale of the enterprise has a good correlation with the risk sources of “registered capital”, the registered capital is used as a measure of risk sources index scale, and the grades are established according to [22] and [26], as shown in Table 2.
A3. Distance from the River I3
In general, the closer the risk sources to the river, the easier the pollutant entering the river and polluting the water environment, and the higher the hazard of the risk sources. Based on Arc-GIS software, a distance measurement (near) tool is used to obtain the nearest distance of each risk sources to the river. The grades are established according to [26], as shown in Table 2.
A4. Downstream Features I4
Referring to the influence of the pollutant produced by the risk sources to the downstream working units through river transportation and diffusion, we consider the two aspects: The distances between the risk source and the downstream working units, and the importance of the downstream working units. The more important the downstream working units (score based on the water function area characteristics of the working units shown in Table 3, and the shorter the distance between the risk source and the downstream working units, the greater the hazard. The following equation is used to describe I4:
I 4 j = ω 1 θ 1 + ω 2 θ 2 + ω n θ n ,
where I 4 j is the score of downstream features; j is a number; n represents the number of risk sources flowing through the working units; ω n stands for the score of the risk sources flowing through the working units, showing in Table 3; θ n is the distance between the risk sources and the working units (10 km).

B Hazard Receptor S2

B1. Sensitive Protection Objectives I5
Two aspects are considered: Whether the water function area to, which the working units belongs are the drinking water source and its population [27], as well as whether there are natural reserves, and their grades are shown in Table 3.
B2. Water Ecological Risk I6
Water ecology means the water effects on biology and the adaptation of biology to various water conditions [28]. When water bodies contain certain characteristic pollutants, they may have different hazards to the water environment. According to the Comprehensive Wastewater Discharge Standards and Sanitary Standards for Drinking Water, the characteristic pollutants are used for grading, as shown in Table 3, and the standard is listed in Table 2.
B3. Exposure Population I7
Use population size of the working units as the exposure population of risk sources. The standard is referred National Emergency Plan for Environmental Emergencies, as shown in Table 2.
B4. River Grade I8
It is generally believed that the higher the river grade that the risk sources pollutant flows into, the greater the possibility of pollutant transport and diffusion, and the greater the influences of the possibility on low-grade rivers. According to the study of [29], river grades are shown in Table 3, and the standard is in Table 2.

C Hazard Control Ability S3

The ability to control the hazard of sudden environmental incidents in a region can be evaluated from two aspects: The hazard control ability of the area and the hazard control ability of the risk sources. The regional hazard control ability is evaluated by the number of beds per 1000 people in the medical institutions where the risk sources lie and the per capita Gross Domestic Product (GDP) in the county. The risk sources ability to control the hazard is evaluated by the superposition of several factors as follows.
C1. Risk Source Ability to Control the Hazard I9
It is superimposed by the following factors, such as Pollution source category, whether to sign the letter of environmental protection undertaking et al. as shown in Table 3.
C2. Per Capita GDP of the County I10
Per capita GDP of the county to which the risk source belongs represents the level of economic development in the region to a certain extent [13]. Generally speaking, the higher the per capita GDP, when the sudden water pollution accidents occur in the region, the more dangerous control capability, the lower corresponding risk grade is, as shown in Table 2.
C3. The Number of Beds per 1000 People in the Medical Institutions
According to the “Standard of Medical Institutions (Trial)”, the more the number of beds per 1000 people in medical institutions means the better health care in the area. Once the sudden water pollution accidents occur, hazard control ability is stronger, and the lower the hazard grade [22], as shown in Table 2.

3.4.3. Determination of Weights of Environmental Risk Sources Hazard Assessment Indexes and Calculation of Results

An improved analytic hierarchy process (AHP) is adopted to determine the weight of the index through expert investigation, and the fuzzy comprehensive evaluation method is used to calculate the hazard grade of the environmental risk source [26].

3.5. Working Unit Hazard Assessment

Most previous work has considered that the working unit hazard grades are equal to the risk sources hazard grade in the unit [19,21]. The present research differs from that. Because once the upstream risk sources release a lot of pollutants in a short period of time and cause the sudden water pollution, leaked dangerous substances will flow through river systems and spread to the downstream, impacting on the downstream working units. Apparently, the downstream working unit hazard grade will be affected by upstream risk sources.
Due to the influence of upstream risk sources on downstream units, the working unit hazard grade is divided into two parts, as shown in Equation (2), one part is the greatest hazard grade of the risk sources in the unit, and the other part is the influence of the upstream risk sources on it.
F = G MAX + i = 1 m C x   ( I = 1 ,   2 ,   3     m ) ,
where F is the total hazard scores of the working unit; GMAX is the largest hazard value of each risk sources score in the working unit (because the sudden water pollution is generally caused by only one risk source although there may be several risk sources in the unit); C x is the effect of one upstream risk source on the working unit through rivers; m is the total number of upstream working units of the working unit under consideration.
The hazard grade of the risk source can be used as a generalized substance or pollutant, such as water quality index BOD or COD. Due to extremely shortage of data, and for simplification purposes, the effect of the upstream risk source on a downstream working unit through rivers is approximated to one-dimensional exponential decay river model as
C x = C 0 exp ( K x μ )
where C x is the influence of upstream risk sources on the working unit through rivers; x is longitudinal distance along the river; C 0 is the upstream risk source grade; μ is the mean velocity (m/s) of the channel under design flow; K is the comprehensive attenuation coefficient of pollutants (1/s).
As indicated above, probabilities that pollutant coincides with all risk sources in a working unit are quite low and usually there exists only one risk source of practical pollution, and hence the highest value of risk grade scores of all risk sources in a working unit is considered. Hazard grades and corresponding scores are shown in Table 4.

3.6. Working Unit Vulnerability Assessment

Vulnerability refers to the extent to which something is vulnerable or injured [30]. Corresponding to the hazard grade, the study divides the vulnerability into four grades, which are higher vulnerability, high vulnerability, moderate vulnerability and low vulnerability. Working unit vulnerability needs to be identified through vulnerability index system construction, quantification, index weight determination and grade calculation.

3.6.1. Establishment of Vulnerability Index System for Working Units

The index system of vulnerability is divided into three layers: Target, topic, and index layers, as in Figure 6.

3.6.2. Working Unit Vulnerability Index System Quantification

D Natural Factors S1

D1. Rainfall I1
The amount of rainfall is related to the amount of biomass in the ecological resources. The more the rainfall in a region, the more complex the biological networking and the stronger stability of the ecosystem, and the less vulnerable the region [26], as shown in Table 5.
D2. Sensitivity I2
Determined by the nature of the load-bearing body itself, the sensitivity of the hazard receptor is mainly reflected by sensitive targets. When the sensitivity is different in the face of the same water pollution incident, the impact is also different. Taking into account the impact of sudden water pollution incidents on residents, and other organisms in the working unit, the two scores are superimposed as a sensitivity score of the working unit as given in Table 6.
D3. Water Quality I3
The better the water quality, the weaker the ability to resist hazard and the higher the vulnerability. According to the Surface Water Environment Quality Standard, and considering that there is no river reach meeting the first grade of the national water quality standard in the study area, the index is defined as the four grades standard, as shown in Table 5.

E Social Factors S2

E1. Per Capita GDP I4
Per capita GDP is an important component of socio-economic factors, and indicates the economic level of an area. Economic factors are closely related to the vulnerability of working units. The more developed the economy, the greater the risk of loss [13], as shown in Table 5.
E2. Per Capita Cultivated Land Area I5
Water resources are vital to cultivated land. The dependence of cultivated land on water resources determines that the per capita cultivated land area is a component of the vulnerability of working units. Per capita cultivated land area is closely related to the vulnerability of regional ecological environment [31]. The more the per capita arable land area, the more vulnerability the ecological environment, and the lower the per capita cultivated land area, the lower the vulnerability is, as shown in Table 5.
E3. Population Density I6
The higher the population density index, the more vulnerable people are affected by the risk, the higher the vulnerability is [15] as shown in Table 5.

3.6.3. Working Unit Vulnerability Index Weight Determination and Result Calculation

An improved analytic hierarchy process (AHP) is adopted through expert investigation to determine the weight of the index, and the fuzzy comprehensive evaluation method is used to calculate the working unit vulnerability [26].

3.7. Risk Assessment

The risk grade of each working unit can be obtained from the risk grade matrix [32], which is based on the evaluation of working units hazard grade and vulnerability grade, as shown in Figure 7.
According to the environmental risk zoning theory [22], the adjacent working units belong to the same water function control unit and the same risk grade are merged, and the watershed risk zoning is obtained.

4. Results Analysis

4.1. Environmental Risk Sources Hazard Assessment Result

Combined with the standard in Table 2 and the fuzzy comprehensive evaluation method, the environmental risk sources risk grade results of the upper-middle reaches of the Luanhe River are obtained. Risk sources evaluation results are shown in Figure 8 and Table 7. There are 43 risk sources that have 15 major risk sources, 14 big risk sources, six general risk sources and eight slight risk sources. From the point of view of the risk level, most risky sources belong to the larger risk. From Figure 8 we can see that in the upper-middle reaches of the Luanhe River, risk sources distribution is concentrated, mainly distributed in the middle and lower segment of the study area.

4.2. Working Units Hazard Assessment Result

Combining the environmental risk sources risk assessment result in Section 4.1 with the calculation method in Section 3.3 the conclusion of hazard assessment of 107 working units in the upper reaches of Luanhe River is obtained, as shown in Figure 9. Most working units in the study area have low grade hazard, i.e., 75 are low hazard working units, eight moderate hazard working units, five high hazard working units, and 19 higher hazard working units, and hence the overall hazard is low.

4.3. Working Units Vulnerability Assessment Result

Combined with the standard in Table 5 and the fuzzy comprehensive evaluation method, the vulnerability assessment results are shown in Figure 10. This shows that there are eight higher vulnerability working units, six high vulnerability working units, 22 moderate vulnerability working units and 71 low vulnerability working units in the upper-middle reaches of Luanhe River Basin, indicating that the study area is overall in lowest vulnerability and the upper-middle reaches of Luanhe River Basin have strong ability to resist risks. High vulnerability areas are mainly located in the upper part of the study area, due to their fragile ecosystem and high sensitivity scores. It is worth noting that the working units for Daheiting reservoir is in high vulnerability, which needs to focus on investigation into protection.

4.4. Working Units Risk Assessment Result

The risk grade of working units can be obtained by evaluating the hazard and vulnerability grade of the working units with the risk matrix. Figure 11 shows that the study area is at low risk (53 units), followed by moderate risk (25 units) and high risk (25 units), least for higher risk (four units), indicating overall low risk in the upper-middle reaches of Luanhe river.

4.5. Risk Zoning Results Based on Water Function Area Control Unit

After obtaining the working units risk grade, according to the environmental risk zoning theory, merging the working units with the same risk grade in a water function area, four risk grades zones are obtained, and they are higher risk zone I, high risk zone II, moderate risk zone III and low risk zone IV, as shown in Figure 12. The risk zoning is scattered distribution in the upper-middle reaches of the Luanhe River, and there are three higher risk zones I, 15 high risk zones II, 14 moderate risk zones III, and 23 low risk zones IV, which illustrate that the risk in the upper reaches of Luanhe River is generally low.

5. Discussion

The risk zoning results for the upper-middle reach of the Luanhe River from sudden water pollution incidents have been recognized by the relevant departments and they apply the methods to practical works. For example, the major risk source “Chengde iron and Steel Group Co., Ltd.”, which is located in Daheiting reservoir and near the drinking water source are identified by the relevant departments and has been shut down.
When studying the downstream features I4 for a very large river basin area, the risk sources hazard upstream to downstream working units generally requires a coupled one-dimensional, two-dimensional or three-dimensional hydrodynamic and water quality model. In particular, in the study area the Panjiakou Reservoir belongs to deep water reservoir and a three-dimensional numerical simulation model is necessary; the Daheiting reservoir is wide with a relatively small water depth and a two-dimensional numerical simulation model is appropriate; and between these two reservoirs, there is a narrow river and a one-dimensional river network model is sufficient. At the same time, the coefficients in those equations need to be determined through actual observations and experiments.
The information about the distribution of water conservancy infrastructures in the study area fails to fully investigate and obtain. For example, a dam may have a hysteresis effect to mitigate the adverse effects of water pollution incidents, including interception, temporary storage, by impounding water when sudden water pollution incidents occur. Therefore, the quantity, quality and distribution of the dam in the area can be used as an index of environmental risk source in hazard control capability. Similarly, the role of other structures can be considered. For instance, tailings dams actually exist in the Luanhe River Basin and it should be taken as risk sources. However, detailed data of tailings dams, such as distribution and quantities in the Luanhe River Basin are not available, and thus the present analysis cannot take them as the environmental risk sources into account.
After identification and obtaining the risk zoning results of sudden water pollution incidents in the upper and middle reaches of Luanhe River Basin, the prevention and control technical measures and schemes of the sudden water pollution events in the upstream of Luanhe River could be proposed based on the national and local policies and regulations, the basic principles of water pollution control, the constituent factors of the risk level of risk zones and the actual situation. They are as follows: (A) Different management measures and schemes should be proposed according to different risk and vulnerability levels in the risk area; (B) mitigate risk levels through controlling risk sources; (C) control risks from multi-aspects, such as legal, social, economic, political and technological factors; and (D) strengthen the protection of risk receptors to reduce vulnerability. For high risk area, for example, because the hazard and vulnerability are in general very high, it is necessary to ban on the construction of new risk sources (e.g., new chemical plants and factories), to prohibit the introduction of new sensitive receptors, to shut down or move major risk sources out in the area, to establish strict and comprehensive environmental pollution monitoring system, to strengthen risk prevention and supervision, as well as to increase the construction of protective barriers for high sensitive receptors.
It is noteworthy to point out that the proposed risk assessment method is based on the water function area control unit, a Chinese and scientific concept and a practical tool in effectively managing water resources and fighting water pollution according to Chinese experience. We strongly believe that this concept could be applied and extended in many other countries and regions, especially with a water shortage and severe water pollution.

6. Conclusions

The risk of water pollution incidents of upper-middle reaches of Luanhe River Basin is assessed in the process of regionalization, leading to the following improvement: (1) Each sub basin within a water function area control unit is identified as a working unit, and taken as the transition unit of the risk assessment regionalization in a sudden water pollution event; (2) in environmental risk source hazard assessment index system, the risk source flows through the downstream different water function areas (working units) and it will affect the environmental risk source hazard grade itself, and thus increase the “Downstream features” of environmental risk sources index; (3) the calculation method of hazard from risk source to working unit is established and the working unit hazard grade is composed of two parts: The risk source hazard grade of the working unit and the influence of the upstream risk source to the downstream working unit hazard; and (4) the risk zoning theory is applied to combine the working units with the same risk grade in the same water function control unit for obtaining the final risk zoning unit.
In the upper-middle reaches of Luanhe River Basin area there are 34 control units. Combining with the sub basin, it is further divided into 107 working units. There are 43 environmental risk sources, among which 15 are the major risk sources, and 14 are the big risk sources which are concentrated and distributed in the middle and lower reaches of the study area. However, there are 75 low hazard working units and the overall hazard is low. Almost working units are in low vulnerability, which means that the upper-middle reaches of Luanhe River basin has the strong ability to resist risks. The risk grade of each working unit is obtained through risk matrix. Working units are combined with those of the same risk grades and the water function area control units, resulting in 55 risk zones, including three higher risk zone I, 15 high risk zone II, 14 moderate risk zone III, and 23 low risk zone IV. This indicates that the risk of the upper-middle reaches of Luanhe River Basin is low, and the corresponding management methods for the different risk zones are suggested. The Daheiting reservoir as one of drinking water sources of Tianjin belongs to the higher risk zone I, and its management is required, which includes separating and cutting off its risk sources, and at the same time increasing the protective barriers of its high-sensitivity receptors.

Author Contributions

Conceptualization, S.H. and Y.X.; Methodology, S.H. and Y.X.; Software, S.H., Y.X. and Y.L.; Validation, S.H., Y.X., F.K. and M.L.; Formal Analysis, S.H. and Y.X.; Investigation, S.H., Y.X. and F.K.; Resources, S.H., Y.X. and F.K.; Data Curation, S.H., Y.X., F.K., M.L. and Y.L.; Writing-Original Draft Preparation, Y.X.; Writing-Review & Editing, S.H. and J.Z.; Visualization, Y.X. and Y.L.; Supervision, S.H.; Project Administration, S.H.; Funding Acquisition, S.H.

Funding

This research was funded by [NSFC] Grant Number 11672139, [NSF of Tianjin] Grant number. 15JCYBJC22500 and the Fundamental Research Funds for the Central Universities of China.

Acknowledgments

The authors are grateful to the editors and the anonymous reviewers for their insightful comments and suggestions, which improved the paper greatly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, X.; Zhu, J.; Li, Z.; Xie, J. An Overview of the Domestic and Foreign Research on Accidental Water Pollution Events. China Rural Water Hydropower 2015, 6, 1–5. [Google Scholar]
  2. Huang, C.; Liu, A.; Wang, Y. A discussion on basic definition of disaster risk. J. Nat. Disasters 2010, 19, 8–16. [Google Scholar]
  3. Orvos, D.R.; Cairns, J. Developing a risk assessment strategy for the Chesapeake Bay. Hydrobiologia 1991, 215, 189–203. [Google Scholar] [CrossRef]
  4. Blair, A.N.; Ayyub, B.M.; Bender, W.J. Fuzzy stochastic risk-based decision analysis with the mobile offshore base as a case study. Mar. Struct. 2001, 14, 69–88. [Google Scholar] [CrossRef]
  5. Mao, H. Discussion on Some Basic Rules and Concepts of Safety Science. China Saf. Sci. J. 2009, 19, 12–16. [Google Scholar]
  6. Xie, Z.; Zhou, J. Re-studying and Understanding the Fundamental Concept of Risk. Insur. Stud. 2013, 2, 101–108. [Google Scholar]
  7. Guo, L.; Yuan, P.; Song, Y.; Wang, L.; Peng, J. Environmental Risk Classification and Management for Chemical Industry Parks. J. Environ. Eng. Technol. 2011, 1, 403–408. [Google Scholar]
  8. Wei, K.; Song, Y.; Peng, J.; Yuan, P.; Wen, L.; Wa, Y. Environmental risk source and its classification. J. Saf. Environ. 2011, 10, 85–89. [Google Scholar]
  9. Hou, D.; Ge, X.; Huang, P.; Zhang, G.; Loáiciga, H. A real-time, dynamic early-warning model based on uncertainty analysis and risk assessment for sudden water pollution accidents. Environ. Sci. Pollut. Res. 2014, 21, 8878–8892. [Google Scholar] [CrossRef] [PubMed]
  10. Qu, C.; Bi, J.; Huang, L.; Li, F.; Yang, J. Dynamic Comprehensive Evaluation on Regional Environmental Risk. Acta Sci. Nat. Univ. Pekin. 2010, 46, 477–482. [Google Scholar]
  11. Xu, Y.; Ma, M.H.; Gao, J.F. The evaluation method of ecological risk assessment based on watershed scale—Take the Taihu watershed as example. China Environ. Sci. 2012, 32, 1693–1701. [Google Scholar]
  12. Salinas-Rodriguez, C.; Gersonius, B.; Zevenbergen, C.; Serrano, D.; Ashley, R. A semi risk-based approach for managing urban drainage systems under extreme rainfall. Water 2018, 10, 384. [Google Scholar] [CrossRef]
  13. Lan, D.; Liu, R.; Zeng, W. Technique to Regional Risk Zonation of Environmental Pollution Accident and Its Application. J. Basic Sci. Eng. 2009, 17, 82–91. [Google Scholar]
  14. Qu, C.; Bi, J.; Ge, Y. Optimization Management of Regional Environmental Risk Based on Risk System. Environ. Sci. Technol. 2009, 32, 169. [Google Scholar]
  15. Xue, P.L.; Zeng, W.H. Shanghai environmental accidents risk regionalization. China Environ. Sci. 2011, 31, 1743–1750. [Google Scholar]
  16. Ma, Y.; Peng, J.F.; Song, Y.; Wang, C. Study on the classification method of environmental risk for sudden accidents in drinking water sources. Acta Sci. Circumst. 2012, 32, 1211–1218. [Google Scholar]
  17. Liu, R.; Liu, J.; Zhang, Z.; Borthwick, A.; Zhang, K. Accidental Water Pollution Risk Analysis of Mine Tailings Ponds in Guanting Reservoir Watershed, Zhangjiakou City, China. Int. J. Environ. Res. Public Health 2015, 12, 15269–15284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Martin, P.H.; LeBoeuf, E.J.; Daniel, E.B.; Dobbins, J.P.; Abkowitz, M.D. Development of a GIS-based Spill Management Information System. J. Hazard. Mater. 2004, 112, 239–252. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, J.; Hu, J.; Chu, J.; Chu, J.; Gao, X. Multi-scale risk assessment of accidental water pollution events in ungauged river basins. J. Tsinghua Univ. Sci. Technol. 2012, 52, 830–835. [Google Scholar]
  20. Golge, M.; Yenilmez, F.; Aksoy, A. Development of pollution indices for the middle section of the Lower Seyhan Basin (Turkey). Ecol. Indic. 2013, 29, 6–17. [Google Scholar] [CrossRef]
  21. Zhang, K.; Liu, R.; Zhang, Z.; Hu, L.; Jia, Z.; Zhuo, Y. A Method of Environmental Risk Assessment for Abrupt Water Pollution Accidents in River Basin. J. Basic Sci. Eng. 2014, 4, 675–684. [Google Scholar]
  22. Zeng, W.; Song, Y.; Yao, X. Multi-Scale Environmental Pollution Emergency Risk Zoning; Science Press: Beijing, China, 2013; pp. 5–17. [Google Scholar]
  23. Li, Y.; Huang, S.; Joseph, D. Temporal and spatial distributions of water and sediment yield in the Luanhe River Basin, China. In Proceedings of the 3rd International Conference on Water Resource and Environment (WRE2017), Qingdao, China, 26–29 June 2017. [Google Scholar]
  24. Li, Y.; Huang, S. Landscape Ecological Risk Response to Land Use Change in the Luanhe River Basin, China. Sustainability 2015, 7, 16631–16652. [Google Scholar] [CrossRef]
  25. General Office of the State Council. National Emergency Plan for Environmental Emergencies. 2014. Available online: http://www.gov.cn/zhengce/content/2015-02/03/content_9450.htm (accessed on 19 December 2014).
  26. Chen, Y.; Qiu, X. Research of risk source recognition of water environment in typical river basin of northwest inland. J. Water Resour. Water Eng. 2014, 25, 188–190. [Google Scholar]
  27. Tang, X.; Liu, B.; You, H.; Li, E. New approach to attacking the unexpected water environmental pollution hazardszone-division. J. Saf. Environ. 2013, 13, 276–278. [Google Scholar]
  28. Zhang, X.; Zhu, J. Problem and Countermeasures of Water Ecology and Water Environment. Environ. Sci. Manag. 2012, 8, 7–9. [Google Scholar]
  29. Gao, P.; Song, Y.; Yang, C. Water function zoning and water environment capacity analysis on surface water in Jiamusi urban area. Procedia Eng. 2012, 28, 458–463. [Google Scholar]
  30. Jiang, Y.; Yuan, D.; Zhang, C.; Kuang, M.; Wang, J.; Xie, S.; Zhang, G.; He, R. Impact of Land Use Change on Soil Proper ties in a Typical Karst Agricultural Region: A Case Study of Xiaojiang Water shed, Yunnan. Acta Geogr. Sin. 2005, 60, 751–760. [Google Scholar]
  31. Seydehmet, J.; Lv, G.-H.; Abliz, A.; Shi, Q.-D.; Turup, A. Irrigation salinity risk assessment and mapping in arid oasis, Northwest China. Water 2018, 10, 966. [Google Scholar] [CrossRef]
  32. Du, J.; He, F.; Shi, P. Integrated flood risk assessment of Xiangjiang River Basin in China. J. Nat. Disasters 2006, 15, 38–44. [Google Scholar]
Figure 1. Upper-middle reaches in the Luanhe River.
Figure 1. Upper-middle reaches in the Luanhe River.
Water 10 01268 g001
Figure 2. Schematic diagram of control unit division.
Figure 2. Schematic diagram of control unit division.
Water 10 01268 g002
Figure 3. The upper-middle reaches of Luanhe River and water function area distribution map.
Figure 3. The upper-middle reaches of Luanhe River and water function area distribution map.
Water 10 01268 g003
Figure 4. Working unit partition map in upper-middle reaches of Luanhe River.
Figure 4. Working unit partition map in upper-middle reaches of Luanhe River.
Water 10 01268 g004
Figure 5. Risk source hazard index system.
Figure 5. Risk source hazard index system.
Water 10 01268 g005
Figure 6. Index system of working unit vulnerability.
Figure 6. Index system of working unit vulnerability.
Water 10 01268 g006
Figure 7. Risk grade matrix diagram.
Figure 7. Risk grade matrix diagram.
Water 10 01268 g007
Figure 8. Risk sources hazard grade map.
Figure 8. Risk sources hazard grade map.
Water 10 01268 g008
Figure 9. Working unit hazard grade distribution map.
Figure 9. Working unit hazard grade distribution map.
Water 10 01268 g009
Figure 10. Working units vulnerability grade distribution map.
Figure 10. Working units vulnerability grade distribution map.
Water 10 01268 g010
Figure 11. Working unit risk grade distribution map.
Figure 11. Working unit risk grade distribution map.
Water 10 01268 g011
Figure 12. Risk zoning in the upper-middle reaches of Luanhe River based on water function area.
Figure 12. Risk zoning in the upper-middle reaches of Luanhe River based on water function area.
Water 10 01268 g012
Table 1. Water function area control unit name.
Table 1. Water function area control unit name.
No.The Name of the Water Function Area Control UnitNo.The Name of the Water Function Area Control Unit
1Daheiting reservoir18Baohe River Water source protection area
2Duolun River Agricultural water area19Sahe river reserve area
3Heifeng river reserve area20Duolun county Agricultural water area
4Heifeng River Agricultural Water area21Shandian River Water source protection area
5Laoniu River drinking water source area22Mongolia Hebei Shandian River buffer area
6Liuhe river buffer area23Shandian river reserve area
7Liuhe river drinking water source area24Shuluun Huh Agricultural water area
8Luanhe River Agricultural water area25Tuligen River industrial water area
9Luanhe river buffer area26Tuligen River Water Conservation area
10Luanhe river reserve area 127Wulie river reserve area
11Luanhe river reserve area 228Wulie River industrial water area
12Luanhe river drinking water source area29Wulie river drinking water source area
13Mongolia Hebei Luanhe River buffer area30Xiaoluanhe River Water source protection area
14Temple palace reservoir31Xingzhou River Water source protection area
15Panjiakou reservoir32Yixun river drinking water source area
16Baohe river buffer area33Yixun River Water source protection area
17Baohe river drinking water source area34Yimatu river reserve area
Table 2. Index classification standard and weight of the hazard index.
Table 2. Index classification standard and weight of the hazard index.
Topic LayerIndex LayerWeight1st2nd3rd4th
Risk of environmental risk sourcesindustry categories0.04850.80.60.40.2
Risk source size (Ten million yuan)0.088210050101
Distance from the river (km)0.05140.10.81.53
Downstream features0.026710631
Hazard receptorSensitive protection objectives0.0797181260
Water ecological risk0.027410.80.60.4
Exposure population (ten thousand people)0.1483151051
River grade0.041610.80.60.4
Risk control ability Risk source ability to control the hazard0.25191234
Per capita Gross Domestic Product (GDP) of the county (dollar/person)0.14984800860012,30016,100
The number of beds per 1000 people in the medical institutions(Bed/thousand people)0.051413610
Table 3. Statistical table of risk sources risk index.
Table 3. Statistical table of risk sources risk index.
Working Unit Water Function Area CharacteristicDrinkingProtection, ReservationBufferIndustry, Agriculture
Score *10.80.50.3
Sensitive protection objectives200 thousand population drinking water source or National Nature Reserve150 thousand population drinking water source or provincial nature reserve100 thousand population drinking water source or municipal Nature ReserveOther or no nature reserves
Score *9630
Characteristic pollutantsCODNutrimentOther chemicalsChromium (six valence)
Score *0.40.40.61
River gradeTrunk streamtributaries of level onetributaries of level twotributaries of level three
Score *10.80.60.4
Risk source ability to control the hazard Pollution source categorySigning the letter of environmental protection undertakingEnvironmental protection illegal informationPollutant discharge permit
Score *1 National gradeYesNoYes
0Other gradeNoYesNo
* Score in the table shows the quantitatively relative importance, and is generally judged subjectively [7,26].
Table 4. Risk source grade and working units hazard grade and score corresponding table.
Table 4. Risk source grade and working units hazard grade and score corresponding table.
Risk Source GradeGWorking Units Hazard Grade F
major risk source10Higher hazardF ≥ 8
big risk source6High hazard4.5 ≤ F < 8
general risk source3Moderate hazard2 ≤ F < 4.5
slight risk source1Low hazard0 ≤ F < 2
Table 5. Vulnerability index classification standard and weight table.
Table 5. Vulnerability index classification standard and weight table.
Vulnerability IndexWeight1st2nd3rd4th
Rainfall (mm)0.12121390490600700
Sensitivity0.180921.81.30.70.2
Water quality0.148132345
Per capita Gross Domestic Product(GDP) (dollar/person)0.1587416,10012,30086004800
Per capita cultivated land area (mu/person)0.169656.23.52.81.1
Population density (person/square kilometer)0.2215639228016552
Table 6. Working unit sensitivity assignment table.
Table 6. Working unit sensitivity assignment table.
Sensitive TargetScoreNature Reserve GradeScore
Drinking1National1
Protection, reservation0.8provincial0.5
buffer0.5Below the provincial0
Industry, agriculture0.3//
Table 7. Environment risk source name.
Table 7. Environment risk source name.
No.Environment Risk Source NameNo.Environment Risk Source Name
1Tianbao cement Ltd.23Chengde city sewage treatment Co., Ltd.
2Jidong Cement Co., Ltd.24Chengde municipal wastewater treatment plant
3Sheng Feng Steel Co., Ltd.25Chengde county sewage treatment plant
4Fuzhenquan Mining Co. Ltd.26Changhong Potato Starch Co. Ltd.
5New brewing Co., Ltd.27Weichang County heating company
6Sihai develops Limited by Share Ltd.28Binfa Potato Starch Co., Ltd.
7Tongxing mining limited liability company29Hengmao Paper Co. Ltd.
8Shuang Fu Mining Co., Ltd.30Jianfeng mining
9Ring energy thermoelectric Co., Ltd.31Jiujiang coal storage and Transportation Co., Ltd.
10Kaixing energy Co. Ltd.32Pingquan sewage treatment plant
11Pingan coal mine33Bishushangzhuang group Wine Co., Ltd.
12Hangyang gas Co. Ltd.34Yingke Fine chemicals Limited by Share Ltd.
13Wang Zhuang collective economy company35China traction group Chengde crankshaft connecting rod Co., Ltd.
14Xinglong County Sewage treatment plant36Heishan Iron Mine
15Liu Yuan sewage treatment Co., Ltd.37Luanping sewage treatment plant
16Datang Duolun Coal Chemical Co., Ltd.38Yushijindan Pharmaceutical Co. Ltd.
17Longxiao heating limited liability company39Jiangyuan Liquor Co., Ltd.
18Jianlong Special Steel Co., Ltd.40Fuze Food Co. Ltd.
19Green Food Co., Ltd.41Hongyuan fruit industry Co., Ltd.
20Qixing Mining Group Co. Ltd.42Longhua County sewage treatment plant
21Guodian Chengde Thermal Power Co., Ltd.43Weichang County wastewater treatment plant
22Chengde iron and Steel Group Co., Ltd.

Share and Cite

MDPI and ACS Style

Xiao, Y.; Huang, S.; Zhou, J.; Kong, F.; Liu, M.; Li, Y. Risk Assessment of Upper-Middle Reaches of Luanhe River Basin in Sudden Water Pollution Incidents Based on Control Units of Water Function Areas. Water 2018, 10, 1268. https://doi.org/10.3390/w10091268

AMA Style

Xiao Y, Huang S, Zhou J, Kong F, Liu M, Li Y. Risk Assessment of Upper-Middle Reaches of Luanhe River Basin in Sudden Water Pollution Incidents Based on Control Units of Water Function Areas. Water. 2018; 10(9):1268. https://doi.org/10.3390/w10091268

Chicago/Turabian Style

Xiao, Yao, Suiliang Huang, Jianguo Zhou, Fanqing Kong, Mingzhe Liu, and Ying Li. 2018. "Risk Assessment of Upper-Middle Reaches of Luanhe River Basin in Sudden Water Pollution Incidents Based on Control Units of Water Function Areas" Water 10, no. 9: 1268. https://doi.org/10.3390/w10091268

APA Style

Xiao, Y., Huang, S., Zhou, J., Kong, F., Liu, M., & Li, Y. (2018). Risk Assessment of Upper-Middle Reaches of Luanhe River Basin in Sudden Water Pollution Incidents Based on Control Units of Water Function Areas. Water, 10(9), 1268. https://doi.org/10.3390/w10091268

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