Analysis of an Ecological Flow Regime during the Ctenopharyngodon Idella Spawning Period Based on Reservoir Operations

: The study of ﬁsh habitats is important for us to better understand the impact of reservoir construction on river ecosystems. Many habitat models have been developed in the past few decades. In this study, a fuzzy logic-based habitat model, which couples fuzzy inference system, two-dimensional laterally averaged hydrodynamic model, and two-dimensional shallow water hydrodynamic model, is proposed to identify the baseline condition of suitable habitat for ﬁsh spawning activities. The proposed model considers the reservoir and the downstream river channel, and explores the comprehensive e ﬀ ects of water temperature, velocity, and water depth on habitat suitability. A real-world case that considers the Ctenopharyngodon idella in the Xuanwei Reservoir of Qingshui River is studied to investigate the e ﬀ ect of in- and outﬂow of reservoir on ﬁsh habitat and the best integrative management measure of the model. There were 64 simulations with di ﬀ erent reservoir in- and outﬂows employed to calculate the weighted usable area and hydraulic habitat suitability. The experimental results show that the ecological ﬂow for Ctenopharyngodon idella spawning can satisfy the basic demand when the reservoir inﬂow is greater than 60 m 3 / s and the reservoir outﬂow is greater than 100 m 3 / s. The habitat ecological suitability is the best when the reservoir outﬂow is 120 m 3 / s. A more reasonable and reliable ecological ﬂow range can be obtained based on the habitat model in this paper, which provides the best scenario for water resources planning and management in the Qingshui River Basin.


Introduction
Reservoirs play key roles in the social and economic development as a water conservancy project that prevents or controls flood and regulates water flow [1], however, the construction of a reservoir changes natural flow regimes of a river and has an impact on the river. These effects are mainly manifested as follows: (1) blocking the migration path of fish to make it impossible to complete its life history process, resulting in a decline in genetic diversity of fish populations [2]; (2) affecting sediment transport, resulting in sedimentation at the bottom of the reservoir and river blockage [3]; (3) affecting nitrogen, phosphorus, and other substances migration, thereby inducing eutrophication [4];

1.
We consider the thermal stratification of reservoirs and analyze the important influence of the hysteresis of the water temperature on the spawning of Ctenopharyngodon idella in habitats.

2.
A fuzzy logic-based habitat suitability model is established to evaluate the influence of hydraulic factors such as water temperature, velocity, and water depth on the suitability of a Ctenopharyngodon idella habitat.

3.
By simulating a variety of different reservoir inflow and outflow conditions, the functional relationship between flow and habitat suitability is established, which provides the best management plan for water resources planning and construction.
The remaining parts of this paper are arranged as follows: In Section 2, we introduce the study area and target ecological species. In Section 3, we describe the construction of hydraulic habitat suitability (HHS) for target ecological species, a method for evaluating the weighted usable area (WUA) based on hydraulic simulation and the fuzzy logic-based habitat model. The results analysis and discussion are presented in Section 4. Finally, in Section 5, we state our conclusion for this study.

Study Area
The Qingshui River, located in the upper reaches of the Yuanjiang River, belongs to the Yangtze River system in China. It is located in the southeast of Guizhou Province and southwest of Hunan Province, with a drainage area of 15,000 km 2 . Hills and mountainous regions compose most of the Qingshui River Basin. The basin has abundant rainfall, warm climate, rich oil, and diverse types of aquatic habitats. The Xuanwei reservoir is located in the upper reaches of the mainstream of the Qingshui River, with a catchment area of 1337 km 2 . The design parameters of the reservoir include a normal water level of 681 m, a dead water level of 657 m, a total storage capacity of 139.80 million m 3 , and an average annual storage flow of 903.39 million m 3 . Among the fish in the study area, the largest number of cypriniformes (34 species) accounted for 65.4% of the total number of fish species by the scene investigation and sampling in the field, followed by 17.3% perciformes (9 species), 15.4% siluriformes (8 species), and 1.9% synbranchiformes (1 species). There is a spawning ground from the reservoir to the downstream river channel with a length of 3 km. The main fish species in the spawning ground is Ctenopharyngodon idella which produces drifting eggs. The bottom of the riverbed is an open rocky beach formed by pebbles and river sand. By considering the 3 km river channel downstream of the reservoir as a study case, the influence of hydraulic factors on fish spawning is simulated and analyzed. The location of the study area is shown in Figure 1.

Ctenopharyngodon idella (C. idella)
Ctenopharyngodon idella (C. idella) is a freshwater species native mainly to China. It plays an important role in aquaculture, being the second highest in world fish production. C. idella, as a food fish or as a biological controller of aquatic vegetation, has been introduced into more than 100 countries in various continental regions, including East Asia, North America, and Europe [36]. The natural reproduction of C. idella has certain demands for water temperature [37], landform, hydrology, and hydrodynamic conditions, especially related to flow and water temperature. C. idella migrates from the middle and lower reaches of the Qingshui River to the upstream for spawning and breeding during the summer flood season. In the autumn, it swims to the deep waters of the lower reaches of the Qingshui River for overwintering. Its spawning occurs during the flood period from mid-May to mid-July each year when the water temperature reaches 18 °C. According to the survey data of C. idella in the Qingshui River, it is most suitable for adult C. idella spawning when the water temperature is between 21 °C and 24.5 °C, the velocity is between 0.3 m/s and 0.9 m/s, and the water depth is between 3 m and 7.5 m. This migratory fish is the target species of reservoir ecological scheduling because the construction of the dam greatly threatens its survival. The ecological scheduling of such reservoirs generally involves releasing a suitable flow for fish migration during the migratory period of the fish and combining certain fish facilities to help the fish reach the spawning ground [38]. Therefore, the natural reproduction of C. idella is required to incorporate the water environment requirements of the spawning habitat into the reservoir scheduling objective. The construction of multi-objective ecological scheduling is one of the ways to effectively protect the river's important economic fish species resources [39]. However, a key step in the implemention of a reservoir ecological scheduling scenario that takes into account the conservation of rare aquatic organisms is to determine the ecological water demand and to clarify the relationship between reservoir scheduling and the suitability of water environment for protecting aquatic habitats.
Practical topographic data of the Xuanwei reservoir and the downstream river channel was acquired from the Hydraulic Exploration Design Institute of the Guizhou Province. Meteorological data such as temperature, relative humidity, wind speed, and cloudage come from the Duyun meteorological station, since 1970. Hydraulic data such as reservoir inflow and water level are provided by the Xiasi hydrological station. The C. idella information of the Qingshui River comes from the sampling data of the past five years in the field. According to the data of the Duyun hydrological station in the study area, the data of reservoir inflow and water temperature of typical years are obtained as follows: wet year (75% guarantee rate of runoff), normal year (50% guarantee rate of runoff), and dry year (25% guarantee rate of runoff). According to the inflow data of the reservoir, the annual large-scale flow is mainly concentrated from mid-May to early-July. The highest water temperature in the year occurred in August, and the lowest water temperature occurred in February. According to the survey, the spawning period of C. idella in the study area is mostly

