Determination of Key Risk Supervision Areas around River-Type Water Sources Affected by Multiple Risk Sources : A Case Study of Water Sources along the Yangtze ’ s Nanjing Section

To provide a reference for risk management of water sources, this study screens the key risk supervision areas around river-type water sources (hereinafter referred to as the water sources) threatened by multiple fixed risk sources (the risk sources), and establishes a comprehensive methodological system. Specifically, it comprises: (1) method of partitioning risk source concentrated sub-regions for screening water source perimeter key risk supervision areas; (2) approach of determining sub-regional risk indexes (SrRI, which characterizes the scale of sub-regional risks) considering factors like risk distribution intensity within sub-regions, risk indexes of risk sources (RIRS, characterizing the risk scale of risk sources) and the number of risk sources; and (3) method of calculating sub-region’s risk threats to the water sources (SrTWS) which considers the positional relationship between water sources and sub-regions as well as SrRI, and the criteria for determining key supervision sub-regions. Favorable effects are achieved by applying this methodological system in determining water source perimeter sub-regions distributed along the Yangtze’s Nanjing section. Results revealed that for water sources, the key sub-regions needing supervision were SD16, SD06, SD21, SD26, SD15, SD03, SD02, SD32, SD10, SD11, SD14, SD05, SD27, etc., in the order of criticality. The sub-region with the greatest risk threats on the water sources was SD16, which was located in the middle reaches of Yangtze River. In general, sub-regions along the upper Yangtze reaches had greater threats to water sources than the lower reach sub-regions other than SD26 and SD21. Upstream water sources were less subject to the threats of sub-regions than the downstream sources other than NJ09B and NJ03.


Introduction
With the progression of urbanization and industrialization, 78% of China's drinking water sources, as a highly sensitive water environment risk receptor, have been distributed with high water pollution risk enterprises in upstream [1].Meanwhile, due to the frequent occurrence of environmental pollution Sustainability 2017, 9, 283; doi:10.3390/su9020283www.mdpi.com/journal/sustainabilityaccidents attributed to the level of enterprise management and other factors, a sharp increase in the probability of water source contamination emerges.Many cases have caused varying degrees of harms to the downstream water sources including Petro China Jilin's benzene leakage accident occurred in the Songhua River in November 2005, Sinopec's phosphoric acid leakage incident at the Changshan oil terminal wharf of Yangtze's Jiangyin section in September 2008, and LCR petrochemical corporation's phenol leakage accident at the wharf of Yangtze's Zhenjiang-Yangzhou section in February 2012.Nevertheless, current China's water resource regulatory authorities generally face problems such as lack of management techniques and poor operability [2].Meanwhile, in the industrially-developed eastern areas, such as the areas along the Yangtze River in Jiangsu Province, water and risk sources are distributed alternately along the river.Moreover, there are numerous risk sources, which present a dense regional distribution.The above water source management problems and risk source distribution will lead to issues such as the difficulty in determining key risk supervision areas by water resource regulatory authorities; improper risk supervision for water sources; and inefficient risk management for water sources [3,4].Therefore, it is necessary to assess the risk sources around the water sources, and screen key risk sources need supervision as well as key risk regulation areas.This will provide a reference for relevant environmental regulatory authorities to improve the management efficiency of drinking water sources.Screening of water source perimeter key risk supervision areas requires assessment of risk scale at the regions where the risk sources lie.Currently, regional risk assessment includes the ecological risk assessment for river basins; regional risk assessment for non-point source distribution [5][6][7][8][9][10][11]; and regional risk assessment for point source distribution.The pollution of non-point sources to water is realized by forming surface runoff and farmland irrigation water through precipitation.Point sources mainly include the industrial sources like industrial enterprises and sewage treatment plants as well as the ports sources distributed along the navigable waterways.Moreover, the industrial sources mainly produce pollution to water by secretly discharging the unprocessed industrial wastewater, while the water pollution accidents like chemical and oil leakage in ports are mainly caused by production accidents.In industrially developed areas, due to low agricultural land use, there is little discharge of farmland irrigation water, and the high urban piped sewage rate also effectively lowers the pollution of early urban rainfall to water.Hence, the point source is the main factor causing environmental risks of water sources in industrially developed areas, and the regional risk assessment in these areas should be carried out targeting at the point source.At present, risk assessment of point source distribution areas is conducted mainly by assessing areas where single point source are located, or by evaluating risk scale of the same industry type industrial zones.For instance, Sadiq and Husain [12] developed a fuzzy-based methodology to assess the environmental risk of drilling waste.Jiang et al. [13] established an environmental risk evaluation index system based on fuzzy membership degree to assess the environmental risk of a given chemical enterprise in Shanghai, China.Liu et al. [14] evaluated the environmental risk of chromium plating enterprise based on the EU (European Union) technical guidance document and EUSES model.Meng et al. [15] applied the information-diffusion theory to carry out a regional environmental risk assessment for the Nanjing Chemical Industry Park, China.Shao et al. [16] developed a simple physical model using existing dispersal patterns and migration models to identify the environmental risks of Chemical Industry Parks in Tianjin Binhai New Area, China.Zhang et al. [17] measured the environmental risk of PBDEs around Industry Park in Huizhou, China.Nadal et al. [18] described the environmental risk of a chemical/petrochemical area in Tarragona, Spain based on GIS technology.Albanese et al. [19] applied GIS technology to conduct a risk assessment of heavy metal pollution in The Zambian Copper belt Province.However, these evaluation approaches either do not consider the distribution of multiple risk sources within a region, or set the multiple risk sources within a region as the same industry type.However, in reality, enterprises around the surface water sources in industrially developed areas, especially those around river-type water sources, are varied in industry type and large in quantity.Moreover, the density of risk distribution within a region will also impact the probability of risk accident occurrence.Therefore, assessment of regional risk sources should specifically consider the risk scale of risk sources within a region; the number of risk sources; and the density of risk distribution in the region.In addition, as the risk sources are non-uniformly distributed in the periphery of water sources, determination of the threats of risk source concentrated local regions (hereinafter referred to as the sub-region) to the water sources will be difficult if the perimeter zone is assessed as a whole.This will lead to difficulty of determining the key supervision sub-regions for water source risk control, thereby reducing the regulatory authorities' efficiency of water source risk supervision.
Based on the above analysis, this paper constructs a methodological system for determination of key risk supervision sub-regions around water sources.The system takes into account the effects of multiple risk sources.Its main contents include: (1) method of partitioning risk source concentrated sub-regions; (2) method of identifying SrRI; and (3) method of calculating SrTWS, and the criteria for determining key supervision sub-regions.Meanwhile, the method is applied to eight water sources along the Yangtze River in Nanjing.

