Next Article in Journal
New Hybrids of ANFIS with Several Optimization Algorithms for Flood Susceptibility Modeling
Previous Article in Journal
Experimental Study of Influence of Karst Aquifer on the Law of Water Inrush in Tunnels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculation of Comprehensive Ecological Flow with Weighted Multiple Methods Considering Hydrological Alteration

1
School of Water Conservancy & Civil Engineering, Northeast Agricultural University, Harbin 150030, China
2
The Key Laboratory of High Efficient Utilization in Agricultural Water Resources of Ministry of Agriculture of China, Harbin 150030, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(9), 1212; https://doi.org/10.3390/w10091212
Submission received: 18 July 2018 / Revised: 16 August 2018 / Accepted: 31 August 2018 / Published: 7 September 2018
(This article belongs to the Section Hydrology)

Abstract

:
Instream ecological flow is an essential determinant of river health. Intra- and interannual distribution characteristics of runoff have been altered to different degrees by dam construction. Historical runoff series with alterations, as basic data for ecological flow calculation, provide minimal instream hydrological process information, which affects the credibility of calculation results. Considering the influence of the alterations in runoff series on ecological flow calculation, the Gini coefficient (GI) is introduced to study the evenness degrees of the intra-annual runoff distribution of four hydrological stations located in the Naolihe basin of the Sanjiang Plain. The hydrological alteration diagnosis system is used to examine the alteration points in the GI series of each hydrological station for selecting reasonable subsequences. Based on the selected subsequences, the ecological flow of each station is calculated using three hydrological methods, and the comprehensive ecological flow is calculated using weighted calculation results from the three hydrological methods. The study results show that ecological flow and natural flow have similar processes with two peaks occurring in the process in May and August, respectively. Also, dams decrease the ecological water requirement damage frequency in dry seasons, but overuse of water resources increases the ecological water requirement damage frequency in flood seasons.

1. Introduction

Riverine ecosystems play a significant role in the survival and development of human society, so protecting riverine ecosystems is vital. Ecological flow is the runoff that maintains the riverine ecosystem health and biodiversity [1,2,3]. Studies on ecological flow have developed rapidly [4]. Many methods have been presented to calculate ecological flow [5,6], of which hydrological methods are widely used due to the simplicity of the data requirements and low calculation costs. A common hypothesis states that the ecological regime becomes well-adjusted to the hydrological regime, and the ecological requirements of river-living organisms have been satisfied from a hydrological standpoint [7]. The hydrological method calculates ecological flow by analyzing the frequency distribution characteristics of historical runoff series. Some representative hydrological methods, such as the Tennant method [8], NGPRP (Northern Great Plains Resource Program) method [9], minimum monthly mean flow method [10], month-by-month minimum ecological flow method [11], flow duration curve method [12], and the 7Q10 (seven-day low flow with a 10-year recurrence interval) method [13]. For example, Tennant defined 10% to 30% of the annual mean flow as the ecological flow range by analyzing the respective relationship between flow and river width, and velocity and depth [14]. Richter [15] classified 33 hydrological parameters from flow magnitude, frequency, duration, timing, predictability, and rate of change to study riverine ecosystem structure and function, and applied quantized hydrologic indicators to study ecological flow. Li [16] proposed a monthly frequency method for ecological flow calculation and evaluated the calculated ecological flow using the Tennant method. The results revealed that monthly ecological flow calculated using the monthly frequency method was higher than that calculated by other ecological methods. The ecological flow calculated by the monthly frequency method allocates more water for environmental use and greater benefits for the river health and habitat conservation in river flood plain and riparian corridors. Using the hydrological method to calculate ecological flow should consider data representativeness and method reasonableness. Namely, users should consider whether the selected historical runoff series represents the natural runoff characteristics and if the method for calculating ecological flow is reasonable. Ecological flow is calculated based on historical runoff series, parts of which can be altered to some extent because of human activities, natural disasters, and so on [17,18]. Therefore, runoff alterations should be evaluated in ecological flow calculations. Analysis on runoff alterations is based on annual runoff in some papers [19,20], but annual runoff poorly reflects the characteristics of intra-annual runoff distribution. Natural monthly flow within a year varies, so diagnosing runoff alterations should consider intra-annual runoff distribution. Each ecological flow calculation method has its own foundational principle and obtains different result for the same river. It is difficult to judge which method is more suitable to the natural situation, due to the real value of ecological flow not being observed in practice [8,9,10,11,12,13]. Furthermore, some methods focus on the maintenance of riverine ecosystem health, while some methods focus on the utilization of water resources. The comprehensive ecological flow took riverine ecosystem health and utilization of water resources into account at the same time by weighting results of different calculation methods.
Above all, a comprehensive ecological flow was defined as an ecological flow obtained by weighting the results of multiple methods to integrate the results of each method. In order to calculate the comprehensive ecological flow under the condition of runoff alteration, the Gini coefficient (GI) was introduced to calculate the degrees of evenness of intra-annual runoff distribution based on all monthly flow in a year, and one GI value is used to holistically reflect intra-annual runoff distribution for a year [21]. Then, a hydrological alteration diagnosis system [20] was used to examine the alteration year based on GI series instead of annual runoff series in each hydrological station to select subsequences for ecological flow calculation [22]. Based on the selected subsequences, the range of variability approach (RVA) method [23], monthly frequency method (MFM) [16], and annual daily mean flow frequency method [24] (ADMFFM) were used to calculate ecological flow. The three indexes (the rate of monthly flow, satisfaction degree of monthly ecological water requirement, and suitability degree of monthly average ecological water requirement) were used to evaluate the calculation results of each method. The weights of each ecological flow calculation result were calculated based on the three evaluation indexes above. Then, a comprehensive ecological flow was calculated based on the weights and results for the three methods.

2. Materials and Methods

2.1. Study Area

