Abstract

The connectivity of the stream network plays an important role in water-mediated transport and river environments, which are threatened by the rugged development process in China. In this study, based on graph theory, a connectivity evaluation index system was built, which includes the Edge Connectivity, Edge-Node rate, Connectivity Reliability, and Edge Weight. The new evaluation standard and calculation method of each index is presented. The river-lake system of Fenhu industrial park in Jiangsu China is simplified to an Edge-Node graph and evaluated by the index system as a case study. The results indicate that the river-lake system of the research area has low Edge Connectivity, a high Edge-Node rate, and high reliability in the current connectivity level. In addition, the Edge Weight index of several channels does not satisfy the standard of the Basic Edge Weight. To solve the connectivity problems, specific project plans include broadening the unqualified channel and building canals linked with the low-connectivity lakes. The results show that, after the planning, the connectivity of the stream network in Fenhu industrial park will increase, and the connectivity evaluation index system is useful in the study area.

1. Introduction

The interconnection of a stream system is notably important, which ensures that the stream network delivers water, nutrients, and sediment through the basin and accelerates the self-cleaning ability of the system [1]. In recent years, with the rapid social and economic development in China, the rugged developing mode makes some towns or villages have serious water connectivity problems [2]. On one hand, the natural sedimentation, block, and drought retard the flow and cause the water exchange ability to decline. On the other hand, because of the pursuit of high utilization rate of land in the urbanization construction process [3], many rivers were filled or truncated to provide more construction land, which artificially intensifies the disconnected situation of the stream network and causes drainage difficulty, water quality deterioration, poor ecological environment, and many types of water problems [46]. Hence, studies about the interconnection of stream systems in China’s urbanization process are emerging.

In recent years, the amount of research regarding stream network connectivity from different perspectives has constantly increased. From the hydrology perspective, connectivity refers to the water-mediated transport of matter, energy, and organisms within or between elements of the hydrologic cycle [7, 8]; from the ecology perspective, river-to-lake connectivity plays an important role in water renewal and aquatic habitat diversity [9, 10]; from the dynamics perspective, connectivity is important to the stream flow and runoff response [11]. To improve the hydraulic relationship between rivers and lakes, many engineering and nonengineering measures such as connecting project and hydrologic dispatch are used in China.

Graph theory is a theoretical model that uses an abstract method to represent the mutual relations among the study objects. The graphs in graph theory do not have characteristics such as size, shape, or mass except objects (represented by nodes) and their relations (represented by edges). The Edge-Node structure in the graph theory system is highly similar to the actual stream network [12]. The rivers in an actual stream network can be simplified to edges, and the lakes and junctions (points where two channels combine into one) can be simplified to nodes [13, 14]. According to the structure characteristic of a stream network and the graph theory method, the stream network can be generalized mathematics [15]. A simplified network can be evaluated by some analysis methods of graph theory, where the connectivity of the stream network can be quantitatively analyzed. Many previous studies examined the stream network analysis using graph theory [16]. Zhao et al. evaluated the connectivity of a channel-beach system using graph theory and provided the Edge Connectivity calculation process [17]. Xu et al. combined graph theory and flow resistance theory, defined a calculation of the Edge Weight, and compared the connectivity before and after channel dredging in Taihu Basin in China [18]; Wang and Xu et al. analyzed the river system connectivity of Wuxi’s central urban area based on graph theory and found that the river system had a decreasing trend [19]; Yang extracted the river network in Jiaodong, China, and evaluated its connectivity based on graph theory [20]. Tejedor et al. presented a framework based on the spectral graph theory and constructed vulnerability maps in the river deltas using the framework [21].

However, in previous studies, most researchers used single indicators to analyze the connectivity of stream networks, which introduced many problems caused by the limitations of individual indices. For example, the Edge Connectivity index can reflect the interrelationship of the nodes but cannot reflect the size of the rivers, and the Edge Weight index is the opposite. Nevertheless, different indicators of graph theory are rarely combined to make a comprehensive analysis of stream systems. In addition, little attention has been paid to the connectivity evaluation for river-lake systems (stream networks with numerous lakes) and the specific project guidance after the evaluation.

