Dam control changes the hydraulic environment and hydrological rhythm of rivers; to explore the causes of fish habitat loss and influence of river regulation in such rivers, the downstream curving reaches of Shaying River were cut off. Grass, silver, and bighead carp were selected as the research targets, and the impact of waterway regulation and the effect of gravel mounds on the target fish habitat were studied using numerical simulation. The results show that the original channel's hydraulic environment is similar to the artificial channel, and the target fish spawning habitat is limited, which is an important reason for the shortage of fish resources in the research area. The reconstructed channel can meet the river's navigation demands and improve the channel's flood discharge capacity. The habitat method was used to estimate the minimum and appropriate ecological flows of the target fish at each growth stage in the original and reconstructed rivers. Under the appropriate ecological flow, the weighted usable area (WUA) of juvenile, migratory, and spawning habitats of the target fish in the reconstructed river increased by 168, 24, and 580%, respectively, compared to the original channel, while the WUA of adult fish that preferred a slow-water environment decreased by 62%, effectively improving the target fish's spawning environment.

  • The habitat distribution of the four major Chinese carps in the whole life cycle was screened and analyzed.

  • The study shows the complete research process from waterway regulation to biological habitat restoration to ecological flow calculation.

  • A biological habitat restoration scheme with gravel mounds was adopted, and the navigation and flood discharge capacity of the waterway was not affected.

Graphical Abstract

Graphical Abstract
Graphical Abstract

Sluice construction and river regulation projects are widely used in hydropower development, flood control, disaster alleviation, urban water supply, etc. (Shukla & Mishra 2019). The curving cutoff is a common river regulation measure (Slowik 2016). Manually cutting and straightening some curved rivers can effectively improve rivers’ navigation conditions and flood discharge capacity. However, while this type of hydraulic engineering provides significant economic benefits, it also alters the natural hydrological rhythm and hydrodynamic conditions of the rivers, resulting in the simplification of aquatic community structures and the destruction of the river ecosystem (Nakamura et al. 2014). This study focuses on considering the ecological function of rivers while meeting their essential functions and the needs for production and existence.

To improve the quality of biological habitat in canalized rivers, scholars have carried out research and practice for many years and accumulated a rich research foundation. The concept of ‘ecological engineering’ was introduced by the American ecologist Odum and Odum in 1962, combining ecological and engineering theories (Odum & Odum 2003). In the mid-1970s, Europe began to carry out ecological engineering practices for river regulation (Seifer 1983). At present, the physical habitat restoration method is used mainly to restore and maintain biodiversity through river geomorphology transformation. Chang et al. (2019) monitored the biological communities around spur dikes used for waterway regulation, indicating that the spur dike group could provide habitat and shelter for fish and benthos to a certain extent, increasing benthos diversity and fish richness. Lee et al. (2010) improved the habitat of endangered fish species by arranging boulders in straight canalized rivers, providing a reference for restoring urban river ecosystem diversity. Mozzaquattro et al. (2020) sampled and analyzed a fish community under different river forms and water retaining structures, and considered that a reasonable arrangement of obstacles in the river had a positive impact on the structure of the fish community.

Physical habitat reconfiguration measures facilitate the maintenance of biological habitat areas under certain river conditions, so it is particularly necessary to determine the optimum conditions for organisms within the reconstructed river. A highly developed river is seriously canalized, and its downstream often suffers from dehydration, causing serious damage to the rivers biological habitat. Maintaining a suitable flow in a river is a prerequisite for ensuring the survival of aquatic organisms. Furthermore, an important approach to resolving this problem is to scientifically determine the ecological flow of such rivers. Chen et al. (2016) proposed the concept of minimum ecological flow: the minimum flow in a river to maintain the river's basic ecological function and shape. Chen et al. (2007) defined appropriate ecological flow as being the water condition of the critical state of ecosystem decline. Currently, the calculation methods of ecological flow mainly include hydrological, hydraulic, habitat, and comprehensive methods (Abdi & Yasi 2015), among which the hydrological and hydraulic methods have no specific ecological objective, making it difficult to accurately reflect ecological requirements. Therefore, these two methods are often combined with water resource planning, which has macroscopic guiding significance for the ecological management of rivers. A comprehensive method for ecological flow considers some of the ecological indicators that have the greatest influence on the protection of the biological habitat. However, this method involves many fields and has a large span of disciplines, making it difficult to be widely adopted. The habitat method, which is extensively used in ecological flow assessment around the world (Liu et al. 2022), is based on the relationship between the quality of the aquatic habitat and the flow in the river and then determines the river environment suitable for the survival of the species to be protected.

The main habitat simulation models include the1-D PHABSIM, EVHA and HABITAT models, the 2-D RIVER2D model and the 3-D DELFT3D software, etc. A 1-D model cannot adequately simulate the hydrodynamic characteristics in habitat assessment, while DELFT3D is mostly used for the analysis of the internal flow field of water bodies, especially in the aspect of pollutant diffusion (Ma et al. 2020). RIVER2D closely combines hydraulic simulation with ecological simulation, which is mainly used for the simulation and evaluation of microhabitats (water depth, velocity, and bottom sediment, etc.). Park et al. (2020) used a RIVER2D model to calculate the ecological flow of endangered fish in the Gam River, and proposed the concept of a drought index based on it, which has certain guidance for the management of upstream dams. Zhou et al. (2019) studied the water requirement to meet the habitat connectivity of Schizothorax during the spawning period, providing a reference for resource regulation and biological habitat conservation in the Yanni wetland. Hou et al. (2020a) compared the ecological flows calculated by the Tennant and habitat methods, and concluded that the habitat method was better to ensure the habitat requirements of Schizopygopsis younghusbandi during spawning in the Yarlung Tsangpo River. However, the above studies on ecological flow mostly focus on evaluating ecological flows of fish at a particular stage in a natural river, and few have looked at the complete process from river regulation to river habitat restoration.

The Shaying River is the largest tributary of the Huaihe River, and it is not only an important channel for flood discharging, water supply, navigation, and irrigation but also a typical sluice and dam control river. This study takes the river regulation measures of curving cutoff and curved river sections downstream to meet the channel's demand for navigation and flood discharge. Additionally, this study investigates the impact of project implementation on aquatic organisms in the river, taking grass, silver, and bighead carp as the research object and combining them with the fish habitat requirements of the juvenile, adult, migration, and spawning stages. The curving cutoff channel is naturalized by arranging gravel mounds, and the change in biological habitat in the reconstructed river is compared with that of the original river. The response relationship between the flow discharge and the weighted usable area (WUA) of the target fish is analyzed, and the ecological flow of the study reach is determined using the habitat method. The results can provide a reference for the naturalization of river biological habitats.