Study Area
Nanjing (31 • 14'N-32 • 36'N, 118 • 22'E-119 • 14'E) lies in the Yangtze River Delta, one of China's three major economic core areas, which is an important central city in the lower reach of the Yangtze River.Nanjing section of the Yangtze River is about 95 km long.As a lower tidal reach of the Yangtze, it is impacted by medium-intensity tides.Every day, there are two tidal peaks and two tidal troughs.Flood tides last about three hours, whereas ebb tides last about nine hours.Reverse flow exists during flood tides.There are a total of 10 county-level or above centralized drinking water sources under the jurisdiction of Nanjing, of which eight (two standby water sources) are distributed along the banks of Yangtze River, namely Jiajiang Water Source (NJ01), Yanziji Water Source (NJ02), Bagua Island Water Source (NJ03), Pukou Water Source (NJ04), Jiangning Zihui Island Water Source (NJ07), Longtan Water Source (NJ08), Bagua Island Standby Water Source (NJ09B) and Qiaolin Standby Water Source (NJ10B).Total designed withdrawal of these water sources is 1.1925 billion t/year, which supply 93.77% of the municipal water in Nanjing.
According to the Pollution Source Census Data 2013-2015 from the Environmental Protection Department of Jiangsu Province, the port wharf data for Yangtze's Jiangsu section from the Jiangsu Maritime Safety Administration, as well as field survey: there are a total of 170 risk sources distributed along the banks of Yangtze's Nanjing section, of which 35 are industrial enterprises' sewage outlets; 120 are wharfs; and 15 are wastewater treatment plants' outlets.Thirteen of 35 industrial enterprises (total wastewater discharge of 3,153,000 t/day) were high-risk industrial enterprises (total wastewater discharge of 220,100 t/day, with six enterprises having a wastewater discharge of over 2000 t/d), which were primarily petrochemical, chemical and pharmaceutical enterprises with certain automatic water quality monitoring abilities.Thirty-one of 120 wharfs (total berthing capacity of 3,395,000 t) were petrochemical and chemical wharfs (total berthing capacity of 619,000 t, with seven wharfs having a berthing capacity of above 40,000 t).Among them, only 18 wharfs were equipped with anti-pollution facilities.All 15 wastewater treatment plants had a wastewater discharge of greater than 2000 t/day (total discharge of 1,338,000 t/day), with complete automatic water quality monitoring abilities.Moreover, as the water treated was domestic sewage, the discharged pollutants were mainly COD, NH 3 -N, TN and TP.
Specific locations of the water sources and risk sources are shown in Figure 1.

Methodology
Risk source concentrated sub-regions were partitioned according to the distribution of risk sources around the water sources.Then, SrRI were calculated.SrTWS were calculated based on the SrRI and the positional relationship between the risk sources within sub-regions and the water sources.Key risk supervision sub-regions around water sources were screened.The specific procedure is shown in Figure 2. Step 1 Step 2 Step 3 Step 4

Methodology
Risk source concentrated sub-regions were partitioned according to the distribution of risk sources around the water sources.Then, SrRI were calculated.SrTWS were calculated based on the SrRI and the positional relationship between the risk sources within sub-regions and the water sources.Key risk supervision sub-regions around water sources were screened.The specific procedure is shown in Figure 2.

Methodology
Risk source concentrated sub-regions were partitioned according to the distribution of risk sources around the water sources.Then, SrRI were calculated.SrTWS were calculated based on the SrRI and the positional relationship between the risk sources within sub-regions and the water sources.Key risk supervision sub-regions around water sources were screened.The specific procedure is shown in Figure 2. Step 1 Step 2 Step 3 Step 4 In addition, according to the relevant provisions of the Law of the People's Republic of China on Prevention and Control of Water Pollution [20], risk sources uncorrelated with water supply and water source protection, such as sewage outlets, direct discharge enterprises and loading and unloading wharfs, are prohibited within primary and secondary water source protection areas.Therefore, if the above-mentioned risk sources appear in the primary and secondary protection areas of water sources, the risk source distribution zones within the areas will be partitioned as separate sub-regions.Moreover, all of these sub-regions are the key risk supervision sub-regions irrespective of their risk indexes.

Minimum Spanning Tree (MST)-Based Method for Partitioning Risk Source Concentrated Sub-Regions
When the risk sources are regarded as the nodes and the lengths of lines linking them are regarded as the weighs, all of the lines linking risk sources constitute a weight map.If we want to connect all the nodes in the weight map with the least number of links, the links selected will have to form a tree [21].By comparing a certain distance threshold to the distance between two adjacent risk sources on the spanning tree (i.e., weight of the two risk sources), we can determine whether these two risk sources are within the same sub-region.Among all spanning trees, the one with the minimum weighted sum is the minimum spanning tree [21].Therefore, it is most accurate to determine whether the adjacent risk sources are within the same sub-region based on the distance between adjacent risk sources on the minimum spanning tree.
If a risk source and its adjacent risk sources are within the same sub-region, the distance between it and its adjacent risk sources should be within a certain threshold; and the minimum distances between risk sources within this sub-region and risk sources within other sub-regions should all be greater than this threshold.A typical example is the risk source f 4 as shown in Figure 3a.Its distances from the risk sources f 5 and f 6 are 1.14 km and 1.62 km, respectively (Table 1), while its distance from the risk source f 2 is 2.74 km (Table 1).Assuming a threshold of 2 km, then the risk sources f 4 -f 6 should be in the same sub-region, whereas the risk source f 2 is in another sub-region.In addition, there is also the case where the distances between a certain risk source and its adjacent risk sources do not differ much, but some of the distances are slightly greater than the threshold.Various risk sources involved in such case can also be partitioned into one sub-region.A typical example is the risk source f 2 as shown in Figure 3a.Its distances from the risk sources f 1 and f 3 are 1.89 km and 2.16 km, respectively (Table 1).If only compared to the threshold of 2 km, the risk source f 3 should not be partitioned into the same sub-region with sources f 1 -f 2 .However, since 1.89 km and 2.16 km are not much different, and 2.16 km is only slightly larger than 2 km, we can also consider that the risk sources f 1 -f 3 are in the sub-region in this case.Therefore, based on the above analysis, we can determine the sub-region partitioning conditions on the basis of constructing the minimum spanning tree of risk sources by comprehensively considering two aspects: distances between risk source and its adjacent risk sources; and deviation between the distances.At present, there are already some classic minimum spanning tree algorithms, such as the Kruskal method [22] and the Prim method [23].In this paper, the Prim algorithm [23] is employed to construct the minimum spanning tree of risk sources.The specific steps are described below.The route between adjacent risk sources on the risk source minimum spanning tree was defined as the node path, whose length was the node distance.With the node distance on the entire risk source minimum spanning tree as the object, if partitioning conditions were satisfied, the corresponding node path would be interrupted.Accordingly, the original single risk source minimum spanning tree would be divided into several minimum spanning trees.Risk sources connected in series by each resulting minimum spanning tree was partitioned into one sub-region.
Sustainability 2017, 9, 283 6 of 24 node path would be interrupted.Accordingly, the original single risk source minimum spanning tree would be divided into several minimum spanning trees.Risk sources connected in series by each resulting minimum spanning tree was partitioned into one sub-region.