Ctenopharyngodon idella (C. idella)
Ctenopharyngodon idella (C. idella) is a freshwater species native mainly to China. It plays an important role in aquaculture, being the second highest in world fish production. C. idella, as a food fish or as a biological controller of aquatic vegetation, has been introduced into more than 100 countries in various continental regions, including East Asia, North America, and Europe [36]. The natural reproduction of C. idella has certain demands for water temperature [37], landform, hydrology, and hydrodynamic conditions, especially related to flow and water temperature. C. idella migrates from the middle and lower reaches of the Qingshui River to the upstream for spawning and breeding during the summer flood season. In the autumn, it swims to the deep waters of the lower reaches of the Qingshui River for overwintering. Its spawning occurs during the flood period from mid-May to mid-July each year when the water temperature reaches 18 • C. According to the survey data of C. idella in the Qingshui River, it is most suitable for adult C. idella spawning when the water temperature is between 21 • C and 24.5 • C, the velocity is between 0.3 m/s and 0.9 m/s, and the water depth is between 3 m and 7.5 m. This migratory fish is the target species of reservoir ecological scheduling because the construction of the dam greatly threatens its survival. The ecological scheduling of such reservoirs generally involves releasing a suitable flow for fish migration during the migratory period of the fish and combining certain fish facilities to help the fish reach the spawning ground [38]. Therefore, the natural reproduction of C. idella is required to incorporate the water environment requirements of the spawning habitat into the reservoir scheduling objective. The construction of multi-objective ecological scheduling is one of the ways to effectively protect the river's important economic fish species resources [39]. However, a key step in the implemention of a reservoir ecological scheduling scenario that takes into account the conservation of rare aquatic organisms is to determine the ecological water demand and to clarify the relationship between reservoir scheduling and the suitability of water environment for protecting aquatic habitats.
Practical topographic data of the Xuanwei reservoir and the downstream river channel was acquired from the Hydraulic Exploration Design Institute of the Guizhou Province. Meteorological data such as temperature, relative humidity, wind speed, and cloudage come from the Duyun meteorological station, since 1970. Hydraulic data such as reservoir inflow and water level are provided by the Xiasi hydrological station. The C. idella information of the Qingshui River comes from the sampling data of the past five years in the field. According to the data of the Duyun hydrological station in the study area, the data of reservoir inflow and water temperature of typical years are obtained as follows: wet year (75% guarantee rate of runoff), normal year (50% guarantee rate of runoff), and dry year (25% guarantee rate of runoff). According to the inflow data of the reservoir, the annual large-scale flow is mainly concentrated from mid-May to early-July. The highest water temperature in the year occurred in August, and the lowest water temperature occurred in February. According to the survey, Water 2019, 11, 2034 5 of 23 the spawning period of C. idella in the study area is mostly concentrated from mid-May to mid-July. The suitable flow range and suitable water temperature range determine the spawning period of the fish. The annual reservoir inflow rate and water temperature of the study area is shown in Figure 2.
concentrated from mid-May to mid-July. The suitable flow range and suitable water temperature range determine the spawning period of the fish. The annual reservoir inflow rate and water temperature of the study area is shown in Figure 2.
C. idella has an important economic value and is the target species in the study area for this study. On the basis of constructing the hydrodynamic mathematical model and habitat suitability evaluation model of the spawning ground in the downstream of the Xuanwei Reservoir, the correlation between the key indicators of reservoir scheduling and the suitability of C. idella spawning habitat have been studied. The research results provide the basis for constructing reservoir ecological scheduling of C. idella protection.

Model Framework
In this study, the C. idella in the study area was used as the target species to study the flow and water temperature effects of C. idella during spawning. This study specifically analyzed the following three main factors affecting the spawning of C. idella: velocity, water depth, and water temperature. Two hydrodynamic models were used to simulate the hydraulic conditions of reservoir and the downstream habitats of reservoir. MIKE21 simulated the spatial distribution of velocity and water depth under different flows in the habitat. CE-QUAL-W2 simulated the water temperature stratification of reservoirs under different inflows and outflows. The simulation results of velocity, water depth, and water temperature were used to construct a fuzzy model [31,40] to comprehensively analyze the effects of the three factors on C. idella spawning. The simulation results of the fuzzy logicbased habitat suitability model were analyzed and evaluated by WUA and HHS, and the best flow management scenario in the study area was finally obtained. The overall framework of the habitat model is shown in Figure 3. C. idella has an important economic value and is the target species in the study area for this study. On the basis of constructing the hydrodynamic mathematical model and habitat suitability evaluation model of the spawning ground in the downstream of the Xuanwei Reservoir, the correlation between the key indicators of reservoir scheduling and the suitability of C. idella spawning habitat have been studied. The research results provide the basis for constructing reservoir ecological scheduling of C. idella protection.

