Next Article in Journal
Multi-Stakeholder Coordinated Operation of Reservoir Considering Irrigation and Ecology
Next Article in Special Issue
Projection of Future Water Resources Carrying Capacity in the Huang-Huai-Hai River Basin under the Impacts of Climate Change and Human Activities
Previous Article in Journal
Metropolis-Hastings Markov Chain Monte Carlo Approach to Simulate van Genuchten Model Parameters for Soil Water Retention Curve
Previous Article in Special Issue
Quantify Runoff Reduction in the Zhang River Due to Water Diversion for Irrigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inverse Trend in Runoff in the Source Regions of the Yangtze and Yellow Rivers under Changing Environments

1
State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, China
2
Research Center for Climate Change, Ministry of Water Resources, Nanjing 210029, China
3
School of Water Resource and Hydropower, Sichuan University, Chengdu 610065, China
4
State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China
5
Yangtze Institute for Conservation and Development, Nanjing 210098, China
6
School of Civil Engineering, Tianjin University, Tianjin 300350, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(12), 1969; https://doi.org/10.3390/w14121969
Submission received: 27 April 2022 / Revised: 14 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022

Abstract

:
The source regions of the Yangtze River (SRYZ) and the Yellow River (SRYR) are sensitive areas of global climate change. Hence, determining the variation characteristics of the runoff and the main influencing factors in this region would be of great significance. In this study, different methods were used to quantify the contributions of climate change and other environmental factors to the runoff variation in the two regions, and the similarities and differences in the driving mechanisms of runoff change in the two regions were explored further. First, the change characteristics of precipitation, potential evapotranspiration, and runoff were analyzed through the observational data of the basin. Then, considering the non-linearity and non-stationarity of the runoff series, a heuristic segmentation algorithm method was used to divide the entire study period into natural and impacted periods. Finally, the effects of climate change and other environmental factors on runoff variation in two regions were evaluated comprehensively using three methods, including the improved double mass curve (IDMC), the slope change ratio of cumulative quantity (SCRCQ), and the Budyko-based elasticity (BBE). Results indicated that the annual precipitation and potential evapotranspiration increased during the study period in the two regions. However, the runoff increased in the SRYZ and decreased in the SRYR. The intra-annual distribution of the runoff in the SRYZ was unimodal during the natural period and bimodal in the SRYR. The mutation test indicated that the change points of annual runoff series in the SRYZ and SRYR occurred in 2004 and 1989, respectively. The attribution analysis methods yielded similar results that climate change had the greatest effect on the runoff variation in the SRYZ, with a contribution of 59.6%~104.6%, and precipitation contributed 65.3%~109.6% of the increase in runoff. In contrast, the runoff variation in the SRYR was mainly controlled by other environmental factors such as permafrost degradation, land desertification, and human water consumption, which contributed 83.7%~96.5% of the decrease in the runoff. The results are meaningful for improving the efficiency of water resources utilization in the SRYZ and SRYR.

1. Introduction

As an essential part of the hydrological cycle, river runoff is one of the forms of water resources and the source of water for human production and living. However, with climate change and the growing impact of human activities, the river runoff of some of the world’s rivers has undergone significant changes [1,2,3]. Statistical results indicated that 22% of the 200 largest rivers in the world showed a significant reduction in annual runoff [4]. Climate change has had significant effect on the water cycle in the basin through factors such as precipitation, temperature, sunshine, wind speed, etc. [5,6], with precipitation being the crucial driving factor [7]. In addition, with rapid economic development, human activities are altering the ecological environment of a basin in the form of greenhouse gas emissions, land-use changes, water conservancy construction, and human water consumption, resulting in an increasingly significant effect on river runoff [8,9]. Therefore, evaluating the influence of climate change and human activities on river runoff is meaningful for rational water resources management, protecting ecological balance, and maintaining sustainable socio-economic development [10].
Three common methods are used for exploring and quantitatively describing the runoff variation in a basin, including statistical methods, hydrological model simulation methods, and Budyko-based elasticity method [11,12]. The empirical statistical methods include the double mass curve and the cumulative slope change methods [13]. Although these methods are simple to calculate, they lack a clear physical basis [14]. The hydrological model simulation methods are based on the actual physical processes, and the physical mechanism is relatively complete [15]. In addition, the simulation accuracy of the hydrological methods is also improved with the innovation and popularization of computer technology. However, the hydrological model contains many parameters [6], and the simulation results depend largely on various assumptions in the model, with a significant degree of uncertainty. In contrast, the Budyko-based elasticity method is based on the hydrothermal balance and has a solid physical basis [1,14]. Therefore, this method is intensively used around the world to quantify the impact of various factors on runoff changes [16,17,18,19], and results indicate that the performance of this method is even comparable to the hydrological models [20,21].
The source regions of the Yangtze River (SRYZ) and Yellow River (SRYR) are located in the Qinghai-Tibet Plateau, which is known as the “China Water Tower” [22,23]. Because of its unique geographical location and climatic conditions, the region is an essential ecological protection barrier for China as well as a “sensitive zone” and “initiation zone” for East Asia and even global climate change [24,25]. However, in recent years, the SRYZ and SRYR have experienced a series of problems such as increasing temperature, decreasing permafrost, retreating glaciers, shrinking lakes, and land desertification [24,26,27,28], which caused significant variation in regional water-cycle processes. More importantly, the runoff in the two regions also affected the water resources in the middle and lower reaches [26,28]. Recent studies have reported the runoff variation characteristics and possible impact factors in the SRYZ and SRYR. Shi et al. [29] analyzed the variation characteristics of precipitation in the Three River Headwaters (TRH) region from 1961 to 2014, and the results indicated an increasing trend in precipitation on the annual scale. Bai et al. [30] suggested that the temperature increase in the SRYZ and SRYR was significantly higher than in most parts of China and the world, and further warming and drying will inevitably have a significant effect on the ecosystems in the source regions. Zhang et al. [26] reported that annual runoff and dry-period runoff in the SRYR decreased significantly during 1958–2005, while the runoff in the SRYZ showed an increasing trend, especially during the flood season. Moreover, some scholars have conducted an attribution analysis of runoff variation in the two regions [31,32,33]. However, few studies have compared and analyzed the similarities and differences in runoff variation and the corresponding drivers in the two regions. Furthermore, previous studies usually used a single evaluation method, which leads to the uncertainty in evaluation results. Glaciers and permafrost are widely distributed in this region, and its change has more obvious influence on runoff than other regions. Hence, glaciers, permafrost, and human activities were categorized as other environmental factors in this study. The impacts of climate change and other environmental factors on runoff variation were distinguished by the Budyko-based elasticity method, and the double mass curve and the cumulative slope change methods were also used to verify the results.
The specific objectives of this paper are to (1) analyze the interannual variation trend of precipitation, potential evapotranspiration, and runoff and determine the mutation points of the annual runoff series; (2) analyze the intra-annual distribution characteristics of runoff and investigate the relationships between runoff versus precipitation and potential evapotranspiration; (3) quantify the contribution of various factors to runoff variation; and (4) comprehensively analyze the internal similarities and differences in the driving mechanism of runoff changes in the SRYZ and SRYR. The specific steps of analysis and calculation in this study are shown in Figure 1. The findings can provide key information for improving the ecological environment of the source regions.

2. Study Area and Datasets

2.1. Study Area

The source regions of the Yangtze River (SRYZ) and the Yellow River (SRYR) are dominated by mountainous landforms with complex topography and an average elevation generally above 4000 m. The climate type of the basin is typical plateau mountain climate, and the vertical climate variation is significant. Glacier resources are abundant, and permafrost is widely distributed in this region.
The area of the SRYZ (32°26′ N~35°45′ N; 90°33′ E~97°15′ E) is about 140,000 km2, accounting for 21.9% of Qinghai province and 8.8% of the Yangtze River basin. There are three important tributaries in the basin, including the Dangqu River, Chumaer River, and Tuotuo River, which converge at the Tongtian River and enter the Yangtze River through Zhimenda station (Figure 2). The average annual precipitation over the SRYZ is about 360 mm during 1966–2013, and the annual average temperature ranges from 3 to 5.5 °C [32,34].
The area of SRYR (32°10′~36°59′ N, 90°54′~103°24′ E) is about 120,000 km2, accounting for about 15% of the Yellow River Basin, and the Tangnaihai station is the export control station in this region (Figure 2). The average annual precipitation in the SRYR varies between 350 mm and 750 mm during 1961–2013, and the mean annual temperature ranges from −4 to −2 °C [35]. This region is a vital runoff generation in the Yellow River Basin, and the mean annual runoff accounts for about 38% of the total runoff in the Yellow River Basin.

