Next Article in Journal
Synthetic Impacts of Internal Climate Variability and Anthropogenic Change on Future Meteorological Droughts over China
Next Article in Special Issue
Evaluating Glacier Volume Changes since the Little Ice Age Maximum and Consequences for Stream Flow by Integrating Models of Glacier Flow and Hydrology in the Cordillera Blanca, Peruvian Andes
Previous Article in Journal
Waters from the Djiboutian Afar: A Review of Strontium Isotopic Composition and a Comparison with Ethiopian Waters and Red Sea Brines
Previous Article in Special Issue
Automated Geospatial Models of Varying Complexities for Pine Forest Evapotranspiration Estimation with Advanced Data Mining
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Barrier-based Longitudinal Connectivity Index for Managing Urban Rivers

1
College of Resource Environment and Tourism, Capital Normal University, State Key Laboratory Incubation Base of Urban Environmental Processes and Digital Simulation, Capital Normal University, Beijing 100048, China
2
College of Environment and Planning, Henan University, Kaifeng 475004, Henan, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(11), 1701; https://doi.org/10.3390/w10111701
Submission received: 12 October 2018 / Revised: 14 November 2018 / Accepted: 16 November 2018 / Published: 21 November 2018
(This article belongs to the Special Issue Applications of Remote Sensing and GIS in Hydrology)

Abstract

:
A large variety of barriers can affect longitudinal connectivity, which leads to shipping blocking and even flood hazard. However, few existing methods can quantify physically the river channel connectivity from the barrier’s details perspective in a watershed. This paper establishes a new model of the River Channel Connectivity Index (RCCI) to quantify the unobstructed degree of river flow in river channels within geographic information system (GIS ) platforms based on the modified concept of time accessibility. A comprehensive classification system of barriers is setup before these barriers are identified by the remote sensing technology. The model is applied to Dashi Watershed in suburban Beijing, China. Results show that submersible bridges and sediment siltation are the main barriers in the watershed. RCCI values in the mountainous areas are generally higher than that of the plains. The assessment results verified by two historical flood events show that the RCCI can reveal where the river channel connectivity is impaired, how serious it is, and what the reason is for managers. Through scenarios’ results, the best restoration measure for each tributary is obtained from the perspective of reducing flood hazards. The new RCCI method not only has methodological significance, but also helps policymakers to enhance river flooding reduction and determine restoration priorities of the river channel.

1. Introduction

Natural watersheds are increasingly being affected by economic and social activities, which can introduce various types of barriers (e.g., dams and culverts) to river channels, reducing hydrologic connectivity, and leading to transport blockage, ecological habitat damage, and heightening flood hazard risks [1,2,3]. Good hydrologic connectivity is a key element for maintaining ecological integrity and it can be used by city managers to reduce flood risks, manage water resources, and restore fluvial ecological systems in practice.
Hydrologic connectivity aims at river ecological protection and mainly refers to water-mediated transferring of matter, energy, and organisms within or between elements of the hydrologic cycle [4,5]. A variety of approaches and metrics have been developed depending on the aim of the research and the target species in hydrological connectivity studies [6].
Hydrological connectivity can be interpreted in terms of static/structural and dynamic/functional connectivity [4,7]. Structural connectivity refers to the spatial patterns in the landscape. It is measureable and easy to understand. A series of metrics on structural connectivity is based on high-resolution digital elevation model (DEM) data. Several topography-based indices have been developed and applied in different contexts [8,9,10,11,12,13,14]. Landscape ecologists are familiar with using graph theory to describe and analyze river networks. Then, a series of assays of landscape connectivity indices are used to measure the structural connectivity [15,16].
In the three dimensions manifested by hydrological connectivity, longitudinal connectivity is defined as connections between upstream and downstream sections of a river network, as opposed to vertical of lateral connections [17,18]. Natural and anthropogenic barriers (e.g., buildings and dams) have impeded longitudinal connectivity severely, especially in urban and suburban watersheds [9,10]. For example, Segurado et al. [19] introduced dams into the analysis of structural connectivity of the Tagus river network in Portugal based on a landscape connectivity metric. When assessing the effect of barriers on aquatic biological communities, the Dendritic Connectivity Index (DCI) is widely used and modified to quantify longitudinal connectivity [1,3,20,21,22,23]. In previous studies, anthropogenic barriers include dams, culvert, and roads, and natural barriers include waterfalls and beaver dams [22,24,25,26]. Combined with field investigation, barriers are defined as anthropogenic or natural deposits that block river flow in the river channels in this paper. Although researches on hydrologic connectivity considered barriers in the calculation of a certain index, they are based on a specific object, such as water, biological species, sediment, or soil-moisture [17,27], and study the effect of barriers on connectivity of these specific objects. The sizes and attributes of barriers are highly variable [3]. Previous researchers may consider the difference of small culverts and large dams. However, culverts are different in size and dams are different in scales. It seems that there is a blank that puts the details’ information of a barrier in the calculation of longitudinal connectivity in the previous studies. Furthermore, the critical item of the river channel is ignored. As the container of runoff or discharge, a river channel directs the movement of river flow, organisms, and materials [28,29]. River channel connectivity as the basis of hydrological connectivity decreased by barriers should be paid more attention.
In this paper, a comprehensive classification system of barriers is developed for a suburban watershed. The blocking degrees of river flow through these barriers are weighted by a method combining subjective and objective weighting strategies. Based on the blocking degrees, a new concept, “river channel connectivity” (RCC), is proposed, which is defined as the unobstructed degree of the river flow in natural or artificial river channels. Moreover, a new assessment model, river channel connectivity index (RCCI), is established to quantify such an unobstructed degree. Then, the RCCI model is applied to Dashi Watershed in the suburb of Beijing, China. The results are validated by using two historical flood events from the potential applications of the RCCI model. Furthermore, scenarios analyses are introduced for scientific and quantitative fluvial environment planning from the perspective of reducing flood hazard.

2. Study Area and Field Data Collection

2.1. Study Area