Model Framework
In this study, the C. idella in the study area was used as the target species to study the flow and water temperature effects of C. idella during spawning. This study specifically analyzed the following three main factors affecting the spawning of C. idella: velocity, water depth, and water temperature. Two hydrodynamic models were used to simulate the hydraulic conditions of reservoir and the downstream habitats of reservoir. MIKE21 simulated the spatial distribution of velocity and water depth under different flows in the habitat. CE-QUAL-W2 simulated the water temperature stratification of reservoirs under different inflows and outflows. The simulation results of velocity, water depth, and water temperature were used to construct a fuzzy model [31,40] to comprehensively analyze the effects of the three factors on C. idella spawning. The simulation results of the fuzzy logic-based habitat suitability model were analyzed and evaluated by WUA and HHS, and the best flow management scenario in the study area was finally obtained. The overall framework of the habitat model is shown in Figure 3.

Hydrodynamic Simulation
The hydrodynamic modeling was divided into two parts. The first part of the reservoir was simulated by CE-QUAL-W2 model, which mainly explored the influence of water temperature change caused by reservoir thermal stratification on downstream habitat. The second part of the downstream channel was simulated by MIKE21 model, which mainly explored the velocity and water depth distribution of the downstream habitat.

CE-QUAL-W2 Model
The CE-QUAL-W2 model is a two-dimensional (longitudinal/vertical) laterally averaged hydrodynamic model for surface water systems [41]. The model is based on the assumption that lateral homogeneity and the flow in a reservoir with a distinct flow direction can be computed from the laterally integrated Navier-Stokes equation [42]. This is particularly suited for relatively long and narrow water bodies exhibiting longitudinal and vertical water temperature gradients. The W2 model was successfully applied to stratified water systems, including reservoirs, lakes, and estuaries.
On the basis of the Boussinesq assumption, the model integrates the homogeneous continuity equation, the momentum equation, and the caloric equation of the three-dimensional turbulent buoyant flow to obtain the equations of the two-dimensional laterally averaged temperature model as follows:

Hydrodynamic Simulation
The hydrodynamic modeling was divided into two parts. The first part of the reservoir was simulated by CE-QUAL-W2 model, which mainly explored the influence of water temperature change caused by reservoir thermal stratification on downstream habitat. The second part of the downstream channel was simulated by MIKE21 model, which mainly explored the velocity and water depth distribution of the downstream habitat.

CE-QUAL-W2 Model
The CE-QUAL-W2 model is a two-dimensional (longitudinal/vertical) laterally averaged hydrodynamic model for surface water systems [41]. The model is based on the assumption that lateral homogeneity and the flow in a reservoir with a distinct flow direction can be computed from the laterally integrated Navier-Stokes equation [42]. This is particularly suited for relatively long and narrow water bodies exhibiting longitudinal and vertical water temperature gradients. The W2 model was successfully applied to stratified water systems, including reservoirs, lakes, and estuaries.
On the basis of the Boussinesq assumption, the model integrates the homogeneous continuity equation, the momentum equation, and the caloric equation of the three-dimensional turbulent buoyant flow to obtain the equations of the two-dimensional laterally averaged temperature model as follows: where B is the river width, u is the longitudinal velocity, w is the vertical velocity, v is the coefficient of kinematic viscosity, v t is the turbulent viscosity coefficient, v e = v + v t is the effective viscosity coefficient, T is the water temperature, T a is the reference temperature of water, ∆T = T − T a , ρ is the water density, ρ a is the reference density of water, β is the coefficient of thermal expansion, and σ T is the Prandtl number of temperature. The Xuanwei Reservoir was represented as 63 longitudinal segments with lengths of 250 m and 48 vertical layers with thicknesses of 1 m according to the topographic data recorded in 2018. The CE-QUAL-W2 turbulence closure algorithm was used. There is only one inflow entrance upstream of the mainstream and there are no other inflow entrances from the tributaries. The boundary data was input with a daily resolution, including reservoir inflow and temperature, outflow, air temperature, wind speed and direction, solar radiation, cloud cover, precipitation, and evaporation. Reservoir grid model is shown in Figure 4. where B is the river width, u is the longitudinal velocity, w is the vertical velocity, v is the coefficient of kinematic viscosity, vt is the turbulent viscosity coefficient, ve = v + vt is the effective viscosity coefficient, T is the water temperature, Ta is the reference temperature of water, ΔT = T -Ta, ρ is the water density, ρa is the reference density of water, β is the coefficient of thermal expansion, and σT is the Prandtl number of temperature.
The Xuanwei Reservoir was represented as 63 longitudinal segments with lengths of 250 m and 48 vertical layers with thicknesses of 1 m according to the topographic data recorded in 2018. The CE-QUAL-W2 turbulence closure algorithm was used. There is only one inflow entrance upstream of the mainstream and there are no other inflow entrances from the tributaries. The boundary data was input with a daily resolution, including reservoir inflow and temperature, outflow, air temperature, wind speed and direction, solar radiation, cloud cover, precipitation, and evaporation.

