With the advancement of hydraulic engineering processes, surface water resources are encountering unprecedented crises that threaten social and ecological sustainability. This study analyzed the response changes and potential correlations between river conditions and climate change in the Dongting Lake (DTL) basin using gray correlation analysis and the LSTM model. Additionally, the InVEST model was employed to further assess the distribution of habitat quality within the basin. The results of the study showed that (1) the range of variation of environmental flows in each tributary became significantly narrower and the frequency of extreme hydrological conditions was reduced in the period after the abrupt hydrological change. (2) The hydrological status had a higher degree of correlation with the rainfall index, and the performance is more prominent on the Xiang River, and the influence of temperature is less. (3) Except in the Xiang River Basin, human activities are the dominant factor of runoff variation, and the annual contribution is generally higher than 50%. (4) Habitat quality in the ‘high’ distribution was the most guarded and always accounted for more than 60% of the basin, while habitat quality in the east-central part of the DTL basin and around the DTL was generally low.

  • Potential linkages between climate change and hydrological regimes.

  • Differences in response to climate change between multiple tributaries and confluences in the basin.

  • Habitat impacts of changes in climate and hydrological regimes.

Necessary hydrological processes within watersheds are one of the most important factors in maintaining the rich diversity of ecosystems, and similarly, hydrological processes within rivers play an essential role in protecting ecological functioning (Sabater et al. 2018; Keith et al. 2022). However, over the last century, extreme climate change and large-scale human-industrial advances have had significant impacts on many rivers and even watersheds globally (Wang et al. 2020). Climate variability and human activities affect the conditions of species in both terrestrial and aquatic ecosystems, not only directly by altering the habitat conditions of species but also indirectly by influencing ecosystem functioning through the biological chain (Zolfagharpour et al. 2022). Therefore, the study of hydrological processes in rivers will be important for the maintenance of water ecology and ecosystem diversity in watersheds (Zhang et al. 2023).

The Brisbane Declaration of 2007 defines environmental flows as the quantity, timing, and quality of water needed to sustain freshwater ecosystems and the livelihoods of the plants and animals that depend on them. Existing studies have shown that river flow regimes play an important role in maintaining the biodiversity, ecological integrity, and restoration of ecosystem services of rivers and the surrounding region, which is more evident in wetlands (Mao et al. 2023). Therefore, many local governments and water management agencies have invested considerable effort in environmental flow releases, such as manipulating dams, reservoirs, or other water management facilities to maintain annual or seasonal flow requirements (Cid et al. 2022; Thompson et al. 2022). Climate change has always been an important factor affecting surface runoff, especially in high-temperature and rainfall zones. In large catchment areas, the hydrological impact of surface runoff has long been a hot research-topic, with researchers worldwide examining the varying factors of surface water from different perspectives (Oliveira & Quaresma 2017). For example, Klotz et al. (2022) simulated changes in surface runoff under different scenarios using hydrological models, providing a visual comparison of the varying effects of various factors. Machine learning has also demonstrated superior performance in simulating natural runoff and restoring runoff changes under natural climatic factors, thereby facilitating comparative analysis of runoff variations (Ayzel & Heistermann 2021). Singh et al. (2023) developed an integrated multi-model framework based on various machine-learning models to reconstruct gridded runoff in Europe over the past 500 years. Previous studies have delved into the mechanisms underlying the interplay between climate and river flow dynamics, yet there remains a gap in understanding how subsequent alterations in watershed habitats correspond to these climate–flow dynamics. Hence, there is a pressing need for further research pertaining to watershed habitats to address the management requirements associated with such changes (Lei & Middleton 2021).

Under the dual impact of environmental changes (surface water, climate, etc.) and human activities, the quality of surface habitats is facing unprecedented challenges (Zhang et al. 2018; Wang et al. 2023). In response to these unknown risks, efforts have been made to assess the ecological and habitat quality of the surface environment through various approaches (Sepulveda et al. 2023). Yuan et al. (2021) evaluated the ecology of the Dongting Lake (DTL) basin using the remote sensing-based ecological index (RSEI) and identified potential influencing factors of ecological changes within the basin. Meanwhile, model-based methods have gained popularity in recent years due to their ease of operation, strong visualization capabilities, and comprehensive theoretical frameworks. The InVEST model is one such example that has been widely adopted (Bastos et al. 2023). Lei et al. (2022) utilized the InVEST model to assess the habitat quality of Hainan Island and analyzed the impact and trend changes of human engineering activities on the island's habitat quality. The DTL basin in China is one of the 200 globally prioritized ecological protection areas and a crucial stopover and wintering ground for migratory birds (Geng et al. 2021a, 2021b; Zhu et al. 2021). Studies by Qi et al. (2023) on the net primary productivity of vegetation in the DTL basin have revealed that some regions within the basin are facing ecological degradation crises. To address these ecological issues, a series of ecological restoration projects have been implemented in the basin. However, due to a lack of regional studies on hydrological and ecological differences, there are significant variations in the effectiveness of ecological restoration across different regions (Wang et al. 2021). Therefore, it is necessary to assess the habitat quality of the DTL basin and investigate the ecological and hydrological responses of its tributary basins. This is crucial for ecological protection and restoration efforts within the basin and holds significant implications for global biodiversity conservation. By conducting such assessments and studies, we can gain insights into the current state of the basin's ecology, identify key areas for restoration, and develop targeted strategies to enhance the overall health and resilience of the DTL basin ecosystem (Cao et al. 2023; Luo et al. 2023).

The aim of this paper is to link the changes in the main rivers in DTL with climate change in the basin to analyze the changes in habitat quality in the basin, and finally to compare and analyze the relationship between the flow regime of the rivers and climate change based on previous studies. To this end, we did the following: (1) analyzed the environmental flow distribution characteristics of the four major rivers in the DTL basin, (2) calculated the distribution of meteorological inputs in the basin based on the RClimDex model, and explored the potential linkages with the four rivers, (3) attributed the changes in runoff from the four rivers to the four rivers on a monthly scale using the long short-term memory (LSTM) model, and (4) assessed the habitat quality of DTL basin based on the InVEST model. The results of the study can provide a reference and theoretical basis for surface water resource management and habitat protection in the DTL basin.

Study area

The DTL basin is located at latitude 24°39′–30°24′ N and longitude 107°16′–114°14′ E (Figure 1). It is situated in the subtropical monsoon humid climate with mild seasons and is a humid zone. The basin has abundant rainfall, with a multi-year average precipitation of 1,437 mm, but the rainfall is unevenly distributed within the year, with rainfall concentrated in March–September, accounting for about 60% (Wu & Zheng 2020). The average temperature over the years is 16.0–17.4 °C, with the highest and lowest temperatures occurring in July and January, and the frost-free period lasts between 255 and 311 days, which is suitable for crop growth. The geomorphology of the basin is dominated by plains, hills, and mountains, with the overall elevation of the central and western regions higher and most of them distributed as hills and mountains, and the northeastern and central-eastern regions lower and more as plains, with the topography of the basin presenting a northwesterly trend, and the overall trend of the elevation decreasing from west to east (Cheng et al. 2018; Zhu et al. 2019). The water system in the DTL basin is dominated by the four rivers of DTL, of which the tributary systems are well developed; and the four rivers, as a sub-basin of DTL, are an important influencing factor affecting the water resources of DTL as well as the lower reaches of the Yangtze River (Geng et al. 2021a, 2021b).
Figure 1

Geographic overview of the study area and distribution of major tributaries.

Figure 1

Geographic overview of the study area and distribution of major tributaries.

Close modal

Data source

For the hydrological data, the measured daily flows of the Yangtze River, the Four Rivers (Li, Yuan, Zi, and Xiang rivers), and the representative hydrological stations of DTL were collected from 1990 to 2020, and were provided by the Yangtze River Water Conservancy Commission (Table 1). Chenglingji (CLJ) Hydrological Station is located at the outlet section of DTL, about 4 km from the confluence of the Yangtze River, and controls the changes in the hydrological situation in DTL. Shimen (SM), Taoyuan (TY), Taojiang (TJ), and Xiangtan (XT) hydrological stations are the last hydrological stations before their respective rivers flow into DTL, representing the hydrological situation of the confluence in the DTL basin. By analyzing the five representative hydrological stations, the evolution of the hydrological situation in the DTL basin can be visualized and accurately reflected.

Table 1

Hydrological station data in the study area