Construction of Risk Source Minimum Spanning Tree
Prim algorithm [23] was used to construct the risk source minimum spanning tree by assuming there were n number of risk sources in the survey area around water source, and the risk source set was F = {f1, f2, … , fn}.The specific steps were as follows: (1) Risk source minimum spanning tree connected the risk sources within the entire survey area in series gradually.The risk source set that had been serially connected by the minimum spanning tree was set as V; the risk source set yet to be serially connected by the minimum spanning tree was set as U; and the node path set on the minimum spanning tree was set as L. Initially, V = { }, U = F, L = { }. (2) fi was added to V starting from any risk source fi.
(3) Risk source fj nearest to all risk sources in V was found from U, and the node path of fj serially connected by minimum spanning tree was added to L. (4) fj was added to V, U = F − V.
(5) Steps (3) and (4) were repeated until U = { }. Specific calculation formula for risk source spacing was as follows: where fix and fiy are the x and y coordinates of risk source fi; and fjx and fjy are the x and y coordinates of risk source fj, respectively.

Judgment Criteria for Sub-Region Partition
From the perspective of the entire survey area, whether two risk sources could be partitioned into the same sub-region depended on the contrast of distance between them to the entire survey area size.If the distance between the two was very small relative to the entire survey area, the two risk sources could be partitioned together.Otherwise, the two risk sources should be subordinate to different sub-regions.Thus, when determining the partitioning distance threshold dmin for judging whether two risk sources were within the same sub-region, the size of the entire survey area should be taken into account.

Construction of Risk Source Minimum Spanning Tree
Prim algorithm [23] was used to construct the risk source minimum spanning tree by assuming there were n number of risk sources in the survey area around water source, and the risk source set was F = {f 1 , f 2 , . . ., f n }.The specific steps were as follows: (1) Risk source minimum spanning tree connected the risk sources within the entire survey area in series gradually.The risk source set that had been serially connected by the minimum spanning tree was set as V; the risk source set yet to be serially connected by the minimum spanning tree was set as U; and the node path set on the minimum spanning tree was set as L. Initially, V = { }, U = F, L = { }. (2) f i was added to V starting from any risk source f i .
(3) Risk source f j nearest to all risk sources in V was found from U, and the node path of f j serially connected by minimum spanning tree was added to L. (4) f j was added to V, U = F − V.
(5) Steps ( 3) and ( 4) were repeated until U = { }. Specific calculation formula for risk source spacing was as follows: where f ix and f iy are the x and y coordinates of risk source f i ; and f jx and f jy are the x and y coordinates of risk source f j , respectively.

Judgment Criteria for Sub-Region Partition
From the perspective of the entire survey area, whether two risk sources could be partitioned into the same sub-region depended on the contrast of distance between them to the entire survey area size.If the distance between the two was very small relative to the entire survey area, the two risk sources could be partitioned together.Otherwise, the two risk sources should be subordinate to different sub-regions.Thus, when determining the partitioning distance threshold d min for judging whether two risk sources were within the same sub-region, the size of the entire survey area should be taken into account.
Risk sources posing threat to water sources were generally distributed along the shoreline.Overall, the spatial size of survey area could be represented by the unilateral shoreline length S of the river in the survey area where the water source was located.Ratio of d min to S was set as 0.1, then: According to the provision of the Ministry of Environmental Protection's Guidelines for Protection of Centralized Drinking Water Source Environment [24], the risk source survey area around water sources covers a 20 km range of upstream of secondary water source protection area.In this paper, S was set as 20 km, so d min was 2 km.
In addition, if the ratio d i /d j of node distances d i and d j corresponding to two connected node paths L i and L j was within a certain range, the ratio range could be set between [0.83, 1.2].When , and one of the node distances was less than or equal to d min , the three risk sources connected in series by L i and L j were partitioned into one sub-region.
Taking into comprehensive consideration the distance between nodes and the deviation between distances of their interconnected nodes, the judgment criteria for sub-region partition were set as follows: where d i is the node distance of node path L i ; L j is the node path with the largest node distance from L i , and d j is its node distance.

Sub-Region Partitioning Procedure
(1) Prim algorithm [23] was used to construct the risk source minimum spanning tree by assuming there were n number of risk sources in the survey area around water sources.Node path set of the minimum spanning tree was defined as L = {L 1 , L 2 , . . ., L n−1 }, while the node distance set corresponding to node paths was defined as D = {d 1 , d 2 , . . ., d n− }. (2) Node distance d i was determined one by one from large to small for whether it satisfied the Equation (3).If satisfied, the node path L i would be interrupted.Accordingly, the minimum spanning tree where L i was located was decomposed into two minimum spanning trees.(3) Risk sources connected in series by the resulting minimum spanning trees were partitioned into one sub-region.

Example of Sub-Region Partitioning
Taking the risk sources f 1 -f 6 in Figure 3a as an example, the risk source minimum spanning tree was constructed from the risk source f 1 .The construction process and results are shown in Figure 3 and Table 1, whereas the minimum spanning tree of risk source constructed is shown in Figure 4.
By comparing the various path distances in Table 1 to the distance threshold 2 km, it was observed that the paths L 2 and L 3 were both greater than the distance threshold.However, the nodal distance ratio of paths L 2 to L 1 (node path with the largest nodal distance from L 2 ) was 1.14, while the nodal distance of L1 was 1.89, which was less than the distance threshold.Thus, L 2 did not satisfy Equation (3), and only path L 3 satisfied the interrupt condition.Through interrupting path L 3 , two sub-regions could be obtained, of which risk sources f 1 -f 3 were partitioned into a sub-region A, and risk sources f 4 -f 6 were partitioned into another sub-region B. The specific partitioning results are shown in Figure 4.

Method for Determining Sub-Regional Risk Indexes
Determination of SrRI should be based on the RIRS of individual risk sources within the subregions.Without considering other influencing factors, the greater the mean value about risk indexes of risk sources (Mean-RIRS) in sub-regions, the greater the value of SrRI theoretically.In addition, considering that risk sources with larger RIRS had far greater risk impact on the sub-regions than those with smaller RIRS, the maximal value about risk indexes of risk sources (Max-RIRS) in subregions should also be taken into account in determining SrRI.However, consideration of the effects of RIRS for various risk sources in the sub-regions only was unable to accurately reflect SrRI.If RIRS for individual risk sources in the sub-regions were the same, the larger the number of risk sources in the sub-regions, the greater the sum of RIRS, the higher the risk threats to the surrounding waters, and the greater the SrRI theoretically.Therefore, besides RIRS for various risk sources in the subregions, the influences of risk source number on SrRI, and the intensive industrial distribution and close distance between risk sources in the economically developed areas need to be considered.Hence, in case of fire, explosion or other accidents at a risk source, sequential fire or explosion at multiple risk sources may be triggered (such as the Tianjin Port explosion that occurred in August 2015).Compared to the fire or explosion accidents at a single risk source, sequential fire or explosion at multiple risk sources will greatly increase the probability of water pollution accidents.Consequently, when determining SrRI, the concentration degree of risk source distribution in subregions should also be considered.
Based on the above analysis, determination of sub-regional risk indexes needs to take into account factors such as Mean-RIRS, Max-RIRS, number of risk sources and degree of risk concentration within a sub-region.The specific calculation formula was as follows: where Kt is the risk index of the t-th sub-region; and m is the number of risk sources in the t-th subregion.t k and kti are the average risk index of risk sources and the risk index of i-th risk source in the t-th sub-region, respectively.They reflected the risk scale of individual risk sources within a subregion.kmax is the maximum value that the risk index of individual risk sources could assume.Gt is the risk distribution density index (RDDI) of the t-th sub-region, which reflected the relative degree of risk concentration in the t-th sub-region.Dt is the risk source quantity index (RSQI) of the t-th subregion, which reflected the risk increment in the t-th sub-region caused by increased number of risk sources.