MIKE21 Model
The fish habitat in the downstream of the Xuanwei Reservoir was simulated by MIKE21, and the spatial distribution of velocity and water depth under different outflows was calculated. The hydrodynamic module is based on the solution of two-dimensional shallow water equations obtained from the Navier-Stokes equations, which are integrated over the water depth including the depth and horizontal momentum equation of the flow just as follows [43]: where t indicates the time; x and y are Cartesian coordinates; T, U, V, and S indicate vectors of the conserved variables, fluxes in the x-and y-direction, and source terms, respectively.
where h is the water depth; g = 9.81 m/s 2 is the gravity acceleration; u and v is the depth-averaged velocity components in the x-direction and y-direction, respectively; S0x and S0y are the bed slopes in the x-direction and y-direction, respectively; Sfx and Sfy are the friction slopes in the x-direction and y-direction, respectively.
Tests of numerical experiments and results of the practical applications showed that this model can simulate the water flow movement over complex terrain accurately [44]. The MIKE21 model was used to divide the downstream habitat of the Xuanwei Reservoir into 796 triangular grids, with a maximum area of 523.49 m 2 and a minimum area of 130.33 m 2 . The measured outflow time curve of the reservoir was taken as the upper boundary of the model. The measured water level time curve at

MIKE21 Model
The fish habitat in the downstream of the Xuanwei Reservoir was simulated by MIKE21, and the spatial distribution of velocity and water depth under different outflows was calculated. The hydrodynamic module is based on the solution of two-dimensional shallow water equations obtained from the Navier-Stokes equations, which are integrated over the water depth including the depth and horizontal momentum equation of the flow just as follows [43]: where t indicates the time; x and y are Cartesian coordinates; T, U, V, and S indicate vectors of the conserved variables, fluxes in the x-and y-direction, and source terms, respectively.
where h is the water depth; g = 9.81 m/s 2 is the gravity acceleration; u and v is the depth-averaged velocity components in the x-direction and y-direction, respectively; S 0x and S 0y are the bed slopes in the x-direction and y-direction, respectively; S fx and S fy are the friction slopes in the x-direction and y-direction, respectively.
Tests of numerical experiments and results of the practical applications showed that this model can simulate the water flow movement over complex terrain accurately [44]. The MIKE21 model was used to divide the downstream habitat of the Xuanwei Reservoir into 796 triangular grids, with a maximum area of 523.49 m 2 and a minimum area of 130.33 m 2 . The measured outflow time curve of the reservoir was taken as the upper boundary of the model. The measured water level time curve at 3 km from the reservoir was taken as the lower boundary of the model. The model simulated the velocity and water depth distribution in the habitat under different outflows. Habitat grid model is shown in Figure 5.

Fuzzy Logic-Based Habitat Modeling
Fuzzy logic can copy the judgment of the uncertainty concept and the way of reasoning thinking of the human brain by using fuzzy sets and fuzzy rules [45,46]. Meanwhile, the model considers the inherent uncertainty of ecological variables and enables expressing nonlinear relations between ecological variables in a transparent way [47,48]. Fuzzy species distribution models transform fuzzy modeling into a technique suitable for species distribution modeling, which can reflect the specific characteristics of ecological problems. This study was based on fuzzy inference for habitat suitability (HS) simulation. The inference process consisted of three steps as follows: (1) fuzzification, in which fuzzifier is used to convert crisp input to fuzzy set using membership functions; (2) fuzzy inference, in which a fuzzy inference system is applied to convert fuzzy input sets to fuzzy output sets; and (3) defuzzification, in which fuzzy output sets are converted to crisp output by a defuzzifier. The modeling for the fuzzy inference process is shown in Figure 6.

Fuzzy Logic-Based Habitat Modeling
Fuzzy logic can copy the judgment of the uncertainty concept and the way of reasoning thinking of the human brain by using fuzzy sets and fuzzy rules [45,46]. Meanwhile, the model considers the inherent uncertainty of ecological variables and enables expressing nonlinear relations between ecological variables in a transparent way [47,48]. Fuzzy species distribution models transform fuzzy modeling into a technique suitable for species distribution modeling, which can reflect the specific characteristics of ecological problems. This study was based on fuzzy inference for habitat suitability (HS) simulation. The inference process consisted of three steps as follows: (1) fuzzification, in which fuzzifier is used to convert crisp input to fuzzy set using membership functions; (2) fuzzy inference, in which a fuzzy inference system is applied to convert fuzzy input sets to fuzzy output sets; and (3) defuzzification, in which fuzzy output sets are converted to crisp output by a defuzzifier.

Fuzzification and Defuzzification
(1) Fuzzification In this study, for the habitat modeling, the membership functions of trapezoidal and triangular shapes are used to convert crisp input values of velocity, water depth, and water temperature lag time to fuzzy sets expressed by one or more linguistic values with membership degrees. Fuzzy sets are usually expressed by linguistic description such as "very low", "low", "high", and "very high". In this study, for the habitat modeling, the membership functions of trapezoidal and triangular shapes are used to convert crisp input values of velocity, water depth, and water temperature lag time to fuzzy sets expressed by one or more linguistic values with membership degrees. Fuzzy sets are usually expressed by linguistic description such as "very low", "low", "high", and "very high". The membership functions of the fuzzy sets have overlapping boundaries and the linguistic statement "the water temperature lag time is quite short but tending to be long" can be translated into a time which has a membership degree of 0.8 to the fuzzy set "high" and of 1 to the fuzzy set "very high".
(2) Fuzzy Inference A fuzzy inference system is composed of a fuzzy rule base and an inference machine. The fuzzy rule consists of if-then fuzzy rules, such as "if x IS A AND y IS B, THEN z IS C", by links the input variables to output. In the same way, the fuzzy rules as a total of 64 (velocity with four linguistic values, water depth with four linguistic values, and water temperature lag time with four linguistic values) describing HS for the spawning period of C. idella are created based on literature reviews and expert experience [49][50][51]. The fuzzy rules for the spawning of C. idella is shown in Table 1. The water temperature lag time is used as an evaluation index of the fuzzy model. The water temperature lag time indicates the time required for the outflow water temperature to reach the same water temperature when the reservoir inflow water temperature (21 • C) is suitable for C. Idella to spawn. By using this fuzzy rule base, the Mamdani-Assilian inference was used in the current application. An example of a fuzzy inference process is described in Figure 7. The first line and second line show three inputs and a single output. The third line depicts how three input values are partial members of three linguistic variables for "velocity", "water depth" and "water temperature lag time". This leads three fuzzy rules to be activated. The implication operator determines the partial membership to the output. These are aggregated into an output surface. Note: Fuzzy sets are usually expressed by linguistic description such as "very low", "low", "high" and "very high". VL, L, H, VH indicate "very low", "low", "high" and "very high", respectively.   (3) Defuzzification The fuzzy inference system uses implication and aggregation operators to scale the membership functions of output linguistic variables before performing defuzzification. Afterwards, a numerical output value of the fuzzy inference process is obtained with defuzzification methods, center of gravity. In this study, the center of gravity defuzzification method is used as follows: where CG is the center gravity of the area, x is the value of the linguistic variable, and x max and x min indicate the range of the linguistic variable. The CG defuzzification method effectively calculates the best compromise between multiple output linguistic terms. This crisp output value calculated by this method denotes HS for C. idella, described as the habitat suitability index (HSI).