Dongting LakeLi RiverYuan RiverZi RiverXiang River
Hydrological stations Chenglingji Shimen Taoyuan Taojiang Xiangtan 
Shorthand CLJ SM TY TJ XT 
Longitude (E°) 113.12 115.41 111.29 112.1 112.9 
Latitude (N°) 29.41 19.10 28.54 28.5 28.16 
Data period 1961–2019 1961–2019 1961–2019 1961–2019 1961–2019 
Dongting LakeLi RiverYuan RiverZi RiverXiang River
Hydrological stations Chenglingji Shimen Taoyuan Taojiang Xiangtan 
Shorthand CLJ SM TY TJ XT 
Longitude (E°) 113.12 115.41 111.29 112.1 112.9 
Latitude (N°) 29.41 19.10 28.54 28.5 28.16 
Data period 1961–2019 1961–2019 1961–2019 1961–2019 1961–2019 

For meteorological data, we selected 38 representative meteorological stations in the DTL basin. We collected day-by-day meteorological data of these weather stations, including precipitation, daily average temperature, maximum temperature, minimum temperature, barometric pressure, corresponding humidity, and solar insolation time; the meteorological data were obtained from the China Meteorological Data Network (http://data.cma.cn/). Due to the rare occurrence of missing or abnormal values for certain dates at a few meteorological stations, we employed the linear interpolation method to fill in the missing parts. Additionally, remote-sensing raster data encompassing meteorological information, digital elevation models (DEMs), and land-use data were sourced from the Geospatial Data Cloud (https://www.gscloud.cn/) and the Resource and Environment Science and Data Center of the Chinese Academy of Sciences (https://www.resdc.cn/) (Figure 2).
Figure 2

Flowchart of the methodology used in this study.

Figure 2

Flowchart of the methodology used in this study.

Close modal

TFPW-MK + MWP trend mutation test

The Mann–Kendall (M-K) trend test method is commonly used to analyze trends in hydrological series and has the advantage that it does not depend on outliers and does not need to satisfy a specific probability distribution. However, the method is affected by the autocorrelation of the original data itself and needs to be improved. In this study, the Trend-Free Pre-Whitening (TFPW) test was used to improve it (Yue et al. 2002). The specific calculation steps of the TFPW-MK method are as follows:
(1)
(2)
(3)
(4)
where ; is the linear trend of the series to be measured; denotes a series without a trend term; is the correlation coefficient; denotes a series with the autocorrelation term removed; is the compensation term; and is a new series without autocorrelation.
(5)
When n ≥ 10, S approximately obeys a positive Tai distribution with mean E(S) and variance var(S):
(6)
where n is the number of data points, ti is the number of connections for the ith value, and m is the number of connection groups.
The standardized test statistic ZMK is:
(7)

Bilateral tests are used here. At the α-significant level, if , the sequence negates the original hypothesis of no trend. The sequence has an increasing or decreasing trend; conversely, the sequence has no increasing or decreasing trend.

For the M-K mutation test, the original time-series statistic UF is calculated at the same time as its inverse sequence column UB (UF = UB), and if an intersection of the two UF and UB curves is observed in the given confidence-level α, it is considered that a mutation occurs there.

Environmental flow components

The environmental flow components (EFCs) are the flow and its processes required to maintain the river's ecosystem and are described by the dry flow, the extra-dry flow, the high flow pulse, the minor flood, and the major flood, a total of 34 ecologically significant indicators in five groups (Table 2). The threshold delineation of flow events for relevant environmental flows is detailed in Peres & Cancelliere (2016).

Table 2

Indicators used in environmental flow

GroupsIndexEcological significance
Low water flow (1–12) Dry water flow by month Maintaining the water table in flood plains, providing drinking water for terrestrial animals, etc. 
Extreme low water flow (13–16) Peak, duration, time of occurrence, and frequency Expansion of plant species in the floodplain, prevention of invasive alien species, etc. 
High flow pulses (17–22) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Shaping the physical characteristics of the river channel, maintaining normal water-quality conditions, ensuring suitable salinity conditions in the estuary, etc. 
Minor flood (23–28) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Provides clues to fish migration and spawning, recharges floodplain levels, controls population structure and distribution of floodplain plants, etc. 
Great flood (29–34) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Maintaining the balance of species in aquatic and riparian communities; facilitating material exchange in channels and floodplains, etc. 
GroupsIndexEcological significance
Low water flow (1–12) Dry water flow by month Maintaining the water table in flood plains, providing drinking water for terrestrial animals, etc. 
Extreme low water flow (13–16) Peak, duration, time of occurrence, and frequency Expansion of plant species in the floodplain, prevention of invasive alien species, etc. 
High flow pulses (17–22) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Shaping the physical characteristics of the river channel, maintaining normal water-quality conditions, ensuring suitable salinity conditions in the estuary, etc. 
Minor flood (23–28) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Provides clues to fish migration and spawning, recharges floodplain levels, controls population structure and distribution of floodplain plants, etc. 
Great flood (29–34) Peak, duration, time of occurrence, frequency of occurrence, and rate of rise and fall Maintaining the balance of species in aquatic and riparian communities; facilitating material exchange in channels and floodplains, etc. 

Note: Numbers 1–34 represent indicators.

Calculation of meteorological index

The RClimDex model is a source of Fortran programs written in R for different indicator calculations. First, for data quality checking, the model checks for invalid data in the dataset, including cases of negative precipitation, minimum temperature greater than or equal to the maximum temperature, and so on. Then a series of statistical calculations are performed to determine 27 meteorological indicators including 11 extreme precipitation indices and 16 extreme temperature indices. In this study, we selected eight extreme precipitation indices (http://etccdi.pacificclimate.org/) and four representative extreme temperature indices as representatives for calculation and analysis (Table 3). More detailed information about meteorological indices can be found in the study by dos Santos et al. (2011).

Table 3

Summary information on meteorological indicators

IndicesDescriptive nameDefineUnit
TN10p Cool nights Number of days with minimum temperature <10% quantile day 
TX10p Cool days Number of days with maximum temperature <10% quantile day 
TN90p Warm nights Number of days with minimum temperature >90% quantile day 
TX90p Warm days Number of days with maximum temperatures >90% quantile day 
Rx1day Monthly maximum 1-day precipitation Annual maximum daily precipitation mm 
Rx5day Monthly maximum consecutive 1-day precipitation Maximum annual precipitation for 5 consecutive days mm 
R95p Very wet days Cumulative annual precipitation with daily precipitation >95% quantile mm 
R10mm Number of heavy precipitation days Total number of days with annual daily precipitation ≥10 mm day 
CDD Consecutive dry days Maximum number of consecutive days with visual precipitation <1 mm day 
CWD Consecutive wet days Maximum number of consecutive days with daily water loss of ≥1 mm day 
PRCPTOT Annual total precipitation in wet days Cumulative annual precipitation with daily precipitation > 1 mm mm 
SDII Simple daily intensity index Ratio of total precipitation ≥1 mm to number of days mm 
IndicesDescriptive nameDefineUnit
TN10p Cool nights Number of days with minimum temperature <10% quantile day 
TX10p Cool days Number of days with maximum temperature <10% quantile day 
TN90p Warm nights Number of days with minimum temperature >90% quantile day 
TX90p Warm days Number of days with maximum temperatures >90% quantile day 
Rx1day Monthly maximum 1-day precipitation Annual maximum daily precipitation mm 
Rx5day Monthly maximum consecutive 1-day precipitation Maximum annual precipitation for 5 consecutive days mm 
R95p Very wet days Cumulative annual precipitation with daily precipitation >95% quantile mm 
R10mm Number of heavy precipitation days Total number of days with annual daily precipitation ≥10 mm day 
CDD Consecutive dry days Maximum number of consecutive days with visual precipitation <1 mm day 
CWD Consecutive wet days Maximum number of consecutive days with daily water loss of ≥1 mm day 
PRCPTOT Annual total precipitation in wet days Cumulative annual precipitation with daily precipitation > 1 mm mm 
SDII Simple daily intensity index Ratio of total precipitation ≥1 mm to number of days mm 

Gray correlation analysis

Gray correlation analysis was first proposed by Professor Deng Julong, an expert in cybernetics, and a series of gray system-related theories were established after that. Gray correlation analysis is one of the important methods in gray system theory, which is a multi-factor analysis method, and research focuses on the main correlation relationship between the factors (Dong et al. 2018). The calculation of gray correlation is analyzed as follows.

Firstly, the reference sequence and comparison sequence in all data are determined, and the reference sequence is set to be:
(8)
Let the comparison sequence be:
(9)
The original sequence is subsequently dimensionless:
(10)
Then find the absolute value of the difference between the reference sequence and the comparison sequence, thus forming the absolute difference sequence:
(11)
Gray correlation coefficients between variables are calculated based on absolute difference series:
(12)
Finally, the mean gray correlation value of the elements between each comparison sequence and the reference sequence is calculated to reflect the correlation between the two sequences of comparison, i.e., the correlation sought:
(13)

LSTM modeling and attribution analysis

The LSTM is a type of deep-learning neural network that has been used by many scholars in recent years for interpolation and prediction of hydrometeorological data. In this study, meteorological data from various watersheds were utilized as inputs to generate five meteorological–flow models. These models were then used to reconstruct the flow processes of the Four Rivers and the DTL region solely under the influence of climate change. This reconstruction allowed for the assessment of the impact of external stressors on the runoff within the watershed. For the development of the watershed models, measured data from 1961 to 1990 were utilized for model calibration, while measured data from 1990 to the year of abrupt change were used to validate the model's accuracy. The performance of the LSTM model was comprehensively evaluated using the Nash–Sutcliffe efficiency (NSE), root mean square error (RMSE), and the coefficient of determination (R2).

Based on the input data characteristics of the LSTM model, this study classifies the influencing factors of flow change variability into climate change and non-climatic factors and assumes that the non-climatic factors are mainly represented by human activities (Wang et al. 2020). Combined with the reconstruction results of the LSTM model, the contributions of climate change and human activities to runoff alteration on annual and monthly scales can be identified, and the effects of model residuals on the results are considered:
(14)
(15)
(16)
(17)
where ΔRA is the difference in runoff volume before and after the mutation year; RAobs,pre-mut and RAobs,post-mut denote the measured flow conditions in the period before and after the mutation year, respectively; RAsim,pre-mut and RAsim,post-mut denote the simulated flow conditions in the period before and after the mutation year, respectively; ΔRAhum and ΔRAcli denote the absolute contribution of human activities and climate change, and their respective impact weights ηhum and ηcli, respectively; ε represents the residuals of the LSTM model in the base period; RAcli denotes the absolute contributions of human activities and climate change, and their respective impact weights ηhum and ηcli, respectively.

InVEST model

The habitat quality module of the InVEST model can calculate habitat quality through the sensitivity of landscape types and the intensity of external threats to landscape types for ecosystem service function assessment and trade-off (Yang et al. 2023). The method is easier for obtaining the required data than other habitat quality evaluation models and has higher reliability than methods such as species surveys, which can evaluate the characteristics of habitat quality changes and spatial and temporal distributions. The calculation formula is as follows:
(18)
where Qxj is the habitat quality index of raster x in land use j; Hj is the habitat suitability of land-use data j; z is the scale constant, z = 2.5; k is the half-saturation constant, which is taken as half of the maximum value obtained from one run of Qxj; and Dxj is the habitat degradation index of raster x in land use j, which indicates the degree of degradation of the habitat after the disturbance, and is calculated by the following formula:
(19)
where r is the threat factor; R is the number of threat factors; y is the number of rasters of threat factor r; x is the number of habitat rasters; ωr is the weight of threat factor r, which indicates the interference of a certain threat factor on all habitat types, and is taken as 0–1; ry is the number of threat factors on the raster cells; βx is the accessibility level of raster x, and is taken as 0–1, which, when closer to 1, indicates that it can be reached more easily; Sjr is the sensitivity of land-use type j to threat factor r, which is taken as 0–1; and irxy is the interference degree of threat factor r on habitat raster x, and is calculated as follows:
(20)
(21)
where dxy is the distance between x and y; and is the maximum influence distance of the threat factor.

The model run should be set up with parameters that are selected considering the InVEST model user's guide and determined considering the results of research in various ecological fields. For threat sources and their weights and the maximum influence distance of threat factors, the parameter settings are shown in Table 4; for the degree of sensitivity of various habitats to threat factors, the parameter settings are shown in Table 5. The higher the sensitivity of natural habitat quality to the threat factors, the greater the influence of this threat source on habitat quality (Guan et al. 2022).

Table 4

Threat sources and their weights and maximum impact distance

Threat factorMaximum impact distance (km)Relative weightType of space recession
Plowed land 0.7 Linear 
Highway 0.5 Linear 
Body of water 0.6 Linear 
Town 10 Exponential 
Rural settlements 0.6 Exponential 
Industrial site 0.6 Exponential 
Threat factorMaximum impact distance (km)Relative weightType of space recession
Plowed land 0.7 Linear 
Highway 0.5 Linear 
Body of water 0.6 Linear 
Town 10 Exponential 
Rural settlements 0.6 Exponential 
Industrial site 0.6 Exponential 
Table 5

Habitat suitability and its relative sensitivity to threat sources

Land cover typeHabitat suitabilitySensitivity of habitat types to stressors
Plowed landHighwayBody of waterTownRural settlementsIndustrial site
Paddy field 0.4 0.3 0.2 0.6 0.5 0.4 0.7 
Arid 0.3 0.3 0.2 0.3 0.5 0.4 0.7 
Woodland 0.6 0.8 0.5 0.9 0.7 0.8 
Low wood 0.8 0.4 0.7 0.6 0.6 0.4 0.7 
Open woodland 0.7 0.6 0.6 0.4 0.8 0.6 0.6 
Other forest land 0.5 0.5 0.5 0.4 0.7 0.5 0.5 
High-cover grassland 0.6 0.5 0.4 0.6 0.7 0.6 0.6 
Medium-cover grassland 0.5 0.6 0.5 0.5 0.6 0.7 0.5 
Low-cover grassland 0.3 0.5 0.3 0.5 0.8 0.5 0.4 
Body of water 0.8 0.5 0.4 0.8 0.6 0.2 
Town 
Rural settlements 
Industrial site 
Unutilized land 0.6 0.5 0.6 0.5 0.8 0.6 0.4 
Land cover typeHabitat suitabilitySensitivity of habitat types to stressors
Plowed landHighwayBody of waterTownRural settlementsIndustrial site
Paddy field 0.4 0.3 0.2 0.6 0.5 0.4 0.7 
Arid 0.3 0.3 0.2 0.3 0.5 0.4 0.7 
Woodland 0.6 0.8 0.5 0.9 0.7 0.8 
Low wood 0.8 0.4 0.7 0.6 0.6 0.4 0.7 
Open woodland 0.7 0.6 0.6 0.4 0.8 0.6 0.6 
Other forest land 0.5 0.5 0.5 0.4 0.7 0.5 0.5 
High-cover grassland 0.6 0.5 0.4 0.6 0.7 0.6 0.6 
Medium-cover grassland 0.5 0.6 0.5 0.5 0.6 0.7 0.5 
Low-cover grassland 0.3 0.5 0.3 0.5 0.8 0.5 0.4 
Body of water 0.8 0.5 0.4 0.8 0.6 0.2 
Town 
Rural settlements 
Industrial site 
Unutilized land 0.6 0.5 0.6 0.5 0.8 0.6 0.4 

Runoff abruptness and intra-annual variability

The M-K mutation test was adopted to analyze the annual runoff of several hydrological stations in the DTL basin, and the years in which mutations occurred were initially obtained, and the results are shown in Figure 3. By analyzing the information in the figure can be seen that the UF and UB values of the annual runoff of each hydrological station intersected many times in the 1970 and 1990s, occurring three times in CLJ station, six times in XT, TY, and TJ stations, and four times in SM station, and the intersection points were basically all within 95% confidence intervals, which reflected that the water situation was complicated and changeable in the DTL basin under the natural conditions before the 1990s, and also revealed that the water conditions in the natural conditions of the DTL basin were not as good as those in the past. After 1990, the UF and UB values of all stations except CLJ station continued to produce intersection, XT station in 2004, TY station in 2003 and 2013, TJ station in 2003, 2005, and 2015, and Shimen station in 2003, which all passed the 95% significance test. The 95% significance test indicates significant changes in annual runoff in each year. Except for the XT station, the M-K statistics of the flow at the other six stations are all negative, and it can be concluded that the flow at the XT station shows an increasing trend on the multi-year scale, while the flow at the other stations shows a decreasing trend.
Figure 3

Schematic diagram of the M-K test mutation.

Figure 3

Schematic diagram of the M-K test mutation.

Close modal
Utilizing the Mann–Whitney–Pettitt (MWP) test, a dual verification of abrupt changes in annual runoff was conducted. As depicted in Figure 4, which presents the analytical results of the MWP test for annual runoff at representative hydrological stations in the DTL basin, two distinct change-points were identified at the CLJ station, occurring in 1983 and 2003, with a |U| value less than 0.5. Consequently, the most probable years for abrupt changes in annual runoff at the CLJ station are 1983 and 2003. The mutation in 1983 may have been triggered by climate change and extreme weather events, while the mutation in 2003 was likely due to the operation of the Three Gorges Dam, which significantly impacted downstream runoff patterns. At the XT station, potential years for abrupt changes in runoff were identified as 1990 and 2004. Similarly, the MWP test revealed a mutation point in annual runoff at the TY station in 2003. At the TJ station, possible mutation years were determined to be 1987 and 2003. Finally, at the SM station, potential mutation years were identified as 1973, 1983, and 2003 (Table 6).
Table 6

Statistical table of mutation test results

Hydrological stationM-K non-parametric statisticsMWP testSynthesize and analyze
CLJ (DTL) 1971/1972/1975 1983/2003 2003 
XT (Xiang River) 1978/1985/1988/2004 1990/2004 2004 
TY (Yuan River) 1977/1988/2003/2013 2003 2003 
TJ (Zi River) 1966/1976/2003/2015 1987/2003 2003 
SM (Li River) 1965/1988/2003/2014 1973/1983/2003 2003 
Hydrological stationM-K non-parametric statisticsMWP testSynthesize and analyze
CLJ (DTL) 1971/1972/1975 1983/2003 2003 
XT (Xiang River) 1978/1985/1988/2004 1990/2004 2004 
TY (Yuan River) 1977/1988/2003/2013 2003 2003 
TJ (Zi River) 1966/1976/2003/2015 1987/2003 2003 
SM (Li River) 1965/1988/2003/2014 1973/1983/2003 2003 
Figure 4

Graph of MWP test results.

Figure 4

Graph of MWP test results.

Close modal
In Figure 5, the changes in annual runoff in the DTL basin before and after the sudden change are shown. As the last controlling hydrological station in the DTL basin, as shown in Figure 5(a), in the period after the mutation, the annual changes except January–March show an increase, and in the rest of the months after the mutation, there is a sharp decrease in the flow, of which the maximum decrease in flow occurs in July. In Figure 5(b)–5(e), annual flow changes in XT, TJ, TY, and SM are shown, and it can be seen that the flow changes in the Four Rivers are relatively small, but the annual trend is that most of the months in the period after the mutation have a decreasing trend, and the changes in the Four Rivers in XT and TJ are the most significant decreases, and this is mostly because not only is the runoff of the Xiang River and Yuan River much higher than that of the Zi River and Li River but also the flows of the Xiang and Yuan rivers are much lower than those of the Zi River and Li River. This is mostly because the runoff of the Xiang and Yuan rivers is not only much higher than that of the Zi and Li rivers, but also the Xiang and Yuan Rivers are close to the flood season of the Yangtze River mainstream, which is the result of the harmonization of water resources in the basin and the Yangtze River.
Figure 5

Schematic diagram of runoff changes before and after the abrupt change at each station.

Figure 5

Schematic diagram of runoff changes before and after the abrupt change at each station.

Close modal

ECF changes

First of all, according to the different rivers before and after the mutation, different hydrological stations of the environmental flow distribution are shown in Figure 6, for the distribution of flow events. As shown in Figure 6(a), in 2003 after the mutation of the runoff of the CLJ Hydrological Station, in the frequency of occurrence of small floods and large floods there is a substantial reduction in time and duration, but the peak value rises. By comparing the extremely low flow and low flow, we found that there were almost no extremely low flow events after 2003. It can be found that the construction of hydraulic engineering has a great influence on the annual flow level of CLJ Station, which has a strong role of ‘thinning the peaks and filling the dry’, and the variation range of extreme flow becomes significantly narrower.
Figure 6

Compositional distribution of environmental flow changes in the main stream of the Yangtze River and within the DTL basin: (a) CLJ station and (b–e) XT, TY, TJ, and SM stations.

Figure 6

Compositional distribution of environmental flow changes in the main stream of the Yangtze River and within the DTL basin: (a) CLJ station and (b–e) XT, TY, TJ, and SM stations.

Close modal

In Figure 6(b)–6(e), it can be found that the high-energy pulses in the four water basins cover most of the time, but in the post-surge period of XT and TJ, it can be found that the frequency of small- and medium-sized floods in the Xiang and Zi rivers is increased, which reflects the influence of the extreme climatic conditions in the DTL basin. The extremely low dry water flow of both TJ and SM reflects the seriousness of the water ecological safety hazard in the basin. In the distribution of the environmental flows of the Four Rivers, it can be seen that the peak values of the environmental flows of the Four Rivers were generally reduced in 2003 and 2004, especially for the small floods, reflecting the reduced probability of floods in the DTL basin as a humid and rainy area.

In general, compared with the period before the mutation, the frequency and magnitude of extra-dry flow events and large and small floods were reduced after the mutation, especially the extra-dry flow events, which seldom occurred after the mutation, and the peak flow of small flood events decreased, especially in the DTL basin where the high-energy pulses took up most of the time intervals, and the extra-dry flow and large floods almost did not occur. The construction of water conservancy projects has played a great role in protecting water security in the basin.

Characteristics of meteorological changes

We calculated the year-to-year changes of eight extreme precipitation indices and four extreme temperature indices by the RClimDex model for the meteorological indices of 38 weather stations in the DTL basin, respectively, and then analyzed the distribution and differences of the climatic states in the DTL basin through inverse-distance weighting interpolation using the mean values of the meteorological indices. The distribution of climate states is shown in Figures 7 and 8. TN10p represents the number of days with minimum temperature <10% quartile and TN90p represents the number of days with minimum temperature >90% quartile. It can be found that TN10p and TN90p are the strongest near the central and western regions, and the two performances are synchronized trends, while in the northwest, the indices of temperature TN10p and TN90p in the northeast direction are always in the lower-range interval.
Figure 7

Map of temperature index distribution in the DTL basin.

Figure 7

Map of temperature index distribution in the DTL basin.

Close modal
Figure 8

Distribution of extreme rainfall indices in the DTL basin.

Figure 8

Distribution of extreme rainfall indices in the DTL basin.

Close modal

Figure 7 shows the regional spatial distribution of temperature indices in the DTL basin, and it can be found that the minimum temperature indices are mostly distributed in the southeast and west, where the mountain ranges are isolated to form an obvious low-temperature holding zone, specifically the east–west trending Nanling Mountain Range that stretches across the border between Hunan and Guangdong, which blocks the cold air from the north from coming down to the south to form a cold-air layer in the sky above Hunan, and the unique mountain range distribution has an impact on the distribution of the minimum temperature. By comparing the distribution of the maximum temperature index, it can be found that in the northeast and southwest of the zone, the climate of high-temperature weather stays higher, under the influence of the mountain range, resulting in the southern zone of the tropical monsoon being blocked. In the northeast of DTL, the TX10p and TX90p indices are significantly higher in the high-temperature areas around the wetlands of DTL than in the other areas due to the plains and fewer hills and mountains.

Figure 8 shows the regional spatial distribution of extreme rainfall indices, and it can be noticed that the high rainfall indices are mostly distributed in the west as well as in the northern belt and in the south, while there is a significant decrease in rainfall closer to the east. In Figure 8(g) of the annual total rainfall index PRCPTOT, it can be found that the high rainfall areas in the basin are less distributed in the southwestern and northern zones, which are in the condition of blocking by the mountain range, and the rainfall events in the whole year are significantly reduced compared with the eastern region. In the Rx1day and Rx5day distributions, the annual maximum daily precipitation and the annual maximum five consecutive days of precipitation in the northern region are significantly smaller than those in the other regions, and the intensity of rainfall over a short period of time represents a magnitude of heavy rainfall that will be much lower than that in the other regions. In the R10mm distribution map (Figure 8(d)), it can only be seen that the rainfall zone in the basin occurs in the eastern part of the basin, and the number of rainfall days in a year is significantly higher than that in the other basins of the Four Rivers due to the coverage of the Xiang River basin in this area. Specifically in the distribution of the Four Rivers watersheds of DTL, it can be found that the rainfall in the Xiang River area is significantly higher than that in the other watersheds. The extreme daily rainfall intensity in the Yuan and Li rivers is significantly higher than that in the other areas, and they are in the same rainstorm area as the Yangtze River, so the impacts on the main stream of the Yangtze River will be much greater than those on the Zi and Xiang Rivers in a combined rainstorm event.

Characterization of the correlation between flow status and meteorological indices

Climatic conditions, as an important influence on surface water resources, change the hydrological process and hydrological-cycle state of surface water, so the analysis of the potential correlation between climate change and flow is of great reference significance for the ecological protection of water resources. Figure 9 shows the distribution of correlation between meteorological indices of each sub-basin in the DTL basin and environmental flow-related indices of each hydrological station. It can be found in Figure 9(a) for CLJ that among all meteorological indices, SDII, PRCPTRT, and R10MM have the highest degree of change on environmental flow, and reach the maximum of 0.886 in R10MM-Aug lowf, and the correlation of TN10P with environmental flow is the highest for all air temperature indices. TN10P has the least correlation to ambient flow and is mostly below 0.8. In terms of environmental flows, the sfldl series are all significantly less correlated with meteorological indicators than the other environmental flow indices, with most of their gray correlation indices below 0.7, and the lowest correlation occurs in the Sfldl Time index.
Figure 9

Distribution of ambient flow–meteorological index correlations by the hydrological station.

Figure 9

Distribution of ambient flow–meteorological index correlations by the hydrological station.

Close modal

Figure 9(b)–9(e) shows the distribution of the correlation of the four sub-basins to the meteorological indices in the basin, and it can be found that most of the correlations between the environmental flow indices and the meteorological indices in the Four Rivers are very similar, and the correlation of the meteorological indices in the Xlowl Dur and the Sfldl Peak are all blank, which may be due to the frequency of large floods and small floods as well as the occurrence of extreme dry water events in the four sub-basins. This may be due to the frequency of major and minor floods as well as extreme dry water events, on the one hand, as the Four Rivers basin is a humid area with sufficient rainfall, and on the other hand, to the difference between the high flow rate in the Four Rivers as Class II and Class IV rivers compared with the Yangtze River. The response to the correlation to the meteorological indices is obvious, even in the river since the sudden change of the river to no more flood events, and the correlation with the meteorological indicators naturally cannot be analyzed. In the subsequent comparison of the four sub-basins, it was found that for the XT and TY hydrographic stations the correlation of the meteorological index of their respective basins is significantly higher than for the TY and SM hydrographic stations analyzed.

LSTM model testing and attribution analysis

To quantify the impacts of climate and human activities on river conditions within a watershed, this study utilized meteorological data as a driver to reconstruct the flow during the variation period, which was considered a natural state solely under the influence of climate change (Figure 10). Under the settings of 200 hidden neurons, an initial learning rate of 0.0016, a learning rate decay period of 40, a learning rate decay factor of 0.12, and a batch size of 4, the LSTM model demonstrated promising performance in both the calibration and validation sets, as presented in Table 7. The NSE values for the reconstructed flows versus the measured flows at each hydrological station were above 0.75, indicating good agreement. Furthermore, the differences in NSE between the baseline period and the calibration period were less than 5%, demonstrating the model's ability to effectively simulate the flow variations at each station under natural conditions.
Table 7

LSTM model validation and natural flow reconstruction results

StationsCalibration period (1961–1990)
Validation period (1990–mutation)
Variation period (mutation–2019)
NSEMAPERMSENSEMAPERMSENSEMAPERMSE
CLJ 0.896 655.3 0.901 0.887 525.6 0.863 0.855 434.6 0.822 
XT 0.908 457.1 0.908 0.889 473.9 0.876 0.841 565.5 0.837 
TY 0.842 565.5 0.837 0.894 489.4 0.903 0.891 517.6 0.887 
TJ 0.755 756.7 0.750 0.808 209.9 0.801 0.797 231.9 0.790 
SM 0.820 189.1 0.815 0.807 212.5 0.802 0.715 169.5 0.709 
StationsCalibration period (1961–1990)
Validation period (1990–mutation)
Variation period (mutation–2019)
NSEMAPERMSENSEMAPERMSENSEMAPERMSE
CLJ 0.896 655.3 0.901 0.887 525.6 0.863 0.855 434.6 0.822 
XT 0.908 457.1 0.908 0.889 473.9 0.876 0.841 565.5 0.837 
TY 0.842 565.5 0.837 0.894 489.4 0.903 0.891 517.6 0.887 
TJ 0.755 756.7 0.750 0.808 209.9 0.801 0.797 231.9 0.790 
SM 0.820 189.1 0.815 0.807 212.5 0.802 0.715 169.5 0.709 
Figure 10

Performance of the LSTM model and reconstruction results of natural flows during the variation period.

Figure 10

Performance of the LSTM model and reconstruction results of natural flows during the variation period.

Close modal
According to the results of the attribution analysis, anthropogenic changes caused a decrease in runoff in the DTL basin in most months of the year (Figure 11). On the annual scale, human activities were the dominant factor for runoff changes, especially at CLJ and TY hydrological stations, where the proportion reached 86% and 89%. On the monthly scale, except for the XT hydrological station on the Xiangjiang River (Figure 11(b)), all of them showed that human activities significantly influenced the flow changes throughout the year. It is especially important to note that at the XT station, the climate contribution at the monthly scale is much larger than the effect of human activities, and climate change is reducing the flow in most of the months. According to Figure 11(b)–11(e), when climate is the main influence on runoff, human activities are reducing the magnitude of flow and thus the flood-control pressure in DTL after the confluence. In Figure 11(b)–11(e), it can be seen that the changes in the XT hydrological station during the year are mainly influenced by climate, and the changes in runoff volume due to human activities are all in the opposite state to the changes in climate. The LSTM model residuals of each hydrological station have less influence on the annual scale, but the monthly-scale residuals of the CJL and TJ hydrological stations are in the larger value, and most of the residuals are between −350 and 175, which have less influence on the overall results.
Figure 11

Impacts of climate change and human activities on runoff from hydrological stations in the DTL basin.

Figure 11

Impacts of climate change and human activities on runoff from hydrological stations in the DTL basin.

Close modal

Changes in habitat quality suitability

The InVEST model uses the habitat quality index to characterize the status of habitat quality. In the model, the habitat quality index in the raster layer changes continuously as a 0–1 value; the closer the value is to 1, the better the indicated quality of the habitat and the habitat is relatively more intact, but also has a certain structure and function, which is conducive to the maintenance of biological diversity. Figure 12 shows the habitat quality of the DTL watershed in 1980, 1990, 2000, 2005, 2010, and 2020 as calculated by the model, which categorizes the habitat quality of the DTL basin into the following four classes: low (0–0.4), medium (0.4–0.6), relatively high (0.6–0.8), and high (0.8–1), and the proportion is calculated of each grade in the basin area (Table 8).
Table 8

Changes in the percentage of habitat quality in the DTL watershed, 1980–2020 (%)

Habitat quality ratingPoint interval (in calculus)1980199020002005201020201980–2020
Low 0–0.2 1.35 1.49 1.6 1.74 1.87 2.86 1.51 
Moderate 0.4–0.6 29.53 28.79 28.77 28.64 28.57 28 −1.53 
Relatively high 0.6–0.8 7.56 6.55 7.75 7.65 7.74 7.66 0.1 
High 0.8–1 61.56 63.17 61.88 61.97 61.82 61.48 −0.08 
Habitat quality ratingPoint interval (in calculus)1980199020002005201020201980–2020
Low 0–0.2 1.35 1.49 1.6 1.74 1.87 2.86 1.51 
Moderate 0.4–0.6 29.53 28.79 28.77 28.64 28.57 28 −1.53 
Relatively high 0.6–0.8 7.56 6.55 7.75 7.65 7.74 7.66 0.1 
High 0.8–1 61.56 63.17 61.88 61.97 61.82 61.48 −0.08 
Figure 12

Distribution of habitat quality within the DTL basin.

Figure 12

Distribution of habitat quality within the DTL basin.

Close modal

Based on the spatial pattern of habitat quality in the DTL basin (Figure 12), it was observed that habitat quality varied considerably spatially, with a wider distribution of higher habitats (0.8–1) in general. Calculations showed that the mean values of different habitat quality in 1980, 1990, 2000, 2005, 2010, and 2020 were 0.8261, 0.82, 0.8188, 0.8183, 0.8172, and 0.8106, respectively, and that the habitat quality in the DTL basin showed a gradually decreasing trend of change, with an overall obvious decline. In terms of the different habitat-quality grades in the DTL basin, the ‘low’ grade showed an increasing trend from 1980 to 2020, with an increase of 1.51%. The ‘medium’ grade shows a decreasing trend, with a decrease of 1.53%, indicating that the rapid economic development under the influence of human activities has led to a clear downward trend in habitat quality. The changes in the ‘relatively high’ and ‘high’ grades were smaller, with the area increasing by 0.1% and decreasing by 0.08%, respectively. In the DTL basin, the ‘high’ grade of habitat quality accounts for more than 60% of the total area, which can indicate that the overall performance of habitat quality in the DTL basin is better, but from 1980 to 2020, the quality of habitat gradually declined, and there is a risk of habitat deterioration in the DTL basin.

Impact of changes in hydrological regimes

In this study, by examining the changes in flow processes recorded at various hydrological stations in the DTL basin, the results show that human activities play a dominant role most of the time. However, the influence played by climatic aspects cannot be ignored; the distribution of the correlation between the RClimDex index and the environmental flow found that the influence of rainfall was significantly higher than that of the temperature index, and rainfall was affected by the geomorphology, and part of the rainfall converged to DTL through the river to form a larger catchment area, and this part of the rainfall played a more pronounced role in evapotranspiration compared with the temperature (Oliveira & Quaresma 2017). Within watersheds with connected lakes, watershed hydrological changes significantly influence the ecological succession of the lakes. Also, within the Poyang Lake watershed in the lower section of the Yangtze River, Li et al. (2015) demonstrated that changes in hydrological regimes within the watershed significantly influence lake connectivity and ecological conditions, and drive the dynamic migration of vegetation and water quality changes in the lake area. Zhang et al. (2021) and others demonstrated that hydrological changes within a watershed lead to changes in ecological status on a larger scale through confluence, and Khatun et al. (2021) accurately identified wetland ecological responses to changes in hydrological regimes while exploring scenarios in which water quality changes and biological responses are leading to declines in the abundance of water bodies in connected lake wetlands. Changes in the hydrological regimes of rivers connected to lakes will significantly alter the ecological status of lakes, and these ecological changes will be more severe if they occur in complex lake wetlands or floodplain wetlands (Lai et al. 2014).

Risk of habitat degradation in watersheds

In recent years, human activities such as urban expansion, lake reclamation for farmland, and industrial pollution have significantly altered the original state of ecological structures, leading to degradation in habitat quality and disruption of ecological balance (Calderon et al. 2019). As evident from the habitat quality results (Figure 11), high-grade habitat quality is primarily distributed in areas far from the DTL water body, where there are numerous mountainous and hilly regions with relatively high forest coverage. These regions benefit from various ecological protection policies for flora and fauna (Yuan et al. 2021). Conversely, human activities centered around the DTL perimeter have led to a decline in ecological conditions in nearby areas, resulting in an inverse relationship between habitat quality and population density (Hu et al. 2015). These spatial differences in habitat degradation suggest that habitats closer to urban and industrial areas are more likely to experience degradation. The habitat quality along the riverbanks within the basin is moderate (0.4–0.6), and there has been a significant increase in low-grade habitat quality near the riverbanks after 2005. This trend is consistent with the high population density along the banks of the DTL basin and the decrease in habitat quality caused by human activities, aligning with research findings in other regions (Mengist et al. 2021).

In previous studies, low-grade habitat quality has similarly been concentrated in industrial areas and densely populated regions, where these areas have been continual subjected to anthropogenic disturbances such as changes in land use, industrial pollution, and habitat fragmentation (Chen et al. 2023). These influencing factors exhibit high inter-correlations, jointly influencing habitat quality changes while also interacting with each other. However, incorporating all these factors into an analysis to identify the drivers of habitat degradation is challenging. Current research indicates that persistent and severe land use and land cover (LULC) changes are the primary factors affecting habitat quality grades (Tang et al. 2021). Figures 13 and 14 depict the distribution of habitat degradation factors and the extent of habitat threats, respectively. Their spatial and temporal distributions closely align with the distribution of low-grade habitat quality. This alignment is primarily due to the rapid development of industrial and urban agglomerations in these areas, which has a significant impact on habitats, leading to a series of crises for the habitat quality of the region (Haddeland et al. 2014). To understand the current status of the ecological environment and incorporate this into nature conservation policies, it is necessary to assess habitat quality. Evaluating the impact of human activities within the basin (such as urbanization, agriculture, and industrial pollution) and the complexity of their interactions is crucial for nature conservation planning (Feng et al. 2021). However, previous research efforts have often overlooked the necessity of comprehensive land-surface resource planning, particularly in regions with limited accessible water resources. Consequently, mapping the spatiotemporal distribution of habitat quality and degradation represents a significant contribution. This mapping can aid in devising practical conservation strategies and provide supportive information for policymakers.
Figure 13

Distribution of habitat degradation in the basin.

Figure 13

Distribution of habitat degradation in the basin.

Close modal
Figure 14

Distribution of habitat threat factors in the watershed.

Figure 14

Distribution of habitat threat factors in the watershed.

Close modal

Implications for water resources management

Within large catchment areas, climate change is a major factor affecting habitats, particularly in terms of rainfall (Astagneau et al. 2022). Precipitation is a direct manifestation of climate change, directly influencing changes in surface water and ecological water demand in the region, with river catchments playing an equally important role. With the rapid development of industry and population increase, the utilization of water resources in the DTL basin is increasing, the contradiction between supply and demand of water resources is becoming more and more serious, and the trade-off relationship between water resources and environmental protection needs to be further studied (Balata et al. 2022). As of 2019, 14,096 reservoirs have been built in the DTL basin, with a total capacity area of 51.412 billion m3, which significantly regulates the hydrological patterns of rivers in the DTL basin and causes serious ecological problems (Zheng et al. 2023). For this reason, the coordination of ecological quality and social development in the region has widely attracted the attention of researchers and managers worldwide (Mao et al. 2019). Due to the dependence of numerous ecological patterns in the watershed on water resources, the ecosystems connected through the vast water system will be significantly affected by the water resources pattern and fluctuate continuously, thus putting more demands on the regulation of water resources.

Limitations and applicability of the study

Uncertainty is a widespread phenomenon in the process of hydrological modeling. The research on climate-driven hydrological changes presented in this article is based on the LSTM deep-learning model. Both model parameters and the quality of input data are significant sources of uncertainty in simulation results (Cui & Yu 2023). To ensure the quality of input data, all data used in the model originate from measured observations, and outlier removal as well as mean normalization have been applied. To minimize the impact of model uncertainty, we employed a manual trial-and-error method to train the optimal hyperparameters. This approach aimed to achieve both minimal model fluctuations and optimal accuracy. Additionally, we utilized a mean-based approach to process the output results. Under these constraints, the flow reconstruction results for all hydrological stations during the baseline period exhibited excellent agreement with the measured values (Table 7). The NSE for both the calibration and validation datasets was higher than 0.75, and the differences in fitting metrics were below 5%. These results indicate that the model's performance during the variation period can be used to evaluate the impact of climate change on hydrology (Frame et al. 2023). Overall, despite the inherent uncertainty in hydrological modeling, the approach taken in this study, including data preprocessing, hyperparameter optimization, and output processing, has effectively minimized the uncertainty and produced reliable results for assessing climate-driven hydrological changes.

This study was conducted to gain insight into the linkages between climate conditions and the hydrology of tributaries in the DTL basin and to analyze the evolution of habitat quality in the basin. We analyzed the evolutionary characteristics of climate indicators and environmental flow indicators of tributaries in the DTL basin, introduced the LSTM model to quantify the contribution of climate change and human activities to runoff changes on short and long time-scales, and introduced the InVEST model to assess the habitat-change status in the basin during 1980–2020. The main conclusions are summarized below:

  • (1) The tributaries in the DTL basin showed a smaller range of environmental flow characteristics after the abrupt change, and the frequency of extreme flow and hydrograph time decreased.

  • (2) Within each tributary, the rainfall index showed a high correlation with the environmental flow index, and the temperature index showed a low correlation except for DTL. The correlation of climate indices was strongest in the Xiang River basin and weakest in the Yuan River basin.

  • (3) The attribution analysis showed that the impacts of human activities on the annual scale were generally above 80%, with the highest being 89% in the Yuan River. On the intra-annual scale, human impacts were mostly dominant, but on the Xiang River, the impacts of climatic factors were more dominant.

  • (4) Higher habitat-quality distribution accounted for the highest percentage, in all years accounting for more than 60% and mainly distributed away from the water body of DTL. Habitat quality is lower in highly populated areas and around DTL.

In the future, further consideration can be given to policy orientation, species invasion and industrial pollution damage effects, and the correlation between them can be analyzed for a more comprehensive assessment.

This study was supported by the Science and Technology Program of Guizhou Province (KT202008) and the Basic Research Project of Key Scientific Research Projects of Colleges and Universities of Henan Province (23ZX012).

X.B. conceived the study and wrote the first draft, the data collection and analysis was done by X.B. and W.Y., the methodology was analyzed with the help of Y.M. and M.Z., literature by 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.

The scope of the work does not require any approval of the ethics committee. The ethical approval is not applicable for the present work.

The scope of the work does not require participation of any candidate. The consent to participate does not apply for the present work.

Consent of all the authors was taken before the communication of the manuscript in its present form.

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

Astagneau
P. C.
,
Bourgin
F.
,
Andréassian
V.
&
Perrin
C.
2022
Catchment response to intense rainfall: evaluating modelling hypotheses
.
Hydrological Processes
36
(
8
),
e14676
.
https://doi.org/10.1002/hyp.14676
.
Balata
E. E.
,
Pinto
H.
&
Moreira da Silva
M.
2022
Latent dimensions between water use and socio-economic development: a global exploratory statistical analysis
.
Regional Sustainability
3
(
3
),
269
280
.
https://doi.org/10.1016/j.regsus.2022.09.004
.
Bastos
M. I.
,
Roebeling
P. C.
,
Alves
F. L.
,
Villasante
S.
&
Magalhães Filho
L.
2023
High risk water pollution hazards affecting Aveiro coastal lagoon (Portugal) – a habitat risk assessment using InVEST
.
Ecological Informatics
76
,
102144
.
https://doi.org/10.1016/j.ecoinf.2023.102144
.
Calderon
M. R.
,
Almeida
C. A.
,
González
P.
&
Jofré
M. B.
2019
Influence of water quality and habitat conditions on amphibian community metrics in rivers affected by urban activity
.
Urban Ecosystems
22
(
4
),
743
755
.
https://doi.org/10.1007/s11252-019-00862-w
.
Cao
S.
,
Hu
X.
,
Wang
Y.
,
Chen
C.
,
Xu
D.
&
Bai
T.
2023
Understanding spatial-temporal interactions of ecosystem services and their drivers in a multi-scale perspective of Miluo using multi-source remote sensing data
.
Remote Sensing
15
,
3479. https://doi.org/10.3390/rs15143479
.
Cheng
J.
,
Xu
L.
,
Wang
X.
,
Jiang
J.
&
You
H.
2018
Assessment of hydrologic alteration induced by the Three Gorges Dam in Dongting Lake, China
.
River Research and Applications
34
(
7
),
686
696
.
https://doi.org/10.1002/rra.3297
.
Cid
N.
,
Erős
T.
,
Heino
J.
,
Singer
G.
,
Jähnig
S. C.
,
Cañedo-Argüelles
M.
,
Bonada
N.
,
Sarremejane
R.
,
Mykrä
H.
,
Sandin
L.
,
Paloniemi
R.
,
Varumo
L.
&
Datry
T.
2022
From meta-system theory to the sustainable management of rivers in the Anthropocene
.
Frontiers in Ecology and the Environment
20
(
1
),
49
57
.
https://doi.org/10.1002/fee.2417
.
Cui
W.
&
Yu
W.
2023
Uncertainty injection: a deep learning method for robust optimization
.
IEEE Transactions on Wireless Communications
22
(
11
),
7201
7213
.
https://doi.org/10.1109/TWC.2023.3249103
.
Dong
W.
,
Liu
S.
&
Fang
Z.
2018
On modeling mechanisms and applicable ranges of grey incidence analysis models
.
Grey Systems: Theory and Application
8
(
4
),
448
461
.
https://doi.org/10.1108/GS-04-2018-0019
.
dos Santos
C. A. C.
,
Neale
C. M. U.
,
Rao
T. V. R.
&
da Silva
B. B.
2011
Trends in indices for extremes in daily temperature and precipitation over Utah, USA
.
International Journal of Climatology
31
(
12
),
1813
1822
.
https://doi.org/10.1002/joc.2205
.
Feng
H.
,
Li
Y.
,
Li
Y.
,
Li
N.
,
Li
Y.
,
Hu
Y.
,
Yu
J.
&
Luo
H.
2021
Identifying and evaluating the ecological network of Siberian roe deer (Capreolus pygargus) in Tieli Forestry Bureau, northeast China
.
Global Ecology and Conservation
26
,
e01477
.
https://doi.org/10.1016/j.gecco.2021.e01477
.
Frame
J. M.
,
Kratzert
F.
,
Gupta
H. V.
,
Ullrich
P.
&
Nearing
G. S.
2023
On strictly enforced mass conservation constraints for modelling the rainfall-runoff process
.
Hydrological Processes
37
(
3
),
e14847
.
https://doi.org/10.1002/hyp.14847
.
Geng
M.
,
Wang
K.
,
Yang
N.
,
Li
F.
,
Zou
Y.
,
Chen
X.
,
Deng
Z.
&
Xie
Y.
2021a
Spatiotemporal water quality variations and their relationship with hydrological conditions in Dongting Lake after the operation of the Three Gorges Dam, China
.
Journal of Cleaner Production
283
,
124644
.
https://doi.org/10.1016/j.jclepro.2020.124644
.
Geng
M.
,
Wang
K.
,
Yang
N.
,
Qian
Z.
,
Li
F.
,
Zou
Y.
,
Chen
X.
,
Deng
Z.
&
Xie
Y.
2021b
Is water quality better in wet years or dry years in river-connected lakes? A case study from Dongting Lake, China
.
Environmental Pollution
290
,
118115
.
https://doi.org/10.1016/j.envpol.2021.118115
.
Guan
D.
,
Wu
L.
,
Cheng
L.
,
Zhang
Y.
&
Zhou
L.
2022
How to measure the ecological compensation threshold in the upper Yangtze River basin, China? An approach for coupling InVEST and grey water footprint
.
Frontiers in Earth Science
10
,
988291
.
https://doi.org/10.3389/feart.2022.988291
.
Haddeland
I.
,
Heinke
J.
,
Biemans
H.
,
Eisner
S.
,
Flörke
M.
,
Hanasaki
N.
,
Konzmann
M.
,
Ludwig
F.
,
Masaki
Y.
,
Schewe
J.
,
Stacke
T.
,
Tessler
Z. D.
,
Wada
Y.
&
Wisser
D.
2014
Global water resources affected by human interventions and climate change
.
Proceedings of the National Academy of Sciences
111
(
9
),
3251
3256
.
https://doi.org/10.1073/pnas.1222475110
.
Hu
C.
,
Fang
C.
&
Cao
W.
2015
Shrinking of Dongting Lake and its weakening connection with the Yangtze River: analysis of the impact on flooding
.
International Journal of Sediment Research
30
(
3
),
256
262
.
https://doi.org/10.1016/j.ijsrc.2014.05.001
.
Keith
D. A.
,
Ferrer-Paris
J. R.
,
Nicholson
E.
,
Bishop
M. J.
,
Polidoro
B. A.
,
Ramirez-Llodra
E.
,
Tozer
M. G.
,
Nel
J. L.
,
Mac Nally
R.
,
Gregr
E. J.
,
Watermeyer
K. E.
,
Essl
F.
,
Faber-Langendoen
D.
,
Franklin
J.
,
Lehmann
C. E. R.
,
Etter
A.
,
Roux
D. J.
,
Stark
J. S.
,
Rowland
J. A.
, Brummitt, N. A., Fernandez-Arcaya, U. C., Suthers, I. M., Wiser, S. K., Donohue, I., Jackson, L. J., Pennington, R. T., Iliffe, T. M., Gerovasileiou, V., Giller, P., Robson, B. J., Pettorelli, N., Andrade, A., Lindgaard, A., Tahvanainen, T., Terauds, A., Chadwick, M. A., Murray, N. J., Moat, J., Pliscoff, P., Zager, I. &
Kingsford
R. T.
2022
A function-based typology for Earth's ecosystems
.
Nature
610
(
7932
),
513–518
.
https://doi.org/10.1038/s41586-022-05318-4
.
Khatun
R.
,
Talukdar
S.
,
Pal
S.
&
Kundu
S.
2021
Measuring dam induced alteration in water richness and eco-hydrological deficit in flood plain wetland
.
Journal of Environmental Management
285
,
112157
.
https://doi.org/10.1016/j.jenvman.2021.112157
.
Klotz
D.
,
Kratzert
F.
,
Gauch
M.
,
Keefe Sampson
A.
,
Brandstetter
J.
,
Klambauer
G.
,
Hochreiter
S.
&
Nearing
G.
2022
Uncertainty estimation with deep learning for rainfall–runoff modeling
.
Hydrology and Earth System Sciences
26
(
6
),
1673
1693
.
https://doi.org/10.5194/hess-26-1673-2022
.
Lai
X.
,
Liang
Q.
,
Jiang
J.
&
Huang
Q.
2014
Impoundment effects of the Three-Gorges-Dam on flow regimes in two China's largest freshwater lakes
.
Water Resources Management
28
(
14
),
5111
5124
.
https://doi.org/10.1007/s11269-014-0797-6
.
Lei
T.
&
Middleton
B.
2021
A US-China EcoPartnership study of disturbed wetland vegetation in West Dongting Lake, China
.
Environmental Progress & Sustainable Energy
40
(
5
),
e13673
.
https://doi.org/10.1002/ep.13673
.
Lei
J.
,
Chen
Y.
,
Li
L.
,
Chen
Z.
,
Chen
X.
,
Wu
T.
&
Li
Y.
2022
Spatiotemporal change of habitat quality in Hainan Island of China based on changes in land use
.
Ecological Indicators
145
,
109707
.
https://doi.org/10.1016/j.ecolind.2022.109707
.
Li
X.
,
Zhang
Q.
,
Xu
C.-Y.
&
Ye
X.
2015
The changing patterns of floods in Poyang Lake, China: characteristics and explanations
.
Natural Hazards
76
(
1
),
651
666
.
https://doi.org/10.1007/s11069-014-1509-5
.
Luo
Y.
,
Liu
C.
,
Wang
Y.
,
Yang
Y.
&
Mishra
S.
2023
Occurrence, distribution and their correlation with different parameters of antibiotics and antibiotic resistance genes in lakes of China: a review
.
Marine Pollution Bulletin
193
,
115189
.
https://doi.org/10.1016/j.marpolbul.2023.115189
.
Mao
X.
,
Wei
X.
,
Jin
X.
,
Tao
Y.
,
Zhang
Z.
&
Wang
W.
2019
Monitoring urban wetlands restoration in Qinghai Plateau: integrated performance from ecological characters, ecological processes to ecosystem services
.
Ecological Indicators
101
,
623
631
.
https://doi.org/10.1016/j.ecolind.2019.01.066
.
Mao
R.
,
Shi
A.
,
Song
J.
,
Xu
W.
,
Tang
B.
&
Li
B.
2023
Response of the runoff process to meteorological drought: baseflow index as an important indicator
.
Journal of Environmental Management
345
,
118843
.
https://doi.org/10.1016/j.jenvman.2023.118843
.
Mengist
W.
,
Soromessa
T.
&
Feyisa
G. L.
2021
Landscape change effects on habitat quality in a forest biosphere reserve: implications for the conservation of native habitats
.
Journal of Cleaner Production
329
,
129778
.
https://doi.org/10.1016/j.jclepro.2021.129778
.
Oliveira
K. S. S.
&
Quaresma
V. d. S.
2017
Temporal variability in the suspended sediment load and streamflow of the Doce River
.
Journal of South American Earth Sciences
78
,
101
115
.
https://doi.org/10.1016/j.jsames.2017.06.009
.
Peres
D. J.
&
Cancelliere
A.
2016
Environmental flow assessment based on different metrics of hydrological alteration
.
Water Resources Management
30
(
15
),
5799
5817
.
https://doi.org/10.1007/s11269-016-1394-7
.
Qi
S.
,
Chen
S.
,
Long
X.
,
An
X.
&
Zhang
M.
2023
Quantitative contribution of climate change and anthropological activities to vegetation carbon storage in the Dongting Lake basin in the last two decades
.
Advances in Space Research
71
(
1
),
845
868
.
https://doi.org/10.1016/j.asr.2022.07.068
.
Sabater
S.
,
Bregoli
F.
,
Acuña
V.
,
Barceló
D.
,
Elosegi
A.
,
Ginebreda
A.
,
Marcé
R.
,
Muñoz
I.
,
Sabater-Liesa
L.
&
Ferreira
V.
2018
Effects of human-driven water stress on river ecosystems: a meta-analysis
.
Scientific Reports
8
(
1
),
11462
.
https://doi.org/10.1038/s41598-018-29807-7
.
Sepulveda
A. J.
,
Gage
J. A.
,
Counihan
T. D.
&
Prisciandaro
A. F.
2023
Can big data inform invasive dreissenid mussel risk assessments of habitat suitability?
.
Hydrobiologia
https://doi.org/10.1007/s10750-023-05156-z
.
Singh
U.
,
Maca
P.
,
Hanel
M.
,
Markonis
Y.
,
Nidamanuri
R. R.
,
Nasreen
S.
,
Blöcher
J. R.
,
Strnad
F.
,
Vorel
J.
,
Riha
L.
&
Raghubanshi
A. S.
2023
Hybrid multi-model ensemble learning for reconstructing gridded runoff of Europe for 500 years
.
Information Fusion
97
,
101807
.
https://doi.org/10.1016/j.inffus.2023.101807
.
Tang
L.
,
Ke
X.
,
Chen
Y.
,
Wang
L.
,
Zhou
Q.
,
Zheng
W.
&
Xiao
B.
2021
Which impacts more seriously on natural habitat loss and degradation? Cropland expansion or urban expansion?
.
Land Degradation & Development
32
(
2
),
946
964
.
https://doi.org/10.1002/ldr.3768
.
Thompson
R. M.
,
Barbour
E. J.
,
Bradshaw
C. J. A.
,
Briggs
S.
,
Byron
N.
,
Grace
M.
,
Hart
B. T.
,
King
A. J.
,
Likens
G. E.
,
Pollino
C. A.
,
Sheldon
F.
,
Stewardson
M. J.
,
Thoms
M.
,
Watts
R. J.
&
Webb
J. A.
2022
Principles for scientists working at the river science–policy interface
.
River Research and Applications
38
(
5
),
819
831
.
https://doi.org/10.1002/rra.3951
.
Wang
Q.
,
Xu
Y.
,
Wang
Y.
,
Zhang
Y.
,
Xiang
J.
,
Xu
Y.
&
Wang
J.
2020
Individual and combined impacts of future land-use and climate conditions on extreme hydrological events in a representative basin of the Yangtze River Delta, China
.
Atmospheric Research
236
,
104805
.
https://doi.org/10.1016/j.atmosres.2019.104805
.
Wang
L.
,
Li
Z.
,
Wang
D.
,
Chen
J.
,
Liu
Y.
,
Nie
X.
,
Zhang
Y.
,
Ning
K.
&
Hu
X.
2021
Unbalanced social–ecological development within the Dongting Lake basin: inspiration from evaluation of ecological restoration projects
.
Journal of Cleaner Production
315
, 128161.
https://doi.org/10.1016/j.jclepro.2021.128161
.
Wang
Y.
,
Xia
J.
,
Cai
W.
,
Liu
Z.
,
Li
J.
,
Yin
J.
,
Zu
J.
&
Dou
C.
2023
Response of fish habitat quality to weir distribution change in mountainous river based on the two-dimensional habitat suitability model
.
Sustainability
15
(
11
),
8698
.
https://doi.org/10.3390/su15118698
.
Wu
H.
&
Zheng
B.
2020
Wetland area identification and waterbird protection management in consideration of lake topography and water level change
.
Global Ecology and Conservation
23
,
e01056
.
https://doi.org/10.1016/j.gecco.2020.e01056
.
Yang
J.
,
Xie
B.
,
Zhang
D.
,
Mak-Mensah
E.
&
Pei
T.
2023
Habitat quality assessment and multi-scenario prediction of the Gansu–Qinghai section of the Yellow River Basin based on the FLUS-InVEST model
.
Frontiers in Ecology and Evolution
11
,
1228558
.
https://doi.org/10.3389/fevo.2023.1228558
.
Yuan
B.
,
Fu
L.
,
Zou
Y.
,
Zhang
S.
,
Chen
X.
,
Li
F.
,
Deng
Z.
&
Xie
Y.
2021
Spatiotemporal change detection of ecological quality and the associated affecting factors in Dongting Lake Basin, based on RSEI
.
Journal of Cleaner Production
302
,
126995
.
https://doi.org/10.1016/j.jclepro.2021.126995
.
Yue
S.
,
Pilon
P.
,
Phinney
B.
&
Cavadias
G.
2002
The influence of autocorrelation on the ability to detect trend in hydrological series
.
Hydrological Processes
16
(
9
),
1807
1829
.
https://doi.org/10.1002/hyp.1095
.
Zhang
Q.
,
Zhang
Z.
,
Shi
P.
,
Singh
V. P.
&
Gu
X.
2018
Evaluation of ecological instream flow considering hydrological alterations in the Yellow River basin, China
.
Global and Planetary Change
160
,
61
74
.
https://doi.org/10.1016/j.gloplacha.2017.11.012
.
Zhang
C.
,
Li
M.
,
Chang
T.
,
Gao
X.
,
Liu
H.
&
Cao
W.
2021
The interaction processes of the fish assemblages between the Yangtze River and Poyang Lake, China
.
Ecology of Freshwater Fish
30
(
4
),
541
550
.
https://doi.org/10.1111/eff.12603
.
Zhang
J.
,
Du
E.
,
Yuan
S.
,
Wu
J.
,
Song
Y.
&
Wang
G.
2023
Initiatives on exploring the mechanism of eco-hydrological response to land surface change and adaptive regulation in the Yellow River Basin
.
River
2
(
2
),
149
157
.
https://doi.org/10.1002/rvr2.50
.
Zheng
F.
,
Gan
Y.
,
Yang
L.
&
Wu
J.
2023
Impact of climate change on the major eco-hydrological parameters in the Dongting Lake Basin
.
Applied Sciences
13
(
17
),
9515
.
https://doi.org/10.3390/app13179515
.
Zhu
L.
,
Liu
X.
,
Wu
L.
,
Tang
Y.
&
Meng
Y.
2019
Long-term monitoring of cropland change near Dongting Lake, China, using the LandTrendr algorithm with Landsat imagery
.
Remote Sensing
11
(
10
),
1234
.
https://doi.org/10.3390/rs11101234
.
Zhu
Y.
,
Wang
H.
&
Guo
W.
2021
The impacts of water level fluctuations of East Dongting Lake on habitat suitability of migratory birds
.
Ecological Indicators
132
,
108277
.
https://doi.org/10.1016/j.ecolind.2021.108277
.
Zolfagharpour
F.
,
Saghafian
B.
&
Delavar
M.
2022
Hydrological alteration and biodiversity change along the river network caused by anthropogenic activities and climate variability
.
Ecological Processes
11
(
1
),
19
.
https://doi.org/10.1186/s13717-022-00362-4
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).