Method for Determining RIRS
RIRS should be determined by comprehensively considering influential factors such as industry type, production scale, technological level, wastewater complexity and risk supervision and emergency response capability [25].However, most of the above factors can only be qualitatively described for their scale of risk threats to the surrounding areas except for the production scale (industrial source can be represented by wastewater discharge capacity, and wharf source can be represented by berthing capacity).Hence, it can hardly be quantitatively determined.In contrast, the

Method for Determining Sub-Regional Risk Indexes
Determination of SrRI should be based on the RIRS of individual risk sources within the sub-regions.Without considering other influencing factors, the greater the mean value about risk indexes of risk sources (Mean-RIRS) in sub-regions, the greater the value of SrRI theoretically.In addition, considering that risk sources with larger RIRS had far greater risk impact on the sub-regions than those with smaller RIRS, the maximal value about risk indexes of risk sources (Max-RIRS) in sub-regions should also be taken into account in determining SrRI.However, consideration of the effects of RIRS for various risk sources in the sub-regions only was unable to accurately reflect SrRI.If RIRS for individual risk sources in the sub-regions were the same, the larger the number of risk sources in the sub-regions, the greater the sum of RIRS, the higher the risk threats to the surrounding waters, and the greater the SrRI theoretically.Therefore, besides RIRS for various risk sources in the sub-regions, the influences of risk source number on SrRI, and the intensive industrial distribution and close distance between risk sources in the economically developed areas need to be considered.Hence, in case of fire, explosion or other accidents at a risk source, sequential fire or explosion at multiple risk sources may be triggered (such as the Tianjin Port explosion that occurred in August 2015).Compared to the fire or explosion accidents at a single risk source, sequential fire or explosion at multiple risk sources will greatly increase the probability of water pollution accidents.Consequently, when determining SrRI, the concentration degree of risk source distribution in sub-regions should also be considered.
Based on the above analysis, determination of sub-regional risk indexes needs to take into account factors such as Mean-RIRS, Max-RIRS, number of risk sources and degree of risk concentration within a sub-region.The specific calculation formula was as follows: where K t is the risk index of the t-th sub-region; and m is the number of risk sources in the t-th sub-region.k t and k ti are the average risk index of risk sources and the risk index of i-th risk source in the t-th sub-region, respectively.They reflected the risk scale of individual risk sources within a sub-region.k max is the maximum value that the risk index of individual risk sources could assume.G t is the risk distribution density index (RDDI) of the t-th sub-region, which reflected the relative degree of risk concentration in the t-th sub-region.D t is the risk source quantity index (RSQI) of the t-th sub-region, which reflected the risk increment in the t-th sub-region caused by increased number of risk sources.

Method for Determining RIRS
RIRS should be determined by comprehensively considering influential factors such as industry type, production scale, technological level, wastewater complexity and risk supervision and emergency response capability [25].However, most of the above factors can only be qualitatively described for their scale of risk threats to the surrounding areas except for the production scale (industrial source can be represented by wastewater discharge capacity, and wharf source can be represented by berthing capacity).Hence, it can hardly be quantitatively determined.In contrast, the semi-quantitative comprehensive index method is based on the calculation of a comprehensive evaluation index that summarizes the indexes of multiple risk elements using weight values [15].This semi-quantitative method has been widely used to assess the risk of chemical and petrochemical areas [18] and mining areas [26].Therefore, in this paper, RIRS can be assessed employing the comprehensive index method.The specific assessment steps were as follows: (1) Each assessment index of individual risk sources was graded according to relevant risk grading criteria and scored.The risk levels were classified into four categories: very low, low, medium and high, which corresponded to fours scores: 1, 2, 3 and 4, respectively; (2) RIRS were determined by weighted summation based on the weight and score values of indexes.Since the sum of index weights was 1, and the maximum score for each index was 4, the maximum value reachable by RIRS was 4, i.e., k max = 4.The risk grading criteria and specific weights of indexes are shown in Table 2 [25].Risk levels of risk sources were determined according to the RIRS based on Table 2 [25].Table 2 could also be used as the sub-regional risk grading criteria.
Table 2. Risk source and sub-regional risk grading criteria.

Risk Level
High Medium Low Very Low 3.2.2.Method for Determining Sub-Regional RDDI RDDI of multiple risk sources within sub-region was determined based on the relative degree of risk concentration between adjacent risk sources on the risk source minimum spanning tree.
Relative degree of risk concentration between adjacent risk sources connected by the i-th node path L ti on the risk source minimum spanning tree in the t-th sub-region could be determined by comparing the node distance d ti corresponding to L ti with the risk distance threshold d max .The specific formula was as follows: where w ti is the index reflecting the relative degree of risk concentration between adjacent risk sources connected by the node path L ti , and w ti ∈ (0,1].When w ti was equal to 1, it indicated that the two risk sources were completely concentrated; otherwise, it indicated relative concentration.The closer the value of w ti to 0 was, the lower the degree of concentration between two risk sources would be.Risk distance threshold d max was determined based on the length of unilateral river shoreline S reflecting the spatial size of survey area.The specific formula was as follows: Ministry of Environmental Protection's Guidelines for Protection of Centralized Drinking Water Source Environment [27] provides that the risk source survey area around water sources covers a 20 km range of upstream of secondary water source protection area.S was considered as 20 km, so d max was considered as 1 km.
The RDDI of the t-th sub-region was calculated as follows: where G t is the RDDI of the t-th sub-region, and G t ∈(0, 1].The closer the value of G t to 0, the more dispersed the risks in the t-th sub-region, and vice versa.k u ti and k d ti denote the risk indexes of risk sources at the upper and lower ends of node path L ti , respectively.Meanings of m and w ti are the same as above.
3.2.3.Method for Determining Sub-Regional RSQI If the impact of risk source risk indexes was not considered, the more the number of risk sources within sub-region, the greater the water environmental risk in the sub-region.RSQI was determined by assuming that the increment of risk source quantity in the t-th sub-region was linearly related to the increment of sub-regional risk.Specific formula was shown below: where D t and m have the same meanings as above.