The Shaying River originates from the Funiu Mountains in Henan Province. It flows through the central part of Henan Province from northwest to southeast and then enters Anhui Province in Jieshou City and finally into the Huai River at the Mohekou in Yingshang County. The Shaying River basin has a high level of development and use of water resources. The basin is a typical artificially transformed river, containing more than 1,800 large- and medium-sized dams and more than 1,400 sluices. The study reach is located in the lower reaches of the Shaying River, 14.2 km from the Mohekou (Figure 1). The length is about 7 km, with an average width of 200 m and a U-shaped cross-section. The depth of the valley is 8–12 m, the channel gradient is 1:3000–1:10000, the roughness coefficient of the river bed is 0.02–0.03. The study area is a warm temperate semi-humid continental monsoon climate with obvious seasonal alternation. According to the records of Yingshangzha Hydrological Station, the maximum annual precipitation is 1,722.5 mm, the minimum is 406.8 mm, and the average annual precipitation is 946.4 mm. The interannual and annual precipitation distribution is very uneven, 40–65% of the annual precipitation is concentrated in June to September, and continuous dry or continuous abundant conditions often occur. The temporal non-uniformity of precipitation distribution leads to large flood peak flow in the flood season, small flow in the dry season, and frequent flow interruption in Shaying River basin.

General situation of fish resources

The aquatic organisms in the study reach are affected by the upstream sluices and dams. In the results of Shu & Wei (2015)'s investigation on the fish composition of the Shaying River, a total of 36 species of fish were collected, including 26 species of Cyprinidae, a decrease of 25% compared to the number of species collected by Li & Wang (1983) from 1973 to 1980. With an average body length of 116.1 mm and average body weight of 62.7 g, the individual fish are small, while the overall fish population density is low and fishery resources are scarce. The downstream reaches of Shaying River are mainly composed of grass carp (Ctenopharyngodon idella), silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) (of the four major Chinese carps), and piscivorous fishes such as culter alburnus, channa argus, and some small omnivorous fish, with grass, silver and bighead carp accounting for over 60% of the total number of fish.

Fish are the climax community in an aquatic ecosystem, and they play a crucial role in the existence and abundance of other communities, particularly in a river ecosystem. Therefore, fish are a better indicator species for ecological flow (Yan et al. 2009). Grass, silver, and bighead carp were chosen as the target species among the four major Chinese carps owing to their high economic value and adaptability to flow conditions. Their spawning types include drifting eggs, with grass carp preferring to live in the lower layer of the water body, silver and bighead carp both inhabiting the upper layer of the water body, which are the main fish cultivated in reservoirs, rivers, and lakes. Moreover, silver carp are filter-feeding fish, which can effectively improve water quality.

General situation of waterway regulation project

The channel in the study section bends in the shape of an ‘S’ due to long-term flow erosion. The Shaying River channel is a natural and canalized grade IV standard river. According to the Ministry of Transport (2014), the water depth of the waterway is 2.3–2.7 m, and the minimum curvature radius is 330 m. The curvature radius of the study river section is just 250 m, and the flood control standard is less than once in 20 years. To enhance the flood carrying capacity of the river and meet navigation conditions, a river regulation scheme of curving cutoff is proposed in this study. As depicted in Figure 2, the starting section pile number of the curving cutoff channel is 2 + 120, the ending section pile number is 7 + 990, and the length is about 4.3 km. The river channel was redesigned, and the mainstream widened to 300 m with floodplains and artificial levees installed outside the mainstream.
Figure 1

River system map of Shaying River Basin; (a) Model verification area, (b) waterway regulation area.

Figure 1

River system map of Shaying River Basin; (a) Model verification area, (b) waterway regulation area.

Close modal
The roadmap of the study is shown in Figure 3. The topographic and fish investigation data of the study area need to be obtained first. The topographic data in the study come from the high-precision topographic survey information of Anhui Provincial Port & Shipping Group Co., Ltd, and the R2D-bed and HEC-RAS river geometric model of the riverbed are established accordingly, where HEC-RAS is used to calculate the water level-discharge relationship for each study condition as the boundary condition of River2D hydrodynamic simulation. The R2D-mash module is used for grid division of the riverbed model, to simulate the flow field in the river. Fish investigation data identifying the target species and establishing their suitability index curves for each physical factor were used, coupled with River2D hydrodynamic simulation, screening and calculation of the combined suitability and habitat area of the target species. As the hydrodynamic environment in the river varies considerably with changing hydrological conditions, which then affects aquatic habitat, it is necessary to set up various flow condition sequences to more accurately and comprehensively explore the habitat changes and determine the ecological flow in the river.
Figure 2

River channel before and after the curving cutoff.

Figure 2

River channel before and after the curving cutoff.

Close modal

Hydrodynamic model

The River2D hydrodynamic module, an average depth finite element model based on the two–dimensional Saint–Venant shallow water equation (Alabyan & Lebedeva 2018), was chosen to simulate the hydrodynamic environment in the study area. It assumes that the vertical direction of the water depth conforms to hydrostatic pressure, and the velocity is constant along the depth. Coriolis and wind forces are assumed negligible. The governing equations of the model are as follows:

  • (1)
    Conservation of mass:
    formula
    (1)
  • (2)
    Conservation of x-direction momentum:
    formula
    (2)
  • (3)
    Conservation of y-direction momentum:
    formula
    (3)
where H is the water depth, u and ʋ are the velocities in the x and y directions respectively. qx, and qy are the flow discharge in the x and y directions respectively; g is the gravitational acceleration and ρ is the density of water. S0x and S0y are the bed slopes in the x and y directions; Sfx and Sfy are the corresponding friction slopes. τxx, τxy, τyx, and τyy are the components of the horizontal turbulent stress tensor.

Bed shear stresses are assumed to be related to the magnitude and direction of the depth–averaged velocity. In the x direction, the friction slope can be calculated by the following equation:
formula
(4)
where τbx is the bed shear stress in the x direction and Cs is a Chezy coefficient.

Mesh subdivision