Dashi Watershed is located in the Fangshan District of Beijing, China. It belongs to the Haihe River Basin. The 126 km long Dashi River drains an area of 1267 km2 with 16 catchments. It is divided into mountains and plains by 100 m contours (Figure 1). It is characterized by a semiarid and semi-humid monsoon climate with an average annual evaporation of 1500 mm and an average annual rainfall of 700 mm. Moreover, the watershed has an average annual runoff of 95.7 million m3, and its rainfalls concentrate in July and August, accounting for 80% of the annual total precipitation. Rainfall usually occurs in high intensity and short duration events, which often leads to severe regional flood events in summer [30].

2.2. Data Sources and Preprocessing

The 30-m shuttle radar topography mission (SRTM) DEM data are used to generate the watershed and catchment boundaries. The river channel and barriers are identified by visual interpretation using the SPOT6 remote sensing image data of 4 August 2015 and multi-period images from Google Earth. Before being used in the model, the above data underwent preprocessing, such as georeferencing, radiation calibration, and image fusion. Local references are used to generate the statistical data, and spatiotemporal information of relevant rainstorm disasters [31,32,33]. These documents are also used as the important basis for the classification and weight assignment of barriers.

3. Methodology

The various barriers are classified and weighted before they are identified from the river channel by the interpretation of the high-resolution image. The RCCI model is first presented in this study by modifying the accessibility method. Additionally, the passability of river flow is the key factor in the model. The RCCI model is applied to calculate the connectivity of different tributaries and river segments, and the results are validated by two historical flood events of 2000 and 2012 in the watershed. Finally, based on the scenarios’ analyses, guides for fluvial environment planning are obtained from the perspective of reducing flood hazard (Figure 2).

3.1. Classification and Weight Assignment of Barriers

All barriers could be classified into four types listed in Table 1 according to the local expertise and the related references [33] in Dashi Watershed. Different barrier types have different blocking weights on the river flow. The blocking weight represents the natural river flow passability, which is the key parameter of the proposed RCCI model. The weight is based on the barrier’s holding capacity of the river flow, which depends on the shape and the type of the barrier [1]. Table 1 presents the subjective blocking weights (in the column of p n i ) of all kind of barriers using the Delphi method [34]. The Delphi method was an anonymous process where ideas were assigned to the participants in the form of a questionnaire [35,36]. There were two questions in the questionnaire. One was whether the classification of the barriers was reasonable; the other was to grade the blocking weights of various barriers. The questionnaires were sent out to 20 experts by E-mail. Then, responses were collected and analyzed. Investigators achieved “group” consensus on the classification, the rank of barriers, and most of the blocking weights, in the second round. Finally, the average values of these different weights were calculated as the final blocking weights if a barrier had several different weights from different experts. Therefore, the anonymous nature of the Delphi process ensures that a single dominant group member does not inordinately influence the group’s outcome [36].

3.2. Identification of River Channels and Barriers

Figure 3 shows the technical flowchart on how to identify the four kinds of barriers (marked in blue) and the river channels. The water bodies are identified firstly by the threshold value (near infrared bands (B4)) and exponential methods (normalized difference water index, NDWI) according to the spectral characteristics before the river channel and reservoirs distinguished according to their shapes. Secondly, all the reservoirs and dams are identified based on the shape characteristics from the image extracted by the river channel and information in the references [31,32]. The submersible bridges and sediment siltation barriers are identified by visual interpretation using multi-period images from Google Earth since the submersible bridges are easily confused with other bridges in the image texture, and the sediment siltation has an irregular shape. Finally, the identification results are verified by the field investigation for ensuring the accuracy of the river channel and barriers (Table 2). Additionally, some other RCCI-related parameters, such as the number and widths of spillways, widths of gates, and heights of submersible bridges and rubber dams, are generated from the local documents and references [31,32] in the case that some key information about river depth and the actual height of barriers cannot be interpreted from remote sensing images (Figure 3).

3.3. Definition and Calculation of RCCI

As the quantitative indicator of river channel connectivity (RCC), RCCI is measured by using an accessibility method defined by the modified traditional accessibility method. Previous researchers defined accessibility as the intensity of the possibility of interaction [37,38]. Here, we define accessibility as the degree of difficulty for a certain volume of river flow to overcome various physical barriers from one particular location to the other one in the river channel. The quantitative assessment indicators of accessibility generally include time accessibility and distance accessibility [39]. The former is adopted since it can be altered by barriers, and the difficulty in the definition refers to the time parameter in this study. In the case of no barriers in a natural river channel, as shown in Figure 4a, the time accessibility required for river flow is the minimum from the upstream point, A, to the downstream point, B, therefore, the connectivity is the maximum. The channel is divided into n + 1 sections if there are n barriers [1], which decreases the river volume in unit time due to the reduced connectivity in the fragmented river channel. Therefore, we assume that the RCCI value equals to the ratio of the time accessibility of river flow in the river channel without any barriers to with some barriers in the condition of the same volume ( V ) of river flow running from one point to another point in a channel by starting from the Equation (1). The range of the RCCI value is [0, 1].
RCCI = t T
where t is the time without any barriers in a river channel, and RCCI equals to 1 in this situation; T is the time with barriers in the same river channel.
The Equation (2) in physics is as follows:
t = L v
where t refers to the time of the river flow from the upstream to the downstream point (s); L is the length of a river channel (m); v is the flow speed (m/s).
Additionally, Equations (3) and (4) can be obtained:
t = V Q
Q = v × S
where t is the time of the river flow going through a cross-section of a channel (s); V is the volume of the river flow (m3); Q refers to the runoff rate (m3/s); and S is the cross-sectional area of a channel (m2).
Therefore, Equation (5) can be referred to from the Equations (2)–(4):
t = L S t V
When there are n barriers with length, l i ( i = 1, 2, 3…n), in the channel, as is the case in Figure 4b, the components of T include the time of the river flow with barriers and without barriers. T is calculated by using Equation (6):
T = ( L 1 n l i ) × S t V + 1 n ( l i S t i V )
where t i is the time of the river flow going through the river channel segment, which is impaired by barrier i (s).
When there is a barrier i in the channel, the cross-sectional area for the river flow changes to S S i . S i refers to the cross-sectional area of barrier i, and it is determined by the maximum width and maximum height in the water of the barrier (m2). In addition, the river depth is set according to the normal flow level. Then, we can refer to Equations (7)–(9) from the established Equations (2)–(6):
t i = S S S i t
T = ( L 1 n l i ) × S t V + 1 n ( l i S S S S i t V )
T = ( L 1 n l i ) × S t V + 1 n ( l i S 1 p s i t V )
The ratio of (S S i ) to S is defined as the river flow passability of the shape ( p s i ) in the channel with barrier i.
Those factors affecting the passability of river flow are related not only to the objective geometry of a barrier, but also to its type, which is treated as the subjective one as introduced in Section 3.1 (Table 1). Both the shape and the nature of a barrier determine the passability of the river flow through a specific barrier ( p i ). The same weights are assigned for the two parameters to calculate p i :
p i = 0.5 × p s i + 0.5 × p n i
The improved T is:
T = ( 1 1 n x i ) L S t V + 1 n x i L S t V × 1 p i
where x i equals the ratio of l i to L .
The calculation of RCCI is done by Equation (12):
RCCI = 1 1 + 1 n ( x i p i x i )
RCCI can not only calculate a tributary connectivity (the connectivity of A to B in Figure 4b), but also characterize a segment connectivity of a specific tributary (the connectivity of segment i in Figure 4b; in which x i equals the ratio of l i to L i ). The RCCI value of each segment depends only on the segment length, barrier length, and the river flow passability of the barrier. However, the RCCI value of each tributary depends on the ratio of each barrier length to the total tributary length, as well as the passability of the river flow of the corresponding barrier.