Method for Determining Water Source Perimeter Key Risk Supervision Sub-Regions
Determination basis of key risk supervision sub-regions around water sources was the sub-region's scale of risk threats to the water sources.In contrast, the aforementioned sub-regional risk index characterized the sub-region's scale of risk threats to its surrounding water environment sensitive receptors.The larger the SrRI, the severer the substandard condition of water for sensitive waters after pollution accidents, and the greater the risk threats to these sensitive waters.In this case, the receptor should be close to the sub-region.If the distance between the two increased, the sub-region's risk threats to the sensitive receptor should be reduced.Similarly, whether the sub-region and sensitive receptor were on the same shoreline; and whether the sub-region was located upstream of sensitive receptor would also impact the sub-region's scale of threats to the receptor.
Based on the above analysis, SrTWS was determined in this paper via the SrRI and the positional relationship between sub-region and water source.The positional relationship between sub-region and water source includes: distance between the sub-region and the water source in the direction of water body's forward flow, i.e., the x directional distance; whether the sub-region and the water source were on the same shoreline, if not, distance between the two in the vertical direction of water body's forward flow, i.e., the y directional distance; and in the case of reciprocating flow of water body, the upstream-downstream positional relationship between the sub-region and the water source.Sub-regional risk supervision level was determined according to SrTWS, thereby identifying the key risk supervision sub-regions around water sources.Risk supervision grading criteria were established based on the risk grading criteria for risk sources [25].The details are shown in Table 3.
Table 3. Sub-regional risk supervision grading criteria.

Sub-Regional Risk Supervision Level Key Supervision Secondary Key Supervision Non-Key Supervision
SrTWS ≥3 [1.5, 3) <1.5 Specific calculation formula for SrTWS was as follows: where WG t is the t-th sub-region's SrTWS; and X t and Y t , are the adjustment coefficients determined by considering the magnitudes of distances between the t-th sub-region and the water source in the x and y directions, respectively.When the sub-region and the water source were located on the same shoreline, Y t was 1. P t was the adjustment coefficient considering the upstream-downstream positional relationship between the t-th sub-region and the water source.G t had the same meaning as above.
According to the provision of the Ministry of Environmental Protection's Guidelines for Protection of Centralized Drinking Water Source Environment [24], the risk source survey area around water sources covers a 20 km range of upstream of secondary water source protection area.In this paper, SrTWS was considered equal to the sub-regional risk index when the x directional distance between sub-region and secondary water source protection area was within 10 km irrespective of other factors.When the x directional distance between sub-region and secondary water source protection area was greater than 10 km, SrTWS decreased linearly with the increasing distance.X t was determined according to this principle, and its calculation formulas were: where S tx is the integrated distance (km) from multiple risk sources within the t-th sub-region to water intake in the x direction.If the t-th sub-region was located upstream of the water source, S tx > 0; otherwise, S tx < 0. LW was the distance (km) from source water intake to the upstream and downstream secondary protection area boundary.When the sub-region was located downstream of the water source, the value of LW was the distance from water intake to downstream secondary protection area boundary.m was the number of risk sources within the t-th sub-region.s tix , k ti were the x directional distance (km) from i-th risk source within t-th sub-region to water intake and the risk index of i-th risk source, respectively.When the risk source was located upstream of the water source, the value of s tix was positive; otherwise, the value was negative.X t had the same meaning as above.
As the lateral flow of the river was small, the sub-regions located on the other shore of water sources had less threat on the water sources.Similarly, reflected in the setting of water source protection areas, the width of protection areas should be significantly less than their length.By referring to the method of determining X t coefficient, Y t in the case where sub-region and water source were not on the same shoreline was determined with protection area width as the standard.The specific calculation formulas were as follows: where S ty and WI p were the y directional integrated distance (km) from the t-th sub-region to the shoreline of water source location and the width (km) of secondary water source protection area at the source water intake; respectively.s tiy was the y directional distance (km) from i-th risk source within t-th sub-region to water intake.Y t and k ti had the same meanings as above.
For reciprocating flow, the sub-regions located downstream of water sources only posed risk threats on the water sources through reverse flow, whose risk threats were smaller compared to the sub-regions upstream of water sources.The influence of the upstream-downstream positional relationship between sub-regions and water sources on the SrTWS could be reflected by the forward and reverse flow durations of water sources.The specific formula was as follows: where T F and T O are the forward and reverse flow durations of water source, respectively; and P t has the same meaning as above.

Sub-Region Partitioning Results
Risk source minimum spanning tree was constructed using Prim algorithm [23] with the risk sources distributed along the Yangtze's Nanjing section as the nodes.Meanwhile, risk source-distributed sub-regions were partitioned.Figure 5 presents the risk source minimum spanning tree and the sub-region partition results.As shown in the figure, risk sources distributed along the Yangtze's Nanjing section can be partitioned into 32 sub-regions from SD01 to SD32.The major risk source-concentrated sub-regions were SD06, SD16, SD19, SD21 and SD26, which contained 14, 14, 21, 22 and 24 risk sources, respectively.In terms of distribution location, except for the sub-region SD06 which was located in the Jiangning Binjiang Industrial Park upstream of Nanjing's main urban district, the other four major sub-regions were located in the Xiaguan Port District, Luhe Chemical Park's riverside portion, Xinshengwei Foreign Trade Harbor District and Longtan Port District downstream of Nanjing's main urban district.Among them, sub-regions SD16 and SD26 were closely neighboring the downstream water source NJ02 and NJ08 protection areas.The remaining sub-regions were smaller than the above-mentioned five major sub-regions, with each containing 1-6 risk sources.Among them, sub-regions SD03, SD05, SD27 and SD32 were distributed within the protection areas of water sources NJ07, NJ10B, NJ08 and NJ09B, respectively.

Sub-Region Partitioning Results
Risk source minimum spanning tree was constructed using Prim algorithm [23] with the risk sources distributed along the Yangtze's Nanjing section as the nodes.Meanwhile, risk sourcedistributed sub-regions were partitioned.Figure 5 presents the risk source minimum spanning tree and the sub-region partition results.As shown in the figure, risk sources distributed along the Yangtze's Nanjing section can be partitioned into 32 sub-regions from SD01 to SD32.The major risk source-concentrated sub-regions were SD06, SD16, SD19, SD21 and SD26, which contained 14, 14, 21, 22 and 24 risk sources, respectively.In terms of distribution location, except for the sub-region SD06 which was located in the Jiangning Binjiang Industrial Park upstream of Nanjing's main urban district, the other four major sub-regions were located in the Xiaguan Port District, Luhe Chemical Park's riverside portion, Xinshengwei Foreign Trade Harbor District and Longtan Port District downstream of Nanjing's main urban district.Among them, sub-regions SD16 and SD26 were closely neighboring the downstream water source NJ02 and NJ08 protection areas.The remaining subregions were smaller than the above-mentioned five major sub-regions, with each containing 1-6 risk sources.Among them, sub-regions SD03, SD05, SD27 and SD32 were distributed within the protection areas of water sources NJ07, NJ10B, NJ08 and NJ09B, respectively.