Criteria for Assessing Habitat Quality
In order to make a better and accurate evaluation for availability and suitability of habitat, some related evaluation methodology such as weighted usable area (WUA), hydraulic habitat suitability (HHS), and habitat suitability maps were calculated in the study. Habitat suitability maps were generated to visualize the habitat suitability within the river under a steady flow condition. The WUA was calculated by integrating habitat quality over the cells of the studied area, which can be widely used to show the habitats available for species. It was assumed that WUA is positively associated with fish biomass in IFIM. The HHS was obtained by dividing the WUA by the total area, which assigns relative values between 0 and 1. It describes the suitability of hydrodynamic variables for the fish species considered. The calculation methods of WUA and HHS are as follows: where A i indicates the i-th grid area, n indicates the total number of grids, and HSI i indicates the i-th habitat suitability index. The best flow management scenario is determined by comparing and analyzing the simulation results under WUA and HHS for each scenario. A larger the WUA value, the larger total habitat weighted usable area for fish spawning. A larger WUA value is more favorable for fish spawning. HHS is a normalized value between 0 and 1. The larger the HHS value, the larger the proportion of habitats area for fish spawning. A larger HHS value is more favorable for fish spawning.

Result of CE-QUAL-W2 Model Simulation
According to the research, the Xuanwei Reservoir has thermal stratification, and therefore it is necessary to analyze the negative effect on the C. idella spawning in downstream of the reservoir caused by the releasing of water with low temperature. In this study, the CE-QUAL-W2 model is used to simulate the time taken to reach an appropriate ecological water temperature under different inflow and outflow conditions. First of all, the initial conditions and boundary conditions of the model are set. The historical runoff data (water level and flow) of the hydrological station and the historical meteorological data (wind speed, wind direction, relative humidity, temperature, and cloudage) of the meteorological station are collected to simulate the daily change process of the water temperature of the Xuanwei Reservoir. The vertical distribution of the water temperature of the reservoir in mid-May is obtained by calculating the simulation results, and the results are shown in Figure 8a. The vertical distribution of water temperature in mid-May was taken as the initial conditions of the model, when the inflow water temperature is 21 • C, and the outflow water temperature is 15.8 • C. According to the sampling survey in the study area, the spawning period of C. idella occurs mostly from mid-May to mid-July of the year, and the optimal water temperature range for spawning is between 21 • C to 24.5 • C. When the inflow water temperature is 21 °C, the time taken for the outflow water temperature to reach 21 °C is shown in Table 2. It can be seen from Table 2 that the larger outflow of the reservoir makes the time to reach the appropriate water temperature shorter for the same reservoir inflow.
Due to the thermal stratification of the reservoir, there appears to exist a time lag phenomenon between inflow water temperature and outflow water temperature. The results of the model simulation show that the larger inflow and outflow of the reservoir will make the water body mix more intensely, resulting in shorter time delay of water temperature. At the same time, the fish reach the appropriate water temperature for spawning in a short time, which is more beneficial to C. idella. Conversely, a long lag time of water temperature delays the spawning of C. idella, which is less favorable for spawning.  Note: / indicates that this scheduling scenario does not meet the requirements. According to this scenario, the water level is higher or lower than the water level limit.

Result of MIKE21 Model Simulation
The MIKE21 model is used to simulate the velocity and water depth distribution of fish habitats under different outflows of the reservoir. The outflow range uses the value simulated by the CE-QUAL-W2 model, which is divided into eight different scenarios. The minimum outflow is 20 m 3 /s and the maximum outflow is 160 m 3 /s with an outflow interval of 20 m 3 /s. The upper boundary of the habitat model uses a flow boundary, and the lower boundary uses a water level boundary. The topographic data is provided by the Hydraulic Exploration Design Institute of the Guizhou Province in China. On the basis of the hydraulic simulation, the spatial distribution of water depth and velocity under eight different modeling scenarios is shown in Figure 9. As can be seen from the figure, the velocity is relatively large near the reservoir location, and the water depth is deeper at the river bend.  Figure 8 and Table 2. Note: / indicates that this scheduling scenario does not meet the requirements. According to this scenario, the water level is higher or lower than the water level limit.
In Figure 8, the vertical distribution of water temperature in the reservoir for one day, three days, and five days is shown by five different scheduling scenarios. The vertical distribution of the reservoir water temperature, respectively, is show in Figure 8b-d, when the reservoir inflow is equal to reservoir outflow. The comparison results show that the large inflow of the reservoir makes the water body mix well, which leads to the decrease of the thermal stratification stability. Therefore, the temperature difference between the inflow and the outflow is relatively small. Similarly, the small inflow of the reservoir makes the water body mix slowly, resulting in enhanced thermal stratification stability. Therefore, the temperature difference between the inflow and the outflow is relatively high. The vertical distributions of the reservoir water temperature are shown in Figure 8c, e, and f, where the reservoir inflow is 100 m 3 /s and reservoir outflows are 100 m 3 /s, 80 m 3 /s, and 120 m 3 /s, respectively. The comparison results show that the larger outflow of the reservoir makes the water body mix more fully for the same reservoir inflow, which leads to the decrease of the thermal stratification stability.
When the inflow water temperature is 21 • C, the time taken for the outflow water temperature to reach 21 • C is shown in Table 2. It can be seen from Table 2 that the larger outflow of the reservoir makes the time to reach the appropriate water temperature shorter for the same reservoir inflow.
Due to the thermal stratification of the reservoir, there appears to exist a time lag phenomenon between inflow water temperature and outflow water temperature. The results of the model simulation show that the larger inflow and outflow of the reservoir will make the water body mix more intensely, resulting in shorter time delay of water temperature. At the same time, the fish reach the appropriate water temperature for spawning in a short time, which is more beneficial to C. idella. Conversely, a long lag time of water temperature delays the spawning of C. idella, which is less favorable for spawning.