2.2. Data Sources

The observational runoff data of Zhimenda station (1966–2013) and Tangnaihai station (1965–2013) were collected from the Hydrological Year Book published by the Hydrology Bureau of the Ministry of Water Resources of China. There are 8 meteorological stations in the SRYZ and 16 stations in the SRYR. The meteorological stations are evenly distributed in space, which can better reflect the spatial distribution characteristics of the climate. A location map of hydrometeorological stations is shown in Figure 2. The meteorological data of the station were collected from the China Meteorological Data Sharing Service Network, including precipitation, temperature, wind speed, and relative humidity, etc. The potential evapotranspiration of the meteorological station was estimated by the Penman–Monteith method. The average potential evapotranspiration and precipitation of the basin area were obtained by the Thiessen polygon method.

3. Methodology

3.1. Abrupt Change Point Detection

Identifying the abrupt change points of runoff series is a necessary step for the attribution analysis of runoff change [36]. The heuristic segmentation algorithm (HSA) method was proposed by Bernaola-Galvan et al. [37] in 2001, which is based on the idea of moving t-test and improved accordingly. The HSA method can not only effectively detect the mean mutation point existing in the original sequence but also greatly enhance its ability to deal with nonlinear and non-stationary sequences [38]. The main calculation steps are as follows:
Suppose a series x(t) composed of N points; the merging deviation SDi at point i can be calculated by:
S D i = N l 1 S l 2 + N r 1 S r 2 N l + N r 2 1 / 2 · 1 N l + 1 N r 1 2             i = 1 , 2 , , N
where Nl and Nr are the length of left and right subsequences of point i, respectively; Sl and Sr are the standard deviations of the left and right subsequences of point i, respectively.
Construct the statistics Ti for the t-test as follows:
T i = u l i u r i S D i
where ul(i) and ur(i) are the mean value of the left and right subsequences of point i, respectively.
The statistical significance P(Tm) corresponding to the maximum value (Tm) in Ti sequence is calculated as:
P T m [ 1 I ν / ν + T m 2 δ ν ,   δ ] γ
γ = 4.19 I n N 11.54
where v = N−2, δ = 0.4, N is the length of the sequence x(t), and Ix(a, b) is an incomplete beta function.
Suppose a critical value of P0; if P(Tm) ≥ P0, the point i is the change point. Repeat the calculation according to the above steps for the two new subsequences, and terminate when the length of the subsequence is less than the minimum segmentation scale t0. The values of P0 range from 0.5 to 0.95, and the value of t0 should not be less than 25 [38].
In order to reduce the uncertainty caused by single-mutation detection method, the Mann–Kendall test (MK) method and ordinal cluster analysis (OC) method were also used in this study. As a parameter-free trend-detection technique, the MK method is widely used in climate change-related research due to its powerful functions and advantages, and the detailed calculation steps have been given in previous studies [11,39]. The OC method is another effective method for mutation test. By estimating the optimal segmentation point, this method makes the sum of squared deviations the smallest among the same class and relatively large between the different classes [32]. After determining the change points of the annual runoff sequence through the above three methods, the study period can be divided into the natural and impacted periods.

3.2. Statistical Methods for Attribution Analysis of Runoff Changes

The slope-change ratio of cumulative quantity (SCRCQ) is a method to analyze the effects of different factors to runoff variation. This method first calculates the change rates of the slope of cumulative annual runoff ( S R ), cumulative annual precipitation ( S P ), and cumulative annual potential evapotranspiration ( S E P ); then calculates the contribution of precipitation and potential evapotranspiration to runoff variation according to the slope change rate; and finally separates the contribution rate of other environmental factors. The calculation details are as follows [11]:
η P = S P S P / S P S R S R / S R × 100 %
η E P = S E P S E P / S E P S R S R / S R × 100 %
η c = η P + η E P
η o = 100 % η c
where and represent data for the natural and impacted periods, respectively; η P and η E P are the contribution rates of precipitation and potential evapotranspiration, respectively; η c and η o are the contribution rates of climate change and other environmental factors, respectively.
The double mass curve method is another popular method besides the SCRCQ method, and the improved double mass curve (IDMC) method used in this study was proposed by Zhang et al. [31]. This method establishes a linear equation of precipitation, potential evapotranspiration, and runoff based on the data of the natural period and then inputs the data of the impacted period into the established relationship to obtain the estimated runoff value, which is compared with the natural runoff to calculate the contribution rate of various factors. The detailed calculation steps can be found in the literature [31].

3.3. Budyko-Based Elasticity Method

The Budyko-based elasticity (BBE) method is based on the principle of water balance. This method calculates the elastic coefficient, which represents the change degree of runoff caused by the change of unit variable, and then quantifies the influence of different factors on runoff. The water balance equation is as follows:
R = P E Δ S
where R , P , and E are the mean annual runoff, precipitation, and actual evapotranspiration during the study period, respectively; Δ S represents the water storage variation in the catchment, which can be neglected in a long study period [32].
Choudhury [40] and Yang et al. [41] further proposed the water–energy balance equation based on Budyko framework, and the Choudhury–Yang equation is as follows:
R = P P × E p P n + E p n 1 / n
where E p is potential evapotranspiration, and n reflects the underlying surface characteristics including land use, topography, and properties of soil, etc.
Schaake [42] first proposed the concept of runoff elasticity, and the sensitivity of runoff to different influencing factor can be calculated by:
ε x = R x × x R
where ε x is the runoff elasticity to P , E p , and n , respectively.
Combined with Equations (10) and (11), the elastic coefficient of ε x can be calculated by:
ε p = 1 + φ n 1 / n + 1 φ n + 1 1 + φ n 1 + φ n 1 / n φ
ε E p = 1 1 + φ n 1 1 + φ n 1 / n
ε n = I n 1 + φ n + φ n I n 1 + φ n n 1 + φ n 1 1 + φ n 1 / n
φ = E P P
Runoff variation caused by different driving factors can be calculated by:
Δ R x = ε x R x Δ x
where Δ R x represents the runoff changes caused by of different factors ( P , E p , and n ); x is the mean annual value ( P , E p , and n ); Δ x is the difference in P , E p , and n between the natural and impacted periods, respectively.
The contributions of P , E p , and n to runoff variation can be calculated by:
η x = Δ R x Δ R × 100 %
where η x is the contributions of P , E p , and n , respectively; Δ R is the change of annual runoff from natural period to impacted period; η c and η o can be calculated by Equations (7) and (8).

4. Result

4.1. Variation Characteristics of Hydrometeorological Variables