The objectives of this study are to (1) select a typical stream network with numerous lakes and abstract it based on graph theory, (2) establish a multidimension evaluation method to estimate the interconnectivity level of the river-lake system, and (3) use specific projects to improve the connectivity of the study area. First, we select the river-lake system of Fenhu industrial park in China as the study area and simplify it to an Edge-Node graph. Second, we build a connectivity evaluation index system that includes the Edge Connectivity, Edge-Node rate, Connectivity Reliability, and Edge Weights, and we explain the meaning and evaluation criterion of each index. Finally, we evaluate the connectivity of the river-lake system in Fenhu using the evaluation system and provide specific project planning to solve the connectivity problems.

2. Study Area and Methods

2.1. Study Area

We use Fenhu industrial park in Jiangsu province, East China, as the research area (Figure 1). Fenhu industrial park belongs to Taihu Basin with an area of 258 km2 and a population of 250 thousand. It has subtropical monsoon climate with plenty of rainfall and four distinct seasons [22]. The area is located on the flood plain with a ground elevation of 2.20-5.20 m and mostly belongs to polder areas. Because of the natural factors, the water level, flow direction, and velocity of the stream often vary.

The study area has a typical stream network structure in flood plains in China with flat, shallow, and narrow channels, a disordered river-lake system, deficient connectivity between rivers and lakes, serious sedimentation, and low draining capability. In addition, because the area is located in the downstream location of Taihu Lake (China’s second largest lake), it is easily threatened by flood disasters. Thus, Fenhu industrial park consists of several polder areas to prevent the inner polder areas from being threatened by flooding. The stream system of each polder area is separated by the dike, sluice, and dam, which make the inner polder area independent from the outside. This setting also introduces many problems such as insufficient flood control and drainage ability, standard-exceeding of water quality, and ecosystem destruction. Because of the relative independence of the inner polder area, in this study, we ignore the inner area and select the outer polder area as the research object (Figure 2).

2.2. River-Lake System Simplification

In this study, we focus on the connectivity of the stream network outside the polder area; we only identify and extract the outer rivers and lakes on the map using labels. As shown in Figure 3(a), following the principle of graph theory, we use numbers to represent lakes and junctions and circled numbers to represent rivers.

After identifying the research objects, we simplify the lakes and river junctions to nodes and the rivers to edges. Because the flow direction in the study area often changes with the climate and rainfall, we define the rivers as undirected edges [23]. As shown in Figure 3(b), there is a simplified, undirected stream network graph that reflects the relationship between rivers and lakes. The research below is based on the simplified stream network graph.

2.3. Connectivity Evaluation Index System

Many indices in graph theory can evaluate the connectivity, for example, the Edge Connectivity, Edge-Node rate, Actual Ring Formation rate, Network Connection degree, and Edge Weights. Currently, there is no unified and complete evaluation index system to perfectly and synthetically evaluate the connectivity of a stream network. In this paper, based on the principles of comprehensiveness, diversity, and operability, four indices (Edge Connectivity, Edge-Node rate, Connectivity Reliability, and Edge Weights) are set as a connectivity evaluation index system (Table 1). In the system, the indices of Edge Connectivity, Edge-Node rate, and Edge Weights are generally used in the graph theory; the Connectivity Reliability is created in this paper for the first time. Through this connectivity evaluation index system, the connectivity of the stream network can be better evaluated. We will introduce the meaning, calculation, and evaluation standard of each index in this system below.

2.3.1. Edge Connectivity

The Edge Connectivity is equal to the minimum number of edges, which must be removed to make an unconnected graph, with notation λ. The Edge Connectivity of an unconnected graph is zero. Graphs with high Edge Connectivity must remove more edges to become an unconnected graph; they have strong connectivity and are difficult to break. In an actual drainage basin, the poor connectivity of a small area always becomes the source of weak flood control ability, deteriorative water quality, and poor ecosystem of the entire basin. The Edge Connectivity index can precisely reflect the extent of the weakest link area in the stream network and use it to represent the connectivity of the entire system. The Edge Connectivity is the determinant index in an evaluation system, by which we classified the networks.

The unconnected stream networks contain isolated lakes, ponds, or independent small subnetworks (e.g., node n1 in Figure 4(a)). The connectivity of these parts is notably poor, and they almost cannot have water exchange with one another. Thus, an unconnected stream network has a notably poor interconnectivity, which is the unacceptable drainage structure form.