4. Results

4.1. Spatial Distribution of Barriers in the Watershed

Table 3 and Figure 5 show the information of the barrier number and their spatial locations for planners. Results show that the number of the submersible bridges accounts for 50% of all barriers as the majority barrier in Dashi Watershed. Deposits in a river channel are the second largest type of barrier, which accounts for 44%. In the mountainous area, there are 158 barriers; the number of submersible bridges accounts for 59% as the majority, and the number of deposits in a river channel accounts for 33%, of which sediment siltation accounts for only 11%. In the plain area, there are 122 barriers; the majority is sediment siltation, representing 47%, followed by 37% of submersible bridges and 11% of illegal buildings. The sediment siltation is the major component of barriers in the plain area due to the slow flow speed and human activities (e.g., dams and levee construction preventing longitudinal and lateral connectivity). Cross-river bridges are basically viaducts in the plain due to the greater number of cities, thus the number of submersible bridges is less than that in the mountain area.

4.2. Assessment Results of RCCI in the Watershed

4.2.1. RCCI Assessment Results for River Segments

The internal RCC differences of one tributary are demonstrated by the analysis of river segments’ RCCI assessment results. Planners can see how a single barrier affects the connectivity of a river channel. There are 341 river segments divided by 280 barriers in Dashi Watershed, and the range of RCCI values is from 0.23 to 1.00, with an average value of 0.89 and a standard deviation of 0.15. The values are divided into three classes by two break points based on the RCCI calculation. The two points are calculated based on the ratio of time accessibility with a barrier to without (T/t). One break point is setup at the ratio of 1.02 referring to more than 2% extra time, while the other break point is setup at the 1.1 referring to more than 10% extra time. Accordingly, the RCCI is defined as 0.98–1.00 for high connectivity, 0.98–0.91 for medium connectivity, and 0–0.91 for low connectivity. The proportions of segment number at various levels of connectivity in each catchment are presented in Table 4 and Figure 5.
The assessment results shown by Figure 5 and Table 4 are as follows: (1) The seven catchments in the plain are all low connected, but there are some high connected segments with no barriers or few barriers of small size, such as submersible bridges and dams. The nine catchments are high connected in the mountain area except the lower part of the upstream Dashi Watershed. However, there are a few low connected segments because of reservoirs or illegal buildings; (2) the RCCI result of segments can reflect the various connectivity characteristics caused by a single barrier in different parts of a tributary. For example, there are 13 segments divided by 12 barriers in Mapaoquan River. In the 12 barriers, there are seven sediment siltations, four submersible bridges, and one reservoir. The lengths of sediment siltation range from 226 m to 880 m. The length ratio ( x i ) of sediment siltation to segment ranges from 0.27 to 1.00, with an average of 0.75. Correspondingly, the RCCI is low with values from 0.65 to 0.85, with an average value of 0.72. However, the average length of the four submersible bridges is 7.25 m, the average x i is 0.06, and the average RCCI is 0.97. It could be concluded that x i (the ratio of the barrier length to the segment) has a negative effect on segment RCCI values. Therefore, RCCI is affected by the barrier type. Three high connected segments are presented in two sites of its downstream area and one site near Xisha River. Thus, we concluded that connectivity is clearly affected by the type and shape of the barriers in the river channel.

4.2.2. RCCI Result for Tributaries

The range of the tributary RCCI value is from 0.74 to 1.00 in Dashi Watershed based on the calculation of the RCCI model. The RCCI values are classified according to the standard in Section 4.2.1. The results show that the number of low connectivity tributaries is 41%, which are mainly in the catchments located in the middle and south of the watershed. The number of medium connectivity tributaries is 11%, and the high connectivity tributaries accounts for 48%, which are mainly within catchments located in the northwest of the watershed (Figure 6).
The RCCI assessment results of each tributary are mainly related to the number and shape of the distributed barriers. The RCCI of Mangniu River is 0.74 as the smallest value due to the large number and various types of barriers. However, Beiyu River shows a high connectivity with no barriers. The RCCI of Liulin River in the mountainous area is much greater than that of Mapaoquan River. The number of barriers in the Mapaoquan River is equal to that of the Liulin River, however, 50% of the barriers in the Liulin River are submersible bridges of a small size while 58% of barriers in the Mapaoquan River are sediment siltation, which are more obstructive to the river channel connectivity. The RCCI of Wajing River and Dongsha River are both about 0.88. The shape attributes of the distributed barriers in the two rivers make the difference. In Wajing River and Dongsha River, the ratio of barrier lengths exceeding 100 m accounts for 28% and 60%, respectively. Furthermore, there is a 775 m sediment siltation as the longest barrier in Dongsha River even though it has fewer barriers than Wajing River.

4.3. Verification of RCCI Assessment Results