Calculation Results of SrRI
A total of 170 risk sources were surveyed.The level and score of assessment indexes for each risk source were determined by referring to Appendix A Table A2.RIRS of various risk sources were determined by weighted summation based on the risk scores and weights of indexes (see Appendix A Table A2).Meanwhile, risk levels of risk sources were determined according to Table 2. Risk assessment results revealed that 29 of 170 risk sources were high risk sources, which were mainly wastewater treatment plants' outlets (15) and wharf risk sources (11).The 11 high-risk wharf loading and unloading articles were oils, bulk toxic liquids, etc.None of these wharfs were equipped with automatic water quality detector or pollution control facility, whose sensitivity and emergency response to sudden pollution accidents were poor.Among them, wharfs njm47, njm83, njm85, njm89 and njm90 all had a berthing capacity above 40,000 t.Therefore, RIRS of these five wharfs were the highest, which were all 3.389 (risk scores for six indexes: wharf type, wharf berthing capacity, technological level of production equipment, management system, emergency prevention system and section monitoring system were 4, 4, 3, 2, 3, and 4, respectively).Jingling Sinopec's first chemical plant (code 271) scored the second highest (3.359) for RIRS due to high industry risk and high complexity of wastewater discharged (for the seven indexes of the industry type, complexity of sewage quality, wastewater discharge, technological level of production equipment, management system, emergency prevention system and section monitoring system, the risk scores were 4, 3, 4, 3, 2, 3, and 4, respectively).As its annual wastewater discharge was nearly 830,000 t, its RIRS was higher than the other two high risk enterprises (both of which were petroleum processing enterprises with RIRS of both 3.197.For the seven indexes of the industry type, sewage quality complexity, wastewater discharge, production equipment level, management system, emergency prevention system and section monitoring system, the risk scores were 4, 3, 1, 4, 4, 3, 4, respectively).All 15 wastewater treatment plants distributed along the Yangtze's Nanjing section were high risk sources.Although annual wastewater discharge of each plant was at least 1,800,000 t, their discharge complexity was lower than the above-mentioned three high risk enterprises since they primarily deal with domestic wastewater.Moreover, they all had certain water quality monitoring capability.Thus, overall, the RIRS of wastewater treatment plants were slightly lower than the high-risk enterprises.Specific locations of high risk sources are shown in Figure 5.The specific RIRS results for high risk sources are listed in Table 4. SrRI were calculated based on the RIRS and locations of various risk sources; meanwhile, sub-regional risk levels (SrRL) were assessed.As shown in Table 5, among the 32 sub-regions, those with high, medium and low risk levels were 21, 10 and 1, respectively.The only low-risk sub-region was SD13, within which there was only one low-risk source, a public dock.
Statistical analysis chart (Figure 6) of SrRI in the high risk sub-regions was plotted.As can be seen, compared to other high risk sub-regions, the sub-regions SD06, SD16, SD19, SD21 and SD26 exhibited significantly larger SrRI owing to their advantages in quantity and RIRS (risk sources with the highest and the second highest RIRS were located in sub-regions SD26 and SD21, respectively).Although the sub-regions SD11, SD12 and SD25 each had only one risk source compared to other high risk sub-regions, they were also listed as high risk sub-regions as the risk sources were large-scale wastewater treatment plants.Nevertheless, their SrRI were slightly smaller than other high risk sub-regions.Although all the risk sources in the sub-regions SD10, SD20, SD27, SD29 and SD31 were medium risks (Max-RIRS of less than 3), SrRI of these sub-regions were greater than 3 due to a certain number of concentrated risk sources within each sub-region, which made them high risk sub-regions.
Statistical analysis chart (Figure 7) of SrRI in the medium risk sub-regions was drawn.Most of the medium risk sub-regions contained only 1 medium risk source.In the medium risk sub-regions SD04, SD05 and SD17, there were slightly more risk sources (3, 5 and 2, respectively).Although the risk sources in the sub-region SD17 were concentrated, the RIRS was small and there were only two risk sources, so the SrRI was also not too large.Since the risk sources were relatively dispersed in the two sub-regions SD04 and SD05, the SrRIof these two sub-regions were not large.Compared to other medium risk sub-regions, sub-region SD05 exhibited the smallest SrRI despite multiple risk sources because the risk sources within it were small-scale shipyards and public docks.

Determination of Key Risk Supervision Sub-Regions
The x and y directional distances from sub-regions to source water intakes were determined based on the risk indexes and locations of risk sources.Appendix A Table A3 lists the specific results.In Appendix A Table A3, the x directional distance from sub-region to water intake would be positive if the sub-region was located upstream of the water source; 0 if the sub-region was located within the water source protection area; and negative if located downstream of the water source.The y directional distance from sub-region to water intake would be 0 if the sub-region and the water source were on the same shoreline; and positive if were on the opposite shorelines.Some sub-regions such as SD14 were located in the Yangtze tributaries, so their x directional distances to water intakes were regarded as the sum of distance from sub-region to tributary estuary and distance from the estuary to water intake.Meanwhile, the y directional distance to water intake was considered as the y directional distance from tributary estuary to water intake.In addition, part of the sub-regions did not pose threats to some water sources as a result of island barrier.In such case, the distance from sub-region to water intake was not considered, such as the distances from the sub-regions SD11-SD14 to the water source NJ01 in Appendix A Table A3.According to the Jiangsu Provincial Program for Partitioning County Level or Above Centralized Drinking Water Source Protection Areas [27], the distances from the source water intake to the secondary protection area upstream and downstream boundaries were 2 km and 1 km, respectively, for eight water sources along the Yangtze's Nanjing section expect

Determination of Key Risk Supervision Sub-Regions
The x and y directional distances from sub-regions to source water intakes were determined based on the risk indexes and locations of risk sources.Appendix A Table A3 lists the specific results.In Appendix A Table A3, the x directional distance from sub-region to water intake would be positive if the sub-region was located upstream of the water source; 0 if the sub-region was located within the water source protection area; and negative if located downstream of the water source.The y directional distance from sub-region to water intake would be 0 if the sub-region and the water source were on the same shoreline; and positive if were on the opposite shorelines.Some sub-regions such as SD14 were located in the Yangtze tributaries, so their x directional distances to water intakes were regarded as the sum of distance from sub-region to tributary estuary and distance from the estuary to water intake.Meanwhile, the y directional distance to water intake was considered as the y directional distance from tributary estuary to water intake.In addition, part of the sub-regions did not pose threats to some water sources as a result of island barrier.In such case, the distance from sub-region to water intake was not considered, such as the distances from the sub-regions SD11-SD14 to the water source NJ01 in Appendix A Table A3.According to the Jiangsu Provincial Program for Partitioning County Level or Above Centralized Drinking Water Source Protection Areas [27], the distances from the source water intake to the secondary protection area upstream and downstream boundaries were 2 km and 1 km, respectively, for eight water sources along the Yangtze's Nanjing section expect