To the stream network, where the Edge Connectivity is 1 (Figure 4(b)), the connectivity of the network remains poor. The lakes whose degree (number of edges that link with the node) is 1 cannot smoothly exchange water with other lakes (e.g., node n1 in Figure 4(b)). It may cause water quality deterioration and allow the deteriorating water to be transmitted to the entire stream network. Thus, a system with an Edge Connectivity of 1 is also unacceptable.

To the stream network, when the Edge Connectivity is 2 (Figure 4(c)), every lake is at least linked with two other lakes as part of a closed loop. The in-out water-delivering channel of the lakes can be build. The pollutant in the stream network will be self-cleaned by a cycle in the loop. Thus, the system with an Edge Connectivity of 2 is acceptable.

To the stream network with Edge Connectivity at 3 or more, the connect abilities between lakes are generally more perfect. The water exchange is smooth, so it is the ideal stream network structure.

2.3.2. Edge-Node Rate

The Edge-Node rate is obtained by dividing the number of rivers by the number of lakes as follows:where is the Edge-Node rate, is the number of rivers, and is the number of lakes or junctions. The Edge-Node rate can reflect the relative density of rivers. This important index can evaluate the connectivity of the entire system. The system with a high Edge-Node rate commonly has a developed stream network, and each lake occupies more rivers.

For the connected network, an Edge-Node rate less than 1 indicates that the stream network is too thin, each lake averagely occupies less than 1 river, and there is no ring structure to circulate water in the network. Hence, it is an unacceptable Edge-Node rate number. When the Edge-Node rate is equal to 1, there is one ring structure in the system (e.g., Figure 4(b)). When the Edge-Node rate is larger than 1, there are two or more ring structures in the system, and the structure of the stream network is more complex (e.g., Figure 4(c)). Thus, the Edge-Node rate should be larger than one to satisfy the connectivity requirement.

2.3.3. Connectivity Reliability

The common connectivity index sometimes cannot reflect the different structures of two similar-size systems, so the Connectivity Reliability index has been created.

The Connectivity Reliability is the possibility of the Edge Connectivity remaining unchanged if one edge randomly disconnects; for a network whose Edge Connectivity is (formula (2)), we havewhere is the Connectivity Reliability of the network with an Edge Connectivity of , is the number of edges that can be removed while keeping the Edge Connectivity invariable, and is the number of all edges in the system.

The Connectivity Reliability has a notably adaptive practical significance: it can simulate the possibility of connectivity reducing by river blocking because of perennial sedimentation or geological disasters. At a certain Edge Connectivity, the stream network structure will become more reliable when the Connectivity Reliability index approaches 1 and more unreliable when it approaches 0. The systems with high Connectivity Reliabilities have a strong ability to resist river breaking; that is, even if a river breaks, the impact will be small. When the Connectivity Reliability is equal to 1, the Edge Connectivity will increase by a level. The Connectivity Reliability cannot be compared when the systems have two different Edge Connectivity levels.

According to statistical analysis and experience, the relationship between the Connectivity Reliability index and the actual stream network is built (Table 2).

2.3.4. Edge Weights

In graph theory, many river characteristics such as the length, width, and depth have been neglected. However, in actual stream networks, some characteristics have obvious effects on the interconnectivity. Thus, the Edge Weight index is developed and can be used to simulate the actual river forms. Many characteristics can affect the river discharge capacity, such as the flow section area, flow velocity, geometrical shape, length, slope, roughness, and evaporation. However, because the research area is flat, the slope tends to be zero, and the flow velocity is notably slow; other characteristics are subordinate and difficult to precisely measure in actual situations. Thus, the flow section of rivers is set as the Edge Weights, which is denoted by . This index has the largest effect on the river discharge capacity.

Rivers with high Edge Weights have large flow section and high discharge capacity to improve the connectivity of the stream network. There are different standards of Edge Weights in different research areas and situations. In this paper, for the lakes concentration areas south of Jiangsu, an evaluation criterion of Edge Weights is created, respectively, from the view of lakes and rivers.

(1) Basic Total River Flow Section to Lakes. The sum of the river flow section areas that are linked with lake is called the Total River Flow Section of lake , which is denoted by . Lakes with larger have greater water exchange capacity and self-cleaning ability. A lake with more water commonly requires larger . Thus, to analyze whether satisfies the lake connectivity requirement, the Basic Total River Flow Section is created and denoted by . Its relationship with the water volumes is given in the following formula:where is the Basic Total River Flow Section linked with lake n (m2); is the Total River Flow Section Coefficient (m2/m3); is the water volume of lake n (m3), calculated by the ordinary water level. is the minimum that the lakes require. When is greater than , the basic connectivity of the lake satisfies the demand. When is less than , the system cannot satisfy the basic connectivity demand and must be improved. Coefficient must be adjusted according to the specific situation of the study area.