Flood hazards might happen due to barriers existing in the river channel. Results of the RCCI model are validated by two historical flood events of 2000 and 2012 in the watershed. The RCCI values are firstly mapped to match spatially to the town locations, with the 2000 flood events reported from the “Report of Flood Control Planning in Fangshan District” [33]. The flood events occurred in Da’anshan Village, Nanjiao Village, Shilou County, Fozizhuang Village, Hebei County, Chengguan Street, Liulihe District, and other towns on 4 July 2000, where the RCCI values of the segments are mostly small in the plain area (Figure 7). In addition, the RCCI average values of each town are counted by the spatial analysis tool in the Arcgis10.0 platform before they are classified according to the standard in Section 4.2.1. The Pearson correlation analysis between the average value and the binary variable whether the town is flooded is conducted by SPSS 19.0 (International Business Machines Corp: New York, NY, USA). The rules are set as follows: High connected is set as 0, medium connectivity is set as 1, and low connected is set as 2. As for the binary variable, it is set as 1 if the town is flooded, and else set as 0. Results show that there is a medium correlation between the flooded town and the RCCI classification (correlation coefficient, r p = 0.405, p < 0.05, n = 18). Therefore, it indicates there is a relationship between the low connectivity and the flooded towns.
The greatest flood hazard in the past 50 years happened due to the heavy rain in Fangshan District on 21 July 2012. Flooded tributaries included the Zhoukoudian River, Xiekuo River, Dongsha River, Mapaoquan River, and Mangniu River, where the RCCI values are small. The Pearson correlation between the RCCI classification of each tributary and the binary variable of whether the tributary is flooded is also carried out by SPSS 19.0. The rules are set the same with the above. Results show that there is a medium correlation between the flooded tributary and the RCCI classification (correlation coefficient, r p = 0.539, p < 0.01, n = 23). Therefore, it also indicates there is a relationship between the low connectivity and the flooded tributaries.
It can be concluded that there is a certain correlation between RCCI and the town/river flooded disaster. The smaller the RCCI value, the more obstructed the river channel, the larger possibility of overtopping, and the more vulnerable the town/river is to flood disasters. This also precisely verifies the rationality of the proposed method in reflecting the river channel connectivity to some extent.

4.4. Scenario Results

The flood events just act as a verification of the results of the proposed model. To some extent, flood disaster might happen when the barriers exist in the river channel because the river flow may pass from the main channel into the floodplain [9]. Through Section 4.2 and Section 4.3, low connected tributaries are in 41%, which may result in a flood risk. To guide fluvial environment planning from the perspective of reducing flood hazard, six scenarios, A-F, are designed to evaluate how the channel connectivity would change under different barrier removal scenarios. Scenario A assumes that all tributaries have no barriers in Dashi Watershed (RCCI of 1.00; baseline); in scenario B-F, submersible bridges, sediment siltation, illegal buildings, reservoirs, and dams are removed to recalculate the RCCI of tributaries, respectively. Using the detailed method in Section 3.3, the RCCI values of scenarios A-F are shown in Table 5. Compared with actual RCCI, an appropriate and quantitative planning scheme for each tributary is obtained. For example, when sediment siltation in Mapaoquan River, Xisha River, Zhoukoudian River, Xiekuo River, and Dongsha River are removed, the RCCI of them would be improved to 0.960, 0.930, 1.000, 0.998, and 0.980, respectively. Therefore, it can be concluded from the scenario-based analysis that removing sediment siltation provides the greatest decrease to flood risk in these rivers, and the most effective planning scheme for these rivers is dredging the river channel. When illegal buildings in Wajing River, Shuangquan River, Baishikou River, Nanjiao River Liulin River, Shibanfang River, and Da’anshan River are removed, the RCCI of them would be improved to 0.960, 1.000, 0.970, 0.949, 0.970, 0.996, and 1.000, respectively. Thus, to improve river connectivity and reduce potential flood risks, the most effective river restoration and planning strategy for these tributaries is to remove illegal buildings.

5. Discussion

5.1. The Application of RCCI in Other Rivers

The key to this work is that the proposed RCCI model is based on the time accessibility method concept, which is applied to the river channel for the first time. The time accessibility is fundamentally reasonable to be applied for assessing the river channel connectivity by considering the barrier issue. The concept of RCC and the definition of RCCI are simple to understand and they could be used in other rivers.
On the other hand, as the study case of this paper, the Dashi Watershed is a representative suburban watershed of northern China. This paper has classified and weighted these barriers in Table 1. It could be used in such small-scale suburban watersheds directly. Moreover, Table 2 in the manuscript is the interpretation sign, which can be used as the standards to extract barriers from remote sensing data in a certain river basin.

5.2. Verification of RCCI Assessment Results by Flood Data

The best method for the RCCI validation might be using the monitoring data of the river flow time in a watershed [40,41]. However, the flow speed data could not be obtained as the river flow in the semiarid areas, like Dashi Watershed, due to its frequent drying up. Therefore, the RCCI validation could be conducted from the potential application of RCCI, such as the relationship to flood disasters. The danger of flooding through levee collapse or overtopping occurred in Zhoukoudian River, Xiekuo River, and Mangniu River in the heavy rain on 21 July 2012 [33]. Barriers, such as sediment siltation in the river channel, decreased the cross-sectional area of a channel. River flow would overtop to the floodplain if there is a large amount of discharge, which might result in a flood hazard. Therefore, results of the RCCI are validated by two historical flood events. As seen from Figure 7, the RCCI in the flooded town indicates a low connected level, and it might be reasonable that the RCCI could explain the flood hazard to some extent. Nevertheless, most segments of the flooded town within a mountainous area are highly connected, and the RCCI is not effective on these places. The RCCI value represents the river channel connectivity directly and it is only one of multiple impacts linked with flood hazards. Some key factors, such as precipitation, hydrological regime, and geomorphology, might potentially cause flood hazards [42,43,44]. It can be concluded that the proposed RCCI has certain limitations in the application of flood hazard prediction. RCCI should be improved by adding other factors according to actual demands in specific applications.

5.3. Reliability Analysis Based on the Time Accessibility Method