The temporal variation trends of different variables in the SRYZ and SRYR are shown in Figure 3. From 1965 to 2013, the precipitation, potential evapotranspiration, and runoff in the SRYZ exhibited an upward trend (Figure 3a–c), with linear trend of 14.38 mm/10a, 2.15 mm/10a, and 6.23 mm/10a, respectively. The interannual runoff variation was evident, with a maximum value of 176.94 mm (2009) and a minimum of 50.52 mm (1979). The MK test results indicated a significant increasing trend for annual runoff and precipitation at 0.05 confidence level. For the SRYR, the annual runoff and precipitation during the study period showed an insignificant decreasing and increasing trend (|Z| < 1.96), respectively (Figure 3d,e), while the potential evapotranspiration increased significantly (|Z| > 1.96) (Figure 3f). The highest mean annual runoff was 267.77 mm (1989) and the lowest was 87 mm (2002).
The statistics of the seasonal runoff variation characteristics in the study area are presented in Table 1. The runoff in the SRYZ indicated an upward trend in all seasons, and the trend in summer was significant, with a linear change rate of 4.377 mm/10a. In the SRYR, obvious differences in seasonal runoff variation trend could be observed, and the summer runoff increased insignificantly, while the runoff in autumn shown a significant decreasing trend (|Z| > 1.96), with the linear trend of −5.1 mm/10a. The reduction in runoff shown in SRYR might be due to the decrease in water volume in autumn, while the runoff variation in the SRYZ might be attributed to the increase in summer runoff. The runoff in the two regions all exhibited an increasing trend in the summer, but the upward trend was greater in the SRYZ.
The intra-annual runoff distribution of the SRYZ and SRYR is shown in Figure 4. A significant change in intra-annual runoff occurred in the SRYZ, indicating a unimodal distribution in the natural and impacted periods, with the peak values occurring in July and August, respectively (Figure 4a). The average monthly runoff during the impacted period increased when compared with the natural period, with the most significant increase in August and October and the smallest in May. For SRYR, the intra-annual runoff showed a bimodal distribution in the natural period, with two peaks in July and September, while the peak value of runoff occurred in July in the impacted period (Figure 4b). The average monthly runoff in the impacted period presented dramatic reductions compared with the natural period, with the greatest reduction in September and the smallest decrease in June.
The scatter plots of runoff versus precipitation and potential evapotranspiration for the SRYZ and SRYR are shown in Figure 5. The runoff was positively correlated with precipitation in the two regions, and the correlation coefficients were all at the 0.01 significance level. In contrast, the runoff and potential evapotranspiration in the two regions were negatively correlated. Moreover, the correlation between runoff and precipitation of the SRYZ during the impacted period was greater than that in the natural period (Figure 5a), and the increase in slope indicated that the same amount of precipitation produced more runoff during the impacted period. However, the runoff generation level of the SRYR during the impacted period was lower than in the natural period (Figure 5c). In terms of the correlation between potential evapotranspiration and runoff, during the impacted period, the correlation level in the SRYZ was higher than that in the SRYR.

4.2. Change Points Analysis of Runoff Series

There are many approaches to detect time-series mutations in climate change research, and using a single method for detection may lead to biased results. Therefore, the heuristic segmentation algorithm (HSA), ordinal cluster analysis, and MK methods were used to determine the mutation points of the runoff series. Based on the above calculation steps, the diagnosed graphs of the HSA method in the SRYZ and SRYR are shown in Figure 6a,d, respectively. The highest T value in the SRYZ occurred in 2004, and the corresponding value of P was 0.991, which was greater than the critical value (0.9); therefore, 2004 was a possible abrupt change point in the SRYZ. Similarly, 1990 was a possible abrupt change point in the SRYR, with the corresponding P(Tm) value of 0.931, which was also greater than the critical value (0.9).
The MK test result of the runoff series in the SRYZ is shown in Figure 6b. The changing trend of the runoff for the past 48 years had obvious fluctuation characteristics. From 1966 to 1980, the UF curve showed a fluctuating state, and the trend of runoff variation was not obvious. From 1985 to 1995, the UF curve declined obviously, indicating that the runoff decreased during this period. After 2005, the UF curve showed a continuous upward trend, indicating a significant upward trend of runoff. MK method detected three possible change points in 2004, 2005, and 2006. Figure 6e represents the MK test results of runoff series in the SRYR. Results indicated that the UF curve increased obviously from 1975 to 1985, indicating that the runoff increased. Conversely, the UF curve exhibited a downward trend from 1985 to 2004 and passed the significance test at the 0.05 level after 2000. The UF and UB curves crossed in 1987 and 1989, indicating that the runoff series of the SRYR may have changed abruptly in 1987 and 1989.
Based on the results of the ordered clustering method, the abrupt change in the runoff series in the SRYZ and SRYR occurred in 2004 and 1989, respectively (Figure 6c,f). By analyzing the test results of the above three methods comprehensively, it was reliable to identify 2004 and 1989 as the change points of the annual runoff series in the SRYZ and SRYR, respectively. Therefore, the period from 1966 to 2003 was the natural period for the SRYZ, and the period of 2004–2013 was regarded as an impacted period. For the SRYR, the natural period was 1965–1988, and the impacted period was 1989–2013.

4.3. Contributions of Different Factors to Runoff Variation

The impacts of different factors on the runoff variation in the SRYZ and SRYR were analyzed quantitatively using the IDMC method, and the parameters estimated using the natural period data are presented in Table 2. The average annual runoff in the SRYZ increased by 27.73 mm during the impacted period, while the mean annual runoff decreased by 24.89 mm in the SRYR. The change in precipitation resulted in an increase of 26.29 mm and 2.14 mm in runoff for the SRYZ and SRYR, respectively. The potential evapotranspiration led to a decrease in runoff of 0.94 mm and 4.78 mm for the SRYZ and SRYR, respectively. The calculated contribution rate of the different factors is shown in Figure 7. The contributions of precipitation and potential evapotranspiration in the SRYZ (SRYR) were 94.8% (−8.6%) and −3.4% (19.2%), respectively. In general, the runoff variation in the SRYZ was attributed to climate change during the impacted period, with a contribution of 91.4%. On the contrary, other environmental factors were responsible for the decrease in the runoff for the SRYR, which accounted for 89.4% of runoff variation.
The cumulative curves of annual precipitation, potential evapotranspiration, and runoff in the SRYZ and SRYR are exhibited in Figure 8, and their linear change slopes in the natural and impacted periods were calculated. All curves had a good-fitting relationship, and the coefficient of determination (R2) was greater than 0.99 (p < 0.01), indicating that more than 99% of the fraction could be reasonably explained. The impacts of climate change and other environmental factors on runoff variation were estimated quantitatively based on the SCRCQ method. The results indicated that precipitation and potential evapotranspiration contributions in the SRYZ were 65.3% and −5.7%, respectively (Figure 7). In the SRYR, the contributions of precipitation and potential evapotranspiration accounted for 8.5% and 7.8%, respectively. Therefore, the contributions of climate change and other environmental factors to runoff variation in the SRYZ (SRYR) were 59.6% (16.3%) and 40.4% (83.7%), respectively.
Based on the BBE method, the calculated elasticity coefficients of runoff to different variables in the SRYZ and SRYR are listed in Table 3. The elasticity coefficients of precipitation, potential evapotranspiration, and surface parameters (n) during the impacted period of the SRYZ were 1.73, −0.73, and −1.4, respectively, indicating that a 1% increase in precipitation, potential evapotranspiration, and n would result in a 1.73% increase, 0.73% decrease, and 1.4% decrease in runoff, respectively. For SRYR, when the potential evapotranspiration and the n increased by 1% during the impacted period, the runoff would decrease by 0.91% and 1.19%, respectively, and the runoff would increase by 1.91% when precipitation increased by 1%. The elastic coefficient of potential evapotranspiration was much lower than that of precipitation, indicating that the sensitivity of runoff to precipitation in the study area was significantly higher than that of potential evapotranspiration.
The BBE method was introduced to further distinguish the impacts of the climatic factors and other environmental factors on runoff variation in the study area, and the results are shown in Table 4 and Figure 7. The relative errors between the estimated runoff changes and the actual runoff changes were within −2%, indicating that this method was reliable and effective for runoff attribution analysis. Precipitation has the most significant effect on runoff variation during the impacted period of the SRYZ, followed by potential evapotranspiration, and n had the least impact, with a contribution of 109.6%, −5.0%, and −4.9%, respectively. In the SRYR, the contributions of precipitation, potential evapotranspiration, and n to runoff variation were −5.3%, 8.4%, and 96.5%, respectively. In general, the contribution rates of climate change and other environmental factors in the SRYZ were 104.6% and −4.9%, respectively. However, in the SRYR, other environmental factors contributed most to the decrease in the runoff, with a contribution of 96.5%.
Although the principles of the IDMC, SCRCQ, and BBE are different, the results of their separation consistently showed that climate change and other environmental factors were the main driving forces for runoff variation in the SRYZ and the SRYR, respectively (Figure 7). The results calculated by the IDMC and BBE were similar, probably because the two methods calculate the contribution based on the relative change in the mean runoff between the impacted and natural periods [31]. However, in the SRYR, the SCRCQ method calculated the contribution value of climate change was lower than the other two methods, probably because the estimated precipitation and potential evapotranspiration data caused errors in the cumulative slope compared with the actual value. Therefore, using multiple methods to separate the effects of different driving forces on runoff variation is more reliable than the results obtained using only one method.