(2) Basic Edge Weight. The Basic Edge Weight is the minimum flow section area that the rivers must have to satisfy the connectivity demand and is denoted by . Lakes connect with other lakes through rivers or channels; because of different water transmission quantities, the connection between large lakes commonly requires a large river flow section. The river flow section and size of lakes linked with the river have a strong positive correlation. Thus, their relationship can be built to find the Basic Edge Weights as follows:where are the numbers of lakes connected by river e; is the Basic Edge Weight of river e calculated from lake a (m2); is the Basic Edge Weight of river e calculated from lake b (m2); is the Basic Total River Flow Section of Lake a (m2); is the Basic Total River Flow Section of lake b (m2); is the volume of lake a (m3); is the volume of lake b (m3); is the total volume of all lakes connected with lake a (m3); is the total volume of all lakes connected with lake b (m3); is the Basic Edge Weight of river e (m2). The formula is calculated, respectively, from two lakes linked by river e. Firstly, We assign and to river e to obtain and ; then we take the maximum of those two as (Figure 5(a)).

When the realistic flow section of a river is larger than the Basic Edge Weight, we consider that the river satisfies the demand of basic connectivity. If every River Flow Section in the system is larger than the Basic Edge Weight, the entire stream system is considered to be smoothly running.

A special drainage structure will sometimes affect the accuracy of the calculation. For example, when two large lakes are connected by a notably small lake, the calculation results may be small and do not affect the actual situation (Figure 5(b)). Thus, we must make some adjustments according to the actual condition in the calculation, such as ignoring the small lakes and directly selecting two large lakes as the research lakes.

2.4. Research Procedure

The research method is summarized in three steps. (1) According to the actual river-lake relationship, abstract the stream network to a nodes-edges graph based on graph theory; (2) bring characteristic data into the connectivity evaluation index system to estimate the connectivity level of the river-lake system; and (3) develop connectivity planning including specific projects to improve the connectivity of the study area (Figure 6).

3. Results and Discussion

3.1. Connectivity Evaluation of Fenhu Industrial Park
3.1.1. Edge Connectivity

In the simplified stream network graph (Figure 3(b)), every node links to others and makes the graph a connected graph. If we remove one specific edge (edge 1 or edge 6), node 1 or node 6 will be separated from the entire network, which makes the system disconnected. Thus, the Edge Connectivity of the current network is 1.

3.1.2. Node-Edge Rate

There are 15 nodes and 19 edges in the graph. According to formula (1), the Node-Edge Rate is 1.27, which indicates that the stream network construction is complex and has more than two cyclic structures.

3.1.3. Connectivity Reliability

After we remove edge 1 or edge 6, the Edge Connectivity will be reduced to 0, and other rivers will not affect the Edge Connectivity if any one of them is removed. According to formula (2), the Connectivity Reliability is 0.89, which is in the range of [0.85, 1.00]; that is, the stream network is notably reliable in the current connectivity level (Table 2).

3.1.4. Edge Weight

The Edge Weight of each river is its river section area, as listed in Table 4. To avoid the influence of river and lake sizes varying and ensure the accuracy of the calculation, all the data we used are up to date. To determine whether the flow sections satisfy the evaluation standard, we should calculate and .

(1) Calculation of the Basic Total River Flow Section of Lakes . The Total River Flow Section of lake n () is calculated by adding all connected river sections of lake n. According to the practical situation in the study area and research experience, the Total River Flow Section Coefficient ω is defined as 0.00006 m2/m3, which implies that every 100000 cubic meters of water need 6 square meters of flow section area to ensure that the lake smoothly runs. Using formula (3), we can obtain of each lake, from which we subtract to obtain the Shortage of Total River Flow Section (Table 3).

When is positive, the lake does not have sufficient drainage channel area, and we must increase to . When MΔ is negative, need not be changed. According to Table 3, lakes 1, 2, 3, 6, 7, and 12 should adjust the flow section around them to satisfy the basic requirement of .