Determination of Key Risk Supervision Sub-Regions
The x and y directional distances from sub-regions to source water intakes were determined based on the risk indexes and locations of risk sources.Appendix A Table A3 lists the specific results.In Appendix A Table A3, the x directional distance from sub-region to water intake would be positive if the sub-region was located upstream of the water source; 0 if the sub-region was located within the water source protection area; and negative if located downstream of the water source.The y directional distance from sub-region to water intake would be 0 if the sub-region and the water source were on the same shoreline; and positive if were on the opposite shorelines.Some sub-regions such as SD14 were located in the Yangtze tributaries, so their x directional distances to water intakes were regarded as the sum of distance from sub-region to tributary estuary and distance from the estuary to water intake.Meanwhile, the y directional distance to water intake was considered as the y directional distance from tributary estuary to water intake.In addition, part of the sub-regions did not pose threats to some water sources as a result of island barrier.In such case, the distance from sub-region to water intake was not considered, such as the distances from the sub-regions SD11-SD14 to the water source NJ01 in Appendix A Table A3.According to the Jiangsu Provincial Program for Partitioning County Level or Above Centralized Drinking Water Source Protection Areas [27], the distances from the source water intake to the secondary protection area upstream and downstream boundaries were 2 km and 1 km, respectively, for eight water sources along the Yangtze's Nanjing section expect for the Jiajiang water source (NJ01), whose corresponding distances were 3.2 km and 3.0 km, respectively.Width of all waters in the secondary protection areas was 0.5 km.
SrTWS were calculated based on the SrRI and the distances from water intakes to secondary protection area boundaries in Table 5 and Appendix A Table A3, respectively.During calculation, the forward and reverse flow durations of water sources were considered according to the durations of ebb and flood in the Yangtze's Nanjing section, respectively.The calculation results are shown in Appendix A Table A4.Appendix A Table A4 also lists the total risk threats from single sub-region to all water sources (hereinafter referred to as the SD-TSrTWS), the total risk threats from all sub-regions to single water source (hereinafter referred to as the WS-TSrTWS), and SrRI.Since SD03 was located within the protection area of water source NJ07, it was a key supervision sub-region for NJ07 regardless of its SrRI.In Appendix A Table A4, the scale of SrTWS on NJ07 from SD03 was 4. Others, such as the scale of SrTWS on NJ10B from SD05 were also taken as 4. In addition, for water sources not affected by some sub-regions due to the barrier of islands in the river, their SrTWS was taken as 0.
Based on Appendix A Table A4, the key and secondary key supervision sub-regions corresponding to each water source were determined according to the sub-region supervision grading criteria established in Table 3.The results are shown in Table 6.Variations of SD-TSrTWS and SrRI for various sub-regions (Figure 8) were plotted in the upstream to downstream order of sub-regions.In the Figure 8, the upstream sub-regions were located to the left of the downstream sub-regions.As shown in Appendix A Table A4 and Figure 8, SD16's SD-TSrTWS was the maximum, 15.873.Although the SrRI of sub-regions SD26, SD21, SD19 and SD06 were all larger than the sub-region SD16, SD26, SD21 and SD19 were located downstream of the majority of water sources.Besides, SD21 and SD19 were located on the north and south watercourse of the Yangtze's Bagua Island, respectively, and did not produce risk threat to the water sources that were not on the same watercourse.As for sub-region SD06, although it was located upstream of most water sources, it was distant from most water sources compared to the sub-region SD16.Thus, among five sub-regions with SrRI greater than 5, SD16's SrTWS to the Nanjing riverside water sources was the greatest.Despite an up to 7.227 SrRI of sub-region SD19, it was only upstream of a water source NJ08 and did not produce risk threat tothe water sources NJ02 and NJ09B that were not on the same watercourse.Moreover, the water sources NJ04 and NJ10B on the same side of shoreline with it were located 20 km upwards of its upstream, so SD19 was only a secondary key supervision sub-region for water source NJ08; and was a non-key supervision sub-region for all other water sources.On the whole, the upstream sub-regions had greater risk threats to the Yangtze riverside water sources in Nanjing.One typical example was the sub-region SD06.Although its SrRI was lower than the downstream sub-regions SD26, SD21 and SD19, its SD-TSrTWS was 15.030 owing to its location upstream of most water sources, which was higher than the above-mentioned three sub-regions.Meanwhile, other upstream sub-regions such as SD01-SD05 and SD07-SD14 had less risk threats to water sources Variations of WS-TSrTWS for various water sources (Figure 9) were plotted in the upstream to downstream order.In the figure, upstream water sources were located to the left of the downstream water sources.As shown in Table A4 and Figure 9, downstream water sources were under greater risk threats from sub-regions than the upstream water sources.For example, WS-TSrTWS was the largest for the most downstream water source NJ08, 42.210.The water source NJ02 exhibited the second largest WS-TSrTWS value, 29.207.However, as it was located on the south watercourse of the Yangtze's Bagua Island, the sub-regions on the north watercourse posed no risk to it.Thus, its WS-TSrTWS was only slightly larger than two upstream water sources NJ01 and NJ04.As for NJ09B and NJ03 located separately on the north and south Bagua Island watercourses, they were also less affected by the risk sources in sub-regions than the two upstream water sources NJ01 and NJ04 since they were not on the same shoreline with most sub-regions apart from the above-mentioned reasons.

Suggestions
According to the judgment results of risk supervision level of risk source and sub-region, it is suggested that water source supervision departments should take the following measures to improve SrRI SD-TSRTWS Sub-regions SD-TSrTWS SrRI Variations of WS-TSrTWS for various water sources (Figure 9) were plotted in the upstream to downstream order.In the figure, upstream water sources were located to the left of the downstream water sources.As shown in Table A4 and Figure 9, downstream water sources were under greater risk threats from sub-regions than the upstream water sources.For example, WS-TSrTWS was the largest for the most downstream water source NJ08, 42.210.The water source NJ02 exhibited the second largest WS-TSrTWS value, 29.207.However, as it was located on the south watercourse of the Yangtze's Bagua Island, the sub-regions on the north watercourse posed no risk to it.Thus, its WS-TSrTWS was only slightly larger than two upstream water sources NJ01 and NJ04.As for NJ09B and NJ03 located separately on the north and south Bagua Island watercourses, they were also less affected by the risk sources in sub-regions than the two upstream water sources NJ01 and NJ04 since they were not on the same shoreline with most sub-regions apart from the above-mentioned reasons.Variations of WS-TSrTWS for various water sources (Figure 9) were plotted in the upstream to downstream order.In the figure, upstream water sources were located to the left of the downstream water sources.As shown in Table A4 and Figure 9, downstream water sources were under greater risk threats from sub-regions than the upstream water sources.For example, WS-TSrTWS was the largest for the most downstream water source NJ08, 42.210.The water source NJ02 exhibited the second largest WS-TSrTWS value, 29.207.However, as it was located on the south watercourse of the Yangtze's Bagua Island, the sub-regions on the north watercourse posed no risk to it.Thus, its WS-TSrTWS was only slightly larger than two upstream water sources NJ01 and NJ04.As for NJ09B and NJ03 located separately on the north and south Bagua Island watercourses, they were also less affected by the risk sources in sub-regions than the two upstream water sources NJ01 and NJ04 since they were not on the same shoreline with most sub-regions apart from the above-mentioned reasons.