Delauney triangulation is adopted in the study model, and the mesh density was set according to the accuracy of the calculation results and the complexity of the riverbed topography. The original river channel model has 56,828 nodes, with a total of 113,600 elements; the mesh size of the mainstream is 3–6 m, and that of the floodplain area is between 5 and 10 m. The curving cutoff river channel is regular, with 24,738 nodes, 84,283 elements, and 8–15 m mesh size in the model. Local refinement was conducted for the area where gravel mounds are arranged and at the junction of the original river channel, to simulate the shape of the gravel mounds, the minimum mesh size was 2 m.

Determination of boundary conditions

Pile numbers 1 + 650 and 8 + 830 were defined as upstream and downstream boundaries, respectively, and their types are flow inlet and water level boundary outlet, respectively. The numerical calculations of the original and curving cutoff river channel were made under a flood discharge of a 20-year return period, to investigate the improvement effect of the curving cutoff project on the flood discharge capacity of the original river. When evaluating the quality of fish habitat in the river channels during the flood season, the flow condition was set at 1,640 m3/s, decreasing to 240 m3/s with 200 m3/s as the flow gradient tolerance. Also, when evaluating the quality of fish habitat during the dry season, the flow condition was set at 100 m3/s, decreasing to 10 m3/s with 10 m3/s as the flow gradient tolerance. Additionally, a 5-year return period flood discharge, average annual discharge of 140 m3/s, and an extremely low discharge of 5 m3/s were set as contrasts, totaling 22 discharge conditions. Using HEC-RAS software, based on Bernoulli's equation (Equation (5)), the upstream and downstream boundary water levels of each operating condition were solved using a step-by-step method, and the discharge-water level relationship under each working condition was obtained, which was used as the boundary condition of numerical simulation:
formula
(5)
where Z1, Z2 are the riverbed elevation; H1, H2 are the water depth of the cross section; α1, α2 are the kinetic energy correction coefficient; ʋ1, ʋ2 are the average velocity of the cross section, and he is the head loss.

Habitat model

The instream flow incremental methodology was adopted to study biological habitats, which uses the WUA to assess habitat quality (Yang et al. 2021). The WUA represents the suitable physical habitat area of the target species, which is as follows:
formula
(6)
where f (Vi, Di, Ci) is the suitability composite function, Vi is the velocity suitability index, Di is the water depth suitability index, Ci is the channel suitability index, and Ai is the area of each element mesh. Because of the lack of suitable data for the target species for the bottom sediment and covering, Ci is 1 by default.
In this study, the water area of the curving cutoff channel was significantly reduced compared with that of the original channel. It is not sufficient to evaluate the habitat quality of the target species using only the WUA value, so the percentage of available area or percent usable area (PUA) (Hung et al. 2022) was introduced as a supplementary basis for evaluating habitat quality. The relationship between the WUA and the PUA is as follows:
formula
(7)
where TA is the total area of the river.

The corresponding flow pattern distribution was obtained from the numerical simulation results of the hydrodynamic environment in the river under different discharges. According to the suitability curve of the target species to the river characteristic conditions, the suitability index of the grid area under each working condition was determined, and the WUA and PUA under each discharge condition were calculated to estimate the ecological flow in the river.

Model verification

As there is readily and publicly available water level monitoring data for the Yingshang sluice and detailed measured cross-sectional velocity data, the Yingshang sluice area was chosen to assess the reliability of the model, and its location is shown in Figure 1(a). To facilitate the establishment of the river model, the model validation flow condition is 1,230m3/s, at this time the Yingshang sluice gates are fully open, and the river inlet and outlet measured water levels are 22.95 m and 22.30 m respectively. The simulation ignores the bridge piers and the sluice piers in the river; the roughness range of the natural channel in the model area is 0.0225–0.0275, that of the floodplain area is 0.0375, that of the artificial slope protection is 0.020–0.025, and that of the sluice pier and sluice flour slab is 0.015. Set the upstream flow inlet and the downstream water level boundary; the water level values at each measurement point and the flow velocity values on the section 0—0 were obtained by hydrodynamic simulation and compared with the measured flow velocity and water level data. The locations of the water level measurement points are shown in Figure 4; water gauges at points 1#, 3#, 4# and 7# provided the actual water level observation data. To fully verify the simulation results of each area, the water level data observed in the model test were compared with those obtained by Ye (2007) (Table 1). The results show that the distribution trend of measured and simulated velocity along the 0—0 section direction is basically consistent (Figure 5), the simulated water level is smaller than the measured level, but the relative error between them is less than 0.10 m, and the relative error with the physical model water level is less than 0.20 m. This indicates that the simulation results are accurate and the model is reasonable.
Table 1

Comparison of water surface elevation (unit: m)