(2) Calculation of Basic Edge Weight . Based on formula (4), we can find two sets of Basic Edge Weight and select the larger one as , which subtract the current river section to obtain the Shortage of River Section (), as shown in Table 4.

When we calculate the Basic Edge Weights of rivers 8, 9, and 10, we select node 7 (Sanbai Lake) and node 9 (Yuan Lake) as lake a and lake b and ignore the small lake 8 (Maiyan) to ensure a correct and reasonable result. When is positive, the River Flow Section is not sufficiently large and we must increase to . When is negative, need not be changed. According to Table 4, rivers 1, 2, 3, 6, 7, 14, 17, and 19 have different insufficient situations in the river section area.

The above calculation and analysis show that the connectivity condition of the river-lake network in Fenhu industrial park is not sufficient. The problems of low Edge Connectivity, narrow channel, and insufficient flow section threaten the local water security and water environment. Therefore, the transformation plan targeted on the river-lake system in the study area to enhance the current connectivity situation is urgently required.

3.2. Transform Plan of the River-Lake System in Fenhu Industrial Park

Based on the current connectivity evaluation in Section 3.1, a connectivity improvement plan is provided with two phases. The first phase aims to make the stream network essentially unimpeded under the limited manpower, material resources, and financial resources. The second phase aims to solve most connected problems to make the river-lake network more unimpeded and smooth. The connectivity condition in the study area will be improved after the two phases.

3.2.1. First-Phase Planning

In Fenhu industrial park, the Edge Connectivity of the river-lake network is 1, which is unacceptable. Because the insufficient connection of nodes 1 and 6 may cause a part of stream network to decrease in water exchange efficacy and deterioration of water quality, which will affect the regional water environment, we must increase the Edge Connectivity. In addition, the Connectivity Reliability is sufficiently large, so a few connection projects can improve the Edge Connectivity, and the plan is feasible.

More specifically, we should focus on linking Xiaoyue Lake (node 1) and Yangwen Lake (node 6) to the entire stream network to increase the Edge Connectivity. According to the design philosophy of low work quantity and respecting the natural river morphology, we attempt to minimize digging new rivers, demolition, and immigration in the project process. Based on the local situation, we connect the water channels between Xiaoyue Lake (node 1) and Mugua Lake (node 15) and between Yangwen Lake (node 6) and Sanbai Lake (node 7). The new water channels are constituted by the surrounding unconnected rivers and lakes in the polder. Figure 7(a) shows the simplified graph after the connection; the new channels are named edge 20 and edge 21 as shown with dotted lines. After the primary improvement, the Edge Connectivity of the system has been increased to 2, which is acceptable.

3.2.2. Second-Phase Planning

After the first phase, the stream network is relatively unobstructed but can be further improved. Because of the difficulty of adding new rivers to increase the Edge Connectivity to 3 and because 2 is an acceptable Edge Connectivity value, we decide to improve the connectivity using other indices.

(1) Improve the Connectivity by Increasing the Connectivity Reliability. After the first phase, the Connectivity Reliability becomes 0.33, which belongs to the unreliable scale and should be improved (Table 2). Following the principle of best comprehensive benefits, we select the edge that can increase the Connectivity Reliability as much as possible and minimize the damage on the local original channel. Thus, we link Bengke Lake (node 3) with Fenhu Lake (node 14) and Hanlang Lake (node 12) with Taipu River (new node 16). Both projects use the existing rivers and lakes in the polder to avoid the disturbance. The simplified graph is shown in Figure 7(b), and the new channels are named edge 22 and edge 24 as shown with dotted lines.

(2) Improve the Connectivity by Increasing the Edge Weight. According to the analysis of Table 3, of some lakes are positive, which implies that the rivers around these lakes must be widened to satisfy the requirement. According to the analysis of Table 4, of some rivers are positive, which implies that we must broaden these rivers to satisfy the requirement.

The narrow rivers mainly appear around Shaobo Lake (node 2) and Sanbai Lake (node 7). Based on the calculation, we must desilt the rivers that link Xiaoyue Lake (node 1), Bengke Lake (node 3), and Mugua Lake (node 15) with Shaobo Lake (node 2) and broaden the rivers that link Yuanhe Lake (node 4) and Fenhu Lake (node 14) with Sanbai Lake (node 7). In addition, the flow section between Yemaowei Lake (node 13) and Hanlang Lake (node 12) is not sufficiently large and must be desilted. Thus, 6 rivers (river 1, river 2, river 7, river 14, river 17, and river 19) must be desilted or broadened. Because the new rivers 21 and 24 reduce the drainage stress of rivers 3 and 6, after the connecting projects, those two narrow rivers can satisfy the standard. After several specific project measures, the flow section of those rivers will be as follows: river 1: 16.52 m2; river 2: 30.24 m2; river 7: 116.34 m2; river 14: 39.05 m2; river 17: 138.18 m2; river 19: 27.18 m2.