The Naolihe River is located in the east of Heilongjiang Province, China. It is the first tributary of the Wusulijiang River. The Naolihe basin has an area of 24,863 km2 with an elongated shape. The average basin length is about 270 km, and the average width is about 90 km. The ratio of length to width is 3:1. The Naolihe basin lies in the temperate continental monsoon climate zone. The multiyear average precipitation is 545 mm, and 70% of the precipitation is mostly concentrated in June, July, August, and September. Moreover, the precipitation in July and August accounts for about 44% of the total. The precipitation in May and June only accounts for 23%.
Four hydrological stations are present in the Naolihe basin: Baoqing, Baoan, Caizuizi, and Hongqiling, as shown in Figure 1. The Baoqing hydrological station is located in the town of Baoqing, with a drainage area of 3689 km2. The Caizuizi hydrological station is located at the trunk stream of Naolihe River, in the town of Sanlitun, and has a control area of 20,896 km2. The Baoan hydrological station is located in the town of Youyi, having a drainage area of 1344 km2. The Hongqiling hydrological station is located on the Hongqiling farm. Its drainage area is 1147 km2. Two leading reservoirs are found in the Naolihe basin: the Longtouqiao Reservoir and the Hamatong Reservoir, shown in Figure 1. The first reclamation in the Naolihe basin began in 1956, and no agricultural activities were performed in the study area prior to 1956. Therefore, 1956 was chosen as the start year for this study. The historical daily flow data, except for December to March in the next year due to the freezing period, is from 1956 to 2012 for the four hydrological stations in the Naolihe Basin for this study. Location of the study area is shown in Figure 1.

2.2. Methodology

The comprehensive ecological flow was calculated using the weighting mean method. The flow chart of the study process is shown in Figure 2.

2.2.1. Evaluation of the Intra-Annual Runoff Distribution Degrees of Evenness

The Gini coefficient is an index used to quantify and evaluate the degrees of distribution evenness. In this paper, it is used to evaluate intra-annual runoff distribution evenness degrees. The GI of runoff series is calculated as follows:
(1)
Group and sort historical monthly flow data. Divide the 12 monthly flows in the ith year into a group, and sort the monthly flows into a group in ascending order. The ascending runoff series of the ith year is { r i 1 , r i 2 , , r i 12 } , where i represents the ith year, r i 1 < r i 2 < < r i 12 .
(2)
Accumulate the ascending runoff data in each group:
R i k = j = 1 k r i j ( k = 1 , 2 , , 12 )
where k is the ranking order.
(3)
Draw a Lorenz curve of the annual distribution for monthly flow. Take k/12 as the abscissa, R i k / R i 12 as the ordinate, and draw a Lorenz curve of the annual distribution for monthly flow. The Lorenz curve is shown in Figure 3.
(4)
The ith year runoff Gini coefficient GIi is computed as:
G I i = S A i / ( S A i + S B i )
where SAi and SBi are the area acreages of Ai and Bi, respectively.
The calculation details of GI series in each station are shown in Appendix A.
The smaller the obtained value, the more even the runoff distribution in a year. Conversely, the larger the obtained value, the more uneven the runoff distribution in a year.

2.2.2. Hydrological Alteration Diagnosis System

A single hydrological alteration diagnosis method occasionally produces unreliable results, and multiple hydrological alteration diagnosis methods produce different results. The hydrological alteration diagnosis system synthetically considers the diagnosis results from multiple diagnosis methods and examines the alteration form and alteration point holistically. The steps we followed for alteration diagnosis are as follows. Firstly, we used the Hurst coefficient [25] method and the moving average method to form a primary diagnosis and judge whether or not the series contains an alteration. If so, then various examination methods were used to conduct a detailed diagnosis, including three trend diagnosis methods (the correlation coefficient method, Spearman rank correlation method, and Kendall rank correlation method) and 11 jump diagnosis methods (the Lee–Heghinan method, rank test, slide F test, R/S method, Mann–Kendall method, Bayesian method, etc.). The diagnosis results were also classified into two types, trend results and jump results, and the results were synthesized. Finally, efficiency coefficients were calculated for identifying the alteration form, and the alteration form was judged to determine if one coefficient was bigger than the other one.

2.2.3. Ecological Flow Calculation Method

  • In the range of the variability approach (RVA) method, ecological flow calculations should consider the hydrological regime. Parameters of the indicators of the hydrologic alteration (IHA) method are all closely related to runoff, and monthly flow influences river-living organisms, soil, etc. Therefore, the RVA analyses the monthly flow frequency distribution of each month and selects the flow corresponding to 25% and 75% in the frequency distribution as the upper and lower limits of the monthly flow, respectively [26]. Ecological flow is calculated as:
    Q e j = W j R V A × Q e j R V A + W j M F M × Q e j M F M + W j A D M F F M × Q e j A D F M M F
    where Q e j R V A is the ecological flow of the jth month, Q ¯ j is monthly mean flow of the jth month, Q j 25 % is the upper limit of RVA, and Q j 75 % is the lower limit of RVA.
  • The steps for ecological flow calculation using the monthly frequency method (MFM) are as follows:
    (1)
    Calculate the monthly flow distribution empirical frequency of the jth month.
    (2)
    Select the probability distribution function and draw the monthly flow distribution theoretical frequency curve of the jth month. The Pearson III distribution curve and generalized extreme value (GEV) distribution curve are commonly used. The GEV distribution curve is better fitting with the runoff data [27].
    (3)
    Select the flow corresponding to the maximum frequency in the GEV distribution curve as the ecological flow of the jth month Q e j M F M [28].
  • The annual daily mean flow frequency method (ADMFFM) assumes that daily flow may appear in a month with a certain probability. The frequency of daily flow in the jth month of the historical runoff series is calculated and the daily flow corresponding to 60% is selected as the ecological flow [24] of the jth month, Q e j A D M F F M . The ecological flow calculated using this method meets the ecological water requirements in different cases.
  • In the comprehensive ecological flow calculation method, ecological flow maintains riverine ecological function. Too much or too little will affect the riverine ecosystem. Therefore, the weights of each ecological flow calculation method are used to comprehensively calculate ecological flow. The steps for this method are as follows:
    (1)
    Evaluate ecological flows from RVA, MFM, and ADMFFM, and the comprehensive evaluation values of RVA, MFM, and ADMFFM of the jth month are W j R V A , W j M F M and W j A D M F F M , respectively.
    (2)
    Calculate the comprehensive ecological flow of the jth month as follows:
    Q e j = W j R V A × Q e j R V A + W j M F M × Q e j M F M + W j A D M F F M × Q e j A D F M M F

