ABSTRACT
The river hydrological regime is a key indicator to ensure the ecological function of rivers. We evaluated the changes in the hydrological regime of the Zi River by combining two methods for assessing the degree of change. Furthermore, we screened the most ecologically relevant hydraulic indicators (ERHIs) using principal component analysis (PCA) and quantified the independent impacts of human activities and climate change using the Whale Optimization Algorithm and Random Forest model (WOA-RF). The results showed that (1) 1987 was the year of the Zi River flow mutation, the overall average annual runoff in the Zi River Basin showed a decreasing trend, and the flow had four cyclic change time scales. (2) The overall hydrological variability in the Zi River Basin was 42%, and after the preferential selection of hydrological indicators, eight indicators were identified as the ERHIs of the basin, and the hydrological variability was 48.5%. The hydrological change in the basin was determined to be moderate. (3) The contribution of climate change is 28.8%, and the contribution of human activities is 71.2%. The construction of reservoirs was the most significant factor contributing to the impact of human activities, with a 65.6% contribution.
HIGHLIGHTS
The study successfully screened out the eight ERHIs and solved the problem of indicator redundancy.
It was found that the natural river flow reconstructed by the WOA-RF model applies to the Zi River Basin.
It was found that the contribution of human activities to flow change was as high as 71.2%, while that of climate change was only 28.8%.
INTRODUCTION
The hydrological regime of a river is a key indicator for ensuring the ecological functioning of the river. It maintains the sustainable development and utilization of water resources and is the basis for many organisms' survival, reproduction, and development (Sun et al. 2014). With increasing climate change and anthropogenic disturbances, globally important rivers such as the Mississippi River (Haider et al. 2025), the Nile (Fathi et al. 2019), the Yellow River (Zhou et al. 2017), and the Yangtze River (Guo et al. 2021) have all experienced different degrees of changes in their hydrological regime. The impacts of climate change and human activities are the two main drivers of river flow (Wei et al. 2020). In recent decades, with socio-economic development and population growth, intense anthropogenic disturbances have altered the hydrological cycle processes in the basin (Guan et al. 2021), which has changed the ecohydrological regime of rivers to varying degrees (Sharifi et al. 2024). Therefore, analyzing changes in hydrological regimes and distinguishing between the impacts of human activities and climate change on hydrological regimes is important for the healthy development of river hydrology.
In recent years, scholars at home and abroad have discussed the quantitative analysis of river hydrological regimes. Richter et al. (1996) proposed the Indicators of Hydrologic Alteration (IHA), which are the most widely used, and the index system contains 33 indicators divided into five groups (flow, event, frequency, duration, and rate of change) and further simplified the common 170 hydrological indicators. It uses changes in river flow and external ecological change to determine a basin's overall degree of change. Later, the IHA method combined with the range of variability approach (RVA) method for univariate and integrated hydrologic evaluation was proposed by Richter et al. (1998). Gunawardana et al. (2021) quantified the pre- and post-impacts of the transboundary Srepok River Basin hydropower development on the hydrologic cycle using the RVA in IHA. Although the IHA system is relatively well developed, there are problems of autocorrelation and overlapping information between IHA variables (Smakhtin et al. 2006). The most ecologically relevant hydrologic indicators (ERHIs) have received more and more attention from researchers. Jiang et al. (2022) screened seven indicators as ERHIs for the Weihe River using principal component analysis (PCA). Li et al. (2022) studied the Han River Basin by selecting strong variability indicators from ERHIs and diagnosing watershed variability with a strong variability index. Replacing the large number of IHA with ERHIs facilitates the study by retaining as much key information as possible from the original data and significantly reducing the redundancy of information from the large number of hydrological variables (Yang et al. 2024).
Meanwhile, for the attribution analysis of river flow changes, scholars in recent years have mainly used various hydrological models and elasticity coefficient methods to analyze the influencing factors of flow changes in different rivers. Mirjalili & Lewis (2016) proposed a new population intelligence optimization algorithm, the Whale Optimization Algorithm (WOA). Wang et al. (2023a, 2023b, 2023b) proposed a method for predicting monthly flow intervals based on the long short-term memory (LSTM) model optimized by WOA and nonparametric kernel density estimation. The Random Forest (RF) model proposed by Wang et al. (2022) for attributing flow changes in the upper Yongding River Basin provides a new method and idea for quantitatively identifying the contribution of flow changes in the basin. Soo et al. (2024) found that the RF model exhibited the best performance in flow simulation and prediction in the Klang River Basin.
Several scholars have studied the hydrological regime of the Zi River Basin from diverse perspectives. Wang et al. (2019) evaluated the hydrological regime using the IHA-RVA method. Meanwhile, Li et al. (2018) employed the cumulative distance level method and the sliding t-test to analyze trends and detect abrupt changes in precipitation and flow series over the past 55 years within the basin. Notably, most hydrological researchers have confined their analysis of flow changes to identifying trends and mutation points (Li et al. 2018; Wang et al. 2019). There is a pressing need for a more in-depth exploration of the variability of various hydrological indicators of flow and the extent of these variations. Furthermore, even when this issue is addressed, the overlap among hydrological indicators must be examined, a particularly acute problem in the Zi River Basin. Additionally, a quantitative analysis of the impacts of human activities and climate change on hydrological changes in the Zi River Basin remains an area requiring further research.
In this study, a strategy based on PCA to select ERHIs combined with the WOA-RF model is proposed to address the above problems, which effectively solves the problem of variable redundancy of traditional IHA method through PCA dimensionality reduction and improves the scientificity of hydrological variability assessment; meanwhile, combining with the WOA-optimized RF model, the generalization ability and reliability of flow attribution analysis are significantly enhanced. Compared with the traditional IHA or single RF method, the optimization and breakthrough of the methodology are achieved.
Therefore, given the limitations of previous studies, the study's main objectives are the following: (1) to analyze the characteristics of long-term series changes of hydrometeorological elements in the Zi River Basin. (2) Access the changes in the hydrological regime of Zi River before and after the sudden change based on IHA and select the indicators of ERHIs of the Zi River Basin by PCA to reduce the influence of autocorrelation information on the assessment results. (3) To reconstruct the natural flow in the Zi River Basin based on the WOA-RF model and to evaluate the impacts of climate change and human activities on the changes in the hydrological regime in the Zi River Basin. This study helps us systematically understand the characteristics of changes in the hydrological regime of the Zi River Basin. The results can provide a new approach to studying hydrological changes in the Zi River Basin and are of great reference value for watershed ecological protection and water resource management. Additionally, they serve as a guide for similar research endeavors.
STUDY AREA AND DATA
Study area
Data
The data used in this study include meteorological data and hydrological data. We collected daily flow data from the Taojiang Hydrological Station in the Zi River Basin, spanning from 1961 to 2017, and this data will be used to analyze the changing pattern of flow in the Zi River Basin. The Taojiang hydrological station is located at the downstream outlet of the Zi River Basin, and its long-term flow data were obtained from the Hunan Provincial Bureau of Hydrology and Water Resources Survey. In addition, the day-by-day meteorological data from 1961 to 2017 were obtained from five meteorological stations: Anhua, Xinhua, Wugang, Shaoyang, and Dongkou. The meteorological data were obtained from the raw data of day-by-day observations of 824 national basic meteorological stations provided by the China Meteorological Network (http://data.cma.cn/), covering indicators such as precipitation, average temperature, relative humidity, wind speed, sunshine duration, and potential evapotranspiration. The geographical locations and altitudes of all the above meteorological stations are shown in Table 1.
Description of the hydrological and meteorological stations in the study area
Station . | Station type . | Longitude (E) . | Latitude (N) . | Altitude (m) . |
---|---|---|---|---|
Taojiang | Hydrological station | 112.12 | 28.53 | 52.10 |
Xinhua | Meteorological station | 111.30 | 27.75 | 211.90 |
Anhua | 111.22 | 28.38 | 128.30 | |
Shaoyang | 111.45 | 27.18 | 311.00 | |
Wugang | 110.63 | 26.73 | 341.00 | |
Dongkou | 110.60 | 27.03 | 339.50 |
Station . | Station type . | Longitude (E) . | Latitude (N) . | Altitude (m) . |
---|---|---|---|---|
Taojiang | Hydrological station | 112.12 | 28.53 | 52.10 |
Xinhua | Meteorological station | 111.30 | 27.75 | 211.90 |
Anhua | 111.22 | 28.38 | 128.30 | |
Shaoyang | 111.45 | 27.18 | 311.00 | |
Wugang | 110.63 | 26.73 | 341.00 | |
Dongkou | 110.60 | 27.03 | 339.50 |
METHODS
Sequence trend, mutation, and periodicity tests
The Mann–Kendall (M–K) test is a typical nonparametric test that can be used to test for trends and significance in hydrological and meteorological time series (Pirnia et al. 2019). The test steps are as follows:
- (3) Arrange the time series x in reverse order, repeat (1) and (2), and produce both:where UB1 = 0, by analyzing statistical sequences UBk and UFk, the time node of mutation of sequence x can be further analyzed to show the mutation region. If the UBk and UFk curves intersect, and the intersection is between the two critical lines, then the moment corresponding to the intersection is when the mutation begins.
The cumulative anomaly method is mainly applied to hydrometeorological studies. It validates the ranking of cumulative anomalies by evaluating their dispersion Xt concerning the mean value of the data. In the case of Xt decreasing, it indicates values below the data mean, while its increasing indicates values above the mean. Fluctuation plots of their long-term series can show increases and decreases in the data. In general, the curve fluctuations reflect the timing of hydrological variations, and the mutation points usually correspond to the turning points of the curve or positions where the slope changes significantly (Wang et al. 2023a, 2023b, 2023b).








Equation (7) follows a t-distribution with . Given a threshold value
, when
, it is considered to be a hydrological mutation point.

In this equation, ‘W’ is the wavelet transform coefficients; ‘f(t)’ is the signal or square productable function; ‘a’ is the scale factor, reflecting the period length of the wavelet; ‘b’ is the translation factor, reflecting the time translation. ‘’ is the complex conjugate of ‘
’.
Evaluation of hydrological change degree
Richter et al. (1996) proposed that the IHA reflect rivers' hydrological regimes. They include indicators of flow, events, frequency, and rate of change, as shown in Table 2.
Indexes of IHA
IHA statistics group . | Characteristics . | Parameters . |
---|---|---|
Monthly water conditions (1–12) | Magnitude, timing | The mean value for each calendar month |
Extreme water conditions (13–23) | Magnitude, duration | Annual 1-, 3-, 7-, 30-, 90-day minimum and maximum; base flow index |
Timing of extreme water conditions (24–25) | Timing | Julian date of each annual 1-day minimum and maximum |
High and low pulses (26–29) | Magnitude, duration, frequency | No. and duration of low/high pulses |
Water condition changes (30–32) | Frequency, variability | No. of rises falls, and reversals |
IHA statistics group . | Characteristics . | Parameters . |
---|---|---|
Monthly water conditions (1–12) | Magnitude, timing | The mean value for each calendar month |
Extreme water conditions (13–23) | Magnitude, duration | Annual 1-, 3-, 7-, 30-, 90-day minimum and maximum; base flow index |
Timing of extreme water conditions (24–25) | Timing | Julian date of each annual 1-day minimum and maximum |
High and low pulses (26–29) | Magnitude, duration, frequency | No. and duration of low/high pulses |
Water condition changes (30–32) | Frequency, variability | No. of rises falls, and reversals |
Po denotes the frequency with which a hydrological parameter falls within the flow management target range following a hydrological change. Pe denotes the expected value, which is the range of the ecohydrological target interval. The expected value for this study is 50%. Di is a positive value, which indicates that the parameter of the hydrological sequence after variation falls within the target interval higher than the expected value; Di is a negative value, which indicates that the parameter of the hydrological sequence after variation falls within the target interval lower than the expected value. When , the hydrological parameters are considered to be highly variable (H); when
, the hydrological parameters are considered to be moderately variable (M); when
, the hydrological parameters are considered to be lowly variable (L).
Black et al. (2005) proposed another broad hydrological alteration index, the Dundee Hydrological Regime Assessment Method (DHRAM), based on the five groups of 32 parameters of IHA, to measure the degree and extent of anthropogenic alteration of hydrological regimes. DHRAM classified the mean and coefficient of variation of the five groups of 32 parameters of IHA into three impact categories: Level 1 (one impact point with a low threshold), Level 2 (two impact points with a medium threshold), and Level 3 (three impact points with a high threshold) (Table 3). Summing the mean and coefficient of variation impact categories for each group of indicators resulted in a score (score range 0–30), by which the degree of impact was assessed for each of the five classes. The higher the score, the more significant the change in the hydrological situation of the river and the greater the likelihood of damage to the river ecosystem.
Hydrological change thresholds used to assign impact points
IHA indicators . | Change in IHA group score . | ||
---|---|---|---|
A point of influence . | Two points of influence . | Three points influence . | |
1a (The first set of means) | 19.9 | 43.7 | 67.5 |
1b (The first group of CVs) | 29.4 | 97.6 | 165.7 |
2a (The second set of means) | 42.9 | 88.2 | 133.4 |
2b (The second group of CVs) | 84.5 | 122.7 | 160.8 |
3a (The third group of means) | 7.0 | 21.2 | 35.5 |
3b (The third group of CVs) | 33.4 | 50.3 | 67.3 |
4a (The fourth group of means) | 36.4 | 65.1 | 93.8 |
4b (The fourth group of CVs) | 30.5 | 76.1 | 121.6 |
5a (The fourth group of means) | 46.0 | 82.7 | 119.4 |
5b (The fifth group of CVs) | 49.1 | 79.9 | 110.6 |
IHA indicators . | Change in IHA group score . | ||
---|---|---|---|
A point of influence . | Two points of influence . | Three points influence . | |
1a (The first set of means) | 19.9 | 43.7 | 67.5 |
1b (The first group of CVs) | 29.4 | 97.6 | 165.7 |
2a (The second set of means) | 42.9 | 88.2 | 133.4 |
2b (The second group of CVs) | 84.5 | 122.7 | 160.8 |
3a (The third group of means) | 7.0 | 21.2 | 35.5 |
3b (The third group of CVs) | 33.4 | 50.3 | 67.3 |
4a (The fourth group of means) | 36.4 | 65.1 | 93.8 |
4b (The fourth group of CVs) | 30.5 | 76.1 | 121.6 |
5a (The fourth group of means) | 46.0 | 82.7 | 119.4 |
5b (The fifth group of CVs) | 49.1 | 79.9 | 110.6 |
Principal component analysis
PCA transforms a large number of correlated variables into a few uncorrelated variables through the idea of dimensionality reduction. It retains as much information as possible in the original data and the principal components (PCs). PCs are the new few composite variables, each a linear combination of the original variables with important applications in fields such as hydrology and ecology (Kholoosi et al. 2016; Yang et al. 2017). After dimensionality reduction, the indicators obtained by the IHA system are called the most ERHIs. These indicators represent many hydrological variables through fewer key hydrological variables and can reflect the most important ecologically relevant information. The ecologically relevant information reflected by the ERHIs screened in this paper is shown in Table 4. The specific process of the PCA implementation and the principle and methodology can be found in the literature (Shlens 2014).
Summary of ERHIs parameters and their ecosystem influences
ERHIs . | Ecosystem influences . |
---|---|
Median flow in March Median flow in May Median flow in September |
|
Minimum 30-day flow Maximum 3-day flow |
|
Number of low-flow pulse occurrences Number of high-flow pulse occurrences |
|
Average rate of increase |
|
ERHIs . | Ecosystem influences . |
---|---|
Median flow in March Median flow in May Median flow in September |
|
Minimum 30-day flow Maximum 3-day flow |
|
Number of low-flow pulse occurrences Number of high-flow pulse occurrences |
|
Average rate of increase |
|
Reconstruction of natural flow
The WOA consists of three phases: encircling predation, bubble net predation, and search predation.
The RF model and the Whale algorithm optimization allowed for a more accurate reconstruction of natural flow. The WOA-RF model in this study includes precipitation, mean air temperature, relative humidity, wind speed, sunshine, and potential evapotranspiration, and the output factor is monthly flow. Based on artificial parameter optimization, the monthly mean natural flow for 1961–2017, considering only the effect of climate change, was recovered using the monthly mean flow before the mutation and the precipitation, mean air temperature, relative humidity, wind speed, insolation, and potential evapotranspiration after the mutation.
In the above equation, and
denote the observed and simulated flow values (m3/s), respectively, and
and
denote the average observed and simulated flow values (m3/s), respectively.
Quantitative evaluation of flow changes







RESULTS
Analysis of hydrometeorological elements
(a) Interannual variability of runoff versus precipitation. (b) Interannual variability of runoff versus potential evaporation.
(a) Interannual variability of runoff versus precipitation. (b) Interannual variability of runoff versus potential evaporation.
Result of the flow mutation test
Hydrological stations . | Mutation point . | |||
---|---|---|---|---|
M–K test . | Cumulative anomaly method . | Sliding t-test . | Mutation year . | |
Taojiang station | 1961, 1963, 1968, 1972, 1975, 1978, 1984, 1987, 2004 | 1977, 1987, 1989, 2006 | 1987 | 1987 |
Hydrological stations . | Mutation point . | |||
---|---|---|---|---|
M–K test . | Cumulative anomaly method . | Sliding t-test . | Mutation year . | |
Taojiang station | 1961, 1963, 1968, 1972, 1975, 1978, 1984, 1987, 2004 | 1977, 1987, 1989, 2006 | 1987 | 1987 |
(a) M–K mutation test. (b) The cumulative anomaly method. (c) Sliding t-test method.
(a) M–K mutation test. (b) The cumulative anomaly method. (c) Sliding t-test method.
Variation of the flow cycle in the Zi River Basin: (a) wavelet contour and (b) wavelet variance.
Variation of the flow cycle in the Zi River Basin: (a) wavelet contour and (b) wavelet variance.
Analysis of overall hydrological changes
Hydrological changes by group and overall
Hydrographic stations . | Degree of change for each group (%) . | Overall hydrological variability (%) . | ||||
---|---|---|---|---|---|---|
Group 1 . | Group 2 . | Group 3 . | Group 4 . | Group 5 . | ||
Taojiang station | 45 (M) | 35 (M) | 26 (L) | 42 (M) | 56 (M) | 42 (M) |
Hydrographic stations . | Degree of change for each group (%) . | Overall hydrological variability (%) . | ||||
---|---|---|---|---|---|---|
Group 1 . | Group 2 . | Group 3 . | Group 4 . | Group 5 . | ||
Taojiang station | 45 (M) | 35 (M) | 26 (L) | 42 (M) | 56 (M) | 42 (M) |
The degree of hydrological change was derived using the DHRAM method (Table 7). From the study results, the total DHRAM score of Taogang station is 4, and the degree of change is 3, which is low. Considering the 32 hydrological change indicators in the IHA, the overall change in the degree group is more moderate than low change, and combining the results of the two methods, it is determined that the degree of change in the hydrological situation of the reservoir is moderate.
Degree of changes in hydrological indicator composites
Hydrographic stations . | Group . | Change ratio (%) . | Influence degree . | Total score . | ||
---|---|---|---|---|---|---|
Mean . | CVs . | Mean . | CVs . | |||
Taojiang station | 1 | 16.6 | 28 | 0 | 0 | 4(3) |
2 | 8.7 | 28.3 | 0 | 0 | ||
3 | 375.6 | 27.9 | 3 | 0 | ||
4 | 34.3 | 30.5 | 0 | 1 | ||
5 | 12.7 | 34.3 | 0 | 0 |
Hydrographic stations . | Group . | Change ratio (%) . | Influence degree . | Total score . | ||
---|---|---|---|---|---|---|
Mean . | CVs . | Mean . | CVs . | |||
Taojiang station | 1 | 16.6 | 28 | 0 | 0 | 4(3) |
2 | 8.7 | 28.3 | 0 | 0 | ||
3 | 375.6 | 27.9 | 3 | 0 | ||
4 | 34.3 | 30.5 | 0 | 1 | ||
5 | 12.7 | 34.3 | 0 | 0 |
Note: CV is the coefficient of variation. Values in parentheses in the ‘total score’ list indicate the level of ecosystem risk due to hydrological change.
Selection of key ERHIs and characterization of their evolution
(a) Eigenvalues and cumulative contribution of each component. (b) Heatmap of factor loadings for the first nine PCs that satisfy the conditions.
(a) Eigenvalues and cumulative contribution of each component. (b) Heatmap of factor loadings for the first nine PCs that satisfy the conditions.
(a) Median flow in March. (b) Median flow in May. (c) Median flow in September. (d) Minimum 30-day flow. (e) Maximum 3-day flow. (f) Number of low-flow pulse occurrences. (g) Number of high-flow pulse occurrences. (h) Average rate of increase. (The red line in the graph indicates a linear increasing trend, the green line indicates a linear decreasing trend, and the yellow line indicates a linear steady trend.)
(a) Median flow in March. (b) Median flow in May. (c) Median flow in September. (d) Minimum 30-day flow. (e) Maximum 3-day flow. (f) Number of low-flow pulse occurrences. (g) Number of high-flow pulse occurrences. (h) Average rate of increase. (The red line in the graph indicates a linear increasing trend, the green line indicates a linear decreasing trend, and the yellow line indicates a linear steady trend.)
RF model and attribution analysis
(a) Comparison of simulated and measured monthly mean flow values at the Taojiang station for 1961–1972. (b) Correlation plot between model and actual for the training period. (c) Comparison of simulated and actual monthly mean flow values at the Taojiang station for 1973–1986. (d) Correlation plots of model versus actual for the validation period.
(a) Comparison of simulated and measured monthly mean flow values at the Taojiang station for 1961–1972. (b) Correlation plot between model and actual for the training period. (c) Comparison of simulated and actual monthly mean flow values at the Taojiang station for 1973–1986. (d) Correlation plots of model versus actual for the validation period.






Simulation results of the annual mean flow at the Taojiang station.
(a) Eigenvalues and cumulative contribution of each component. (b) Heatmap of factor loadings for the first two PCs that satisfy the conditions.
(a) Eigenvalues and cumulative contribution of each component. (b) Heatmap of factor loadings for the first two PCs that satisfy the conditions.
DISCUSSION
Ecological impacts of changes in hydrological regimes
This study used PCA to screen eight ERHIs in the Zi River Basin. The eight metrics are median flow in March, median flow in May, median flow in September, minimum 30-day flow, maximum 3-day flow, number of low-flow pulse occurrences, number of high-flow pulse occurrences, and average rate of increase. Median monthly flows (median flows in March, May, and September) mainly affect habitat availability for riverine aquatic organisms as well as water affecting water temperature, oxygen levels, photosynthesis in the water column (CONSERVANCY T N 2009); minimum 30-day flow, maximum 3-day flow are related to channel morphology and structure and natural habitat conditions (CONSERVANCY T N 2009); the number of low-flow pulse occurrences, and the number of high-flow pulse occurrences mainly affects aquatic organisms habitat availability in the river rookery as well as the relationship between the river and the nutrient and organic matter exchanges between the river and the rambler (CONSERVANCY T N 2009); The average rate of increase, in turn, correlates with the trapping (rising numbers) of organisms on islands and flood plains (CONSERVANCY T N 2009). Gao et al. (2009) found some variation in the results of analyses across watersheds and datasets. However, it tended to include three common elements: at least one monthly mean flow statistic, at least two extreme event statistics representing high and low extremes, and at least one statistic correlating the frequency of low and high flows. It agrees with the eight ERHIs initially selected for this paper.
Impact of human activities on changes in the hydrological regime
Previous studies have established that human activities are the primary factor influencing changes in the hydrological regime of the basin. These activities directly or indirectly affect the quality, quantity, and flow processes through the construction of hydraulic infrastructure and modifications to subsurface conditions (Fan et al. 2023). Li et al. (2018) excluding evapotranspiration from their analyses, quantitatively assessed the impacts of precipitation and human activities on flow changes during the variability period in the Zi River Basin. They concluded that human activities were the predominant contributors to these changes. Further research has underscored that the alterations in flow within the Zi River Basin are primarily attributable to human interventions, including reservoir construction, urbanization, artificial water consumption, and other integrated factors (Li 2018). While the present study's findings align with these conclusions, notable differences may stem from the exclusion of evapotranspiration in prior studies, variations in the definition of base and variant periods, and differences in the data series and methodologies employed. Nonetheless, the overarching consensus remains that human activities are the dominant factor driving flow variability. The downstream of the Zi River Basin in Yiyang City has constructed terrace power stations, such as Zhexi and Majitang (Li et al. 2018). Under the effect of the reservoirs and dams, the average flow in May of the flood season was reduced, and the average flow in March of the dry season was slightly increased. The analyses in the ERHIs section of this study are in agreement. After the construction of a reservoir on a river, there is a tendency for the flow volume to decrease significantly, especially when it is first put into operation. However, this tendency will gradually slow down or rebound after regular operation. Therefore, the direct effect of dam construction on total flow is limited (Chen et al. 2021). However, in the 1990s, the Zhexi Reservoir needed more flood control capacity and high discharge. The Zi River Basin often suffered from floods, which, to some extent, increased Zi River flow (Xiong & Liu 1999).
(a) Land use changes over time in the Zi River Basin. (b) Land use changes before and after reservoir construction.
(a) Land use changes over time in the Zi River Basin. (b) Land use changes before and after reservoir construction.
Applicability and limitations
In this study, we propose a structured and robust framework for assessing the hydrological regime change and its ecological impacts. It can be adapted to a wide range of watersheds with different hydrological, climatic, and anthropogenic disturbances. PCA can be applied across watersheds to screen for ERHIs. Still, initial hydrological indicators need to be selected in accordance with the characteristics of the watersheds, e.g., in arid zones, priority should be given to the extreme flow indicators. The WOA-RF model can effectively differentiate the effects of human activities and climate change on flow and is applicable to basins with different climatic conditions. However, the implementation of the WOA-RF model needs to pay attention to the availability of data, the main drivers, and the model's calibration. Alternative methods or small-scale datasets can be explored for data-scarce basins, and basins dominated by climate change should incorporate more climatic variables. The model should be recalibrated so that each basin can adapt to its hydrological processes.
Sobol sensitivity analysis (Features 1–6 represent input factors as precipitation, mean air temperature, relative humidity, wind speed, sunshine, and potential evapotranspiration, respectively).
Sobol sensitivity analysis (Features 1–6 represent input factors as precipitation, mean air temperature, relative humidity, wind speed, sunshine, and potential evapotranspiration, respectively).
CONCLUSION
In order to quantitatively analyze the characteristics of hydrological regime changes in the Zi River Basin and the degree of influence of climate change and human activities on the hydrological regime changes in the Zi River Basin, we started with the inter-annual characteristics of hydrological regime changes. On this basis, the key ERHIs were selected to reduce the influence of redundant information on the hydrological regime change. The WOA-RF model was also used to quantitatively separate the degree of changes in the flow of the Zi River due to human activities and climate change. The main conclusions are as follows:
(1) The mutation results show that the flow in the Zi River Basin underwent a mutation in 1987. The overall mean annual flow and mean annual precipitation in the basin showed a decreasing trend, and the potential evaporation showed an increasing trend. The results of wavelet analyses indicated significant cycle changes on time scales of 29 and 5 years, with average cycles of 25 and 4 years.
(2) Combining the IHA-RVA and DHRAM methods, the 32 ecohydrological indicators in the Zi River Basin were mainly medium and low changes, with an overall change degree of 42%, and the total DHRAM score was 4. After the hydrological indicators were preferred, eight ERHIs were identified, with a combined change degree of 48.5%, a medium change. Compared with the traditional RVA method and DHRAM method, the degree of change of the hydrological regime of Zi River is more intuitively shown.
(3) The results of the WOA-RF model showed that climate change caused a downward trend in flow, while human activities caused an upward trend in flow. Meanwhile, the results of attribution analysis showed that the contribution of human activities and climate change to the change of flow of the capital water was 71.2 and 28.8%, respectively, and that human activities were the main factor affecting the change of flow of the capital water. The construction of reservoirs was the most significant factor contributing to the impact of human activities, with a 65.6% contribution.
ACKNOWLEDGEMENTS
The authors thank their brothers at North China University of Water Resources and Electric Power for their comments and help with this study.
FUNDING
This study was supported by Basic Research Project of Key Scientific Research Projects of Colleges and Universities of Henan Province (24ZX007).
AUTHOR CONTRIBUTIONS
X.Z. conceived the study and wrote the first draft, X.Z. and X.B. collected the data and performed the analysis, L.H. and X.J. analyzed the methodology, W.G. and H.W. supervised the paper, and all authors provided comments and assistance on the first few versions of the manuscript. All authors read and approved the final manuscript.
DATA AVAILABILITY STATEMENT
Data cannot be made publicly available; readers should contact the corresponding author for details.
CONFLICT OF INTEREST
The authors declare there is no conflict.