3.3. Discussion

After two phases, the connectivity indices of the study area are changed. The Edge Connectivity is changed to 2, which implies that every lake is at least linked with two other lakes to compose a closed loop. The in-out water-delivering channel of each lake can be built, and the pollutant in the stream network will be self-cleaned by the cycle in the loop; there are 24 edges and 16 nodes in the graph, so the Edge-Node rate is 1.5; that is, the structure becomes more complex and more conducive to water exchange. Among all 24 edges, if we remove any edge in a group of 12 edges (edges 2, 3, 7, 8, 9, 16, 17, 18, 19, 22, 23, and 24), the Edge Connectivity will not change. Hence, in the situation where Edge Connectivity is 2, the index of Connectivity Reliability P2 is 0.5, which is in the reliable range. Two new edges (edges 20 and 21) are added to the Edge Weight index. The calculation results show that both and are negative after the projects, which indicates that all rivers satisfy the demand of minimum flow section area. All projects in planning are shown in Figure 8, including 4 connecting projects and 6 widening projects. All indices satisfy the standard requirements according to the connectivity evaluation index system. The connectivity between rivers and lakes is unobstructed, and the water security and water environment goals can be basically achieved.

The results show that the connectivity planning is effective. However, if the size of the river-lake system changes due to erosion and so forth in the future, we should use the evaluation system and the new data to reevaluate its connectivity. We still need to further validate the effectiveness of the evaluation system by investigating the actual connectivity of study area after the completion of the specific projects. In addition, more statistical analysis is needed to cross-validate the evaluation in the future to verify the universal applicability.

4. Conclusion

In this paper, based on the graph theory, we have built a connectivity evaluation index system of the river-lake system and introduced the calculation method of each index. The standard of a smooth river-lake network is as follows: the Edge Connectivity is larger than 2, the Edge-Node rate is larger than 1, the Connectivity Reliability is larger than 0.5, and the Edge Weights are larger than the Basic Edge Weight. According to different evaluation results, we can provide different planning schemes. For the stream networks with Edge Connectivity below 2, we add new channel to the systems to enhance the Edge Connectivity. For the river networks whose Edge Connectivity are larger than 2, we can enhance the Edge-Node rate, Connectivity Reliability, and Edge Weights, which is an effective method to improve the connection level.

We used the river-lake network of Fenhu industrial park in Suzhou China as an example and analyzed the connection situation using the established evaluation methodology. The result indicates that the river-lake network in Fenhu industrial park has several connectivity problems: low Edge Connectivity, narrow flow section, lack of Total River Flow Section of some lakes, and so forth. Hence, specific planning advice was provided, which includes designing new channels to link the isolated or low-node-degree lakes (such as Xiaoyue Lake and Yangwen Lake) and broadening the present narrow rivers (including 6 specific rivers). According to the calculation, after two phases of connecting plan are implemented, all connectivity indices will satisfy the standard requirements, and the connectivity of the river-lake system in Fenhu industrial park will be efficiently promoted.

The evaluation index system and planning method are notably useful in our study area. However, a limitation is worth noting. The Total River Flow Section Coefficient ω reflects the relationship between the lake volume and the Total River Flow Section of the lake; in our case study, it is defined by the interview of the front-line governor of Fenhu and our research experience. This index can be accurately calculated in future research.

Data Availability

The river data of Fenhu river-lake system were provided by Wujiang district water conservancy bureaus, which included the width, average surface elevation, bed elevation, and side slope ratio of each river. The flow section of each river was calculated from the basic river data. The lake volume data were obtained from the planning document of Wujiang Lakes Protection Planning and Basic Statistics of Lake over 50 mu. Requests for access to these data should be made to Chenguang Xiao ([email protected]).

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This study has been supported by the National Key Research and Development Program of China (2017YFC040320502), a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and Australian Research Council for Discovery Project (Grant DP170104138).