The proposed RCCI model is established based on the time accessibility method. Although the flow speed variable is introduced in Equation (2), the river flow speed is affected by various variables, such as elevation, slope, and roughness of a river channel. There are two hypotheses in the proposed model: (1) The time of the river flow going through the cross-section of a river channel increases when the cross-sectional area decreases due to an existing barrier if a barrier is the key issue for the river flow to go across; (2) river flow speed is assumed to be the same with or without a barrier in a channel. The RCCI value equals to the ratio of the time of a river channel without any barriers to the same river channel with some barriers, hence the river flow speed variable would be eliminated from the equation due to the above two hypotheses. Therefore, the only time difference between a river channel with barriers and without barriers is the increased time due to the presence of barriers. The flow speed might change after the river flow passes through the existing barriers. A constant, K, can be set theoretically as the ratio of the flow speed before and after passing through a barrier no matter how the field flow speed changes, though the K value varies largely within the different rivers and different kinds of barriers. As the RCCI value is a ratio value, the K value is assumed to be 1 for the convenience of calculation in this paper. Therefore, RCCI of different rivers in the study area is compared and analyzed under the equal standard. A future work will focus on how to get the actual K value for each barrier by using the field survey method for the improvement of the current RCCI model.
The methodology proposed in this paper is an original metric of river channel connectivity. It measures a structural connectivity of the river channel from the longitudinal dimension. The river channel is set as the object; and river flow is introduced into the construction of the RCCI model. However, flow speed is eliminated in the calculation. The proposed model is a pure structural connectivity metric based on the structures of the river channel and barriers. In contrast, when considering barriers, previous hydrological longitudinal connectivity researchers prefer to set the river flow and aquatic biological communities as the study object [1,21,26]. The Dendritic Connectivity Index (DCI) is a popular assessment index when setting aquatic biological communities as the study object. When setting river flow as the study object, Grill et al. modified the DCI model by combing the HydroROUT model [21]. The hydrological model needs a discharge database, which limits the application in the semiarid areas, like Dashi Watershed. The data for the RCCI model is easily obtained and the assessment method based on time accessibility method is understandable. The proposed method in this paper is suitable in the semiarid areas and may be modified to be used in other areas.

5.4. The Weakness of the RCCI Model

It seems that the use of the Delphi method makes it difficult to apply the RCCI model to other rivers because of the expert knowledge on barrier impacts. In Section 5.1, we have discussed that the Dashi Watershed is a representative suburban watershed of northern China; Table 1 could be used in other rivers directly. However, factors affecting the unobstructed degree of a river channel are not limited to the four categories in this paper. The channel morphology, such as lithology, elevation, and slope, are all key factors to the channel connectivity [10]. The above factors are not included in the RCCI calculation as this paper is focused on analysis of the various physical barriers on RCC. Nonetheless, sediment siltation, submersible bridges, and other artificial blocking factors are common barriers in the suburb of northern China. Considering hydrologic features, such as hydrologic topography, gradient, and rainfall, could be taken as further study on RCCI.

6. Conclusions

Reduced river channel connectivity can result in wide impacts on ecology and hydrology along a river or within a watershed. The proposed RCCI model provides a simple and convenient method for river channel connectivity assessment without any hydrological database only from the perspective of barriers. Furthermore, it keeps wide potential applications for the future impact assessments linking with hydrology and ecology, which gives rise to more new assessment indices modified from the core idea of RCCI. In this paper, RCCI presents a new vision to address fluvial planning problems from the perspective of reducing flood hazards, and it can guide river regulation and determine restoration priorities of a river channel for the government. This paper applies the RCCI model to Dashi Watershed in suburban Beijing, China. Results show that submersible bridges and sediment siltation are the main barriers in the watershed, while the sediment siltation is more common in the plain area. The segment RCCI indicates the spatial variety of river connectivity along each river channel and reveals the corresponding barrier structure. The range of tributary RCCI is from 0.74 to 1.00 in Dashi Watershed, which decreases from the Beiyu River to the Mangniu River. Assessment results from both tributaries and segments indicate that the RCC in the northwest mountainous areas are generally better than that of the south middle plains. The validation of the RCCI model shows that the RCCI can reflect the connectivity of the river channel. Based on scenario-based evaluation, effective restoration measures for each tributary to improve its connectivity is obtained, which provides useful information to support the decision making of planners in general and flood managers. For example, the most effective planning scheme for Mapaoquan River, Xisha River, Zhoukoudian River, Xiekuo River, and Dongsha River is dredging the river channel. For Wajing River, Shuangquan River, Baishikou River, Nanjiao River Liulin River, Shibanfang River, and Da’anshan River, removing illegal buildings will improve their connectivity effectively, and reduce flood risks.

Author Contributions

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

Funding

This research is funded by National Key R&D Program of China (2017YFC0505903) and the National Natural Science Foundation of China (NSFC41671509).

Acknowledgments