5. Discussion

As the carrier of water resources utilization, runoff variation is crucial to sustainable social and economic development [43], especially in the paramos regions, such as the SRYZ and SRYR [25]. The results of this study showed that climate change had a positive effect on the increase in runoff in the SRYZ, which is consistent with the results of Liu et al. [32] and Zhang et al. [31]. In the SRYR, other environmental factors played a leading role in runoff variation, which is similar to the conclusions obtained by Wang et al. [33].

5.1. Influencing Factors of Runoff in the SRYZ

From 1966 to 2013, the runoff and precipitation in the SRYZ exhibited a significant increasing trend. Precipitation in the SRYZ was mainly concentrated from May to September, which accounted for 88.8% of the annual precipitation, and the runoff accounted for about 72.5% of the total annual runoff [44]. The correlation analysis results showed that the runoff and precipitation were highly correlated during the impacted period, and the summer precipitation showed a significant increase trend based on the MK test. As shown in Figure 3, annual precipitation increased significantly after 2003, but potential evapotranspiration showed a decreasing trend. Precipitation contributed the most (65.3∓112.7%) to the increase in runoff in the SRYZ. The simulation results from the cold region hydrological model (CRHM) also proved that runoff was mainly affected by precipitation in the SRYZ [45], and annual runoff was related to precipitation (62.96%) [46]. In addition, the area of the lake increased by 15.7% in 2010 compared to 1976 due to increased precipitation [47]. Attribution analysis showed that other environmental factors contributed between −5% and 40%, and the contribution of this part comes from glaciers permafrost and human activities. The SRYZ has the highest concentration of glaciers and permafrost in China, accounting for 0.95% and 75% of the whole basin area, respectively [6]. In recent decades, the temperature in the SRYZ increased significantly, with an increase of 0.6 °C/10a, which was about four times of the global mean [48]. Under the background of increasing temperature, the melting of snow and glaciers in the basin increased, and the river runoff also increased [34]. According to the statistics of the Chinese Glacier Inventory (CGI), the glaciers and snow meltwater in the SRYZ accounted for about 9.2% of the total discharge [49]. Han et al. [50] further quantified the effects of snow and glacier meltwater on runoff variation during 2003–2014, which accounted for about 7% and 5% of total runoff, respectively. The degradation of permafrost is another consequence of the rising temperatures. Based on the statistics of the second Tibet Plateau Examination, the area of permafrost decreased by 16% during the past 50 years [34]. From 1981 to 2015, the mean decline rate of the maximum frozen depth of seasonally frozen ground in the SRYZ was 12.3 cm/10a, and the permafrost degradation contributed −11.6 mm to the runoff changes [51]. The decrease in permafrost thickness may increase the soil water storage capacity, resulting in a corresponding reduction in surface runoff. The parameter n in the BBE method exhibited a decreasing trend in the impacted period compared with the natural period, which was related to the vegetation cover of the basin and possibly related to the change of glaciers and permafrost. Despite the absence of the large-scale water-conservancy facilities and a small population in the study area, it has a large area of grassland, and the vegetation degradation caused by overgrazing also affects runoff changes.

5.2. Influencing Factors of Runoff in the SRYR

The SRYR is an important water conservation area in the Yellow River Basin. Under the combined influence of climate change and other environmental factors, the runoff in the SRYR showed an insignificant decreasing trend from 1965 to 2013. On the one hand, the increase in temperature caused significant changes in the glaciers and permafrost in this region. The glacier area in this region only accounts for 0.11% of the basin area, and the contribution of glacier meltwater to the runoff was small (less than 2% of the annual flow) [52]. However, the permafrost area in the SRYR accounts for more than 80%, the permafrost depth of this region exhibited a significant downward trend driven by the increase in temperature, the regional permafrost area ratio decreased by 1.1% per year, and the mean permafrost depth decreased by 0.012 m per year between 1981 and 2015 [53]. Wang et al. [54] further found that, compared with the period from 1990 to 2003, permafrost degradation contributed 38.5% to annual runoff changes during 2004–2015. On the other hand, with the development of the social economy, human activities in this region intensified, which had a greater effect on the runoff change. Ji et al. [55] suggested that the spring, summer, and winter runoff variation were mainly controlled by human activities. The area of arable land increased by about 10% from 1980 to 2010. Irrigation of crops (wheat and rapeseed) in spring requires large amounts of water, but under the influence of global climate change, the water deficit of crops is increasing further, and the demand for irrigation water will further increase [55,56]. Between 1961 and 2004, residential water consumption increased obviously as the population increased nearly twofold, which also contributed to the runoff decrease. The domestic water and livestock drinking water accounted for 1.7% and 9.3% of the natural runoff, respectively [57]. The degradation of grasslands and desertification of land because of overgrazing and cultivation of arable land, and the land desertification in the eastern and southern parts of the SRYR showed an upward trend [58]. However, since the implementation of the Three-River Headwaters Region ecological protection project in 2005, the vegetation degradation and land desertification had improved significantly, and the water storage capacity of the soil improved [58]. In the past decades, the number of lakes in the SRYR decreased significantly, and the two major lakes (Eling Lake and Zhaling Lake) experienced a decline in water levels, which also can be seen as a response to the decrease in the runoff [57]. Additionally, the Huangheyuan hydropower station was completed in 2000, and the total reservoir capacity is 25.01 × 108 m3, accounting for 4.31% of the annual runoff of the entire Yellow River Basin [59]. Considering the need for flood control and power generation, it contributed to the decrease in the runoff. From 1965 to 2013, the contribution rate of climate change to the reduction of runoff in the SRYR ranged from 3.1% to 16.3%. Runoff in this area mainly mainly comes from precipitation, accounting for about 64.1% of all runoff sources [28,60]. As shown in Figure 3, the runoff severely decreased in the 1990s. The high number of light precipitation events and the long duration of precipitation in the 1990s caused runoff decreased significantly during the flood season [61,62]. In addition, with the increase in temperature, the effect of evapotranspiration increased and the soil surface moisture decreased quickly [63], and evapotranspiration played an essential role in affecting runoff variation in the SRYR [35]. Hu et al. [64] also suggested that the decrease in precipitation and the increase in temperature during 1959–2008 were important reasons for the decreased runoff in the SRYR.

5.3. Comparison of Runoff Variation

The monthly runoff in the SRYR has two peak values during the natural period but only one peak value in the SRYZ. The difference in the intra-annual distribution of precipitation between the two regions was responsible for the different monthly runoff processes [65]. In addition, because the temperature in the SRYZ was generally lower than that in the SRYR, the difference in the formation time of glacier and snowmelt runoff was also one of the reasons for the different intra-annual runoff distribution. In spring and winter, the flow of the SRYR was twice that of the SRYZ, the difference in flow was relatively small in summer, and the difference in flow between the two regions mainly occurred in autumn [66]. The increase in summer runoff in the SRYZ may lead to more severe summer floods in downstream areas, while the decrease in autumn runoff in the SRYR may lead to more severe water shortages in the downstream area in autumn and winter [66]. In summary, the precipitation in the SRYR increased slowly during the study period, but the evaporation increased significantly. Combined with the effects of regional grassland degradation, permafrost ablation, the retreat of glaciers and lakes, and human water consumption, the runoff generation conditions on the underlying surface of the basin deteriorated. Under the same precipitation conditions, the amount of runoff generated will be significantly decreased. This finding was also proven by the results in Figure 5a, where the runoff–precipitation relationship curve for the impacted period was lower than that in the natural period. However, in the SRYZ, runoff is more sensitive to climate change due to its higher elevation, fewer changes in the underlying surface, and fewer human activities. After 2004, annual precipitation increased significantly, and potential evapotranspiration decreased (Figure 3), coupled with the fact that it has more glacial meltwater than the SRYR, thereby leading to a significant increase in runoff. The runoff in the SRYR was more affected by human activities due to its low altitude: the closer to the river source, the greater the impact of climate change, and the closer to the middle river, the greater influence of human activities [33].