Staff gauge point1#2#3#4#5#6#7#
Actual measurement 22.91  22.62 22.41   22.31 
Physical model (Ye 200722.88 22.80 22.70 22.48 22.39 22.35 22.30 
Numerical simulation 22.81 22.74 22.55 22.36 22.36 22.33 22.30 
Staff gauge point1#2#3#4#5#6#7#
Actual measurement 22.91  22.62 22.41   22.31 
Physical model (Ye 200722.88 22.80 22.70 22.48 22.39 22.35 22.30 
Numerical simulation 22.81 22.74 22.55 22.36 22.36 22.33 22.30 
Figure 3

Research framework.

Figure 3

Research framework.

Close modal
Figure 4

Model verification area.

Figure 4

Model verification area.

Close modal
Figure 5

Comparison of flow velocity at 0—0 cross section.

Figure 5

Comparison of flow velocity at 0—0 cross section.

Close modal

Waterway regulation influence and biological habitat naturalization

To more accurately and comprehensively evaluate the impact of river regulation on the habitat and life cycle of many types of fish, through a field investigation and literature research of the target fish habitat, combined with the relevant experience of fishery workers, the flow velocity and water depth indexes were selected as the main environmental factors affecting the survival of fish, and the suitability index curves of the target species in different life stages were constructed accordingly (Figure 6). A habitat suitability index (HSI) was used to quantify the suitability of the target fish habitats, and the higher the HSI value, the better the suitability (Zhou et al. 2022).
Figure 6

Water velocity and depth suitability index curve for target fish; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Figure 6

Water velocity and depth suitability index curve for target fish; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Close modal

The growth stage of the target fish can be classified into juvenile and adult stages, and the habitat environment in these stages requires the enrichment of nutrients, most of which are slow-water environments. In the feeding stage of juvenile fish (July to August), the optimum velocity range is 0.12–0.2 m/s, and the optimum water depth range is 0.7–1.6 m (Bai et al. 2013). The optimum velocity range of adult fish in the feeding stage (July to November) is 0.1–0.25 m/s, and the optimum water depth range is 2–4 m. The four major Chinese carps are all spawning migration fishes. The optimum range to meet the migratory passage of the target fish is 0.4–0.6 m/s, and the optimum water depth range is 2–4 m (Dang et al. 2018).

The spawning stage is the most fragile and important in the fish life cycle. The hydrodynamic habitat conditions in this stage are of great significance for protecting fish species resources and habitat restoration. The four major Chinese carps are typical drifting egg fish, and the fish eggs need to drift in the water for some time before they can develop into juvenile fish with swimming ability. In this process, if the eggs sink to the bottom and break or are buried in the sand, hypoxia and other phenomena will cause the eggs to die. Hence, it is necessary to meet certain hydrodynamic conditions to keep the eggs suspended in water. With economic development and urbanization, examples of hydraulic engineering, such as artificial transformation and sluice dams, are common in natural rivers. The hydraulic environment in a river changes greatly, which may lead to flow velocity, pattern, and other conditions that cannot maintain the drifting and hatching of fish eggs, destroying the complete life cycle of fish, affecting the number of communities, and even endangering the survival of species. The optimum velocity range of the target fish spawning habitat is between 1.4 and 1.6 m/s, and the optimum water depth range is between 0.75 and 1.5 m (Zhang 2017; Dang et al. 2018).

Influence of the curving cutoff

Figure 7(a) shows the Froude number (Fr) of the river before and after the curving cutoff under a 20-year flood discharge design return period (3,760 m3/s). The water flow pattern in the original channel is mainly subcritical flow, and the area is evenly distributed in the mainstream. The Fr value ranges between 0.13 and 0.27. However, when the floodplain area is submerged, the Fr value ranges between 0 and 0.18, with the standing water area being widely distributed. The Fr value of the curving cutoff channel is mainly concentrated in the range of 0.15–0.32, and the Fr in the floodplain area is larger than that of the original channel. The Fr value alternately changes in the primary flow direction of the curving cutoff channel, and the Fr value in the vertical primary flow direction decreases toward both banks. At the same time, compared with the shear velocity magnitude of the river before and after the curving cutoff in Figure 7(b), the shear velocity magnitude in the mainstream of the original channel remains almost unchanged, and there is an obvious grading phenomenon in the floodplain area. There is a sudden change in the intersection of the curving cutoff channel and the original river, and the shear velocity magnitude in the curving cutoff channel is relatively small, indicating that the flood discharge in the curving cutoff channel is unobstructed. Combined with the analysis of water surface elevation, there is a slight hydraulic drop at the beginning of the curving cutoff channel, and then the flow tends to be stable along the course. At the end of the curving cutoff channel, local backwater will occur because of the narrowing of the river.
Figure 7

Comparison of the river before and after diversion; (a) Fr, (b) shear velocity magnitude.

Figure 7

Comparison of the river before and after diversion; (a) Fr, (b) shear velocity magnitude.

Close modal

After the curving cutoff of the original channel, the shape of the riverbed changed significantly. The analysis of the Fr value and shear velocity magnitude in the channel shows that the flow velocity in the channel increases while the viscous force between the flow layers decreases, which is beneficial for improving the flood discharge and sediment transport capacity of the river channel. Under 3,760 m3/s, the water depth of the original channel is mainly concentrated at 11–14 m, and the water depth of the curving cutoff channel is about 8.7 m. The water level of sections 1 + 650 and 8 + 830 dropped by 2.53 and 1.93 m, respectively. At the same time, the flow velocity of the original channel is mainly concentrated at 1.5–2.6 m/s, and the flow velocity is evenly distributed in the mainstream. The flow velocity in the curving cutoff channel has increased but the flow velocity distribution is similar to that of the original channel, and the velocity in the mainstream is 1.6–2.9 m/s.

Habitat classification is an important aspect of river ecological assessment. The target fish are screened for a suitable environment at each growth stage and classified into different habitat types according to flow patterns. Hydraulic gradients or Froude numbers are widely used to identify habitat types worldwide. Based on the Fr value, the habitats of the target fish are classified into three types (Jowett 1993): supercritical flow (0.18 ≦ Fr < 0.41), riffle (Fr ≧ 0.41), and deep pools (0 < Fr < 0.18). The suitable habitat types of the target fish at each growth stage are classified according to the above methods, as shown in Table 2.

Table 2

Distribution of target fish habitat types (unit: m2)

River typeOriginal channel
Curving cutoff channel
Habitat typesRiffleDeep poolsSupercritical flowRiffleDeep poolsSupercritical flow
Juvenile 10 228,828 1,941 68,440 885 
Adult 228,430 187 22,085 60 
Migration 22 215,923 31,559 16 50,735 15,479 
Spawning 10,130 15,402 56,097 2,850 7,668 17,452 
River typeOriginal channel
Curving cutoff channel
Habitat typesRiffleDeep poolsSupercritical flowRiffleDeep poolsSupercritical flow
Juvenile 10 228,828 1,941 68,440 885 
Adult 228,430 187 22,085 60 
Migration 22 215,923 31,559 16 50,735 15,479 
Spawning 10,130 15,402 56,097 2,850 7,668 17,452 

Under 3,760 m3/s discharge, the deep pool habitats of the target fish in the original river are the most widely distributed in the entire life cycle. The area of all types of habitats decreases to varying degrees after the curving cutoff, with the adult fish deep pool habitat showing the greatest reduction, at only 10% of that of the original river. The riffle and supercritical flow habitats in the original and cutoff rivers are quite limited, and the hydraulic environment in the river channel is relatively simple. At this point, the original river is quite similar to the flow pattern of the artificial canalized river. According to relevant research, the habitat stability of aquatic organisms in this type of river environment is poor, which is also an important reason for the lack of fish resources in the lower reaches of the Shaying River.

Naturalization

Referring to the restoration measures in previous research (Nan et al. 2021), it was proposed that the ecological restoration scheme of arranging gravel mounds in the river channel be adopted. As the flow on the concave bank of the river was obviously affected by the gravel mounds, an arrangement scheme of several groups of gravel mounds had the most obvious effect on the restoration of the fish habitat. At the same time, to ensure the navigation requirements of the channel, the sites selected for the gravel mounds were positioned at five bisectrixes of the channel concave bank. Each group of gravel mounds was composed of three single gravel mounds. Considering the influence of the gravel mounds on the navigation and flood discharge capacity of the river, it was determined that the edge length of a single gravel mound was 20 m, the distance between adjacent mounds was 9 m, and the height was 5 m. The inside of the mounds was filled with blocks, and the outside was fixed using ecological cages (Figure 8).
Figure 8

Naturalized reconstruction scheme of the river channel.

Figure 8

Naturalized reconstruction scheme of the river channel.

Close modal

Analysis of flood control in reconstructed reach

The arrangement of the gravel mounds is conducive to the formation of a certain deep pool–riffle pattern in the canalized channel, enriching the flow pattern in the channel and increasing the habitat stability of biological communities. Its essence is to use the retaining water effect of the gravel mounds, so that a local backwater area will be formed near the gravel mound group, and a certain range of low-speed flow area and vortex flow will appear downstream. The water level of each section 20 m upstream of the gravel mound groups was monitored under a flood discharge of 3,760 m3/s to evaluate the influence of the gravel mound groups' arrangement on the flood discharge capacity of the curving cutoff river. An average of 25 measuring points were set up in each section, and an average value of 25 points was taken as the water surface elevation of this section, which was then compared to the river channel without the arrangement (Table 3). Under this flow condition, the water overflowed the main channel and flooded the gravel mound groups. For the convenience of comparison, the locations of the gravel mound groups are marked in Figure 9.
Table 3

Average water surface elevation of a typical section of the curving cutoff river channel

Section1–12–23–34–4
Water surface elevation/m Channel with gravel mound groups 23.414 23.364 23.339 23.295 
Unarranged channel 23.401 23.356 23.327 23.287 
Section1–12–23–34–4
Water surface elevation/m Channel with gravel mound groups 23.414 23.364 23.339 23.295 
Unarranged channel 23.401 23.356 23.327 23.287 
Figure 9

Water surface elevation of the curving cutoff river channel before and after the arrangement of gravel mound groups; (a) Channel with gravel mound groups, (b) unarranged channel.

Figure 9

Water surface elevation of the curving cutoff river channel before and after the arrangement of gravel mound groups; (a) Channel with gravel mound groups, (b) unarranged channel.

Close modal

Generally, there is little change in the water surface elevation in the river before and after the arrangement of the gravel mound groups. Local backwater will be formed on the concave bank of the channel without arrangement, and the flow is most obvious when it enters the curving cutoff channel from the original river channel. The maximum water level elevation difference of each measuring point in section 1–1 is 0.046 m. After arranging the gravel mound groups, the backwater area on the concave bank side increased, in which the range of water level elevation at each measuring point of section 1–1 increased to 0.085 m. The average water level of section 1–1 is 0.013 m higher than that of the river channel without arrangement, and the water level increase in other sections is less than that of section 1–1. The average water level elevation difference of each section is less than 0.02 m, and the water level range of each section's measuring point is less than 0.05 m. To sum up, the arrangement of the gravel mound groups has no obvious influence on the overflowing capacity of the curving cutoff channel, which can meet certain flood discharge requirements.

Ecological flow research

This study's purpose is to reconstruct a complete biological habitat; therefore, determining ecological flow is crucial. Referring to the scheme of the habitat method in related research (Fu et al. 2021), the flow corresponding to the obvious turning point of the Q-WUA curve is taken as the minimum ecological flow of the corresponding stage of the target fish in the original river, and the flow corresponding to the highest point is the appropriate ecological flow. The biological habitats of the juvenile, adult, migration, and spawning stages of the target fish were numerically simulated to estimate the ecological flow in the river more comprehensively and accurately. The flow conditions were classified into two categories: the flow condition in wet (240–1,840 m3/s) and dry periods (10–100 m3/s). Additionally, flood discharges of 5 and 140 m3/s for 5- and 20-year return periods, respectively, are specified as the characteristic flow conditions as references and supplements for ecological flow study. To make the ecological flow determination process more intuitive, only the flow conditions with relatively drastic changes are displayed when drawing the Q-WUA curves (Figure 10). The calculation results under other large flow conditions will not influence the determination of ecological flow and are briefly analyzed.
Figure 10

Relationship between flow discharge and fish WUA; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Figure 10

Relationship between flow discharge and fish WUA; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Close modal

Juvenile fish

The habitat area of juvenile fish feeding in the original river is less affected by the change in flow discharge. The obvious turning point of the Q-WUA curve is 20 m3/s, and the flow discharge is determined to be the minimum ecological flow of the original river. When the flow discharge is 40 m3/s, the WUA of the original river reaches the maximum value of 192,961 m2, the PUA is 6.44%, so 40 m3/s is determined to be the appropriate ecological flow of the original river. As the discharge gradually increases, the WUA gradually decreases and stabilizes between 21,000 and 26,000 m2. When the discharge is greater than 1,640 m3/s, the floodplain area begins to submerge, the WUA increases again, reaching 230,779 m2 under a flood discharge of 3,760 m3/s. Although the increased water area brings uncertainty to flood control safety and biological habitat, it cannot be used as a stable habitat for fish; thus, it is not considered an ecological flow.

In reconstructing the river channel, the juvenile feeding habitat is made significantly larger than the original river. This is because compared to the large-scale deep water and slow-water environment of the original river, the velocity in the river increases, and the water depth decreases after the curving cutoff, which can better meet the feeding demand of juvenile fish. The flow discharge of 40 m3/s is the obvious turning point of the Q-WUA curve, indicating that this flow is the minimum ecological flow of the reconstructed river. When the flow discharge is 60 m3/s, the WUA reaches the maximum value of 523,207 m2, the PUA is 27.05%, and indicates that 60 m3/s is the appropriate ecological flow for the reconstructed river. With the increase in discharge, the velocity and water depth further increases, and the WUA decreases rapidly. After 240 m3/s, the WUA became smaller than the original river, with the smallest being at 840 m3/s, which is 7,000 m2. In the process of a subsequent increase in discharge, WUA increases slightly. The WUA increases to 69,835 m2 under a flood discharge of 3,760 m3/s, indicating that the juvenile fish habitat in the reconstructed river is greatly affected by the flow discharge. Therefore, the flow discharge in the reconstructed river should not be greater than 240 m3/s at the juvenile stage of the target fish.

Adult fish

The flow pattern of the original river is mainly in a slow or still water state, and the target fish prefer a slow-water environment during adult fish feeding, so the habitat of adult fish in the original river is widely distributed. The WUA of adult fish increases rapidly in the flow discharge range of 5–40 m3/s, where the obvious turning point of the Q-WUA curve is 30 m3/s, which is determined to be the minimum ecological flow of the original river, and WUA = 480,985 m2 is the peak value of the curve at 40 m3/s, and the PUA is 16.05%, which is determined to be the appropriate ecological flow of the original river.

Because the river is widened by 100 m when the curve channel is cut and straightened, and the water depth in the reconstructed channel is smaller under low discharge, the adult fish habitat obviously decreases, and the flow velocity gradually increases when the flow increases. Therefore, the adult fish habitat does not obviously increase in the high-flow environment. The turning point of the Q-WUA curve of the reconstructed river is also 30 m3/s, which is determined to be the minimum ecological flow of the reconstructed river. The WUA = 183,040 m2 is the peak value of the curve at 60 m3/s, and it is determined to be the appropriate ecological flow for the reconstructed river. Within the flow discharge range of 240–1,040 m3/s, the shrinking trend of WUA in the reconstructed river is similar to the original river. In the flow discharge range of 1,040–1,900 m3/s, the WUA in the original river will be stable at 60,000 m2, but the WUA in the reconstructed river still decreases. Combined with the curve change trend of the juvenile fish stage, it shows that the biological habitat adjustment ability of the reconstructed river is weaker than that of the original river.

Migration

The target fish migrate in groups before and after the spawning period. Figure 10(c) shows that the changing trend of the Q-WUA curve in the original river and the reconstructed river is basically the same, and 140 m3/s is the obvious turning point of the curve, which was determined as the minimum ecological flow of the river. The maximum point 240 m3/s of WUA is determined to be the appropriate ecological flow; at this time, the WUA for the original river is 583,300 m2, and that of the reconstructed river is 725,621 m2, indicating that the reconstructed river will not destroy the target fish's migratory passage.

Spawning

Within the flow range shown in Figure 10(d), the rapid growth point and peak flow discharge corresponding to the Q-WUA curve of the original and reconstructed rivers are the same, being 30 m3/s and 140 m3/s, respectively. The peak value of WUA in the reconstructed river is 87,237 m2, which is 581% higher than that in the original river, and the quality of the target fish spawning habitat in the river is greatly improved. The WUA is the minimum when the flow discharge is 840 m3/s, and it increases slightly after this flow discharge. Under a flood discharge of 3,760 m3/s, the WUA of the original river reaches the maximum, which is 81,629 m2, but the WUA in the reconstructed river is only 30,084 m2 at this time, and the WUA is a theoretical value considering only the factors of water depth and velocity. In practice, the flood environment has high sediment content, poor water quality and unstable flow environment, all of which are unsuitable for fish survival. To sum up, the minimum and appropriate ecological flows of the study river are 30 and 140 m3/s, respectively.

Discussion

To meet the navigation conditions of the curved reach of the lower reaches of the Shaying River, the curving river was cut off, gravel mound groups were arranged in the cut-off channel, and the influence of the waterway regulation project on the living environment of fish has been discussed. Prior studies have noted the importance of flow to habitat (Nan et al. 2022); the annual precipitation distribution in the Shaying River basin is extremely uneven, the flood flow is much higher than the annual average flow, water reduction and dehydration occur frequently in the middle and lower reaches of the river in the dry season, which seriously destroys the habitats of aquatic organisms and blocks the migratory passage of fish resulting in a decline in the population of aquatic organisms (Ge et al. 2019). In this study, the ecological flow of the reconstructed river and the original river were estimated by the biological habitat method, the minimum and appropriate ecological flows of the target fish in the juvenile, adult, migration, and spawning stages were obtained respectively, and the results were summarized in Table 4. Li et al. (2012) used the 7Q10 method to estimate the minimum ecological flow of the Shaying River as 30.86 m3/s, which is in the ‘middle’ ecological water demand grade of the Tennant method and is close to the minimum ecological flow for adult fish and the spawning stage obtained in this study. Hou et al. (2020b) pointed out that even if the minimum ecological flow is met, biodiversity and ecological integrity are reduced. The appropriate ecological flow obtained by the habitat method takes into account the needs of organisms, so it can better ensure the health of a river ecosystem.

Table 4

Estimation table of the ecological flow of target fish in the river channel (unit: m3/s)

Growth stageOriginal river
Reconstructed river
MinimumAppropriateMinimumAppropriate
Juvenile (July to August) 20 40 40 60 
Adult (July to November) 30 40 30 60 
Migration 140 240 140 240 
Spawning (March to June) 30 140 30 140 
Growth stageOriginal river
Reconstructed river
MinimumAppropriateMinimumAppropriate
Juvenile (July to August) 20 40 40 60 
Adult (July to November) 30 40 30 60 
Migration 140 240 140 240 
Spawning (March to June) 30 140 30 140 

Habitat change under appropriate flow

Figure 11 shows the distribution of suitable habitats of target fish in the original river and reconstructed river under the appropriately ecological flow. In the juvenile, adult, and migration stages of the target fish, continuous habitats are formed in the lower reaches of the gravel mound groups in the reconstructed river, and the HSI is obviously superior to the convex bank of the river channel, in which the WUA of juvenile fish and migration habitats increased by about 168 and 24% respectively compared with the original river; the suitability and continuity of the adult fish habitat in the lower reaches of the gravel mound groups are improved compared to those without arrangement. The spawning stage requires higher water depth and velocity conditions in the river. Under the long-term sluice control and artificial transformation of the Shaying River, the mainstream is close to the artificial canalized channel, the flow pattern in the channel is single, and most of the water is in a slow water or still water state (Zhang 2017). At this time, the large-scale slow water environment in the channel can no longer meet the spawning demand of the target fish, the spawning WUA is only 12,817 m2 and extremely discontinuous, while the spawning stage is the most fragile and important stage in the fish life cycle (King et al. 2016). As a result, the individual populations of the four major Chinese carps in the Shaying river have greatly reduced. As shown in Figure 11(d), with the increase of velocity in the reconstructed river, the area of spawning habitat is increased and the continuity is improved in the channel, and the WUA is increased by 581% of that of the original river. Although the HSI value in the low-speed flow area behind the gravel mound groups is small, the high-velocity waters formed on the left side of the gravel mound groups has increased, and the habitat suitability has improved. To sum up, it is concluded that the arrangement of gravel mound groups enriches the habitat structure in the curving cut-off channel, which is suitable for the habitat environment of the target fish during the whole life cycle.
Figure 11

Habitat distribution maps of different rivers; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Figure 11

Habitat distribution maps of different rivers; (a) Juvenile, (b) Adult, (c) Migration, (d) Spawning.

Close modal

Deficiencies and prospects

In the study, only two hydraulic factors, flow velocity and water depth, are considered, while the factors affecting fish survival are various. Studies show that temperature and the life behavior of fish are inextricably linked (King et al. 2016), and chemical factors such as ammonia nitrogen concentration and dissolved oxygen also have certain effects on fish habitat quality (Boets et al. 2021). In future studies, combining other habitat influencing factors should be considered to have a more scientific evaluation of habitat.

Based on the fish habitat changes in the river channel before and after waterway regulation, the minimum and appropriate ecological flows of grass, silver, and bighead carp in the study reach were estimated. However, the study area river length is about 19% of the length from Yingshang Sluice to Mohekou, and it is not certain that the ecological flows obtained in this study are optimal in other sections of the river. In the next phase of the study, the research area should be expanded, combining the life history of fish with the hydrological conditions of the basin, to estimate the ecological flow values for all life stages of fish in the entire area, and make reasonable recommendations for upstream reservoir dispatching.

  • (1)

    The river channel can meet navigation requirements after the curving cutoff, the flow velocity increases, and the viscous force between the flow layers decreases. The water levels of the inlet and outlet decrease by 0.64 and 2.23 m, respectively, under a 20-year return period flood discharge, and the flood discharge capacity of the channel is improved. However, the flow pattern in the curving cutoff channel is single, and the habitat type of the target fish is mainly deep pools.

  • (2)

    The minimum and appropriate ecological flows of the original and reconstructed rivers were estimated, and it was established that the minimum ecological flow for the grass, silver, and bighead carp spawning stage (March – June) was 30 m3/s and that the appropriate ecological flow was 140 m3/s. The minimum and appropriate ecological flows for the migration stage were 140 and 240 m3/s, respectively. In the growth stage (July – November), during which there are juvenile and adult fish stages, the appropriate ecological flow of the original and reconstruction rivers were 40 and 60 m3/s, respectively. From July to August, the ecological flow regulation should be dominated by the needs of the juvenile fish. Furthermore, the minimum ecological flow of the original and reconstructed rivers were 20 and 40 m3/s, respectively. From September to November, adult fish dominate, and the minimum ecological flow of the river was 30 m3/s.

  • (3)

    The natural reconstruction of the curving cutoff river is conducted using an arrangement of gravel mound groups. Under the condition of appropriate flow, the peak area of juvenile, migratory, and spawning habitats of the target fish in the reconstructed river increased by 168, 24, and 581%, respectively, compared to the original river channel. Meanwhile, the habitat area of adult fish which prefer a slow-water environment decreased by 62%, and the habitat environment was more balanced.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict of interest.

Alabyan
A. M.
&
Lebedeva
S. V.
2018
Flow dynamics in large tidal delta of the Northern Dvina River: 2D simulation
.
J. Hydroinf.
20
(
4
),
798
814
.
doi:10.2166/hydro.2018.051
.
Bai
Y.
,
Lu
B.
,
Luo
J.
,
Liang
Y.
,
Liu
D.
,
Chen
Q.
&
Shi
X.
2013
Induction velocity of juvenile grass carp, silver carp, and darkbarbel catfish
.
Chin. J. Ecol.
32
(
08
),
2085
2089
.
doi:10.13292/j.1000-4890.201804.027
.
Boets
P.
,
Dillen
A.
,
Mertens
J.
,
Vervaeke
B.
,
Van Thuyne
G.
,
Breine
J.
,
Goethals
P.
&
Poelman
E.
2021
Do investments in water quality and habitat restoration programs pay off? An analysis of the chemical and biological water quality of a lowland stream in the Zwalm River basin (Belgium)
.
Environ. Sci. Policy
124
,
115
124
.
doi:10.1016/j.envsci.2021.06.017
.
Chang
L.
,
Xu
B.
,
Zhang
P.
&
Tang
W.
2019
The influence of spur dikes for deep-water channel regulation on fish habitat
.
J. Hydraul. Eng.
50
(
09
),
1086
1094
.
doi:10.13243/j.cnki.slxb.20190207
.
Chen
A.
,
Sui
X.
,
Liao
W.
&
Chen
K.
2016
Review study on instream ecological base flow in China
.
J. China Inst. Water Resou. Hydropower. Res.
14
(
06
),
401
411
.
doi:10.13244/j.cnki.jiwhr.2016.06.001
.
Chen
M.
,
Feng
H.
,
Wang
L.
&
Chen
Q.
2007
Calculation methods for appropriate ecological flow
.
Adv. Water. Sci.
5
(
05
),
745
750
.
10.3321/j.issn:1001-6791.2007.05.018
.
Dang
L.
,
Ma
C.
&
Lian
J.
2018
Effects of reservoir regulation on downstream fish habitat suitability
.
J. Tianjin Univ. (Sci. Technol)
51
(
6
),
566
574
.
doi:10.11784/tdxbz201705046
.
Fu
Y.
,
Leng
J.
,
Zhao
J.
&
Wu
W.
2021
Quantitative calculation and optimized applications of ecological flow based on nature-based solutions
.
J. Hydrol.
598
,
126216
.
doi:10.1016/j.jhydrol.2021.126216
.
Ge
J.
,
Peng
W.
,
Zhang
W.
&
Qu
X.
2019
Comparison of several hydrological methods to determine the ecological flow for instream river
.
South-to-North Water Transfers and Water Sci. Technol.
17
(
02
),
75
80
.
doi:10.13476/j.cnki.nsbdqk.2019.0036
.
Hou
J.
,
Huang
Y.
,
Miao
L.
,
Shao
G.
&
Wu
M.
2020a
Calculating environmental flows in middle reach of Yarlung Tsangpo River based on fish habitat requirements
.
Water. Resour. Protect.
36
(
04
),
8
12 + 20
.
Hou
J.
,
Pei
J.
,
Huang
Y.
,
Liang
X.
,
Gong
X.
,
Ding
W.
&
Wu
M.
2020b
Study on environmental flow of lower reach of Xi County hub project based on fish demand
.
Water. Resour. Protect.
36
(
02
),
8
12 + 20
.
doi:10.3880/j.issn.1004-6933.2020.02.002
.
Hung
H. J.
,
Lo
W. C.
,
Chen
C. N.
&
Tsai
C. H.
2022
Fish habitat area and habitat transition in a river under ordinary and flood flow
.
Ecol. Eng.
179
,
106606
.
doi:10.1016/j.ecoleng.2022.106606
.
Jowett
I. G.
1993
A method for objectively identifying pool, run, and riffle habitats from physical measurements
.
New. Zeal. J. Mar. Fresh.
27
(
2
),
241
248
.
doi:10.1080/00288330.1993.9516563
.
King
A. J.
,
Gwinn
D. C.
,
Tonkin
Z.
,
Mahoney
J.
,
Raymond
S.
&
Beesley
L.
2016
Using abiotic drivers of fish spawning to inform environmental flow management
.
J. Appl. Ecol.
53
(
1
),
34
43
.
doi:10.1111/1365-2664.12542
.
Lee
J. H.
,
Kil
J. T.
&
Jeong
S.
2010
Evaluation of physical fish habitat quality enhancement designs in urban streams using a 2D hydrodynamic model
.
Ecol. Eng.
36
(
10
),
1251
1259
.
doi:10.1016/j.ecoleng.2010.05.004
.
Li
Z.
&
Wang
C.
1983
The distribution of fishes in Henan
.
J. Henan Norm. Univ. (Nat. Sci. Technol)
1983
(
04
),
65
75
.
doi:10.16366/j.cnki.1000-2367.1983.04.010
.
Li
Y.
,
Xu
D.
,
Liu
T.
&
An
Y.
2012
Estimation of minimum ecological water demand in Shaying River
.
Eng. Constr.
26
(
02
),
159
161 + 167
.
Liu
D.
,
Xie
J.
,
Yang
X.
,
Wei
N.
,
Gao
S.
&
Yang
Z.
2022
Adaptive ecological baseflow computation based on the ecological requirements: a case study in Weihe Basin, China
.
Water Supply
22
(
9
),
7111
7124
.
doi:10.2166/ws.2022.316
.
Ma
B.
,
Dong
F.
,
Peng
W. Q.
,
Liu
X. B.
,
Huang
A. P.
,
Zhang
X. H.
&
Liu
J. Z.
2020
Evaluation of impact of spur dike designs on enhancement of aquatic habitats in urban streams using 2D habitat numerical simulations
.
Global Ecol. Conserv.
24
,
e01288
.
doi:10.1016/j.gecco.2020.e01288
.
Ministry of Transport
2014
Navigation Standard of Inland Waterway. GB 50139-2014. Ministry of Transport of the People’s Republic of China, Beijing, China
.
Mozzaquattro
L. B.
,
Dala-Corte
R. B.
,
Becker
F. G.
&
Melo
A. S.
2020
Effects of spatial distance, physical barriers, and habitat on a stream fish metacommunity
.
Hydrobiologia
847
,
3039
3054
.
doi:10.1007/s10750-020-04309-8
.
Nakamura
F.
,
Ishiyama
N.
,
Sueyoshi
M.
,
Negishi
J. N.
&
Akasaka
T.
2014
The significance of meander restoration for the hydrogeomorphology and recovery of wetland organisms in the Kushiro River, a lowland river in Japan
.
Restor. Ecol.
22
(
4
),
544
554
.
doi:10.1111/rec.12101
.
Nan
J.
,
Liu
Y.
,
Chen
Y.
&
Zhang
S.
2022
Reconstruction of biological habitat and ecological flow estimation in curving cut-off river channel
.
Water Resour. Protect.
38
(
03
),
189
197
.
doi:10.3880/j.issn.1004 6933.2022.03.026
.
Nan
J.
,
Zhang
S.
,
Yuan
F.
&
Yang
X.
2021
Naturalization of biological habitat based on arrangement of gravel groups
.
Adv. Water Sci.
32
(
04
),
608
617
.
doi:10.14042/j.cnki.32.1309.2021.04.012
.
Odum
H. T.
&
Odum
B.
2003
Concepts and methods of ecological engineering
.
Ecol. Eng.
20
(
5
),
339
361
.
doi:10.1016/j.ecoleng.2003.08.008
.
Park
S. Y.
,
Sur
C.
,
Lee
J. H.
&
Kim
J. S.
2020
Ecological drought monitoring through fish habitat-based flow assessment in the gam river basin of Korea
.
Ecol. Indic.
109
(
Feb.
),
105830.1
105830.10
.
doi:10.1016/j.ecolind.2019.105830
.
Seifer
A.
1983
Naturnäherer Wasserbau
.
Deutsche Wasser Wirtschaft.
33
(
12
),
361
366
. (In German).
Shu
W.
&
Wei
C.
2015
Analysis on species composition and characteristics of fish in Shaying River
.
Harnessing the Huaihe River
2015
(
01
),
27
28
.
doi:10.16747/j.cnki.cn61-1109/tv.2022.07.008
.
Shukla
S. S.
&
Mishra
M.
2019
Tracing of palaeochannels of Baku-lahi river system in Uttar Pradesh, India
.
Arabian J. Geosci.
12
(
9
),
1
9
.
doi:10.1007/s12517-019-4429-6
.
Yang
L.
,
Hou
J.
,
Cheng
L.
,
Wang
P.
,
Pan
Z.
,
Wang
T.
,
Ma
Y.
,
Gao
X.
&
Liu
N.
2021
Application of habitat suitability model coupling with high-precision hydrodynamic processes
.
Ecol. Modell.
462
,
109792
.
doi:10.1016/j.ecolmodel.2021.109792
.
Ye
H.
2007
Model Research of Flow Pattern at Lock Approach Based on Yingshang Lock of Fuyang, Anhui. Master's Degree
.
Zhang
Y.
2017
Ecohydraulics Study on Spawning Behavior of Grass Carp of the Four Major Chinese Carps in the Middle Reach of the Yangtze River. Master's Degree. Chongqing Jiaotong University, Chongqing, China
.
Zhou
X.
,
Xu
M.
,
Lei
F.
,
Zhang
J.
,
Wang
Z.
&
Luo
Y.
2022
Responses of macroinvertebrate assemblages to flow in the Qinghai-Tibet Plateau: establishment and application of a multi-metric habitat suitability model
.
Water Resour. Res
e2021WR030909
.
doi:10.1029/2021WR030909
.
Zhou
Z.
,
Deng
Y.
,
Li
Y.
&
An
R.
2019
The ecological water demand of Schizothorax in Tibet based on habitat area and connectivity
.
Int. J. Environ. Res. Public Health
16
(
17
),
3045
.
doi:10.3390/ijerph16173045
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).