Thanks go to National Key R&D Program of China (2017YFC0505903) and the National Natural Science Foundation of China (NSFC41671509), as well as to Bowen Qiao for the algorithm discussion. I also thank Jun Yang from Tsinghua University for his comments on the article and Hong Zhang for providing administrative boundaries, DEM and vector data of flooding victim sites.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cote, D.; Kehler, D.G.; Bourne, C.; Wiersma, Y.F. A new measure of longitudinal connectivity for stream networks. Landsc. Ecol. 2009, 24, 101–113. [Google Scholar] [CrossRef] [Green Version]
  2. Jackson, C.R.; Pringle, C.M. Ecological benefits of reduced hydrologic connectivity in intensively developed landscapes. BioScience 2010, 60, 37–46. [Google Scholar] [CrossRef]
  3. Buddendorf, W.B.; Malcolm, I.A.; Geris, J.; Wilkinson, M.E.; Soulsby, C. Metrics to assess how longitudinal channel network connectivity and in-stream Atlantic salmon habitats are impacted by hydropower regulation. Hydrol. Process. 2017, 31, 2132–2142. [Google Scholar] [CrossRef] [Green Version]
  4. Bracken, L.J.; Wainwright, J.; Ali, G.A.; Tetzlaff, D.; Smith, M.W.; Reaney, S.M.; Roy, A.G. Concepts of hydrological connectivity: Research approaches, pathways and future agendas. Earth Sci. Rev. 2013, 119, 17–34. [Google Scholar] [CrossRef] [Green Version]
  5. Pringle, C. What is hydrologic connectivity and why is it ecologically important? Hydrol. Process. 2003, 17, 2685–2689. [Google Scholar] [CrossRef]
  6. Nadeau, T.L.; Rains, M.C. Hydrological connectivity between headwater streams and downstream waters: How science can inform policy. J. Am. Water Resour. Assoc. 2007, 43, 118–133. [Google Scholar] [CrossRef]
  7. Hohl, A.; Václavík, T.; Meentemeyer, R.K. Go with the flow: Geospatial analytics to quantify hydrologic landscape connectivity for passively dispersed microorganisms. Int. J. Geogr. Inf. Sci. 2014, 28, 1626–1641. [Google Scholar] [CrossRef]
  8. Jacobson, R.B.; Janke, T.P.; Skold, J.J. Hydrologic and geomorphic considerations in restoration of river-floodplain connectivity in a highly altered river system, Lower Missouri River, USA. Wetl. Ecol. Manag. 2011, 19, 295–316. [Google Scholar] [CrossRef]
  9. Trigg, M.A.; Michaelides, K.; Neal, J.C.; Bates, P.D. Surface water connectivity dynamics of a large scale extreme flood. J. Hydrol. 2013, 505, 138–149. [Google Scholar] [CrossRef] [Green Version]
  10. Wohl, E.; Rathburn, S.; Chignell, S.; Garrett, K.; Laurel, D.; Livers, B.; Patton, A.; Records, R.; Richards, M.; Schook, D.M.; et al. Mapping longitudinal stream connectivity in the North St. Vrain Creek watershed of Colorado. Geomorphology 2016, 277, 171–181. [Google Scholar] [CrossRef]
  11. Lane, S.N.; Reaney, S.M.; Heathwaite, A.L. Representation of landscape hydrological connectivity using a topographically driven surface flow index. Water Resour. Res. 2009, 45, 2263–2289. [Google Scholar] [CrossRef]
  12. Borselli, L.; Cassi, P.; Torri, D. Prolegomena to sediment and flow connectivity in the landscape: A GIS and field numerical assessment. Catena 2008, 75, 268–277. [Google Scholar] [CrossRef]
  13. Shore, M.; Mechan, S.; Cushen, M.; Jordan, P.; Mellander, P.E.; Kelly-Quinn, M.; Melland, A. Extent and role of ditches in affecting hydrological connectivity in agricultural landscapes. In Proceedings of the EGU General Assembly Conference, Vienna, Austria, 22–27 April 2012; p. 4956. [Google Scholar]
  14. Cavalli, M.; Trevisani, S.; Comiti, F.; Marchi, L. Geomorphometric assessment of spatial sediment connectivity in small Alpine catchments. Geomorphology 2013, 188, 31–41. [Google Scholar] [CrossRef]
  15. Cook, E.A. Landscape structure indices for assessing urban ecological networks. Landsc. Urban Plan. 2002, 58, 269–280. [Google Scholar] [CrossRef]
  16. Pascual-Hortal, L.; Saura, S. Comparison and development of new graph-based connectivity indices: Towards the prioritization of habitat patches and corridors for conservation. Landsc. Ecol. 2006, 21, 959–967. [Google Scholar] [CrossRef]
  17. Thompson, C.J.; Fryirs, K.; Croke, J. The disconnected sediment conveyor belt: Patterns of longitudinal and lateral erosion and deposition during a catastrophic flood in the Lockyer Valley, South East Queensland, Australia. River Res. Appl. 2016, 32, 540–551. [Google Scholar] [CrossRef]
  18. Ward, J.V.; Stanford, J.A. The four-dimensional nature of lotic ecosystems. J. N. Am. Benthol. Soc. 1989, 8, 2–8. [Google Scholar] [CrossRef]
  19. Segurado, P.; Branco, P.; Ferreira, M.T. Prioritizing restoration of structural connectivity in rivers: A graph based approach. Landsc. Ecol. 2013, 28, 1231–1238. [Google Scholar] [CrossRef]
  20. Mahlum, S.; Kehler, D.; Cote, D.; Wiersma, Y.F.; Stanfield, L. Assessing the biological relevance of aquatic connectivity to stream fish communities. Can. J. Fish. Aquat. Sci. 2014, 71, 1852–1863. [Google Scholar] [CrossRef]
  21. Grill, G.; Lehner, B.; Lumsdon, A.E.; MacDonald, G.K.; Zarfl, C.; Liermann, C.R. An index-based framework for assessing patterns and trends in river fragmentation and flow regulation by global dams at multiple scales. Environ. Res. Lett. 2015, 10, 55–60. [Google Scholar] [CrossRef]
  22. Fuller, M.R.; Doyle, M.W.; Strayer, D.L. Causes and consequences of habitat fragmentation in river networks. Ann. N. Y. Acad. Sci. 2015, 1355, 31–51. [Google Scholar] [CrossRef] [PubMed]
  23. Bourne, C.M.; Kehler, D.G.; Wiersma, Y.F.; Cote, D. Barriers to fish passage and barriers to fish passage assessments: The impact of assessment methods and assumptions on barrier identification and quantification of watershed connectivity. Aquat. Ecol. 2011, 4, 389–403. [Google Scholar] [CrossRef] [Green Version]
  24. Favaro, C.; Moore, J.W. Fish assemblages and barriers in an urban stream network. Freshw. Sci. 2015, 34, 991–1005. [Google Scholar] [CrossRef]
  25. Cheng, S.T.; Herricks, E.E.; Tsai, W.P.; Chang, F.J. Assessing the natural and anthropogenic influences on basin-wide fish species richness. Sci. Total Environ. 2016, 572, 825–836. [Google Scholar] [CrossRef] [PubMed]
  26. Mahlum, S.; Cote, D.; Wiersma, Y.F.; Kehler, D.; Clarke, K.D. Evaluating the barrier assessment technique derived from FishXing software and the upstream movement of brook trout through road culverts. Trans. Am. Fish. Soc. 2014, 143, 39–48. [Google Scholar] [CrossRef]
  27. Kim, J.; Mohanty, B.P. A physically based hydrological connectivity algorithm for describing spatial patterns of soil moisture in the unsaturated zone. J. Geophys. Res. Atmos. 2017, 122, 2096–2114. [Google Scholar] [CrossRef]
  28. Forman, R.T.T. Land Mosaics: The Ecology of Landscapes and Regions; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  29. Sánchez-Montoya, M.M.; Moleón, M.; Sánchez-Zapata, J.A.; Tockner, K. Dry riverbeds: Corridors for terrestrial vertebrates. Ecosphere 2016, 7, 10–24. [Google Scholar] [CrossRef]
  30. Liu, T.G.; Li, F. Study on optimization of water resources in Dashi River. Beijing Water Resour. 2004, 6, 45–47. [Google Scholar]
  31. Water Resources Authority of Fangshan District. Flood and Drought Disasters in Fangshan District; China Water & Power Press: Beijing, China, 2003. [Google Scholar]
  32. Liu, Y.K. Dictionary of Water Affairs in Beijing; China Water & Power Press: Beijing, China, 2008. [Google Scholar]
  33. Beijing Institute of Water; Fangshan Water Authority. Report of Flood Control Planning in Fangshan District; China Water & Power Press: Beijing, China, 2015. [Google Scholar]
  34. Wang, D.Q.; Tang, L.H.; Gu, Z.C.; Zheng, X.; Jiang, W.K. Malaria Transmission Potential in the Three Gorges Reservoir of the Yangtze River, China. Biomed. Environ. Sci. 2013, 26, 54–62. [Google Scholar] [PubMed]
  35. Williams, P.L.; Webb, C. The Delphi technique: A methodological discussion. J. Adv. Nurs. 1994, 19, 180–186. [Google Scholar] [CrossRef] [PubMed]
  36. Graham, B.; Regehr, G.; Wright, J.G. Delphi as a method to establish consensus for diagnostic criteria. J. Clin. Epidemiol. 2003, 56, 1150–1156. [Google Scholar] [CrossRef]
  37. Hansen, W.G. How accessibility shapes land-use. J. Am. Inst. Plan. 1959, 25, 73–76. [Google Scholar] [CrossRef]
  38. Zhang, L.; Lu, Y.Q. Regional accessibility of land traffic network in the Yangtze River Delta. J. Geogr. Sci. 2007, 17, 351–364. [Google Scholar] [CrossRef]
  39. Jiang, B.; Chu, N.C.; Xiu, C.L. Comprehensive evaluation and comparative analysis of accessibility in the four vertical and four horizontal HSR networks in China. Acta Geogr. Sin. 2016, 71, 591–604. [Google Scholar]
  40. Rushlow, C.R.; Godsey, S.E. Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost, North Slope, Alaska, USA. Hydrol. Process. 2017, 31, 4092–4106. [Google Scholar] [CrossRef]
  41. Golosov, V.; Collins, A.L.; Tang, Q.; Zhang, X.; Zhou, P.; He, X.; Wen, A. Sediment transfer at different spatial and temporal scales in the Sichuan Hilly Basin, China: Synthesizing data from multiple approaches and preliminary interpretation in the context of climatic and anthropogenic drivers. Sci. Total Environ. 2017, 598, 319–329. [Google Scholar] [CrossRef] [PubMed]
  42. Peñuela, A.; Javaux, M.; Bielders, C.L. How do slope and surface roughness affect plot-scale overland flow connectivity? J. Hydrol. 2015, 528, 192–205. [Google Scholar] [CrossRef]
  43. Zhang, Z.T.; Gao, C.; Liu, Q.; Zhai, J.; Wang, Y.; Su, B.; Tian, H. Risk assessment on storm flood disasters of different return periods in Huaihe River basin. Geogr. Res. 2014, 33, 1361–1372. [Google Scholar]
  44. Słapińska, M.; Chormański, J.; Glińska-Lewczuk, K. Relation between inundation frequency and habitat Conditions of floodplain lakes—A case study of the lowland Biebrza River (ne Poland). Environ. Eng. Manag. J. 2016, 15, 1311–1321. [Google Scholar]