5.4. Uncertainties and Limitations

The calculation results in this study are based on the measured data of the gauge station, and the updating errors in observational instruments and measurement rules may lead to non-homogeneity in climate data series, which will affect the accurate description of climate change. However, existing studies have shown that the non-homogeneity of climate data has less effect on large-scale climate change research [67]. This study analyzes the overall climate change characteristics of the study area from a regional perspective, which reduces the possible uncertainty of the climate sequence. Then, the distribution of meteorological stations at high altitudes is sparse, and the errors in measuring instruments and manual operations also lead to uncertainties in the collected hydrometeorological data. It is also difficult to accurately estimate potential evapotranspiration. The relative contribution of the influencing factors was calculated by the BBE method and verified by the IDMC and SCRCQ methods. These methods all assume that climatic and human influences are independent of each other, thereby isolating their contribution to runoff variation. In fact, human activities and climate change are interactive, and ignoring the effects of human activities during the natural period will introduce errors [68]. Additionally, the impact of the spatial and temporal distribution of precipitation on runoff was not considered in the model; it is possibly another error [69]. Finally, only precipitation and potential evapotranspiration were used to represent the impact of climate change on runoff, and the influence of other factors on runoff variation was regarded as the change of parameter n. Thus, it is still a challenge to further distinguish the contribution of various factors on runoff.

6. Conclusions

This study selected the SRYZ and SRYR as the study areas. The annual and seasonal variation characteristics of runoff were analyzed through statistical methods, and the three methods (HSA, MK, and OC) were used to further detect the change points of the runoff series. Then, the IDMC, SCRCQ, and BBE methods were used to quantify the impacts of climate change (precipitation and potential evapotranspiration) and other environmental factors on runoff variation. Finally, the driving mechanism of runoff change in the two regions was analyzed. The main findings are as follows:
(1)
The annual runoff and precipitation in the SRYZ indicated a significant increasing trend during 1966–2013, while the potential evapotranspiration increased insignificantly. From 1966 to 2013, the annual runoff for SRYR exhibited an insignificant decreasing trend, with the rate of −6.45 mm/10a, whereas potential evapotranspiration increased significantly, with the rate of 7.89 mm/10a.
(2)
Runoff in the SRYZ showed an upward trend in all seasons, with the summer runoff increasing significantly. The summer runoff in the SRYR exhibited an insignificant upward trend, and the runoff in the autumn decreased significantly. During the natural period, the intra-annual distribution of runoff in the SRYZ was unimodal and was bimodal in the SRYR.
(3)
By evaluating the three mutation detection results comprehensively, it was determined that the annual runoff series in the SRYZ and SRYR changed abruptly in 2004 and 1989, respectively.
(4)
Climate change was the dominant factor for runoff variation in the SRYZ, with a contribution of 59.6%~104.6%, while the runoff in the SRYR was mainly controlled by other environmental factors, contributing 83.7%~96.5%. Moreover, precipitation changes accounted for 65.3%~109.6% of the increase in annual runoff in the SRYZ, and the potential evapotranspiration contributed 7.8%−19.2% of runoff variation in the SRYR. Generally, the significant increase of runoff in the SRYZ could be attributed to the precipitation and glacial meltwater, while the decrease of runoff in the SRYR was the combined effects of the permafrost degradation, land desertification and human water consumption.

Author Contributions