Result of MIKE21 Model Simulation
The MIKE21 model is used to simulate the velocity and water depth distribution of fish habitats under different outflows of the reservoir. The outflow range uses the value simulated by the CE-QUAL-W2 model, which is divided into eight different scenarios. The minimum outflow is 20 m 3 /s and the maximum outflow is 160 m 3 /s with an outflow interval of 20 m 3 /s. The upper boundary of the habitat model uses a flow boundary, and the lower boundary uses a water level boundary. The topographic data is provided by the Hydraulic Exploration Design Institute of the Guizhou Province in China. On the basis of the hydraulic simulation, the spatial distribution of water depth and velocity under eight different modeling scenarios is shown in Figure 9. As can be seen from the figure, the velocity is relatively large near the reservoir location, and the water depth is deeper at the river bend.

Suitable Ecological Flow Calculation
The simulation results of the CE-QUAL-W2 model and the MIKE21 model are applied to the suitability evaluation of representative species in the study area. The effects of flow, water depth, and water temperature hysteresis on habitats are reflected in the suitability evaluation. The implementation of the fuzzy inference system represented in the previous section gives the outputs for all possible inputs, as shown in Figure 10. The combined effects of each of the two indicators on habitat suitability are shown in Figure 10, with a total of three sets of analysis results. The results

Suitable Ecological Flow Calculation
The simulation results of the CE-QUAL-W2 model and the MIKE21 model are applied to the suitability evaluation of representative species in the study area. The effects of flow, water depth, and water temperature hysteresis on habitats are reflected in the suitability evaluation. The implementation of the fuzzy inference system represented in the previous section gives the outputs for all possible inputs, as shown in Figure 10. The combined effects of each of the two indicators on habitat suitability are shown in Figure 10, with a total of three sets of analysis results. The results show that when velocity and water depth are excessively large or small, habitat suitability is so poor that the hydraulic condition cannot meet the spawning requirements of C. idella. According to the survey data of C. idella in the Qingshui River (Section 2, 2.2. Ctenopharyngodon idella) it is most suitable for adult C. idella spawning when the water temperature is between 21 • C and 24.5 • C, velocity is between 0.3 m/s and 0.9 m/s, and water depth is between 3 m and 7.5 m. The results derived from the 3D fuzzy surface graph are basically the same as the field measured data, so the fuzzy logic-based habitat model should be reliable and adoptable. On the basis of the hydraulic simulation and the fuzzy logic-based habitat model, the weighted usable area (WUA) and hydraulic habitat suitability (HHS) are evaluated. The specific operational steps are to simulate the distribution of velocity, water depth, and water temperature lag time under different inflows and outflows through two hydraulic models. Then, the fuzzy-based habitat model is used to construct the relationship between hydraulic factors and habitat suitability index. Finally, the WUA and the HHS under different reservoir inflows and reservoir outflows are obtained by fuzzy calculation. This study simulates a combination of eight different reservoir inflows and eight different reservoir outflows, including a total of 64 scenarios. The results of 64 scenarios are analyzed, but 19 of them do not meet the scheduling requirements. According to these scenarios, the water level is higher or lower than the water level limit. The simulation results under 45 scenarios are shown in Figure 11 and Table 3.
The WUA and HHS are calculated from the habitat quality maps by using Equations (8) and (9). As can be seen from Figure 11 and Table 3, the change tendencies of the WUA and HHS are almost On the basis of the hydraulic simulation and the fuzzy logic-based habitat model, the weighted usable area (WUA) and hydraulic habitat suitability (HHS) are evaluated. The specific operational steps are to simulate the distribution of velocity, water depth, and water temperature lag time under different inflows and outflows through two hydraulic models. Then, the fuzzy-based habitat model is used to construct the relationship between hydraulic factors and habitat suitability index. Finally, the WUA and the HHS under different reservoir inflows and reservoir outflows are obtained by fuzzy calculation. This study simulates a combination of eight different reservoir inflows and eight different reservoir outflows, including a total of 64 scenarios. The results of 64 scenarios are analyzed, but 19 of them do not meet the scheduling requirements. According to these scenarios, the water level is higher or lower than the water level limit. The simulation results under 45 scenarios are shown in Figure 11 and Table 3. 3. Determination of the suitable ecological flow range based on 64 scenarios. According to the results of 64 scenarios, it shows that the ecological flow for C. idella spawning can satisfy the basic demand when the reservoir inflow is greater than 60 m 3 /s and the reservoir outflow is greater than 100 m 3 /s. The habitat ecological suitability is the best when the reservoir outflow is 120 m 3 /s. Therefore, compared with historical data, the suitable ecological flow from 100 m 3 /s to 160 m 3 /s for C. idella propagation is reliable, and the ecological flow can meet the spawning requirements of the fish species. Figure 11. Weighted usable area (WUA) curves for C. idella in the study area. Note: / indicates that this scheduling scenario does not meet the requirements. According to this scenario, the water level is higher or lower than the water level limit.