Figure 1. Location of study area and Dashi Watershed.
Figure 1. Location of study area and Dashi Watershed.
Water 10 01701 g001
Figure 2. The technical flowchart of the river channel connectivity index (RCCI ) model.
Figure 2. The technical flowchart of the river channel connectivity index (RCCI ) model.
Water 10 01701 g002
Figure 3. The technical flowchart for the identification of river channels and barriers.
Figure 3. The technical flowchart for the identification of river channels and barriers.
Water 10 01701 g003
Figure 4. Conceptualized diagram of different barriers and parameters in the RCCI model. (a) is a natural river channel with no barriers; and (b) is the same river channel with some barriers.
Figure 4. Conceptualized diagram of different barriers and parameters in the RCCI model. (a) is a natural river channel with no barriers; and (b) is the same river channel with some barriers.
Water 10 01701 g004
Figure 5. RCCI assessment results for segments.
Figure 5. RCCI assessment results for segments.
Water 10 01701 g005
Figure 6. RCCI assessment result for each tributary.
Figure 6. RCCI assessment result for each tributary.
Water 10 01701 g006
Figure 7. RCCI model validation by the flooded town data of 4 July 2000.
Figure 7. RCCI model validation by the flooded town data of 4 July 2000.
Water 10 01701 g007
Table 1. Classification, grade, and blocking weight assignment of barriers.
Table 1. Classification, grade, and blocking weight assignment of barriers.
TypesNamesDescription p n i
ReservoirsMedium-sized Storage capacity is more than 10 million m3.0.30
Small-(I)-sizedStorage capacity is of 1–10 million m3.0.40
Small-(II)-sizedStorage capacity is of 0.1–1 million m3.0.50
PondThe pond is a small water storage built in mountainous or hilly areas, and its storage capacity of local runoff is less than 100,000 m3.0.60
DamsSluicesThe grading standards are the same as the reservoir. The storage capacities of the two sluices are equivalent to the medium-sized reservoir in Dashi Watershed.0.30
Rubber DamAlso known as a rubber sluice, crest can overflow. The storage capacities of the two rubber dams are equivalent to the small-(II)-sized reservoir in Dashi Watershed.0.50
Submersible Bridges Submersible bridges are simple ordinary bridges, which are constructed across the river channel. When water rises slightly, river flow will go through above the bridge.0.60
Deposits in a river channelIllegal BuildingsIllegal buildings are built in the river channel partly or wholly, such as village houses or cemeteries, etc. 0.60
Sediment siltationThe sediment siltation in a river channel is produced naturally or man-made, such as sand mining activity.0.70
Table 2. The barriers in the field landscape and remote sensing imagery.
Table 2. The barriers in the field landscape and remote sensing imagery.
NamePhotoImage
Reservoir Water 10 01701 i001
Rubber Dam Water 10 01701 i002 Water 10 01701 i003
Submersible Bridge Water 10 01701 i004 Water 10 01701 i005
Illegal Buildings Water 10 01701 i006
Sediment siltation Water 10 01701 i007 Water 10 01701 i008
Table 3. Number and ratio of various types of barriers in Dashi Watershed.
Table 3. Number and ratio of various types of barriers in Dashi Watershed.
TypesThe Whole RegionMountainous AreaPlain Area
NumberRatioNumberRatioNumberRatio
Reservoirs155%117%43%
Dams41%21%22%
Submersible Bridges13850%9359%4537%
Deposits in a river channelSediment siltation7527%1811%5747%
Illegal Buildings4817%3422%1411%
Total280100%158100%122100%
Table 4. Percentage of river segments with different levels of connectivity in each catchment.
Table 4. Percentage of river segments with different levels of connectivity in each catchment.
Name of Sub-BasinSegmentsThe Percentage of Each Level Segments in NumberThe Length of Segment/mThe Percentage of Each Level Segments in LengthClassificationBasic Morphological Types
High/%Medium/%Low/%High/%Medium/%Low/%
Upper Part of the Downstream Dashi Watershed701486861701882LowPlain
Middle Part of the Downstream Dashi Watershed702872657123077LowPlain
Zhoukoudian Watershed584085235,90746153LowPlain
Dongsha Watershed650050632830070LowPlain
Lower Part of the Downstream Dashi Watershed933224524,808414316LowPlain
Xiekuo Watershed7341154451,763531037LowPlain
Dingjiawa Watershed1833283914,983324127LowPlain
Lower Part of the Upstream Dashi Watershed520404013,914145036LowMountain
Middle of Dashi Watershed2857113235,714661618HighMountain
Shijiaying Watershed2463162133,96972199HighMountain
Nanjiao Watershed2365181731,13374224HighMountain
Baishikou Watershed2370131722,4619118HighMountain
Shibanfang Watershed12503317674948475HighMountain
Upper Part of the Upstream Dashi Watershed358361169,9879028HighMountain
Da’anshan Watershed12848823,6309451HighMountain
Beiyu Watershed11000010,00110000HighMountain
Note: River channel connectivity level defined for each sub-basin is based on majority segments of a certain level; that is, the proportion of certain segments exceeded the average of 33%.
Table 5. RCCI values at different scenarios.
Table 5. RCCI values at different scenarios.
NameScenario AScenario BScenario CScenario DScenario EScenario FActual RCCI
Beiyu River1.0001.0001.0001.0001.0001.0001.000
North River1.0001.0001.0001.0001.0001.0001.000
Baikoumen River1.0001.0001.0001.0001.0001.0001.000
Shangshuiyu River1.0001.0001.0001.0001.0001.0001.000
Yanglin Rvier1.0001.0001.0001.0001.0001.0001.000
Qiulinpu River1.0001.0001.0001.0001.0001.0001.000
Shijiaying1 River1.0001.0001.0001.0001.0001.0001.000
Jinjitai River1.0000.9971.0000.9991.0001.0000.996
Yuzi River1.0001.0000.9900.9900.9900.9900.990
Zhongjiao River1.0000.9880.9900.9980.9900.9900.987
Beijiao River1.0001.0000.9900.9900.9900.9900.990
Simatai River1.0000.9900.9901.0000.9900.9900.990
Da’anshan River1.0000.9900.9801.0000.9800.9800.980
Shijiaying2 River1.0000.9890.9800.9980.9800.9800.980
Shibanfang River1.0000.9720.9700.9960.9700.9700.970
Mao’ershan River1.0000.9300.9300.9500.9800.9300.930
Liulin River1.0000.9210.9200.9700.9450.9200.920
Nanjiao River1.0000.9040.9020.9490.9450.9020.902
Baishikou River1.0000.9030.9300.9700.9000.9000.900
Dashi River1.0000.9000.9300.9200.9400.8940.890
Dongsha River1.0000.8810.9800.8900.8800.8800.880
Shuangquan River1.0000.8810.8801.0000.8800.8800.880
Wajing River1.0000.8780.9100.9600.8760.8800.876
Xiekuo River1.0000.8710.9980.8700.8700.8700.869
Zhoukoudian River1.0000.8601.0000.8600.8600.8600.858
Dingjiawa River1.0000.8540.9000.8520.9400.8520.852
Xisha River1.0000.7930.9300.8100.8100.7890.789
Mapaoquan River1.0000.7790.9600.7780.8000.7780.778
Mangniu River1.0000.7440.8100.7600.8700.7430.743

Share and Cite

MDPI and ACS Style

Li, H.; Zhou, D.; Hu, S.; Zhang, J.; Jiang, Y.; Zhang, Y. Barrier-based Longitudinal Connectivity Index for Managing Urban Rivers. Water 2018, 10, 1701. https://doi.org/10.3390/w10111701

AMA Style

Li H, Zhou D, Hu S, Zhang J, Jiang Y, Zhang Y. Barrier-based Longitudinal Connectivity Index for Managing Urban Rivers. Water. 2018; 10(11):1701. https://doi.org/10.3390/w10111701

Chicago/Turabian Style

Li, Heying, Demin Zhou, Shanshan Hu, Jianchen Zhang, Yuemei Jiang, and Yue Zhang. 2018. "Barrier-based Longitudinal Connectivity Index for Managing Urban Rivers" Water 10, no. 11: 1701. https://doi.org/10.3390/w10111701

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

Article Metrics

Back to TopTop