H.W. wrote the paper with contributions from all co-authors; Z.B. and G.W. conceived and designed the research; J.W. and C.L. supervised the research; Y.Y., C.Z., D.Z. and S.L. analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (grant no. 2017YFA0605002, 2017YFA0605004), the National Natural Science Foundation of China (grant no. 41961124007, 52121006, 41830863, 51879164), and the “Six top talents” in Jiangsu province (grant no. RJFW-031).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author or the first author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yan, D.; Lai, Z.Z.; Ji, G.X. Using Budyko-Type Equations for Separating the Impacts of Climate and Vegetation Change on Runoff in the Source Area of the Yellow River. Water 2020, 12, 3418. [Google Scholar] [CrossRef]
  2. Milly, P.C.D.; Dunne, K.A.; Vecchia, A.V. Global pattern of trends in streamflow and water availability in a changing climate. Nature 2005, 438, 347–350. [Google Scholar] [CrossRef] [PubMed]
  3. Saifullah, M.; Li, Z.J.; Li, Q.L.; Zaman, M.; Hashim, S. Quantitative Estimation of the Impact of Precipitation and Land Surface Change on Hydrological Processes through Statistical Modeling. Adv. Meteorol. 2016, 2016, 6130179. [Google Scholar] [CrossRef] [Green Version]
  4. Walling, D.E.; Fang, D. Recent trends in the suspended sediment loads of the world’s rivers. Glob. Planet. Chang. 2003, 39, 111–126. [Google Scholar] [CrossRef]
  5. Mlynski, D.; Cebulska, M.; Walega, A. Trends, Variability, and Seasonality of Maximum Annual Daily Precipitation in the Upper Vistula Basin, Poland. Atmosphere 2018, 9, 313. [Google Scholar] [CrossRef] [Green Version]
  6. Guo, M.J.; Li, J.; Wang, Y.S.; Bai, P.; Wang, J.W. Distinguishing the Relative Contribution of Environmental Factors to Runoff Change in the Headwaters of the Yangtze River. Water 2019, 11, 1432. [Google Scholar] [CrossRef] [Green Version]
  7. Abbas, S.A.; Xuan, Y.Q. Impact of Precipitation Pre-Processing Methods on Hydrological Model Performance using High-Resolution Gridded Dataset. Water 2020, 12, 840. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, P.P.; Cai, Y.P.; Yang, W.; Yi, Y.J.; Yang, Z.F. Climatic and anthropogenic impacts on water and sediment generation in the middle reach of the Jinsha River Basin. River Res. Appl. 2020, 36, 338–350. [Google Scholar] [CrossRef]
  9. Zeng, F.; Ma, M.G.; Di, D.R.; Shi, W.Y. Separating the Impacts of Climate Change and Human Activities on Runoff: A Review of Method and Application. Water 2020, 12, 2201. [Google Scholar] [CrossRef]
  10. Wang, G.Q.; Yan, X.L.; Zhang, J.Y.; Liu, C.S.; Jin, J.L.; Liu, Y.L.; Bao, Z.X. Detecting evolution trends in the recorded runoffs from the major rivers in China during 1950-2010. J. Water Clim. Chang. 2013, 4, 252–264. [Google Scholar] [CrossRef]
  11. Liang, S.Q.; Wang, W.S.; Zhang, D.; Li, Y.Q.; Wang, G.Q. Quantifying the Impacts of Climate Change and Human Activities on Runoff Variation: Case Study of the Upstream of Minjiang River, China. J. Hydrol. Eng. 2020, 25, 05020025. [Google Scholar] [CrossRef]
  12. Zhang, K.; Ruben, G.B.; Li, X.; Li, Z.J.; Yu, Z.B.; Xia, J.; Dong, Z.C. A comprehensive assessment framework for quantifying climatic and anthropogenic contributions to streamflow changes: A case study in a typical semi -arid North China basin. Environ. Modell. Softw. 2020, 128, 104704. [Google Scholar] [CrossRef]
  13. Mokhtar, A.; He, H.M.; Zhao, H.F.; Keo, S.; Bai, C.Y.; Zhang, C.J.; Ma, Y.; Ibrahim, A.; Li, Y.; Li, F.; et al. Risks to water resources and development of a management strategy in the river basins of the Hengduan Mountains, Southwest China. Environ. Sci. Water Res. Technol. 2020, 6, 656–678. [Google Scholar] [CrossRef]
  14. Li, L.J.; Li, B.; Liang, L.Q.; Li, J.Y.; Liu, Y.M. Effect of climate change and land use on stream flow in the upper and middle reaches of the Taoer River, northeastern China. For. Stud. China 2010, 12, 107–115. [Google Scholar] [CrossRef]
  15. Jiang, S.H.; Ren, L.L.; Yong, B.; Fu, C.B.; Yang, X.L. Analyzing the effects of climate variability and human activities on runoff from the Laohahe basin in northern China. Hydrol. Res. 2012, 43, 3–13. [Google Scholar] [CrossRef]
  16. Li, H.J.; Shi, C.X.; Zhang, Y.S.; Ning, T.T.; Sun, P.C.; Liu, X.F.; Ma, X.Q.; Liu, W.; Collins, A.L. Using the Budyko hypothesis for detecting and attributing changes in runoff to climate and vegetation change in the soft sandstone area of the middle Yellow River basin, China. Sci. Total Environ. 2020, 703, 135588. [Google Scholar] [CrossRef] [PubMed]
  17. Gan, G.J.; Liu, Y.B.; Sun, G. Understanding interactions among climate, water, and vegetation with the Budyko framework. Earth-Sci. Rev. 2020, 212, 103451. [Google Scholar] [CrossRef]
  18. Jiang, C.; Xiong, L.H.; Wang, D.B.; Liu, P.; Guo, S.L.; Xu, C.Y. Separating the impacts of climate change and human activities on runoff using the Budyko-type equations with time-varying parameters. J. Hydrol. 2015, 522, 326–338. [Google Scholar] [CrossRef]
  19. Liang, W.; Bai, D.; Wang, F.Y.; Fu, B.J.; Yan, J.P.; Wang, S.; Yang, Y.; Long, D.; Feng, M. Quantifying the impacts of climate change and ecological restoration on streamflow changes based on a Budyko hydrological model in China’s Loess Plateau. Water Resour. Res. 2015, 51, 6500–6519. [Google Scholar] [CrossRef]
  20. Zuo, D.P.; Xu, Z.X.; Wu, W.; Zhao, J.; Zhao, F.F. Identification of Streamflow Response to Climate Change and Human Activities in the Wei River Basin, China. Water Resour. Manag. 2014, 28, 833–851. [Google Scholar] [CrossRef]
  21. Liang, W.; Bai, D.; Jin, Z.; You, Y.C.; Li, J.X.; Yang, Y.T. A Study on the Streamflow Change and its Relationship with Climate Change and Ecological Restoration Measures in a Sediment Concentrated Region in the Loess Plateau, China. Water Resour. Manag. 2015, 29, 4045–4060. [Google Scholar] [CrossRef]
  22. Ma, Y.M.; Han, C.B.; Zhong, L.; Wang, B.B.; Zhu, Z.K.; Wang, Y.J.; Zhang, L.; Meng, C.; Xu, C.; Amatya, P.M. Using MODIS and AVHRR data to determine regional surface heating field and heat flux distributions over the heterogeneous landscape of the Tibetan Plateau. Theor. Appl. Climatol. 2014, 117, 643–652. [Google Scholar] [CrossRef] [Green Version]
  23. Xu, M.; Kang, S.C.; Chen, X.L.; Wu, H.; Wang, X.Y.; Su, Z.B. Detection of hydrological variations and their impacts on vegetation from multiple satellite observations in the Three-River Source Region of the Tibetan Plateau. Sci. Total Environ. 2018, 639, 1220–1232. [Google Scholar] [CrossRef]
  24. Chu, H.B.; Wei, J.H.; Qiu, J.; Li, Q.; Wang, G.Q. Identification of the impact of climate change and human activities on rainfall-runoff relationship variation in the Three-River Headwaters region. Ecol. Indic. 2019, 106, 105516. [Google Scholar] [CrossRef]
  25. Meng, X.H.; Chen, H.; Li, Z.G.; Zhao, L.; Zhou, B.R.; Lv, S.H.; Deng, M.S.; Liu, Y.M.; Li, G.W. Review of Climate Change and Its Environmental Influence on the Three-River Regions. Plateau Meteorol. 2020, 39, 1133–1143. (In Chinese) [Google Scholar]
  26. Zhang, Y.Y.; Zhang, S.F.; Zhai, X.Y.; Xia, J. Runoff variation and its response to climate change in the Three Rivers Source Region. J. Geogr. Sci. 2012, 22, 781–794. [Google Scholar] [CrossRef]
  27. Luo, D.L.; Jin, H.J.; Lu, L.Z.; Zhou, J. Spatiotemporal changes in extreme ground surface temperatures and the relationship with air temperatures in the Three-River Source Regions during 1980-2013. Theor. Appl. Climatol. 2016, 123, 885–897. [Google Scholar] [CrossRef]
  28. Yang, J.P.; Ding, Y.J.; Chen, R.S. Climatic causes of ecological and environmental variations in the source regions of the Yangtze and Yellow Rivers of China. Environ. Geol. 2007, 53, 113–121. [Google Scholar] [CrossRef]
  29. Shi, H.Y.; Li, T.J.; Wei, J.H.; Fu, W.; Wang, G.Q. Spatial and temporal characteristics of precipitation over the Three-River Headwaters region during 1961-2014. J. Hydrol. Reg. Stud. 2016, 6, 52–65. [Google Scholar] [CrossRef] [Green Version]
  30. Bai, X.L.; Wei, J.H.; Xie, H.W. Characteristics of wetness/dryness variation and their influences in the Three River Headwaters region. Acta Ecol. Sin. 2017, 37, 8397–8410. (In Chinese) [Google Scholar]
  31. Zhang, D.; Wang, W.S.; Yu, S.Y.; Liang, S.Q.; Hu, Q.F. Assessment of the Contributions of Climate Change and Human Activities to Runoff Variation: Case Study in Four Subregions of the Jinsha River Basin, China. J. Hydrol. Eng. 2021, 26, 05021024. [Google Scholar] [CrossRef]
  32. Liu, J.L.; Chen, J.; Xu, J.J.; Lin, Y.R.; Yuan, Z.; Zhou, M.Y. Attribution of Runoff Variation in the Headwaters of the Yangtze River Based on the Budyko Hypothesis. Int. J. Environ. Res. Public Health 2019, 16, 2506. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, Y.Q.; Yuan, Z.; Xu, J.J.; Yan, B.; Hong, X.F. Research on the attribution identification of source runoff variation in the Yellow River Source Region based on water and energy balance model. In Proceedings of the 5th International Conference on Water Resource and Environment (WRE), Macao, China, 16–19 July 2019; p. 344. [Google Scholar]
  34. Li, Z.X.; Li, Z.J.; Feng, Q.; Zhang, B.J.; Gui, J.; Xue, J.; Gao, W.D. Runoff dominated by supra-permafrost water in the source region of the Yangtze River using environmental isotopes. J. Hydrol. 2020, 582, 124506. [Google Scholar] [CrossRef]
  35. Meng, F.C.; Su, F.G.; Yang, D.Q.; Tong, K.; Hao, Z.C. Impacts of recent climate change on the hydrology in the source region of the Yellow River basin. J. Hydrol. Reg. Stud. 2016, 6, 66–81. [Google Scholar] [CrossRef] [Green Version]
  36. Jiang, C.; Li, D.Q.; Gao, Y.N.; Liu, W.F.; Zhang, L.B. Impact of climate variability and anthropogenic activity on streamflow in the Three Rivers Headwater Region, Tibetan Plateau, China. Theor. Appl. Climatol. 2017, 129, 667–681. [Google Scholar] [CrossRef]
  37. Bernaola-Galvan, P.; Ivanov, P.C.; Amaral, L.A.N.; Stanley, H.E. Scale invariance in the nonstationarity of human heart rate. Phys. Rev. Lett. 2001, 87, 168105. [Google Scholar] [CrossRef] [Green Version]
  38. Huang, S.Z.; Huang, Q.; Leng, G.Y.; Liu, S.Y. A nonparametric multivariate standardized drought index for characterizing socioeconomic drought: A case study in the Heihe River Basin. J. Hydrol. 2016, 542, 875–883. [Google Scholar] [CrossRef]
  39. Zhang, H.W.; Xu, W.; Xu, X.T.; Lu, B.H. Responses of Streamflow to Climate Change and Human Activities in a River Basin, Northeast China. Adv. Meteorol. 2017, 2017, 1023821. [Google Scholar] [CrossRef]
  40. Choudhury, B.J. Evaluation of an empirical equation for annual evaporation using field observations and results from a biophysical model. J. Hydrol. 1999, 216, 99–110. [Google Scholar] [CrossRef]
  41. Yang, H.B.; Yang, D.W.; Lei, Z.D.; Sun, F.B. New analytical derivation of the mean annual water-energy balance equation. Water Resour. Res. 2008, 44, 893–897. [Google Scholar] [CrossRef]
  42. Schaake, J.C. Climate change and U.S. water resources: From climate to flow. Wiley 1990, 89, 129–130. [Google Scholar]
  43. Fowler, K.J.A.; Peel, M.C.; Western, A.W.; Zhang, L.; Peterson, T.J. Simulating runoff under changing climatic conditions: Revisiting an apparent deficiency of conceptual rainfall-runoff models. Water Resour. Res. 2016, 52, 1820–1846. [Google Scholar] [CrossRef] [Green Version]
  44. Jiang, Y.C. Analysis of runoff evolution and driving factors in the source region of Yangtze River under climate change. Master’s Thesis, Xihua University, Chengdu, China, 2021. (In Chinese). [Google Scholar]
  45. Han, L.; Song, K.C.; Zhang, W.J.; Liu, L.; Jiang, H.R. Temporal and Spatial Variations of Hydrological Factors in the Source Area of the Yangtze River and Its Responses to Climate Change. Mt. Res. 2017, 35, 129–141. (In Chinese) [Google Scholar]
  46. Qian, K.Z.; Wang, X.S.; Lv, J.J.; Wan, L. The wavelet correlative analysis of climatic impacts on runoff in the source region of Yangtze River, in China. Int. J. Climatol. 2014, 34, 2019–2032. [Google Scholar] [CrossRef]
  47. Duan, S.Q.; Liu, T.; Cao, G.C.; Wu, Q.; Li, Y.; Liu, X.S. Expansion of the lakes and its causes in the source region of the Yangtze River. Arid Zone Res. 2015, 32, 15–22. (In Chinese) [Google Scholar]
  48. Zhang, Y.; Kang, S.; Grigholm, B.; Zhang, Y.; Kaspari, S.; Morgenstern, U.; Ren, J.; Qin, D.; Mayewski, P.A.; Zhang, Q.; et al. Twentieth-century warming preserved in a Geladaindong mountain ice core, central Tibetan Plateau. Ann. Glaciol. 2016, 57, 70–80. [Google Scholar] [CrossRef] [Green Version]
  49. Wu, S.S.; Yao, Z.J.; Huang, H.Q.; Liu, Z.F.; Chen, Y.S. Glacier retreat and its effect on stream flow in the source region of the Yangtze River. J. Geogr. Sci. 2013, 23, 849–859. [Google Scholar] [CrossRef]
  50. Han, P.F.; Long, D.; Han, Z.Y.; Du, M.D.; Dai, L.Y.; Hao, X.H. Improved understanding of snowmelt runoff from the headwaters of China’s Yangtze River using remotely sensed snow products and hydrological modeling. Remote Sens. Environ. 2019, 224, 44–59. [Google Scholar] [CrossRef]
  51. Shi, R.J.; Yang, H.B.; Yang, D.W. Spatiotemporal variations in frozen ground and their impacts on hydrological components in the source region of the Yangtze River. J. Hydrol. 2020, 590, 125237. [Google Scholar] [CrossRef]
  52. Zhang, L.L.; Su, F.G.; Yang, D.Q.; Hao, Z.C.; Tong, K. Discharge regime and simulation for the upstream of major rivers over Tibetan Plateau. J. Geophys. Res. Atmos. 2013, 118, 8500–8518. [Google Scholar] [CrossRef]
  53. Qin, Y.; Yang, D.W.; Gao, B.; Wang, T.H.; Chen, J.S.; Chen, Y.; Wang, Y.H.; Zheng, G.H. Impacts of climate warming on the frozen ground and eco-hydrology in the Yellow River source region, China. Sci. Total Environ. 2017, 605, 830–841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, T.H.; Yang, H.B.; Yang, D.W.; Qin, Y.; Wang, Y.H. Quantifying the streamflow response to frozen ground degradation in the source region of the Yellow River within the Budyko framework. J. Hydrol. 2018, 558, 301–313. [Google Scholar] [CrossRef]
  55. Ji, G.X.; Wu, L.Y.; Wang, L.D.; Yan, D.; Lai, Z.Z. Attribution Analysis of Seasonal Runoff in the Source Region of the Yellow River Using Seasonal Budyko Hypothesis. Land 2021, 10, 542. [Google Scholar] [CrossRef]
  56. Haddeland, I.; Heinke, J.; Biemans, H.; Eisner, S.; Florke, M.; Hanasaki, N.; Konzmann, M.; Ludwig, F.; Masaki, Y.; Schewe, J.; et al. Global water resources affected by human interventions and climate change. Proc. Natl. Acad. Sci. USA 2014, 111, 3251–3256. [Google Scholar] [CrossRef] [Green Version]
  57. Zhang, S.F.; Jia, S.F.; Liu, C.M.; Cao, W.B.; Hao, F.H.; Liu, J.Y.; Yan, H.Y. Study on the changes of water cycle and its impacts in the source region of the Yellow River. Sci. China Ser. E Technol. Sci. 2004, 47, 142–151. [Google Scholar] [CrossRef]
  58. Guo, B.; Wei, C.X.; Yu, Y.; Liu, Y.F.; Li, J.L.; Meng, C.; Cai, Y.M. The dominant influencing factors of desertification changes in the source region of Yellow River: Climate change or human activity? Sci. Total Environ. 2022, 813, 152512. [Google Scholar] [CrossRef]
  59. Ma, L.J.; Liu, Z.; Zhao, B.F.; Lyu, J.W.; Zheng, F.M.; Xu, W.; Gan, X.B. Variations of runoff and sediment and their response to human activities in the source region of the Yellow River, China. Environ. Earth Sci. 2021, 80, 552. [Google Scholar] [CrossRef]
  60. Lan, Y.C.; Kang, E.S.; Ma, Q.J.; Yang, W.H.; Yao, Z.Z. Runoff Forecast Model for Inflow to the Longyangxia Reservoir in the Upper Yellow River Basin During Spring. J. Glaciol. Geocryol. 1999, 21, 391–395. (In Chinese) [Google Scholar]
  61. Zhao, E.R.; Chen, H.C.; Zhu, S.L.; Zhao, L.Q. Analysis on Variations of Runoff and Reasons of Headwater Region of the Yellow River. Yellow River 2007, 29, 15–16. (In Chinese) [Google Scholar]
  62. Zhou, D.G.; Huang, R.H. Exploration of Reason of Runoff Decrease in the Source Regions of the Yellow River. Clim. Environ. Res. 2006, 11, 302–309. (In Chinese) [Google Scholar]
  63. Wu, X.L.; Zhang, X.; Xiang, X.H.; Zhang, K.; Jin, H.J.; Chen, X.; Wang, C.H.; Shao, Q.Q.; Hua, W.J. Changing runoff generation in the source area of the Yellow River: Mechanisms, seasonal patterns and trends. Cold Reg. Sci. Technol. 2018, 155, 58–68. [Google Scholar] [CrossRef]
  64. Hu, Y.R.; Maskey, S.; Uhlenbrook, S.; Zhao, H.L. Streamflow trends and climate linkages in the source region of the Yellow River, China. Hydrol. Process. 2011, 25, 3399–3411. [Google Scholar] [CrossRef]
  65. Xie, C.W.; Ding, Y.J.; Liu, S.Y.; Wang, G.X. Comparison Analysis of Runoff Change in the Source Regions of the Yangtze and Yellow Rivers. J. Glaciol. Geocryol. 2003, 25, 414–422. (In Chinese) [Google Scholar]
  66. Wu, Y.H.; Qin, N.S.; Peng, X.B. Similarities and Differences of Runoff Changes in the Source Region of the Yangtze River and the Yellow River. Plateau Mt. Meteorol. Res. 2021, 41, 43–47. (In Chinese) [Google Scholar]
  67. Liang, X.W.; Yang, M.X.; Wan, G.N.; Wang, X.J.; Li, Q. Research on the homogeneity of air temperature series over Tibetan Plateau. J. Glaciol. Geocryol. 2015, 37, 275–285. [Google Scholar]
  68. Hu, J.; Ma, J.; Nie, C.; Xue, L.; Zhang, Y.; Ni, F.; Deng, Y.; Liu, J.; Zhou, D.; Li, L.; et al. Attribution Analysis of Runoff Change in Min-Tuo River Basin based on SWAT model simulations, China. Sci. Rep. 2020, 10, 2900. [Google Scholar] [CrossRef] [Green Version]
  69. Xu, X.Y.; Yang, D.W.; Yang, H.B.; Lei, H.M. Attribution analysis based on the Budyko hypothesis for detecting the dominant cause of runoff decline in Haihe basin. J. Hydrol. 2014, 510, 530–540. [Google Scholar] [CrossRef]