Suitable Ecological Flow Discussion
For further analysis of the simulation results, several scenarios are selected for analysis. The spatial distribution of habitat suitability for several scenarios is shown in Figures 12 and 13. According to the simulation results of velocity and water depth in Figures 12 and 13, the velocity is  Figure 11. Weighted usable area (WUA) curves for C. idella in the study area. Note: / indicates that this scheduling scenario does not meet the requirements. According to this scenario, the water level is higher or lower than the water level limit.
The WUA and HHS are calculated from the habitat quality maps by using Equations (8) and (9). As can be seen from Figure 11 and Table 3, the change tendencies of the WUA and HHS are almost consistent. The suitability of hydraulic habitats under different reservoir inflows and outflows are shown in Table 3. The best flow management scenario is determined by comparing and analyzing the simulation results under the WUA and HHS for each scenario as follows: 1. The effects of different outflow on the WUA and HHS for the same reservoir inflow. It can be seen from the table that an increase of the reservoir outflow will make the HHS value gradually increase for the same reservoir inflow. When the reservoir outflow reaches 120 m 3 /s, the HHS value is the largest. However, the HHS value gradually decreases when the reservoir outflow is greater than 120 m 3 /s. At the same time, it can be seen from the table that when the reservoir outflow reaches 100 m 3 /s or more, the comprehensive index of HHS can generally reach 0.6 or more.
2. The effect of the WUA and HHS when the reservoir inflow is equal to outflow. A total of eight scenarios (20,40,60,80,100,120,140, and 160 m 3 /s) were simulated for the same reservoir inflow and outflow. The habitat ecological suitability is the best when the reservoir outflow is 120 m 3 /s. When the flow is greater than 100 m 3 /s, the HHS index can reach 0.6 or more, which is beneficial to C. idella spawning.
3. Determination of the suitable ecological flow range based on 64 scenarios. According to the results of 64 scenarios, it shows that the ecological flow for C. idella spawning can satisfy the basic demand when the reservoir inflow is greater than 60 m 3 /s and the reservoir outflow is greater than 100 m 3 /s. The habitat ecological suitability is the best when the reservoir outflow is 120 m 3 /s. Therefore, compared with historical data, the suitable ecological flow from 100 m 3 /s to 160 m 3 /s for C. idella propagation is reliable, and the ecological flow can meet the spawning requirements of the fish species.

Suitable Ecological Flow Discussion
For further analysis of the simulation results, several scenarios are selected for analysis. The spatial distribution of habitat suitability for several scenarios is shown in Figures 12 and 13. According to the simulation results of velocity and water depth in Figures 12 and 13, the velocity is larger at the dam site, but the excessive velocity leads to poor habitat suitability. The water depth at the bend of the river is deep, and the deepwater depth leads to poor habitat suitability.
The spatial distribution of habitat suitability indicators is shown in Figure 12, where the reservoir inflow is equal to outflow. It can be seen from the results that the large inflow of the reservoir makes the water body mix well, which leads to a decrease of the thermal stratification stability. Therefore, the temperature difference between the inflow and the outflow is relatively small, which is beneficial to C. idella spawning. Similarly, the small inflow of the reservoir makes the water body mix slowly, resulting in an enhanced thermal stratification stability. Therefore, the temperature difference between the inflow and the outflow is relatively high, which is unfavorable for C. idella spawning. Comparing the spatial distribution characteristics of the four different scenarios in Figure 12, the overall suitability is optimal at 120 m 3 /s. Habitat suitability is extremely poor when the reservoir outflow is below 40 m 3 /s. This shows that too much or too little flow will be detrimental to C. idella spawning. The results show that when reservoir inflow and reservoir outflow are excessively large or small, habitat suitability is so poor that the hydraulic condition cannot meet the spawning requirements of C. idella. Figure 13 shows the spatial distribution of habitat suitability index when the reservoir inflow is 120 m 3 /s and reservoir outflows are 80 m 3 /s, 100 m 3 /s, 120 m 3 /s, 140 m 3 /s, and 160 m 3 /s, respectively. It can be seen from the results that the larger outflow of the reservoir makes the water body mix well for the same reservoir inflow, which leads to the decrease of the thermal stratification stability. Therefore, the temperature difference between the reservoir inflow and the outflow is relatively small, which is beneficial to C. idella spawning, however, increasing the outflow of the reservoir leads to excessive velocity in the local area of the habitat, and also leads to deeper water depth. A too high velocity and too deepwater depth can negatively affect C. idella spawning. On the basis of the spatial distribution characteristics of habitats under five different reservoir outflow scenarios, the overall suitability is optimal at 120 m 3 /s. This further validates the results in the previous section. It is most beneficial to the reproduction of C. idella when the reservoir outflow reaches 120 m 3 /s. Figure 14 shows the spatial distribution of habitat suitability index when the reservoir outflow is 120 m 3 /s and the reservoir outflows are 40 m 3 /s, 60 m 3 /s, 80 m 3 /s, 100 m 3 /s, 120 m 3 /s, 140 m 3 /s, and 160 m 3 /s, respectively. From the analysis in Figure 13, 120 m 3 /s is a suitable reservoir outflow. From the results of Figure 14, it can be seen that under this suitable outflow, the larger reservoir inflow will make the overall habitat suitability index better. The larger reservoir inflow makes the water body mix well for the same reservoir outflow, which leads to a decrease of the thermal stratification stability. Therefore, the temperature difference between the reservoir inflow and the outflow is relatively small, which is beneficial to C. idella spawning. On the basis of the spatial distribution characteristics of habitats under seven different reservoir inflow scenarios, the overall suitability is optimal at 160 m 3 /s. This shows that reservoir inflow has a certain impact on C. idella spawning. The reservoir inflow provides a space for rational allocation of water resources.