Suggestions
According to the judgment results of risk supervision level of risk source and sub-region, it is suggested that water source supervision departments should take the following measures to improve SrRI SD-TSRTWS Sub-regions SD-TSrTWS SrRI

Suggestions
According to the judgment results of risk supervision level of risk source and sub-region, it is suggested that water source supervision departments should take the following measures to improve the anti-risk capability of risk source and sub-region as well as prevent the occurrence of water pollution accidents.
(1) Compared with sewage treatment plants, the water source supervision departments should strengthen the risk supervision for high-risk port sources listed in Table 4, and urge them to equip with pollution prevention devices and automatic water quality monitor so as to improve their sensitivity to liquid chemical and oil leakage accidents.Meanwhile, considering that the follow-up processing difficulty and pollution degree to water of chemical and oil leakage accidents are higher than those of other types of pollution accidents, the water source supervision departments should also enhance inspection frequency to high-risk ports, and urge responsible units for the corresponding ports to strengthen the daily safety management of the ports so as to eliminate the hidden risks of safety accident and prevent the pollution accidents caused by safety accidents.(2) Similarly, for the three high-risk enterprises (petrochemical and petroleum processing enterprises) listed in Table 4, the water source supervision departments also should urge them to install automatic water quality monitors at the wastewater discharge ports.Like the wastewater processing plants listed in Table 4, the water quality monitor data need to simultaneously transfer to online monitoring platform of supervision departments so that they can conduct real-time monitor over the water quality of wastewater discharge ports in industrial sources and avoid the occurrence of secretly discharging pollutants.(3) Water source supervision departments should strengthen daily inspection for key supervision sub-regions and secondary key supervision sub-regions surrounding the various water sources, especially the high-risk sources in sub-regions.Meanwhile, considering the concentrated distribution of risk sources in sub-regions, accidents like explosions and fire may cause continuous safety accidents and increase the occurrence rate of water pollution accidents.Therefore, water source supervision departments should also strengthen inspection for the safe production of risk sources in key and secondary key supervision sub-regions.

Limitations
Based on the high degree of industrial development and urbanization of areas along Nanjing of Yangzi River, the methodology proposed in this paper only considers the threats of point sources to the water sources.For water sources in some industrially underdeveloped areas (such as Northern Jiangsu) whose threats are primarily from non-point sources, this methodology is not appropriate.Hence, further studies are needed to determine the key risk sources and key risk supervision sub-regions around water sources by comprehensively considering the threats of point and non-point sources.
In addition, after a pollution accident at the risk sources, its threat to the water sources will be affected by the flow of water sources.Larger flow will reduce the quality substandard multiple of water sources resulting from pollution accidents.Accordingly, the threats imposed to water sources will also diminish.In this paper, water source flow has little influence on identifying the key risk sources and key risk supervision sub-regions since the flow rates of the eight water sources are basically the same.Nonetheless, if the water sources in a sub-region are not in the same river, flow difference between water sources will be large.From the overall regional perspective, the influence of water source flow should be considered in the identification of key risk sources and key risk supervision sub-regions.Therefore, it is necessary to consider the influence of water source flow in the follow-up study.

Conclusions
In this paper, a feasible methodological system for determining the key risk supervision areas is put forward for river-type water sources surrounded by multiple risks sources.The system comprises: (1) method of partitioning risk source concentrated sub-regions around water sources proposed based on Prim's minimum spanning tree algorithm [23] that considers the distribution of risk sources around water sources; (2) approach of determining sub-regional RDDI proposed based on the distance between adjacent risk sources on the minimum spanning tree of sub-regional risk sources and their respective RIRS; method of determining RSQI proposed by considering the number of risk sources in sub-regions; and the approach of determining SrRI reflecting the risk scale of sub-regions proposed based on the RDDI and RSQI combined with the RIRS within sub-regions; and (3) method of calculating SrTWS based on the positional relationship between sub-regions and water sources as well as SrRI, and the criteria for determining key supervision sub-regions.By applying this methodological system for determining water source perimeter key risk supervision areas along the Yangtze's Nanjing section, it is observed that the system can effectively identify the key and secondary key risk supervision sub-regions for each water source, and the high risk level risk sources within the entire region.Thus, it provides a basis for water source regulatory authorities to improve the efficiency of drinking water source supervision.

Figure 1 .
Figure 1.Distribution locations of water sources and risk sources in the study area.

Figure 2 .
Figure 2. Flow chart for determination of key risk supervision sub-regions around water sources.

Figure 1 .
Figure 1.Distribution locations of water sources and risk sources in the study area.

Figure 1 .
Figure 1.Distribution locations of water sources and risk sources in the study area.

Figure 2 .
Figure 2. Flow chart for determination of key risk supervision sub-regions around water sources.

Figure 2 .
Figure 2. Flow chart for determination of key risk supervision sub-regions around water sources.

Figure 3 .
Figure 3. Construction process of risk source minimum spanning tree.

Figure 3 .
Figure 3. Construction process of risk source minimum spanning tree.

Figure 6 .
Figure 6.Statistical analysis chart of risk indexes for high risk sub-regions.

Figure 7 .
Figure 7. Statistical analysis chart of risk indexes for medium risk sub-regions.

Figure 6 . 24 Figure 6 .
Figure 6.Statistical analysis chart of risk indexes for high risk sub-regions.

Figure 7 .
Figure 7. Statistical analysis chart of risk indexes for medium risk sub-regions.

Figure 7 .
Figure 7. Statistical analysis chart of risk indexes for medium risk sub-regions.

Figure 8 .
Figure 8. Variations of SD-TSrTWS and SrRI for various upstream and downstream sub-regions.

Figure 9 .
Figure 9. Variations of WS-TSrTWS for various upstream and downstream water sources.

Figure 8 .
Figure 8. Variations of SD-TSrTWS and SrRI for various upstream and downstream sub-regions.

Figure 9 .
Figure 9. Variations of WS-TSrTWS for various upstream and downstream water sources.

Figure 9 .
Figure 9. Variations of WS-TSrTWS for various upstream and downstream water sources.

Table 1 .
Construction results of risk source minimum spanning tree.

Table 4 .
RIRS calculation results for high risk sources.

Table 5 .
Sub-regional risk index calculation results.

Table 6 .
Key and secondary key supervision sub-regions corresponding to various water sources.

Table A3 .
Integrated distances from sub-regions to source water intakes.Unit: km.