Figure 1. Flowchart of quantifying the contribution of various factors to runoff variation in the SRYZ and SRYR.
Figure 1. Flowchart of quantifying the contribution of various factors to runoff variation in the SRYZ and SRYR.
Water 14 01969 g001
Figure 2. Geographical location of study area and hydrometeorological stations.
Figure 2. Geographical location of study area and hydrometeorological stations.
Water 14 01969 g002
Figure 3. Annual runoff (R), precipitation (P) and potential evapotranspiration (Ep) variation characteristics in the SRYZ (ac) and SRYR (df) (Z represents the statistics of MK test).
Figure 3. Annual runoff (R), precipitation (P) and potential evapotranspiration (Ep) variation characteristics in the SRYZ (ac) and SRYR (df) (Z represents the statistics of MK test).
Water 14 01969 g003
Figure 4. Average monthly runoff between the natural and impacted periods in the SRYZ (a) and SRYR (b).
Figure 4. Average monthly runoff between the natural and impacted periods in the SRYZ (a) and SRYR (b).
Water 14 01969 g004
Figure 5. Correlation analysis of precipitation and potential evapotranspiration versus runoff during the natural and impacted periods in the SRYZ (a,b) and SRYR (c,d) (r represents the correlation coefficient; p is the corresponding significance level; black and red represent the natural period and impacted period, respectively).
Figure 5. Correlation analysis of precipitation and potential evapotranspiration versus runoff during the natural and impacted periods in the SRYZ (a,b) and SRYR (c,d) (r represents the correlation coefficient; p is the corresponding significance level; black and red represent the natural period and impacted period, respectively).
Water 14 01969 g005
Figure 6. Result of the heuristic segmentation algorithm (HSA), Mann–Kendall (MK), and ordinal cluster (OC) mutation test of the annual runoff series in the SRYZ (ac) and SRYR (df) (T represents the test statistics in the HSA method; Sn(τ) represents the squared deviations in the OC method).
Figure 6. Result of the heuristic segmentation algorithm (HSA), Mann–Kendall (MK), and ordinal cluster (OC) mutation test of the annual runoff series in the SRYZ (ac) and SRYR (df) (T represents the test statistics in the HSA method; Sn(τ) represents the squared deviations in the OC method).
Water 14 01969 g006
Figure 7. Contributions of climate change (precipitation and potential evapotranspiration (b)) and other environmental factors (a) to runoff variation in the SRYZ and SRYR (BBE, Budyko-based elasticity; IDMC, improved double mass curve; SCRCQ, slope-change ratio of cumulative quantity).
Figure 7. Contributions of climate change (precipitation and potential evapotranspiration (b)) and other environmental factors (a) to runoff variation in the SRYZ and SRYR (BBE, Budyko-based elasticity; IDMC, improved double mass curve; SCRCQ, slope-change ratio of cumulative quantity).
Water 14 01969 g007
Figure 8. Cumulative curves of annual runoff, precipitation, and potential evapotranspiration for the SRYZ (a) and SRYR (b).
Figure 8. Cumulative curves of annual runoff, precipitation, and potential evapotranspiration for the SRYZ (a) and SRYR (b).
Water 14 01969 g008
Table 1. The MK test results of annual and seasonal runoff series in the SRYZ and SRYR.
Table 1. The MK test results of annual and seasonal runoff series in the SRYZ and SRYR.
RegionResultSpringSummerAutumnWinter
SRYZTrend (mm/10a)0.2774.3771.4820.126
Z1.2892.142 *1.5731.484
SRYRTrend (mm/10a)−1.2260.057−5.100−0.177
Z−1.845−0.259−1.974 *−0.948
Note: Z represents the statistics of MK test; * indicates a significance level of 0.05.
Table 2. Parameters and runoff variation calculated by the improved double mass curve method.
Table 2. Parameters and runoff variation calculated by the improved double mass curve method.
RegionNatural
Period
Impacted
Period
k1k2c Δ R P   ( mm ) Δ R Ep   ( mm ) Δ R   ( mm )
SRYZ1966–20032004–20130.41−0.06−27.7726.29−0.9427.73
SRYR1965–19881989–20130.93−0.40−11.932.14−4.78−24.89
Table 3. Multi-year elastic coefficient of runoff on various impact factors in the SRYZ and SRYR.
Table 3. Multi-year elastic coefficient of runoff on various impact factors in the SRYZ and SRYR.
RegionPeriodR/mmP/mmEp/mmnEp/PεPεEpεn
SRYZ1966–200387.99328.81827.731.042.521.76−0.76−1.55
2004–2013115.72392.20844.281.052.151.73−0.73−1.40
SRYR1965–1988178.47525.28769.521.191.461.76−0.76−1.10
1989–2013153.58527.58781.391.351.481.91−0.91−1.19
Table 4. Contribution of various impact factors to runoff variation in the SRYZ and SRYR based on the BBE method.
Table 4. Contribution of various impact factors to runoff variation in the SRYZ and SRYR based on the BBE method.
RegionNatural
Period
Impacted
Period
ΔRPΔREpΔRnRe
(%)
ηP
(%)
ηEp
(%)
ηo
(%)
ηc
(%)
SRYZ1966–20032004–201330.40−1.40−1.36−0.3109.6−5.0−4.90.3
SRYR1965–19881989–20131.32−2.09−24.02−0.4−5.38.496.50.4
Note: Re represents the relative deviation between simulated runoff variation and actual variation; ηc is the contribution of the coupled other factors.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, H.; Bao, Z.; Wang, J.; Wang, G.; Liu, C.; Yang, Y.; Zhang, D.; Liang, S.; Zhang, C. Inverse Trend in Runoff in the Source Regions of the Yangtze and Yellow Rivers under Changing Environments. Water 2022, 14, 1969. https://doi.org/10.3390/w14121969

AMA Style

Wu H, Bao Z, Wang J, Wang G, Liu C, Yang Y, Zhang D, Liang S, Zhang C. Inverse Trend in Runoff in the Source Regions of the Yangtze and Yellow Rivers under Changing Environments. Water. 2022; 14(12):1969. https://doi.org/10.3390/w14121969

Chicago/Turabian Style

Wu, Houfa, Zhenxin Bao, Jie Wang, Guoqing Wang, Cuishan Liu, Yanqing Yang, Dan Zhang, Shuqi Liang, and Chengfeng Zhang. 2022. "Inverse Trend in Runoff in the Source Regions of the Yangtze and Yellow Rivers under Changing Environments" Water 14, no. 12: 1969. https://doi.org/10.3390/w14121969

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

Article Metrics

Back to TopTop