Suitable Spawning Areas Discussion
In order to discuss suitable spawning activity areas for C. Idella, this study calculates and analyzes the weighted usable area with a suitability of 0.6 or more. As shown in Table 4, the following results can be seen from the table: (1) under the seven scenarios with low flow, the weighted usable area with a suitability of 0.6 or more is 0; (2) when the reservoir inflow is 100 m 3 /s and outflow is 160 m 3 /s, the weighted usable area with a suitability of 0.6 or more is the largest. The maximum value is 135,396 m3, accounting for 62.6% of the total habitat area of the study area; and (3) when the reservoir inflow is greater than 60 m 3 /s and outflow is greater than 100 m 3 /s, the C. idella has more than 100,000 m 3 area suitable for spawning activity.

Suitable Spawning Areas Discussion
In order to discuss suitable spawning activity areas for C. Idella, this study calculates and analyzes the weighted usable area with a suitability of 0.6 or more. As shown in Table 4, the following results can be seen from the table: (1) under the seven scenarios with low flow, the weighted usable area with a suitability of 0.6 or more is 0; (2) when the reservoir inflow is 100 m 3 /s and outflow is 160 m 3 /s, the weighted usable area with a suitability of 0.6 or more is the largest. The maximum value is 135,396 m 3 , accounting for 62.6% of the total habitat area of the study area; and (3) when the reservoir inflow is greater than 60 m 3 /s and outflow is greater than 100 m 3 /s, the C. idella has more than 100,000 m 3 area suitable for spawning activity.
The main spawning area of C. idella spawning is shown in Figure 15 (select a scenario where the reservoir inflow is 100 m 3 /s and the outflow is 160 m 3 /s). According to the survey data, the C. idella historical sampling points are mainly concentrated in three areas, which is close to the best suitability area of the simulation results. The model has a certain reliability under less data samples and is a better choice.
The main spawning area of C. idella spawning is shown in Figure 15 (select a scenario where the reservoir inflow is 100 m 3 /s and the outflow is 160 m 3 /s). According to the survey data, the C. idella historical sampling points are mainly concentrated in three areas, which is close to the best suitability area of the simulation results. The model has a certain reliability under less data samples and is a better choice.

Future Work
There are still some shortcomings in the current research and further research should be carried out which includes the following: (1) Fish habitats not only provide living space for fish, but also provide all environmental factors such as water temperature, substrate, flow rate, pH value, and dissolved oxygen that satisfy the survival, growth, and reproduction of fish. The impact of various hydraulic factors on fish habitats should be considered in the future. (2) In addition to C. idella in the study area, other economically valuable local fish should be considered as research protection targets. Habitat models that consider multiple species should be established in the future. (3) The current research is only for the fish spawning period, and the whole life cycle of the fish should be considered in the future, including the hydraulic conditions of the feeding ground, the spawning ground, and the wintering ground.

Conclusion
Construction of water conservation projects changes natural flow regimes of rivers and causes adverse impacts on rivers such as obstruction of fish migration paths, as well as altering water quality and water temperature, and thereby causing profound impacts on the aquatic ecosystem. By coupling the fuzzy inference system, a two-dimensional laterally-averaged hydrodynamic model (CE-QUAL-W2), and a two-dimensional shallow water hydrodynamic model (MIKE21), the fuzzy logic-based habitat model is developed for evaluating the effect of inflow and outflow of the Xuanwei Reservoir on the spawning habitat suitability of Ctenopharyngodon idella (C. idella) in the Qingshui River. The

Future Work
There are still some shortcomings in the current research and further research should be carried out which includes the following: (1) Fish habitats not only provide living space for fish, but also provide all environmental factors such as water temperature, substrate, flow rate, pH value, and dissolved oxygen that satisfy the survival, growth, and reproduction of fish. The impact of various hydraulic factors on fish habitats should be considered in the future. (2) In addition to C. idella in the study area, other economically valuable local fish should be considered as research protection targets. Habitat models that consider multiple species should be established in the future. (3) The current research is only for the fish spawning period, and the whole life cycle of the fish should be considered in the future, including the hydraulic conditions of the feeding ground, the spawning ground, and the wintering ground.

Conclusion
Construction of water conservation projects changes natural flow regimes of rivers and causes adverse impacts on rivers such as obstruction of fish migration paths, as well as altering water quality and water temperature, and thereby causing profound impacts on the aquatic ecosystem. By coupling the fuzzy inference system, a two-dimensional laterally-averaged hydrodynamic model (CE-QUAL-W2), and a two-dimensional shallow water hydrodynamic model (MIKE21), the fuzzy logic-based habitat model is developed for evaluating the effect of inflow and outflow of the Xuanwei Reservoir on the spawning habitat suitability of Ctenopharyngodon idella (C. idella) in the Qingshui River. The model is constructed considering the reservoir and the downstream river channel, and explores the comprehensive effects of water temperature, velocity, and water depth on habitat suitability.
On the basis of the historical runoff data of the hydrological station, 64 different reservoir inflow and outflow scenarios are designed and simulated. The weighted usable area (WUA) and hydraulic habitat suitability (HHS) under each scenario are calculated and analyzed to obtain a suitable ecological flow range. The results show that the ecological flow for C. idella spawning can satisfy the basic demand when the reservoir inflow is greater than 60 m 3 /s and the reservoir outflow is greater than 100 m 3 /s. The habitat ecological suitability is the best when the reservoir outflow is 120 m 3 /s. The suitable ecological flow range provides a reference and guidance for the ecological scheduling of reservoirs and plays an important role in conserving fish habitats. Further works should focus on data collection to consider other important physical variables such water quality, dissolved oxygen, and sediment in order to model a complex river ecosystem with more accuracy. Habitat simulation can evaluate river health and provide support for river management, and therefore it is needs to be quickly put into practice to improve local river ecology and environment.