1. Introduction
Rivers exert many kinds of ecosystem service functions and thus are receiving increasing attention. Rivers’ ecosystem service functions include water supply and related services, such as transportation and hydroelectric generation, ecological supporting functions, aesthetic and cultural services, and so on. Although rivers’ ecological significance is well known, they are suffering worldwide ecological degradation. A key cause is water shortage in rivers. Over half of the accessible surface water in the world has been regulated by humans, and experts estimate that this proportion will increase to about 70% by 2025 [
1,
2]. Water resource exploitation activities, such as water impoundment and withdrawal, within-basin and interbasin water transfers, aquifer water exploitation for the water supply to agriculture, industry, and households and for hydropower generation, have obviously altered the natural hydrological regime of rivers [
3]. Nilsson et al. [
4] investigated 225 basins throughout the world and found that 83 and 54 rivers in these basins were highly or moderately regulated, respectively.
The concept of environmental flows (e-flows) has been advanced to meet ecosystem demands for water [
5]. E-flows are defined as the volume of water that should remain in a river and the variation of this provision over time, to maintain specific indicators of ecosystem health [
6]. A considerable number of methods have been developed to determine e-flows. E-flow determination methods can be classified into four major categories, i.e., hydrological, hydraulic rating, habitat simulation, and holistic methods. Hydrological methods are the simplest, relying primarily on hydrological data, usually in the form of naturalized, historical monthly, or daily flow records. They often set one or several proportions of flows, usually termed the minimum e-flow, to maintain the survival of target species or other highlighted ecological features at some acceptable level. Due to the simplicity and low data requirement, hydrological methods are widely used at the planning stage of water resource development or for a preliminary estimate of e-flows [
7]. From the 1970s onwards, rapid development of e-flow determination methods has taken place to try to develop a quantitative relationship between water discharges and the quality of instream resources, such as fishery habitats [
8]. These methods examine the effects of discharge increments on instream habitats [
9]. Hydraulic rating methodologies use changes in simple hydraulic variables, such as wetted perimeter and maximum depth, as surrogates for habitat-sensitive factors that are known or considered to be important to target biota. The implicit assumption underpinning these methods is that sustaining a habitat-sensitive hydraulic parameter lower than threshold values will maintain the target species and/or the biota. The physical habitat change against discharge changes is modeled. Commonly, the physical habitat is described by one or more hydraulic variables, including depth, velocity, substratum composition, cover, and complex hydraulic indices (e.g., benthic shear stress). The simulated habitat conditions are linked with the preferred habitat conditions for target species, often depicted by habitat suitability index curves. The resultant outputs are expressed in the form of habitat–discharge curves for one or several target species. Based on the curves, the optimum e-flows are recommended. These holistic methodologies have been advocated by freshwater ecologists for over a decade [
10,
11], as they try to determine the e-flows for the entire riverine ecosystem rather than for one or several target species. Holistic methodologies usually apply the tools for hydrological, hydraulic, and physical habitat analysis used in the other three types of e-flow determination methods to establish the e-flows of the riverine ecosystem [
12]. Nowadays, they also seek to incorporate some quantitative flow-ecology models, especially if they are expected to predict the ecological results after the implementation of the designed e-flows [
13].
Urban rivers (i.e., the river stretches in urban areas), as an important river type, have their own unique characteristics. During urbanization, a significant proportion of the world’s urban rivers have been channelized [
14]. This situation is not confined to a particular geographic region, but rather is a global problem [
15]. Habitat provision, pollutant dilution, and recreation are three common requirements considered in e-flow assessment of urban rivers. The habitat provision requirement could be assessed by the e-flow assessment methods of natural rivers mentioned in the previous paragraph. Mass balance is the key model to determine the e-flow requirement for pollutant dilution [
16,
17]. For the recreation requirement, the river bed should not be bare, and the water depth is usually set to no lower than 0.2 m [
18]. Besides the three common requirements, Yin et al. [
19] suggested that the requirement of hydrological connectivity should also be incorporated in e-flow assessment of urban rivers. Three hydrological connectivity scenarios are proposed, i.e., high, medium, and low hydrological connectivity, corresponding to different environmental protection targets.
However, previous studies on urban rivers have not considered the influence of channel morphological changes on e-flows [
5,
20]. Weirs are very common in urban rivers, which could increase the water elevation and consequently enhance the aesthetic effects. Meanwhile, weirs reduce flow velocity, leading to siltation in urban rivers. Thus, the morphology in urban rivers is dynamic. Previous research has assumed that the river morphology is fixed. Siltation is inevitable under the influence of weirs. It is necessary to explore the influence of morphological changes on e-flows.
This research will explore the influence of morphological changes on e-flows in urban rivers. In the following sections, we will introduce the e-flow assessment method established by Yin et al. [
19] This method is suitable for urban rivers with weirs. Then, a method is proposed to explore the influence of morphological changes on e-flows. The Shiwuli River, an urban river impacted by several weirs, is adopted as the case study.
3. Results
As sediments are trapped by weirs, the channel morphology will change. In this paper we assume that siltation initially takes place near weirs (upstream), and then gradually expands upward. To explore the effects of morphological changes on e-flow requirements, the siltation depths are set to 0.1–0.5 m at increments of 0.1 m. The new e-flow requirements under the high hydrological connectivity scenario are listed in
Table 1 (flow velocity = 0.1 m/s) and
Table 2 (flow velocity = 0.2 m/s). The tables show that the e-flow requirements are significantly impacted by the channel morphological changes for all river segments. With the increased depth of siltation, the e-flow requirements under the high hydrological connectivity scenario will decrease. This indicates that if the water resources available for e-flow supply are too limited to meet the e-flow requirements, modifying the channel morphology is a potential solution.
In addition, the change ratio of e-flow requirements under the high hydrological connectivity scenario is listed in
Table 3 and
Table 4. The change ratio of e-flow requirements varies for different river segments. Among them, the change ratio for river segment 6 is the largest, while the change ratio for segment 2 is the smallest. The e-flow change ratio is closely related to the weir height. The height of weir 6 is the lowest, at only 0.7 m. Thus, a siltation depth of 0.1 m is significant compared with the weir height.
The change ratio of e-flow requirements under a flow velocity of 0.1 m/s or 0.2 m/s is also different. The change ratio at a low flow velocity is obviously greater than that under a high flow velocity. For example, under a flow velocity of 0.1 m/s, the e-flow decreases by 6.75% and 16.72% for river segments 1 and 6 with a siltation depth of 0.1 m, respectively, while under a flow velocity of 0.2 m/s, the e-flow only decreases by 1.05% and 7.84% for river segments 1 and 6 with a siltation depth of 0.1 m, respectively. Thus, modifying channel morphology to aid e-flow supply is more effective under low flow velocity requirements.
In terms of medium hydrological connectivity, siltation also has obvious effects on e-flow requirements (
Table 5 and
Table 6). With an increase in siltation depth, the e-flow requirements decrease. For some river segments, the e-flow requirements even become 0. For example, when the siltation depth exceeds 0.3 m, the e-flow requirements for river segment 6 become 0. This means that no continuous water input is required to maintain hydrological connectivity (the evaporation and seepage are not considered) after the water level exceeds 0.2 m. This is because the channel slope for river segment 6 is very low. When the siltation depth increases by 0.3 m, the channel slope becomes nearly 0 and thus the flow velocity decreases to nearly 0.
Besides the significant changes in e-flow requirements in river segments 6 and 1, the change ratio in segments 7 and 8 also exceeds 20%. However, the change ratio of e-flow requirements for river segment 2 is very limited. It is only 6.57% when the siltation depth increases by 0.5 m. This is due to the highest channel slope being in river segment 2.
Under the low hydrological connectivity regime, e-flow requirements also decrease with an increase in siltation depth. In this research the water storing time
T0 is set to three days. The e-flow reduction ratio for segments 6‒8 is the greatest, higher than 70%. The reduction ratio for segments 1 and 4 is the lowest, no more than 30% (
Table 7 and
Table 8).
Table 1,
Table 2,
Table 3,
Table 4,
Table 5,
Table 6,
Table 7 and
Table 8 also show that the sensitivity of e-flow requirements to morphological changes is different under high, medium, and low hydrological connectivity. Under the high hydrological connectivity scenario, river segment 1 is the least sensitive and the sensitivity sequence is segment 6 > segment 8 > segment 3 > segment 7 > segment 5 > segment 4 > segment 2 > segment 1 for a planned flow velocity of 0.2 m/s. while under the medium hydrological connectivity scenario, river segment 2 is the least sensitive, and the sensitivity sequence is segment 6 > segment 1 > segment 7 > segment 8 > segment 3 > segment 4 > segment 5 > segment 2; under the low hydrological segment, the sensitivity sequence is segment 6 > segment 7 > segment 8 > segment 3 > segment 5 > segment 2 > segment 4 > segment 1. The difference in the sensitivity sequences is due to different primary influence factors. For high hydrological connectivity, the relative weir height is the primary influence factor, while for medium hydrological connectivity the channel slope is the primary influence factor, and for low hydrological connectivity all the channel morphological characteristics, including channel width, length, and slope, exert an influence.
4. Discussion
4.1. Tradeoff between Flood Control and E-Flow Supply
Although siltation could reduce the e-flow requirements, which is positive in terms of sufficient e-flow, siltation is negative for flood control. For the Shiwuli River, the design standard is a 100-year flood.
Table 9 lists the flow magnitude that the river can tolerate under different siltation depths. With an increase in siltation depth, the maximum flows allowed that will not lead to flooding decrease. When the siltation depth increases to 0.5 m, the reduction ratio of flood magnitude allowed exceeds 20% for all river segments (
Table 10). Especially for segment 7, the reduction ratio is 37.06%. Segment 7 should be given the highest priority for siltation control to avoid flood disasters.
The benefits of siltation for e-flow supply and flood control are compared in
Figure 2,
Figure 3,
Figure 4 and
Figure 5. When the planned flow velocity for e-flow requirements is 0.1, the e-flow reduction ratio in river segments 3, 5, and 6 is greater than the flood reduction ratio; in segments 3 and 6 only, the e-flow reduction ratio is greater than the flood reduction ratio. Thus, if siltation is used to satisfy the e-flow requirements, segments 3 and 6 should be the primary area. The flood reduction ratios for segments 1, 2, 4, and 7 are greater than the e-flow reduction ratio under a planned flow velocity 0.1 m/s or 0.2 m/s. In these four river segments, dredging has a more significant effect.
Under the medium hydrological connectivity scenario, the comparison results between e-flow supply and flood control are different from those under the high hydrological connectivity scenario. The e-flow reduction ratio in segments 1 and 6 is obviously greater than the flood reduction ratio, while in the other river segments the e-flow reduction ratio is lower. Thus, if siltation is used to satisfy the e-flow requirements, segments 1 and 6 should be the primary area to allow siltation for e-flow satisfaction under the medium hydrological scenario.
Under the low hydrological connectivity scenario, the comparison results between e-flow supply and flood control are obviously different from those under the high and medium hydrological connectivity scenarios. The e-flow reduction ratio in segments 2, 3, 5, 6, 7, and 9 is obviously greater than the flood reduction ratio. In these segments, siltation may be a potential solution to satisfy the e-flow requirements. In comparison, the flood reduction ratio is greater than the e-flow reduction ratio in segment 1, and thus, dredging should be adopted in this segment.
Due to the positive effects of siltation on e-flow supply, the timing, frequency, and volume of dredging should be further optimized. Dredging could begin before the flood season. A 100-year flood is the design standard for the Shiwuli River. The standard is very high and such large floods represent a very low possibility. A reliable flood prediction model is valuable to balance the needs of flood control and e-flow supply. With a reliable flood prediction model, river managers could determine the frequency and volume of dredging. If a specific flood risk is allowed, the dredging frequency and volume could be reduced to some degree.
4.2. E-Flow Comparison under Different E-Flow Assessment Methods
In this research, following the new e-flow assessment method established by Yin et al. [
19], we explore the influence of channel morphological changes on e-flow requirements. However, as mentioned in the introduction, besides the method by Yin et al. [
19], hydrological, hydraulic rating, habitat simulation, and holistic methods are four commonly used method types to determine e-flows. These four methods usually assume that the channel morphology does not change. It is necessary to further explore the influence of channel morphological changes on e-flow requirements under the four methods.
The hydrological methods are based on hydrological data. The underlying assumption is that hydrological data, i.e., historical flow data, can effectively reflect the characteristics of a river. However, the historical flow data cannot reflect continuous changes in channel morphology. If hydrological methods are to be adopted in future research (when only hydrological data are available), the latest flow data, which correspond to the new channel morphology, should be used to determine e-flows.
Morphological data are necessary for the hydraulic rating, habitat simulation, and holistic methods. In hydraulic rating methodologies, hydraulic variables, such as wetted perimeter, will be determined. These hydraulic variables are closely related to river channel morphology. The e-flow results under the hydraulic rating methodologies can reflect the influence of river morphological changes. Although the e-flow results under the habitat simulation and holistic methods are influenced by changes in river morphology, in real-world research the river morphology is usually assumed to be fixed to reduce the complexity of e-flow assessment and management. To reduce the deviation of e-flow results, it is suggested that we update the e-flow results every few years, especially for rivers with a high sediment concentration.