2.2.4. Weight Calculation of Comprehensive Ecological Flow

The three evaluation indexes listed in Section 2.2.4 were used to evaluate the calculation results following the above ecological flow calculation methods. Three single evaluation results were synthesized into a comprehensive evaluation result. The weights of each single ecological flow calculation result were obtained. The detailed calculation steps are as follows.
  • We first determined the deviation rate of monthly flow. The median of the natural daily mean flow series in the jth month of N years was called the standard value. This index is a ratio of ecological flow in the jth month to the standard value. It reflects the deviation degree between the ecological flow and the natural flow in the same period [28,29]. It is computed as follows:
    (1)
    Calculate the natural daily mean flow in the jth month of N years Q i j , where i is the year, i = 1, 2, …, N; N is equal to the length of the runoff series; and j is the month, j = 1, 2, …, 12;
    (2)
    Sort Q i j in the ascending order Q 1 j , Q 2 j , , Q n j , ( Q 1 j < Q 2 j < < Q n j , n is the serial order), where the standard value of the jth month is Q m i d , j ;
    (3)
    The calculation formula of the deviation rate of the monthly flow is:
    C j = Q e j Q m i d , j
    where C j is the deviation rate of the monthly flow of the jth month and Q e j is the ecological flow of the jth month, in m3/s. If the value is close to 1, the calculated flow is approaching the natural flow and the riverine ecosystem is healthy.
  • Satisfaction degree of the monthly ecological flow. This a ratio of the days where the natural flow is equal to or greater than the ecological flow to the total number of days in the same month. The formula to calculate the satisfaction of the monthly ecological water requirement is:
    P j = D e j D j = i = 1 N k = 1 K n u m sgn i j k D j
    sgn i j k = { 1 , Q i j k > Q e j 0 , Q i j k < Q e j
    D j = N × K n u m , j
    where P j is the satisfaction degree of the monthly ecological flow of the jth month, D e j is the number of days where the natural flow was equal to or greater than the ecological flow in the jth month, D j is the total number of days in the jth month in N years, K n u m , j is the number of days in the jth month, Q i j k is the natural daily flow on the kth day of the jth month in the ith year in m3/s, and Q e j is the ecological flow of the jth month in m3/s. The greater the satisfaction degree of the monthly ecological flow, the healthier the riverine ecosystem.
  • Suitability degree of the monthly ecological flow. The monthly ecological flow discrete coefficient is the sum of the discrete degree between the median and the characteristic extreme value of the ecological flow and natural flow. It reflects the suitability degree between the ecological flow and natural flow in a month.
    F j = ( Q e m i d , j Q m i d , j Q m i d , j ) 2 + ( Q e max , j Q max , j Q max , j ) 2
    where F j is the discrete coefficient of the monthly ecological flow in the jth month, Q e m i d , j is the median of the jth month’s ecological flow of over N years, Q m i d , j is the median of the jth month’s natural flow of N years, Q e max , j is the maximum of the jth month’s ecological flow of N years, and Q max , j is the maximum of the jth month’s natural flow of N years. In the case of F j being greater than 10, it is considered as completely discrete and the value of F j is taken as 10.
    The suitability degree of the monthly ecological flow is computed as:
    F j = 1 F j / 10
    When F j equals 1, a complete suitableness exists between the ecological flow and natural flow in the jth month; when it equals 0, absolutely no suitableness exists between the ecological flow and natural flow in the jth month.
    This index reflects the degree of suitability between the total ecological flow requirement and the natural flow in a month. When F j is 1, the total ecological flow requirement and the natural flow completely match. When F j is 0, the flow is completely unsuitable. A discrete coefficient greater than 10 can be considered completely discrete and the value of the discrete coefficient is equal to 10.
  • Ecological flow comprehensive evaluation. The above indexes compare the ecological flow and natural flow from different aspects. The deviation rate of the monthly flow evaluates the degree and magnitude of deviation between the ecological flow and natural flow. The satisfaction degree of the monthly ecological flow temporally analyses the degree of satisfaction between the ecological flow and natural flow. The suitability degree of the monthly ecological flow evaluates the degree of suitability between the ecological flow and natural flow in discrete degrees. The above three indexes were synthesized into a comprehensive index α j to evaluate ecological flow in the jth month. The formula for calculating α j is:
    α j = C j P j F j 3
    where α j [ 0 , 1 ] and a value of α j close to 1 means the calculated ecological flow is approaching the natural flow and the riverine ecosystem is healthy.
  • Weight calculation of the comprehensive ecological flow. The geometric mean method was used to calculate the weights for each ecological flow calculation method result.
    W j R V A = α j R V A α j R V A + α j M F M + α j A D M F F M
    W j M F M = α j M F M α j R V A + α j M F M + α j A D M F F M
    W j A D M F F M = α j A D M F F M α j R V A + α j M F M + α j A D M F F M

3. Results

3.1. Hydrological Alteration Diagnosis System

The results of the change in trend in the GI series obtained using the moving average method are shown in Figure 4; the moving average curve for the Hongqiling and Caizuizi stations showed a rising trend. In addition, slight fluctuations were observed in the moving average curves for the Baoqing and Baoan stations. Also, a downward trend exists in the mean annual flow curves shown in Figure 5. Additionally, GI is the index used to quantify and evaluate the evenness degrees of mean annual flow distribution, so there were different trends in GI curves and mean annual flow curves for the same station.

3.2. Results of Ecological Flow Calculation

Based on the selected subsequence above, the RVA method, monthly frequency method, and annual daily mean flow frequency method were used to calculate the ecological flows of the four hydrological stations. Comprehensive evaluation results for each method were determined by comprehensive evaluation in Section 2.2.4. The weights of the ecological flow calculated using each method were calculated from comprehensive evaluation results, and the comprehensive ecological flow based on the above methods was obtained. The ecological flows determined using the above methods and their weights are exhibited in Table 2 and Table 3, respectively. The comprehensive ecological flow calculation results are shown in Figure 6.
As shown in Figure 6, comparing the ecological flow values determined by different calculation methods, the ecological flow values from the RVA method were generally low and the ecological flow values from the monthly frequency method were relatively high. However, the comprehensive ecological flow values were all within the extreme values of the three calculation results. The comprehensive ecological flow and the ecological flow obtained from the three methods were all less than the monthly average flow, which is consistent with the characteristics of the ecological flow. In addition, the ecological flow process exactly coincides with the natural flow process, which is shown by the identical timing of extreme ecological flow and extreme natural flow.
The monthly flow meets the ecological water requirements when the monthly flow is greater than the ecological flow, otherwise the ecological water requirement is not satisfied, which means the riverine ecosystem may be damaged. The ecological water requirement damage frequency is a ratio of the number of times the monthly flow was less than the ecological flow to the series length. The calculation results of the ecological water requirement damage frequency of the subsequences and the whole series are shown in Table 4.

4. Discussion

The ecological flow determined using the hydrological method, in which there is no consideration of the effect of some special situations such as overbank flow on the riverine ecosystem, was calculated from historical runoff data. Therefore, the calculated ecological flow regime is in concordance with the natural runoff regime. As the Naolihe basin is located in the cold region of China, precipitation from December to March in the next year is mostly dominated by sleet and solid precipitation involving the formation of snow. The sleet and solid precipitation flow into rivers in the form of snowmelt runoff in April and May in the following year. Flow in April and May accounts for a large proportion of annual flow. Therefore, rivers in the Naolihe basin generally have two flood processes: the spring freshet in April and May and the summer freshet in July and August. As shown in Figure 5, the ecological flow process, which corresponds to the natural flow process, also has two peaks. The maximum peak appears in May, which is consistent with organism reproduction needs in the alpine region.
In Figure 6, the comprehensive ecological flow values were between the values calculated by the three methods. Compared to the bigger values of the three methods, it is can be considered that more water in the river could be utilized for irrigation, drinking, power generation, and so on based on values of the comprehensive ecological flow. At the same time, the water utilization, according to the comprehensive ecological flow, is enough to guarantee the water requirement for the river ecosystem, because they exceed the smaller values calculated by the three methods.
In the Baoqing and Baoan stations, located in the upstream region of the Naolihe river, human activities have relatively little impact on the riverine ecosystem. Therefore, differences in the ecological water requirement damage frequency for the subsequence series and the whole series were not significant. However, the ecological water requirement damage frequency increased obviously in July, August, and September. This is because the study area is located in an agricultural area, and crops need more water in this period. Human use of water resources for irrigation tends to enhance the ecological water requirement damage frequency. However, ecological water requirement damage frequency averages for the subsequence series and the whole series of the Hongqiling station and Caizuizi station were basically the same. This is because the Hongqiling station and Caizuizi station are located in the middle and downstream portions of the Naolihe river, respectively, and the construction of the Longtouqiao and Hamatong Reservoirs effectively solves the problem of instream flow decreasing during the non-flood season. Similarly, during the flood season, due to the impoundment process of reservoirs, part of the natural flow is stored. The instream flow decreases and the ecological water requirement damage frequency increases. In addition, the average ecological water requirement damage frequencies in the Hongqiling station and Caizuizi station were more than in the Baoqing station and Baoan station. This finding is also because of human development.

Author Contributions

Conceptualization, Z.X. and X.G.; Methodology, Y.W., Y.J., J.W.; Validation, Z.X., Y.W., Y.J.; Formal Analysis, Y.W., Z.X.; Investigation, Z.X., Y.W.; Resources, Q.F., Z.X.; Data Curation, Y.W.; Z.X.; Writing-Original Draft, Z.X.; Writing-Review & Editing, Z.X., Y.W.; Visualization, Y.W.; Supervision, X.G., Q.F.; Projection Administration, Z.X., X.G.; Funding Acquisition, X.G.

Funding

This research was funded by the National Key R&D Program of China (No. 2017YFC0406004), the Natural Science Foundation of Heilongjiang province of China (No. E2015024), the Projects for Science and Technology Development of Water Conservancy Bureau in Heilongjiang Province of China (No. 201402, No. 201404, No. 201501), and the Academic Backbones Foundation of Northeast Agricultural University (No. 16XG11).

Acknowledgments

The present study was supported by the Administrative Center for China’s Agenda 21, the Council of Nature Science Funding of Heilongjiang Province, and Water Conservancy Bureau of Heilongjiang Province and Northeast Agricultural University.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. An Example for GI Value Calculation

An example for the GI of the Baoqing station in 1956 is calculated as follows:
The monthly mean flow (monthly flow) series was listed in Table A1.
Table A1. Historical monthly flow of Baoqing station in 1956.
Table A1. Historical monthly flow of Baoqing station in 1956.
Month123456789101112
Monthly mean flow0.080.0090.02718.437.727.489.560.147.518.66.60.57
First, sort the monthly flows into a group in ascending order; the result is listed in column (3) in Table A2.
Second, accumulate the ascending monthly flow series by Equation (1); the result is listed in column (4).
Third, calculate the runoff cumulative frequency; for example, 0.000117461 = 0.036/306.486; and calculate the time cumulative frequency; for example, 0.166667 = 2/12.
Fourth, draw a Lorenz curve with the series in column (5) and (6) in Table A2.
Table A2. Detailed calculation results.
Table A2. Detailed calculation results.
MonthOrderMonthly FlowAccumulated ValueRunoff Cumulative FrequencyTime Cumulative Frequency
(1)(2)(3)(4)(5)(6)
210.0090.0090.0000029650.083333
320.0270.0360.0001174610.166667
130.080.1160.0003784840.25
1240.570.6860.0022382750.333333
1156.67.2860.0237727010.416667
4618.425.6860.083808070.5
10718.644.2860.1444959970.583333
6827.471.6860.2338964910.666667
5937.7109.3860.3569037410.75
91047.5156.8860.511886350.833333
81160.1216.9860.7079801360.916667
71289.5306.48611
Figure A1. Lorenz curve of the annual distribution for monthly flow in 1956. Note: A i is a region bounded by Lorenz curve and dashed diagonal, B i is a region bounded by Lorenz curve and vertical diagonal.
Figure A1. Lorenz curve of the annual distribution for monthly flow in 1956. Note: A i is a region bounded by Lorenz curve and dashed diagonal, B i is a region bounded by Lorenz curve and vertical diagonal.
Water 10 01212 g0a1
Then, calculate the areas of A i and B i ; the area of A i is 0.2876 and the sum of the areas of A i and B i is 0.5. So, the GI value was calculated by Equation (2), and the result was 0.5752.
The calculation results for each year at the four stations are shown in Table A3.
Table A3. GI of the four stations.
Table A3. GI of the four stations.
BaoqingBaoanCaizuiziHongqiling
19560.5751670.5428040.4655290.428152
19570.6637210.5994170.5625830.477699
19580.6548110.6214350.59550.494503
19590.5535470.5300330.654460.524602
19600.5721320.5534990.4379410.414069
19610.6654260.6464240.433190.411643
19620.5592090.5531480.4965310.443979
19630.6165810.5150540.3523640.370382
19640.68560.656230.6309190.512584
19650.6031780.5808380.3776960.383314
19660.6501940.5455970.433940.412027
19670.6102270.6060120.7131390.554558
19680.6348820.4694680.548390.470453
19690.429670.4419010.3487850.368555
19700.7385350.6158520.6861960.540803
19710.6073470.5883130.4694130.430135
19720.6179860.5432920.5643690.47861
19730.6073210.6074210.3755520.382219
19740.5550220.5823960.2874870.337262
19750.7134960.6613350.712160.554058
19760.6844770.5705560.7138340.554912
19770.5952980.5381280.836570.617569
19780.5202680.5458120.6366220.515496
19790.6111830.6324460.610840.470432
19800.5153410.5340830.6270680.44741
19810.6130040.5816390.5977520.603207
19820.6003860.5565390.7593670.617718
19830.6977380.6280780.6553490.477994
19840.5243240.4108190.4803680.389041
19850.6053470.5138350.5847530.416934
19860.5476960.4931030.6817880.551018
19870.5990180.5580520.4769330.445092
19880.6323030.5334060.6357090.492935
19890.6306920.4123650.7148850.527368
19900.5374020.5072990.5877280.455828
19910.7058460.6923140.5611830.567586
19920.7085550.630150.7442960.472671
19930.5770320.4457920.6248470.508258
19940.6319190.5827620.6002890.546941
19950.4874170.4520520.6683520.472467
19960.7170550.5905420.7193610.550217
19970.6268370.5586180.65470.407772
19980.7122930.4657540.7214650.589609
19990.6542390.6512770.6975950.650854
20000.7232930.6375650.6670180.640227
20010.7104920.6293170.7043670.595746
20020.5642980.5358560.5546640.500066
20030.5815970.5469150.6301940.512214
20040.5859740.5497130.6305970.51242
20050.5859710.5491780.7366040.566536
20060.6318460.5589460.6377880.516091
20070.6219540.6237160.7228090.559494
20080.5869230.5974980.7386380.567575
20090.6045150.6390230.711460.622316
20100.6800250.6516260.7933930.595616
20110.6520220.5466810.6922610.555792
20120.5458630.5673070.6830580.539201

References

  1. Zhang, Q.; Zhang, Z.J.; Shi, P.J.; Singh, V.P.; Gu, X.H. Evaluation of ecological instream flow considering hydrological alterations in the Yellow River basin, China. Glob. Planet. Chang. 2018, 160, 61–74. [Google Scholar] [CrossRef]
  2. Xia, Z.Q.; Li, Q.F.; Guo, L.D.; Li, J.; Zhang, W.; Wang, H.Q. Computation of minimum and optimal instream ecological flow for the Yiluohe River. IAHS-AISH Publ. 2007, 315, 142–148. [Google Scholar]
  3. Yang, F.C.; Xia, Z.Q.; Yu, L.L.; Guo, L.D. Calculation and analysis of the instream ecological flow for the Irtysh River. Procedia Eng. 2012, 28, 438–441. [Google Scholar] [CrossRef]
  4. Xu, Z.X.; Wu, W.; Yu, S.Y. Ecological baseflow: Progress and challenge. J. Hydroelectr. Eng. 2016, 35, 1–11. [Google Scholar] [CrossRef]
  5. Yang, Z.F.; Cui, B.S.; Liu, J.L. Estimation methods of eco-environmental water requirements: Case study. Sci. China Ser. D Earth Sci. 2005, 48, 1280–1292. [Google Scholar] [CrossRef]
  6. Chen, M.J.; Feng, H.L.; Wang, L.Q.; Chen, Q.Y. Calculation methods for appropriate ecological flow. Adv. Water Sci. 2007, 18, 745–750. [Google Scholar] [CrossRef]
  7. Wang, X.Q.; Liu, C.M.; Yang, Z.F. Method of resolving lowest environmental water demands in river course (I)—Theory. Acta Sci. Circumstantiae 2001, 21, 544–547. [Google Scholar] [CrossRef]
  8. Geoffrey, E.P. Instream flow science for sustainable river management. J. Am. Water Resour. Assoc. 2009, 45, 1071–1086. [Google Scholar] [CrossRef]
  9. Northern Great Plains Resource Program (NGPRP). Instream Needs Sub-Group Report; Work Group C, Water; US Fish and Wildlife Service, Office of Biological Service: Washington, DC, USA, 1974. [Google Scholar]
  10. Li, L.J.; Zheng, H.X. Environmental and ecological water consumption of river systems in Haihe-Luanhe basins. Acta Geogr. Sin. 2000, 55, 495–500. [Google Scholar] [CrossRef]
  11. Yu, L.J.; Xia, Z.Q.; Du, X.S. Connotation of minimum ecological runoff and its calculation method. J. Hehai Univ. (Nat. Sci.) 2004, 32, 18–22. [Google Scholar]
  12. Ravindra, K.V.; Shankar, M.; Sangeeta, V.; Surendra, K.M. Design flow duration curves for environmental flows estimation in Damodar River Basin, India. Appl. Water Sci. 2017, 7, 1283–1293. [Google Scholar] [CrossRef]
  13. Danial, C.; Nassir, E.J. Comparison and regionalization of hydrologically based instream flow. Can. J. Civ. Eng. 1995, 5, 235–246. [Google Scholar] [CrossRef]
  14. Donald, L.T. Instream flow regimens for fish, wildlife, recreation and related environmental resources. Fish. Mag. 1976, 1, 6–10. [Google Scholar] [CrossRef]
  15. Brian, D.R.; Jeffrey, B.; Robert, W.; David, P.B. How much water does a river need? Freshw. Biol. 1997, 37, 231–249. [Google Scholar] [CrossRef]
  16. Li, J.; Xia, Z.Q.; Ma, G.H. A new monthly frequency computation method for instream ecological flow. Acta Ecol. Sin. 2007, 27, 2916–2921. [Google Scholar]
  17. Zhan, C.S.; Wang, H.X.; Xu, Z.X.; Li, L. Quantitative study on the climate variability and human activities to hydrological response in the Weihe River Basin. In Proceedings of the China Academy of Natural Resources Academic Annual Meeting 2011, Wulumuqi, China, 25 July 2011. [Google Scholar]
  18. Poff, N.L.; Bledsoe, B.P.; Cuhaciyan, C.O. Hydrologic variation with land use across the contiguous United States: Geomorphic and ecological consequences for stream ecosystems. Geomorphology 2006, 79, 264–285. [Google Scholar] [CrossRef]
  19. Huang, S.Z.; Chang, J.X.; Huang, Q.; Wang, Y.M.; Chen, Y.T. Calculation of the instream ecological flow of the Wei River based on Hydrological Variation. J. Appl. Math. 2014, 11, 1–9. [Google Scholar] [CrossRef]
  20. Cai, S.Y.; Lei, X.H.; Meng, X.Y. Ecological flow analysis method based on the comprehensive variation diagnosis of Gini coefficient. J. Intell. Fuzzy Syst. 2018, 34, 1025–1031. [Google Scholar] [CrossRef]
  21. Xie, P.; Chen, G.C.; Lei, H.F.; Wu, F.Y. Hydrological alteration diagnosis system. J. Hydroelectr. Eng. 2010, 29, 85–91. [Google Scholar]
  22. Hu, C.X.; Xie, P.; Xu, B.; Chen, G.C.; Liu, X.Y.; Tang, Y.S. Variation analysis method for hydrologic annual distribution homogeneity based on Gini coefficient. A case study of runoff series at Longchun station in Dongjiang river basin. J. Hydroelectr. Eng. 2012, 31, 7–13. [Google Scholar]
  23. Yang, T.; Zhang, Q.; Chen, Y.D. A spatial assessment of hydrologic alteration caused by dam construction in the middle and lower Yellow River, China. Hydrol. Process. 2008, 24, 262–267. [Google Scholar] [CrossRef]
  24. Li, Z.Y.; Liu, D.F.; Huang, Q.; Zhang, S.A. Study on the ecological flux of the Ziwu River in Hanjiang River based on several hydrological methods. J. North China Univ. Water Resour. Electr. Power (Nat. Sci. Ed.) 2017, 38, 8–12. [Google Scholar] [CrossRef]
  25. Xie, P.; Chen, G.C.; Lei, H.F. Hydrological alteration analysis method based on hurst coefficient. J. Basic Sci. Eng. 2009, 17, 32–39. [Google Scholar]
  26. Liu, G.H.; Zhu, J.X.; Xiong, M.Y.; Wang, D.; Qi, S.H. Assessment of hydrological regime alteration and ecological flow at Meigang station of Xinjiang River. J. China Hydrol. 2016, 26, 51–57. [Google Scholar]
  27. Li, J.F.; Zhang, Q.; Chen, X.H.; Jiang, T. Study of ecological instream flow in Yellow River, considering the hydrological change. Acta Geogr. Sin. 2011, 66, 99–110. [Google Scholar]
  28. Zhang, Q.; Cui, Y.; Chen, Y.Q. Evaluation of ecological instream flow of the Pearl River basin. Ecol. Environ. Sci. 2010, 19, 1828–1837. [Google Scholar]
  29. Halwatura, D.; Najim, M.M.M. Environmental flow assessment—A analysis. J. Environ. Proffessionals Sri Lanka 2014, 3, 1–11. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Water 10 01212 g001
Figure 2. Flow chart of the study process. Note: RVA is the range of variability approach method.
Figure 2. Flow chart of the study process. Note: RVA is the range of variability approach method.
Water 10 01212 g002
Figure 3. Lorenz curve of annual distribution for monthly flow. Note: A i is a region bounded by Lorenz curve and dashed diagonal, B i is a region bounded by Lorenz curve and vertical diagonal.
Figure 3. Lorenz curve of annual distribution for monthly flow. Note: A i is a region bounded by Lorenz curve and dashed diagonal, B i is a region bounded by Lorenz curve and vertical diagonal.
Water 10 01212 g003
Figure 4. The values and their moving average of the Gini coefficient (GI) from 1965 to 2012: (a) Baoqing station, (b) Baoan station, (c) Hongqiling station, and (d) Caizuizi station. Note: 5a moving average is a moving average of 5 years, and 10a moving average is a moving average of 10 years.
Figure 4. The values and their moving average of the Gini coefficient (GI) from 1965 to 2012: (a) Baoqing station, (b) Baoan station, (c) Hongqiling station, and (d) Caizuizi station. Note: 5a moving average is a moving average of 5 years, and 10a moving average is a moving average of 10 years.
Water 10 01212 g004
Figure 5. The mean annual flow from 1965 to 2012: (a) Baoqing station, (b) Baoan station, (c) Hongqiling station, and (d) Caizuizi station. The results of the hydrological alteration diagnosis system show a jump alteration in the GI series in 1966 at the Baoqing station, a trend alteration in 1990 in the GI series at the Baoan station, and jump variations in the GI series in 1987 and 1990 at the Hongqiling and Caizuizi stations, respectively. The alteration in the GI series of the Baoqing station occurred in 1966, and no obvious alteration was perceived in the GI series in the subsequent 46 years. Based on the alteration analysis above, runoff data from 1966 to 2012 were selected as the subsequence to calculate ecological flow at the Baoqing station. Runoff data for the periods 1956–1990, 1956–1987, and 1956–1990 were selected as the subsequence to calculate ecological flow at the Baoan, Hongqiling, and Caizuizi stations, respectively. Results of the hydrological alteration diagnosis system are shown in Table 1.
Figure 5. The mean annual flow from 1965 to 2012: (a) Baoqing station, (b) Baoan station, (c) Hongqiling station, and (d) Caizuizi station. The results of the hydrological alteration diagnosis system show a jump alteration in the GI series in 1966 at the Baoqing station, a trend alteration in 1990 in the GI series at the Baoan station, and jump variations in the GI series in 1987 and 1990 at the Hongqiling and Caizuizi stations, respectively. The alteration in the GI series of the Baoqing station occurred in 1966, and no obvious alteration was perceived in the GI series in the subsequent 46 years. Based on the alteration analysis above, runoff data from 1966 to 2012 were selected as the subsequence to calculate ecological flow at the Baoqing station. Runoff data for the periods 1956–1990, 1956–1987, and 1956–1990 were selected as the subsequence to calculate ecological flow at the Baoan, Hongqiling, and Caizuizi stations, respectively. Results of the hydrological alteration diagnosis system are shown in Table 1.
Water 10 01212 g005
Figure 6. Ecological flow calculation results: (a) Baoqing, (b) Baoan, (c) Hongqiling, and (d) Caizuizi stations.
Figure 6. Ecological flow calculation results: (a) Baoqing, (b) Baoan, (c) Hongqiling, and (d) Caizuizi stations.
Water 10 01212 g006
Table 1. Results of the hydrological variation diagnosis system.
Table 1. Results of the hydrological variation diagnosis system.
Diagnosis MethodIntra-Annual Runoff Distribution of GI (1956–2012)
Station
BaoqingBaoanHongqilingCaizuizi
Primary diagnosis Hurst coefficient0.7410.8960.7270.615
Alteration existenceYesYesYesYes
Detailed diagnosisJump diagnosisSliding F test1966
(+)
1990
(+)
1987
(−)
1990
(−)
Sliding T test1966
(−)
1998
(−)
2009
(−)
1981
(+)
Lee–Heghinan method1966
(+)
1990
(+)
1989
(+)
1990
(+)
Orderly cluster method1966
(+)
1974
(−)
1987
(+)
1990
(−)
R/S1989
(+)
1990
(−)
1979
(−)
1990
(+)
Brown–Forsythe1964(−)1990(+)1987(+)1992
(−)
Sliding run test method1992
(+)
1990
(−)
1992
(+)
1985(+)
Sliding rank test method1985
(+)
1992
(+)
1987
(+)
1991
(+)
Optimum information dichotomy1966
(−)
1989
(+)
1989
(−)
1990
(+)
Mann–Kendall1992
(+)
1992
(−)
1987
(−)
1991
(−)
BSYES1966
(+)
1989
(−)
1987
(+)
1989
(+)
Trend diagnosisTrend alteration degreeSignificant alterationSignificant alterationSignificant alterationSignificant alteration
Correlation coefficient method(+)(+)(+)(+)
Spearman(+)(+)(−)(+)
Kendall(+)(+)(+)(+)
Jump point1966199019871990
JumpComprehensive weight0.550.670.420.34
Comprehensive significance4 (+)5 (+)4 (+)4 (+)
TrendComprehensive significance3 (+)3 (+)2 (−)3 (+)
Comprehensive diagnosisAlteration form selectionJump efficiency coefficient/%45.335.152.550.3
Trend efficiency coefficient/%38.346.239.137.5
Diagnosis result1966199019871990
Note: + in the table represents a strong alteration, and − represents a weak alteration. The runoff series was altered when the Hurst coefficient was not equal to 0.5 [25].
Table 2. Ecological flow of the selected series and the whole series in the four stations in the Naolihe basin (m3/s).
Table 2. Ecological flow of the selected series and the whole series in the four stations in the Naolihe basin (m3/s).
StationMethodAprilMayJuneJuly
(1)(2)(1)(2)(1)(2)(1)(2)
BaoqingRVA10.056.8915.8111.485.676.802.291.28
MFM11.196.8517.9912.8511.957.727.394.95
ADMFFM7.645.9016.8412.558.356.556.524.35
BaoanRVA6.954.926.984.942.611.122.681.15
MFM4.323.638.436.026.034.444.423.36
ADMFFM5.324.926.985.442.612.942.682.54
HongqilingRVA8.268.6115.1112.536.613.355.690.75
MFM10.179.1816.5912.945.625.395.043.82
ADMFFM7.617.6014.7211.78.216.356.214.70
CaizuiziRVA57.0837.1830.1338.1830.1032.4717.4429.62
MFM35.2834.9462.3955.3043.6939.4433.7630.13
ADMFFM42.3138.6567.895635.7431.8529.3228.3
StationMethodAugustSeptemberOctoberNovember
(1)(2)(1)(2)(1)(2)(1)(2)
BaoqingRVA12.215.147.023.251.320.492.220.98
MFM15.4211.667.445.094.443.012.261.63
ADMFFM14.8211.659.425.53.222.033.822.5
BaoanRVA5.892.975.818.32.161.462.131.45
MFM8.376.514.564.223.913.512.162.08
ADMFFM5.895.045.815.602.162.122.131.86
HongqilingRVA10.265.275.491.374.352.672.461.03
MFM9.769.334.785.494.664.082.222.41
ADMFFM12.6311.756.345.753.923.903.122.55
CaizuiziRVA17.3635.6615.9126.4115.1819.0933.4219.10
MFM33.9933.0531.4128.6029.1126.6119.4920.74
ADMFFM39.4537.821.2222.0519.4122.3523.3322
Note: (1) is the ecological flow calculated by subsequence series selected by the hydrological alteration diagnosis system. (2) is the ecological flow calculated by the whole series. RVA is the range of variability approach method, MFM is monthly frequency method and ADMFFM is annual daily mean flow frequency method.
Table 3. Weights of ecological flow calculated by RVA, MFM, and ADMFFM based on the subsequence series and the entire series, respectively.
Table 3. Weights of ecological flow calculated by RVA, MFM, and ADMFFM based on the subsequence series and the entire series, respectively.
StationMethodAprilMayJuneJuly
(1)(2)(1)(2)(1)(2)(1)(2)
BaoqingRVA0.310.290.440.320.440.310.420.32
MFM0.250.350.320.430.300.400.240.42
ADMFFM0.440.290.240.320.260.310.340.31
BaoanRVA0.320.340.410.380.390.410.310.36
MFM0.290.340.350.390.170.430.410.39
ADMFFM0.390.350.240.390.440.430.280.35
HongqilingRVA0.310.340.340.420.290.390.310.36
MFM0.270.330.330.420.290.420.290.38
ADMFFM0.420.290.330.320.420.310.400.30
CaizuiziRVA0.340.340.420.390.310.460.380.38
MFM0.370.350.220.340.240.440.240.37
ADMFFM0.290.360.360.350.450.460.380.40
StationMethodAugustSeptemberOctoberNovember
(1)(2)(1)(2)(1)(2)(1)(2)
BaoqingRVA0.340.290.410.300.320.310.430.32
MFM0.420.400.230.390.250.420.270.43
ADMFFM0.240.300.360.300.430.310.300.33
BaoanRVA0.430.410.180.400.290.370.290.41
MFM0.310.370.410.420.410.410.420.43
ADMFFM0.260.380.410.410.300.390.290.43
HongqilingRVA0.250.420.200.390.300.420.250.33
MFM0.420.420.430.420.410.430.330.34
ADMFFM0.330.310.370.310.290.310.420.30
CaizuiziRVA0.400.390.240.330.290.380.230.31
MFM0.400.360.410.370.300.370.370.33
ADMFFM0.200.370.350.360.410.400.400.38
Note: (1) is the ecological flow calculated by subsequence series selected by the hydrological alteration diagnosis system. (2) is the ecological flow calculated by the whole series. RVA is the range of variability approach method, MFM is monthly frequency method and ADMFFM is annual daily mean flow frequency method.
Table 4. Ecological water requirement damage frequency in each station (in percent).
Table 4. Ecological water requirement damage frequency in each station (in percent).
MonthBaoqingBaoanHongqilingCaizuizi
(1)(2)(1)(2)(1)(2)(1)(2)
427.2616.8916.1414.6725.8614.5718.1214.61
514.6319.3512.2316.5732.4328.0721.8612.87
616.1517.2914.9716.1224.7934.6523.4435.32
715.8916.2314.0323.8724.6833.7424.3433.47
819.5524.5618.1022.4827.6532.8022.8023.49
912.5618.9619.5723.6330.0836.3621.4226.31
1028.6516.2312.1814.3526.7913.5722.6214.06
1128.6315.9912.7014.3016.3614.7315.9813.04
Average20.4218.1914.9918.2526.0826.0621.3221.65
Note: (1) is the ecological water requirement damage frequency by subsequence series selected by the hydrological alteration diagnosis system. (2) is the ecological water requirement damage frequency according to the whole series.

Share and Cite

MDPI and ACS Style

Xing, Z.; Wang, Y.; Gong, X.; Wu, J.; Ji, Y.; Fu, Q. Calculation of Comprehensive Ecological Flow with Weighted Multiple Methods Considering Hydrological Alteration. Water 2018, 10, 1212. https://doi.org/10.3390/w10091212

AMA Style

Xing Z, Wang Y, Gong X, Wu J, Ji Y, Fu Q. Calculation of Comprehensive Ecological Flow with Weighted Multiple Methods Considering Hydrological Alteration. Water. 2018; 10(9):1212. https://doi.org/10.3390/w10091212

Chicago/Turabian Style

Xing, Zhenxiang, Yinan Wang, Xinglong Gong, Jingyan Wu, Yi Ji, and Qiang Fu. 2018. "Calculation of Comprehensive Ecological Flow with Weighted Multiple Methods Considering Hydrological Alteration" Water 10, no. 9: 1212. https://doi.org/10.3390/w10091212

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

Article Metrics

Back to TopTop