Next Article in Journal
Discrete Analysis of Local Scour Morphology of Bridge Piers Affected by Sediment Storage Dam
Previous Article in Journal
Indication of Long-Term Changes of Algae Communities in a Hydrologically Transformed Estuary Sasyk, Black Sea, Ukraine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Runoff Changes in the Wei River Basin, China: Confronting Climate Change and Human Activities

1
State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Xianyang 712100, China
2
Institute of Soil and Water Conservation, Northwest A & F University, Xianyang 712100, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
Yangtze River Basin Monitoring Center Station for Soil and Water Conservation, Changjiang Water Resources Commission, Wuhan 430012, China
*
Authors to whom correspondence should be addressed.
Water 2023, 15(11), 2081; https://doi.org/10.3390/w15112081
Submission received: 29 March 2023 / Revised: 23 May 2023 / Accepted: 25 May 2023 / Published: 30 May 2023

Abstract

:
Abrupt runoff reduction in the Wei River Basin (WRB) has attracted extensive attention owing to climate change and human activities. Nevertheless, previous studies have inadequately assessed the respective contributions of climate variability and human activities to runoff change on different spatial scales. Using Mann–Kendall and Pettitt’s methods, this study identified long-term (1970–2018) changes in hydro-meteorological variables. Furthermore, the Budyko-based method was used to quantify the influence of climate change and human activities on runoff change at different spatial scales of the WRB, including the whole WRB, three sub-basins, and sixteen catchments. The results show that a significant decrease trend was identified in runoff at different spatial scales within the WRB. Runoff in almost all catchments showed a significant downward trend. Temperature, potential evapotranspiration, and the parameter n showed significant increases, whereas no significant trend in precipitation was observed. The change in runoff was mainly concentrated in the mid-1990s and early 2000s. Anthropogenic activities produced a larger impact on runoff decrease in the WRB (62.8%), three sub-basins (53.9% to 65.8%), and most catchments (–47.0% to 147.3%) than climate change. Dramatic catchment characteristic changes caused by large-scale human activities were the predominant reason of runoff reduction in the WRB. Our findings provide a comprehensive understanding of the dominate factors causing runoff change and contribute to water resource management and ecosystem health conservation in the WRB.

1. Introduction

Runoff, a key component of the hydrological cycle, is strongly influenced by climate variability and anthropogenic activities, and therefore, changes dramatically at various spatial and temporal scales [1,2]. Climate change includes the changes in precipitation, temperature, and potential evaporation, which largely affect the future runoff of a basin [3]. In addition to climate change, anthropogenic factors such as ecological restoration programs, urbanization, industrialization, reservoir operation, large-scale irrigation, and drainage have strongly shaped the runoff regime [4]. Previous studies have shown that global runoff has increased over the 20th century [5]. However, river flow has decreased significantly in many basins worldwide, such as the Spanish basins [6], the Lake Issyk-Kul Basin [7], East African basins [8], as well as some basins in China, such as the Yellow River Basin [9] and the Shiyang River Basin [10], especially in arid and semi-arid regions [11,12]. The continuing decline in runoff has caused a severe water crisis and ecosystem imbalance, threatening the high state and sustainable development of the economy and society [13].
According to the IPCC report [14], the global climate has been undergoing extensive and complicated change in the past century, which aggravates the uneven distribution of precipitation (P) [15,16] and potential evapotranspiration (PET) [17], and consequently influences water balance. For rivers in tropical areas at low latitudes, the long-term tendency of runoff change is generally consistent with that of P [18,19]; however, for rivers in middle latitudes and some alpine regions, the correlation between runoff and P is weak [20,21]. The direct or indirect interference of runoff by human activities, including deforestation, afforestation, reservoir operation, large-scale irrigation, and urbanization, is an important reason for this regional difference [9,22,23]. Quantitative assessment of the impacts of climate variability and anthropogenic activities on runoff regimes at the regional scale is therefore vital to water resource management.
Several methods have been applied to analyze the influence of climate variability and human intervention on alterations in runoff [24,25,26]. Empirical statistical methods can achieve results quickly and easily; however, it is difficult to interpret complex and non-linear hydrological behaviors. Hydrological models, although powerful, can be constrained by the inherent variability and uncertainty in the impact of human activities on hydrological systems [27], and the difficulty in obtaining numerous input data related to human activities, such as water abstraction and irrigation data [28,29]. With the Budyko [30] hypothesis, various Budyko-based approaches have been widely used to explore the non-linear interactions among climate, landscape, and hydrology due to its straightforward calculations and physical mechanisms [31]. It suggests that the ratio of precipitation to potential evapotranspiration (P/PET) is the primary factor determining the hydrological balance at basin level and thereby dictates the dryness or wetness of a basin. Early empirical formulations [30,32,33] represented the general hydrological reactions of catchments to global climate forcing, as they were based on continental or global runoff information. In 1981, Fu [34] incorporated the catchment characteristic parameter (parameter n) into catchment water balance. Parameter n comprehensively reflects catchment characteristics, including the characteristics of soil, topography, geologic properties, and vegetation cover [31,35]. In 2008, Yang et al. [36] derived the Choudhury–Yang formula based on Choudhury’s [37] proposed empirical theoretical equation. The early consistency of the rich theory and the later improvements to the methodology increased the credibility of the method. Therefore, coupled with the elasticity coefficient method, the Choudhury–Yang method has been successfully employed to analyze and separate the impacts of climate change and human activities on changes in runoff [3,27,38,39].
As the largest tributary of the Yellow River, the Wei River discharges approximately 19.7% of water into the Yellow River every year, which plays a significant role in the management and sustainable development of water resources [40]. It provides ~85% of the water supply for irrigation, domestic, and industrial use, and nourishes more than 22 million people in this basin, making it an essential factor for regional economic development [41]. However, significant runoff reductions have been reported in many studies [40,42]. With the growth in population and development of the economy, the demand for water resources tends to increase, and water scarcity becomes more severe [43,44]. Moreover, the northern part of the Wei River Basin (WRB) is located in the Loess Plateau, which suffers from serious soil erosion. Catchment characteristics and runoff patterns have been significantly affected by the execution of ecological restoration projects [44,45,46]. Despite the fact that many studies have separated the influence of climate change and human intervention on WRB runoff [24,40,46], different studies have arrived at disparate conclusions due to their varying methods and periods. Moreover, most studies only focused on alterations in runoff of the entire basin or a single hydrological station but ignored the runoff change in smaller catchments. Therefore, we conducted a comprehensive attribution analysis of runoff change at different catchment scales of the WRB from 1970 to 2018 by using the Budyko-based method. The main objectives of this study were to (1) depict the trend of runoff, P, PET, temperature, and parameter n in different catchments and the change for abrupt year; (2) assess the spatial-temporal variability of the runoff elasticity coefficients to P, PET, and parameter n in different catchments; and (3) quantitatively estimate the changes in runoff attributable to climate variables and land cover and clarify the possible influencing factors at different spatial scales.

2. Study Area and Data

2.1. Study Area

The WRB (34–38° N; 104–111° E) is located in the lower midstream of the Yellow River, covering an area of approximately 1.35 × 105 km2 (Figure 1). The Wei River rises in the Niaoshu Mountain of Gansu Province, through three provinces (Gansu, Shaanxi, and Ningxia), and across the Loess Plateau before entering the Yellow River. Its two primary tributaries are the Jing River and Beiluo River, which divide the Wei River Basin into three sub-basins: the Jing River (JR) basin, Beiluo River (BLR) basin, and the upstream of the Wei River (UWR) basin. Traditionally, the region is named the Jing–Luo–Wei Region, and the catchment characteristics of the three tributaries are very different in space. The WRB is situated in the continental monsoon climate zone, and experiences high temperatures and abundant P during summer, whereas winter brings low temperatures and scarce P. The average annual P is approximately 536 mm. The distribution of P is uneven across space, with the majority of it (over 60%) being concentrated within the months from June to September (Figure 2).

2.2. Data

Annual runoff data for 16 hydrological stations from 1970 to 2018 were made available by the Yellow River Water Resources Commission. Daily meteorological data (1970–2018) from 22 meteorological stations were collected from the National Meteorological Information Centre of China (NMIC) (https://www.nmic.cn/ (accessed on 1 January 2020)), and consisted of precipitation, wind speed, relative humidity, sunshine duration, and average, maximum, and minimum temperature. The information of the hydrological stations is shown in Table S1 in the Supplementary Materials. Land cover datasets from 1980, 1995, and 2018 were available at the Chinese Academy of Sciences Resource and Environmental Science Data Center (http://www.resdc.cn/ (accessed on 1 January 2020)). The land cover was divided into six classifications: cropland, forest, settlement, grassland, water, and others, with a spatial resolution of 1 km. The 15 day Normalized Difference Vegetation Index (NDVI) dataset was acquired from the Global Inventory Modeling and Mapping Studies (GIMMS) NDVI3g vegetation index dataset designed by NASA (https://ladsweb.modaps.eosdis.nasa.gov/ (accessed on 1 January 2020)), with a spatial resolution of 500 m from 1982 to 2015. The annual mean and maximum NDVI values were used in this study.
The entire basin was divided into 16 catchments, considering the distribution of hydrological stations, to examine the spatial features of runoff change. The catchment controlled by a single hydrological station was divided into one catchment, and the runoff in that catchment was equal to the observed runoff of this hydrological station. For instance, the runoff value of catchment 1 was equal to the runoff observed at the Wuqi station. The catchment controlled by adjacent hydrological stations was divided into one catchment, and the runoff value of this catchment was calculated using the runoff difference between the upstream and downstream hydrological stations. For example, the runoff of catchment 2 was the observed runoff of Liujiahe minus that of Wuqi. Table S2 in the Supplementary Materials presented comprehensive details about the catchments.

3. Methodology

3.1. Trend Detection

The non-parametric Mann–Kendall (M–K) method was implemented to detect the trends for annual runoff, P, PET, and parameter n in the WRB [47,48]. This method was selected for trend detection in hydro-meteorological records because its robustness demands non-normality and censoring data.
For a given dataset X (i = 1, 2, 3, …, n), the test statistic (S) is given by
S = i = 1 n 1 i + 1 i sgn ( X i X j )
where
sgn ( x j x i ) = { 1 0 1 x j x i > 0 x j x i = 0 x j x i < 0
A standardized statistic Z is obtained by
Z = { ( S 1 ) / n ( n 1 ) ( 2 n + 5 ) / 18 0 ( S + 1 ) / n ( n 1 ) ( 2 n + 5 ) / 18 S > 0 S = 0 S < 0
To reduce the influence of the serial correlation, the trend-free pre-whitening test was applied to identify any temporal trends in hydro-meteorological variables [49], which is defined as
Y t = X t r X t 1
where Yt is the de-trended and pre-whitened series, referred to as the residual series, and Xt represents the original time series. If the lag-one serial coefficient r of the de-trended series is statistically significant at the 5% level, it means that the serial correlation has been eliminated. Additionally, a Mann–Kendall test was conducted on the whitened series to obtain a Z-statistic. A positive Z-value indicates an upward trend, while a negative value indicates a downward trend [40]. |Z| < 1.96 means no significant trend was detected (p > 5%). |Z| ≥ 1.96 and |Z| ≥ 2.58 indicate that the significance test has passed with a confidence level of 95% and 99%.

3.2. Breakpoint Analysis

In this study [49], Pettitt’s test was applied to identify the abrupt year in the hydrological data by utilizing a significance level of 5% through a non-parametric approach [50]. This test detects the occurrence of breakpoints using a non-parametric method at a significance level of 5%. This test applies a function of the Mann–Whitney statistic (Ut,n), which has been mentioned in the Mann–Kendall test. When a significant breakpoint is detected, it indicates that the two subseries (x1, …, xt and xt+1, xt+2, …, xn) are from different populations. Ut,n is expressed as follows:
U t , n = i = 1 t j = 1 n sgn ( X t X j )   If   t = 2 ,   ,   n
The test statistic Kn is defined as maximum value of |Ut,n|
K n = m a x | U t , n |
The associated probability (p) with Kn is given by:
p = exp ( 6 ( K n ) 2 n 3 + n 2 )
The null hypothesis of Pettitt’s test is that there is no change point.

3.3. Assessing the Impacts of Climatic and Anthropogenic Factors on Runoff Change

3.3.1. Budyko Framework

For a given catchment, the long-term water balance equation can be written as
R = P A E T Δ S
where R, P, AET, and ΔS are denoted as the runoff depth (mm), precipitation (mm), actual evapotranspiration (mm), and water storage change, respectively. Here, the variation in water storage can be ignored for the given hydrological year and is assigned a value of zero.
Budyko [30] proposed that the ratio of PET to P over a long-term average can be used to describe the hydrological drought conditions of a region. This basis for the hypothesis is that the water and energy balance processes are regulated by climate and topography, and thus any changes in these factors will lead to corresponding variations in runoff. Over the years, this hypothesis has been widely utilized in hydrology and hydrological modelling to understand the impacts of climate and land surface characteristics on runoff processes [30,36,37]. To separate the influence on runoff independently, this study adopted the elasticity coefficient method coupled with the Budyko frameworks. AET can be derived from the Choudhury–Yang equation:
A E T = P × P E T ( P n + P E T n ) 1 / n
where PET represents the potential evapotranspiration, and parameter n comprehensively reflects the catchment characteristics, including the characteristics of soil, topography, geologic properties, and vegetation cover [35].

3.3.2. Elasticity Coefficient

Combining Equations (8) and (9), we obtain an equation R = f (P, PET, n). P is the observed annual mean precipitation. PET can be estimated using the FAO Penman–Monteith [51], which is currently recognized as an effective method for estimating PET. The method takes into account variations in meteorological factors such as temperature, humidity, wind speed, and solar radiation, as well as the complexity of vegetation cover, and is used to estimate PET in different regions/basins to help design irrigation systems and improve crop management. The parameter n can be derived from this equation using MATLAB R2022b software.
The total differential equation of R with different factors can be written as
d R = f P d P + f P E T d P E T + f n d n
The elasticity of runoff to a specific independent variable change is the degree of the response of runoff to the changes in variables, e.g., P, PET, or n, which can be estimated by the proportional change in simulated runoff. The elasticity of runoff is given as
ε x i = R x i × x i R
where ε x i is the elasticity of runoff to P, PET, or n. Therefore, Equation (10) is written as
d R = ε x P R x P d x P + ε x P E T R x P E T d x P E T + ε x n R x n d x n
The elasticity of runoff is expressed as [3]:
ε p = 1 [ ( P E T / P ) n 1 + ( P E T / P ) n ] 1 / n + 1 1 [ ( P E T / P ) n 1 + ( P E T / P ) n ] 1 / n
ε P E T = 1 1 + ( P E T / P ) n 1 1 [ 1 + ( P E T / P ) n ( P E T / P ) n ] 1 / n
ε n = A B [ 1 + ( P E T / P ) n ] 1 / n 1
A = P n ln ( P ) + P E T ln ( P E T ) P n + P E T
B = ln ( P n + P E T n ) n

3.3.3. Attribution Analysis of Runoff Changes

According to the breakpoint, the whole time series can be divided into two periods: “period 1” which represents the base period before the abrupt year, and “period 2” which represents the human impact period after the abrupt year. Thus, the change in annual runoff from period 1 to period 2 is given by
Δ R = R 2 R 1
where ΔR is the difference between the mean annual observed runoff in the two different periods, R1 is the average runoff in the base period, and R2 is the average runoff in the human impact period.
The total runoff change over a long period can be estimated as follows:
Δ R total = Δ R clima + Δ R h u m a n
where ∆Rclima is the change in runoff depth due to climatic variables, including the change in runoff caused by P (∆RP) and PET (∆RPET). ∆Rhuman is the runoff change caused by human activities, equal to the change in runoff caused by the catchment characteristics of parameter n. The simulated runoff change can be calculated as follows:
Δ R P = ε P R P Δ P Δ R P E T = ε P E T R P E T Δ P E T Δ R n = ε n R n Δ n
∆P, ∆PET, and ∆n represent the difference in P, PET, and parameter n from period 1 to period 2.
The relative contribution is the percentage of P, PET, and parameter n to the runoff change, which is expressed as
C p = Δ R P Δ R × 100 % C P E T = Δ R P E T Δ R × 100 % C n = Δ R n Δ R × 100 %
where CP, CPET, and Cn are the contribution proportions of climate variation (P and PET) and anthropogenic activities to runoff change, respectively.
The performance of the Budyko-based method can be evaluated by Nash–Sutcliffe efficiency values (NSE) and the absolute values of the relative error (Re).
N S E = 1 i = 1 n ( Δ R 1 Δ R 2 ) 2 i = 1 n ( Δ R 1 Δ R 1 a ) 2
R e = Δ R 1 Δ R 2 Δ R 1 × 100 %
∆R1 is the observed runoff change, ∆R2 is the runoff change estimated by Budyko-based method. ∆R1a is the average observed runoff change.

4. Results

4.1. Trends in Hydro-Meteorological Factors

To better understand the variation in runoff and its influencing factors at different scales, we analyzed the trends in the annual runoff, P, PET, temperature, and parameter n within the three sub-basins (UWR, JR, and BLR) of the WRB from 1970 to 2018 (Figure 3). Significantly decreasing trends were observed for runoff (p < 0.01), with average rates of 0.87 mm/y, 0.60 mm/y, and 0.19 mm/y in the UWR, JR, and BLR, respectively. P presented downward trends in the UWR and BLR and an upward trend in the BLR, though all trends were found to be non-significant. PET in the UWR and JR showed upward trends (no significance at p > 0.05), and both increased by 1.15 mm/y, while the PET of the BLR presented downward trends. Temperature showed a significant increase trend (p < 0.01) in all three regions with average rates of 0.03 °C/y, 0.04 °C/y, and 0.02 °C/y in the UWR, JR, and BLR, respectively. In terms of parameter n, the UWR and JR displayed a significant increasing trend (p < 0.01), with an average rate of 0.01 mm/y and 0.02 mm/y, respectively, while no significant trend was found for parameter n in the BLR.
The trends in annual runoff, P, PET, and parameter n from 1970 to 2018 are shown in Figure 4. Out of all the catchments, 12 catchments (which make up 75% of the catchments) exhibited a notable decrease trend (p < 0.05) in annual runoff. P showed an increase trend in most catchments, but all of them were statistically insignificant. PET in half of the catchments showed significant increasing trends. A significant increase trend (p < 0.05) in temperature was observed in 15 catchments (accounting for 94% of catchments). In addition, parameter n showed a significant increase trend (p < 0.05) in 75% of the catchments. The significant increase trend of parameter n indicates that the disturbance of catchment characteristics is gradually intensified.

4.2. Abrupt Years of Runoff Change

The abrupt years of annual runoff were determined using Pettitt’s method in the WRB and three sub-basins (UWR, JR, and BLR) from 1970 to 2018 (Table 1). In the WRB, the runoff underwent an abrupt change in 1993 (p < 0.05). In the three sub-basins, significant abrupt years (p < 0.01) were detected in the UWR (1993), JR (1996), and BLR (1994). It was found that the change in runoff and P decreased after the abrupt year, while that in PET and parameter n increased in the WRB and the three sub-basins. 62.5% of the catchments exhibited significant abrupt changes (p < 0.05) around the early 1990s and around 2000. The change in hydro-meteorological characteristics and the contribution of climatic and anthropogenic factors to runoff changes in 16 catchments are shown in Table S3 in the Supplementary Materials.

4.3. The Elasticity Coefficient of Runoff to Climate Change and Anthropogenic Factors

Table 1 presents the elasticity coefficients of runoff with respect to P, PET, and parameter n for the WRB and the three sub-basins. Compared to the base period, there was a noticeable increase in the absolute values of the runoff elasticity coefficient compared to the base period, implying an increasing sensitivity of runoff change to climate variables and parameter n. For the WRB, the elasticity coefficients after the abrupt year were 3.28 for P, –2.28 for PET, and –2.38 for parameter n. These values suggest that a 10% reduction in P would cause a 32.8% decline in runoff, whereas a corresponding decrease in PET or parameter n would result in an increase of 22.8% and 23.8%, respectively. The results indicated that the highest and positive response of runoff was to P, while it exhibited a negative response to PET and parameter n.
The elasticity coefficients of runoff in 16 catchments concerning P, PET, and parameter n are shown in Figure 5. The spatial distribution patterns of the runoff elasticity coefficient in response to P exhibit similarities with those of PET, but are different in polarity. The elasticity coefficient of runoff to P ranged from 2.6 to 4.2 (Figure 5a). Contrary to P, the elasticity coefficients (–3.2 to –1.6) of runoff were negative for PET (Figure 5b). In the middle and downstream regions of the BLR, the absolute elasticity coefficients of runoff with respect to precipitation (P) and potential evapotranspiration (PET) were found to be higher compared to other regions across the WRB. This suggests that changes in P and PET have a greater impact on the amount of runoff in these regions. The elasticity coefficients of runoff to parameter n ranged from –3.4 to –1.6, and higher absolute values were found in the JR and BLR (Figure 5c). This finding suggested that alterations in catchment characteristics have a greater impact on the runoff variations in these specific regions as compared to others.

4.4. Evaluating the Impacts of Climatic Anthropogenic Factors on Runoff

The performance of the Budyko-based method is shown in Figure 6. Upon comparing the observed and simulated changes in runoff, it was found that the regression slopes of both close to 1. The NSE values were greater than 0.99, and the absolute values of the Re were less than 5%. This indicates that the method has a strong modeling capability for runoff changes in most catchments.
Table 2 demonstrates the runoff alterations caused by hydro-meteorological variables (P and PET) and the catchment characteristics parameter n in the WRB and the three sub-basins. The findings indicate that the decreased runoff, induced by P, PET, and parameter n in the WRB and the three sub-basins, may be affected by the decrease in P and the increase in PET and parameter n. The values obtained for parameter n are related to integrated catchment characteristics, including vegetation cover, soil properties, and topography [52,53]. Regarding soil properties and topography as static variables, land cover change is an important factor affecting catchment hydrological behavior [3]. Thus, the impact of various human activities such as afforestation, deforestation, soil and water conservation projects, and surface water consumption will ultimately be indicated by parameter n. The study found that climate variability, which includes P and PET, contributed to 37.2% of the changes in runoff, while human activities accounted for 62.8%. Furthermore, the impact of human activities on the decrease in runoff was particularly significant in each sub-basin, with a contribution rate exceeding 50%.
∆RP, ∆RPET, and ∆Rn are the runoff changes induced by the change in P, PET, parameter n after the abrupt year, respectively. CP, CPET, and Cn are the relative contributions of P, PET, parameter n to runoff change.
It can be observed from Figure 7 that the runoff changes before and after the abrupt year induced by hydro-meteorological variables and parameter n is spatially inconsistent. The change in P positively impacted runoff (runoff increase) in most catchments, while changes in PET and parameter n had a negative influence on runoff (runoff decrease) in all catchments. The runoff change induced by P, PET, and parameter n between two periods ranged from –8.5 to 15.6 mm, –7.5 to –1.2 mm, and –38.1 to 12.5 mm, respectively. Figure 8 shows the contribution proportions of climate variables (P and PET) and catchment characteristics (parameter n) to runoff change. Except for catchments 3 and 4, all catchments showed that catchment characteristic changes induced by human activities contributed more than climate variables to runoff change. Therefore, catchment characteristic changes were the dominant factors shaping runoff in most catchments, while climate change contributed the most to runoff in two catchments (3 and 4) in the middle of the BLR. For sub-region 3, the increased PET (38.5 mm) from the base period to the human impact period had the most significant impact on runoff reduction, accounting for 113%. For catchment 4, P (52.2 mm) increased by 52 mm after the abrupt year, and the runoff change induced by P was 8.3 mm, which accounted for 107%.

5. Discussion

5.1. Contribution of the Climatic Factors to Runoff Change

Water resource systems are dramatically influenced by climatic variation, especially in arid and semi-arid areas [7,54]. Previous studies have reported that variations in temperature, P, and PET are manifestations of climate change, which has a significant impact on runoff in the WRB [55,56]. It has been shown that P and temperature are positively associated with terrestrial water storage, and that decreased terrestrial water storage is driven by decreased P and increased PET [57]. Our study indicates that the temperature of the WRB has shown a notable rise. Pang et al. [58] observed that the temperature in the WRB has increased by 0.42 °C/10a. The high temperature accelerates evapotranspiration, reducing the magnitude of water entering the river channel and peak runoff [38]. The drastic runoff reduction caused by increasing PET in the WRB was also observed by Zhao et al. [55]. However, as the main source of runoff, our study showed that the annual P declined 12.3 mm after 1993, but there was no noticeable decreasing trend in the WRB from 1970 to 2018. Similar decreasing trends in P were also detected in other studies in the WRB [46,59]. Our results indicated that P showed the highest and most positive correlation with runoff change, but the change in P (20.3%) contributing to the decrease in runoff was slightly greater than that for PET (16.9%) in this region. This is in agreement with the previous results, which identified that changes in PET have a weaker impact on runoff change than those in P [52]. In addition, climate change has led to the reduction in runoff in the WRB from 1970 to 2018, although its impact is far less than that of human activities. In contrast, Luo et al. [60] demonstrated that climate change has increased the annual runoff in the source area of the Yangtze River, which has offset a large part of the annual runoff reduced by the ecological protection plan. This may be attributed to the decrease in P and the increase in PET in the WRB.

5.2. Contribution of the Catchment Characteristics’ Change to Runoff Change and Potential Anthropogenic Activities

To investigate runoff changes caused by human activities in the WRB, we evaluated three satellite images that recorded land cover. Croplands, forests, and grasslands constituted the major types of land cover in the WRB. Our results indicate the change point in runoff regime in the WRB around 1995, leading to this study focusing on land cover change after the transition (from 1995 to 2018). During this period, the main land cover conversion in the WRB was concentrated from cropland to grassland, settlement, and forest; from grassland to cropland and forest; and from forest to grassland, with a total conversion area of 54,851 km2 (Figure 9a). Consequently, we selected the six main conversion types and computed the proportion of each land conversion type in the overall conversion area based on the total conversion area (Figure 9b). Conversion between cropland and grassland had the highest percentage (around 30%) followed by conversion between forest and cropland (about 11%), while transformations of some cropland to settlement and forest were also presented, which account for 6.4% and 6.5%. Meanwhile, the net change in four main land cover types underwent greater changes from 1995 to 2018 were analyzed: cropland, settlement, forest, and grassland. During this period, cropland decreased by 4203 km2, while forest and grassland increased by 1162 km2 and 1114 km2, respectively. These changes can be partly attributed to large-scale ecological conservation programs, such as the Grain for Green (GFG) projects implemented in 1999, which led to the reconversion of substantial areas of cropland to forest and grassland, substantially increasing the vegetation cover in the WRB [61]. Chen et al. [62] also observed a dramatic change in land cover owing to an increase in vegetation restoration. Moreover, we observed a significant increase in settlement area, with an increase of 2012 km2, which may be related to the urbanization of the local area.
Given the execution of large-scale ecological conservation programs, we assessed vegetation coverage variations using the annual mean NDVI and annual maximum NDVI of the WRB from 1982 to 2015 (Figure 10). The results show that vegetation cover increased significantly (p < 0.05) after 2000, which may be attributed to the GFG project executed after 1999. Research has suggested that large-scale vegetation restoration will locally decrease the annual mean water availability and runoff [63,64]. Increased vegetation cover allows more precipitation to be intercepted by the vegetation canopy and root system [65], thereby diminishing terrestrial runoff [66]. Yang et al. [67] found that increased vegetation coverage was the main reason for runoff reduction in the middle reaches of the Yellow River Basin, which decreased runoff by 35% from 2000 to 2005. Similar results were obtained by Liu et al. [68] and Zuo et al. [42], who concluded that vegetation coverage was the dominant factor responsible for runoff reduction. In addition, large-scale soil and water conservation measures (e.g., afforestation, terraces, fish-scale pits, and grazing ban) have played a vital role in vegetation restoration in the WRB. Since the implementation of these measures in the 1980s, runoff has declined by approximately 9.9 × 108 m3/year (15.5% of runoff at Huaxian station) in the WRB [55]. Moreover, interesting results were observed showing that 62.5% of the catchments have significant abrupt changes (p < 0.05), which may attribute to the soil and water conservations measures implemented after the 1990s and Grain for Green project since 2000 [9]. However, for the runoff regime in the WRB and three sub-basins, there was no significant abrupt change detected around 2000. It may be that vegetation restoration is a slow and long-term process, and the hydrological effect of vegetation on runoff in large catchment has a hysteresis effect [69]. Zhao et al. [70] also found that the response of runoff change to vegetation cover change in large catchment (>1000 km2) weakened with the increase in catchment scale.
The impacts of human activities on catchment characteristics are complex, especially in the context of the extensive implementation of ecological conservation programs. In addition to these ecological conservation programs, the WRB has also witnessed many other intensive human activities in the past 30 years, such as urbanization, and large-scale infrastructure construction (water infrastructures). A large amount of cropland has been converted into settlement (3503 km2) from 1995 to 2018; this change mainly occurred in the Guanzhong Plain owing to urbanization. Intensified urbanization would have also resulted in changes in catchment characteristics by increasing impervious areas, consequently leading to increases in runoff and decreases in infiltration [22]. Additionally, the reservoir/dams construction has a direct impact on the catchment characteristics of parameter n, mainly by affecting the vegetation coverage in the irrigation area to affect the runoff [38]. Our results showed that human activities play a dominant role in runoff decreasing, which is consistent with the findings in other catchments such as Loess Plateau [22], Yellow River [71], and Yangtze River [60].

5.3. Limitations of the Study

This study has some notable limitations. Equation (10) is a first-order approximation, which may lead to errors in the calculated impacts of climatic variations and anthropogenic activities. In addition, this framework of contribution estimation supposes that anthropogenic activities and climatic variation are completely independent, but it is presently difficult to easily separate and decompose them [71]. Both climate change and human activities will affect the change in catchment characteristic parameter n. Attributing land cover change to human activities will overestimate the impact of local human activities on runoff change [3].

6. Conclusions

This study quantified the impact of climate change and human-induced land cover variations on runoff change in the WRB using the elasticity coefficient method coupled with decomposition methods based on the Budyko framework. The results show the following:
(1)
The annual runoff and P presented a statistically downward trend and a significant decrease in runoff was detected in all three sub-basins (p < 0.01) and 75% of the catchments (p < 0.05). However, PET and parameter n showed an increasing trend in the WRB, three sub-basins, and most catchments.
(2)
The significant abrupt years of runoff change were observed in the mid-1990s at the WRB and the three sub-basins (p < 0.01); 62.5% of the catchments exhibited significant abrupt changes (p < 0.05) in the early 1990s and around 2000.
(3)
The results of the elasticity coefficient revealed that P had the highest positive correlation with runoff change, but PET and parameter n had a negative correlation. During the human impact period, a 10% decrease in P would result in a 32.8% drop in runoff, whereas a 10% decrease in PET and parameter n would induce a 22.8% and 23.8% increase in runoff, respectively.
(4)
Land cover change induced by human activities had a greater impact on runoff decline than climatic factors in the entire basin, three sub-basins, and most catchments. The impacts of P, PET, and parameter n resulted in average runoff reductions of 20.3, 16.9, and 62.8%, respectively, in the WRB.
These findings highlight the importance of systematically investigating the impacts of climate change and land cover on runoff reduction using the Budyko-based method at different catchment scales. This study could contribute to better quantitative clarification of the driving factors of runoff variation and provide a reference foundation for sustainable water management and ecosystem health conservation across the WRB and midstream of the Yellow River. However, considering the complex interactions among diverse human activities, accurately separating the contributions of other human activities remains challenging. More human activity metrics (e.g., GDP and population) should be considered when interpreting runoff changes, and more methods should be applied to verify the accuracy of the results in future studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15112081/s1, Table S1: Information on the hydrological stations; Table S2: Information on the hydro-meteorological characteristics of the 16 sub-regions in the Wei River Basin; Table S3: The change in hydro-meteorological characteristics and contribution of climatic and human activity factors to runoff changes in catchments.

Author Contributions

Conceptualization, R.X., X.M. and P.G.; Methodology, R.X., C.G., D.Q. and C.W.; Resources, C.G., D.Q., C.W., X.M. and P.G.; Writing—original draft, R.X.; Writing—review & editing, D.Q. and P.G.; Supervision, X.M.; Funding acquisition, P.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program of China, Grant/Award Number: 2016YFC0501707.

Data Availability Statement

The daily series of meteorological data used in this study are available on National Meteorological Information Centre of China (NMIC) (https://www.nmic.cn/), while the data sharing of runoff are not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. PAGES2k Consortium. A global multiproxy database for temperature reconstructions of the Common Era. Sci. Data 2017, 4, 170088. [Google Scholar] [CrossRef] [PubMed]
  2. Martin, J.T.; Pederson, G.T.; Woodhouse, C.A.; Cook, E.R.; McCabe, G.J.; Anchukaitis, K.J.; Wise, E.K.; Erger, P.J.; Dolan, L.; McGuire, M.; et al. Increased drought severity tracks warming in the United States’ largest river basin. Sci. Data 2020, 117, 11328–11336. [Google Scholar] [CrossRef]
  3. 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]
  4. Wagener, T.; Sivapalan, M.; Troch, P.A.; McGlynn, B.L.; Harman, C.J.; Gupta, H.V.; Kumar, P.; Rao, P.S.C.; Basu, N.B.; Wilson, J.S. The future of hydrology: An evolving science for a changing world. Water Resour. Res. 2010, 46, W05301. [Google Scholar] [CrossRef]
  5. Labat, D.; Goddéris, Y.; Probst, J.L.; Guyot, J.L. Evidence for global runoff increase related to climate warming. Adv. Water Resour. 2004, 27, 631–642. [Google Scholar] [CrossRef]
  6. Pulido-Velazquez, D.; Collados-Lara, A.J.; Pérez-Sánchez, J.; Segura-Méndez, F.J.; Senent-Aparicio, J. Climate change impacts on the streamflow in Spanish basins monitored under near-natural conditions. J. Hydrol. Reg. Stud. 2021, 38, 100937. [Google Scholar] [CrossRef]
  7. Alifujiang, Y.; Abuduwaili, J.; Groll, M.; Issanova, G.; Maihemuti, B. Changes in intra-annual runoff and its response to climate variability and anthropogenic activity in the Lake Issyk-Kul Basin, Kyrgyzstan. Catena 2021, 198, 104974. [Google Scholar] [CrossRef]
  8. Guzha, A.C.; Rufino, M.C.; Okoth, S.; Jacobs, S.; Nóbrega, R.L.B. Impacts of land use and land cover change on surface runoff, discharge and low flows: Evidence from East Africa. J. Hydrol. Reg. Stud. 2018, 15, 49–67. [Google Scholar] [CrossRef]
  9. Zhao, G.; Mu, X.; Wen, Z.; Wang, F.; Gao, P. Soil Erosion, Conservation, and Eco-Environment Changes in the Loess Plateau of China. Land Degrad. Dev. 2013, 24, 499–510. [Google Scholar] [CrossRef]
  10. Xue, D.X.; Zhou, J.J.; Zhao, X.; Liu, C.F.; Wei, W.; Yang, X.M.; Li, Q.Q.; Zhao, Y.R. Impacts of climate change and human activities on runoff change in a typical arid watershed, NW China. Ecol. Indic. 2021, 121, 107013. [Google Scholar] [CrossRef]
  11. Wang, W.; Shao, Q.; Yang, T.; Peng, S.; Xing, W.; Sun, F.; Luo, Y. Quantitative assessment of the impact of climate variability and human activities on runoff changes: A case study in four catchments of the Haihe River basin, China. Hydrol. Process. 2013, 27, 1158–1174. [Google Scholar] [CrossRef]
  12. Giuntoli, I.; Renard, B.; Vidal, J.P.; Bard, A. Low flows in France and their relationship to large-scale climate indices. J. Hydrol. 2013, 482, 105–118. [Google Scholar] [CrossRef]
  13. Han, H.; Hou, J.M.; Huang, M.S.; Li, Z.B.; Xu, K.Y.; Zhang, D.W.; Bai, G.G.; Wang, C. Impact of soil and water conservation measures and precipitation on streamflow in the middle and lower reaches of the Hulu River Basin, China. Catena 2020, 195, 104792. [Google Scholar] [CrossRef]
  14. Dong, X.; Yao, Z.; Chen, C. Runoff Variation and Responses to Precipitation in the Source Regions of the Yellow River. Resour. Sci. 2007, 29, 67–73. [Google Scholar]
  15. Berg, P.; Moseley, C.; Haerter, J.O. Strong increase in convective precipitation in response to higher temperatures. Nat. Geosci. 2013, 6, 181–185. [Google Scholar] [CrossRef]
  16. Abbott, B.W.; Bishop, K.; Zarnetske, J.P.; Minaudo, C.; Chapin, F.S.; Krause, S.; Hannah, D.M.; Conner, L.; Ellison, D.; Godsey, S.E.; et al. Human domination of the global water cycle absent from depictions and perceptions. Nat. Geosci. 2019, 12, 533–540. [Google Scholar] [CrossRef]
  17. Yin, Y.H.; Wu, S.H.; Chen, G.; Dai, E.F. Attribution analyses of potential evapotranspiration changes in China since the 1960s. Theor. Appl. Climatol. 2009, 101, 19–28. [Google Scholar] [CrossRef]
  18. Mengistu, D.; Bewket, W.; Dosio, A.; Panitz, H.-J. Climate change impacts on water resources in the Upper Blue Nile (Abay) River Basin, Ethiopia. J. Hydrol. 2021, 592, 125614. [Google Scholar] [CrossRef]
  19. Jahfer, S.; Vinayachandran, P.N.; Nanjundiah, R.S. Long-term impact of Amazon river runoff on northern hemispheric climate. Sci. Rep. 2017, 7, 10989. [Google Scholar] [CrossRef]
  20. Zhang, D.; Hong, H.Y.; Zhang, Q.; Li, X.H. Attribution of the changes in annual streamflow in the Yangtze River Basin over the past 146 years. Theor. Appl. Climatol. 2014, 119, 323–332. [Google Scholar] [CrossRef]
  21. Liu, X.W.; Xu, Z.X.; Liu, W.F.; Liu, L. Responses of hydrological processes to climate change in the Yarlung Zangbo River basin. Hydrol. Sci. J. 2019, 64, 2057–2067. [Google Scholar] [CrossRef]
  22. Liang, W.; Bai, D.; Wang, F.; Fu, B.; Yan, J.; 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]
  23. Potter, N.J.; Chiew, F.H.S. An investigation into changes in climate characteristics causing the recent very low runoff in the southern Murray-Darling Basin using rainfall-runoff models. Water Resour. Res. 2011, 47, W00G10. [Google Scholar] [CrossRef]
  24. Huang, S.Z.; Liu, D.F.; Huang, Q.; Chen, Y.T. Contributions of climate variability and human activities to the variation of runoff in the Wei River Basin, China. Hydrol. Sci. J. 2016, 61, 1026–1039. [Google Scholar] [CrossRef]
  25. Li, C.Y.; Fang, H.Y. Assessment of climate change impacts on the streamflow for the Mun River in the Mekong Basin, Southeast Asia: Using SWAT model. Catena 2021, 201, 105199. [Google Scholar] [CrossRef]
  26. Wang, G.S.; Xia, J.; Chen, J. Quantification of effects of climate variations and human activities on runoff by a monthly water balance model: A case study of the Chaobai River basin in northern China. Water Resour. Res. 2009, 45, W00A11. [Google Scholar] [CrossRef]
  27. Song, X.Y.; Sun, W.Y.; Zhang, Y.Q.; Song, S.B.; Li, J.Y.; Gao, Y.J. Using hydrological modelling and data-driven approaches to quantify mining activities impacts on centennial streamflow. J. Hydrol. 2020, 585, 124764. [Google Scholar] [CrossRef]
  28. Cong, Z.; Yang, D.; Gao, B.; Yang, H.; Hu, H. Hydrological trend analysis in the Yellow River basin using a distributed hydrological model. Water Resour. Res. 2009, 45, W00A13. [Google Scholar] [CrossRef]
  29. Liu, Q.; Yang, Z.; Cui, B.; Sun, T. The temporal trends of reference evapotranspiration and its sensitivity to key meteorological variables in the Yellow River Basin, China. Hydrol. Process. 2010, 24, 2171–2181. [Google Scholar] [CrossRef]
  30. Budyko, M.I. Climate and Life; Academic Press: Cambridge, MA, USA, 1974. [Google Scholar]
  31. Gan, G.J.; Liu, Y.B.; Sun, G. Understanding interactions among climate, water, and vegetation with the Budyko framework. Earth-Sci. Rev. 2021, 212, 103451. [Google Scholar] [CrossRef]
  32. Pike, J.G. The estimation of annual run-off from meteorological data in a tropical climate. J. Hydrol. 1964, 2, 116–123. [Google Scholar] [CrossRef]
  33. Turc, L. The water balance of soils. Relation between precipitation, evaporation and flow. Ann. Agron 1954, 5, 491–569. [Google Scholar]
  34. Baw-puh, F. ON the calculation of the evaporation from land surface. Chin. J. Atmos. Sci. 1981, 5, 23–31. [Google Scholar]
  35. Roderick, M.L.; Farquhar, G.D. A simple framework for relating variations in runoff to variations in climatic conditions and catchment properties. Water Resour. Res. 2011, 47, W00G07. [Google Scholar] [CrossRef]
  36. 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, W03410. [Google Scholar] [CrossRef]
  37. 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]
  38. Hu, Y.N.; Duan, W.L.; Chen, Y.N.; Zou, S.; Kayumba, P.M.; Sahu, N. An integrated assessment of runoff dynamics in the Amu Darya River Basin: Confronting climate change and multiple human activities, 1960–2017. J. Hydrol. 2021, 593, 126905. [Google Scholar] [CrossRef]
  39. He, Y.; Qiu, H.; Song, J.; Zhao, Y.; Zhang, L.; Hu, S.; Hu, Y. Quantitative contribution of climate change and human activities to runoff changes in the Bahe River watershed of the Qinling Mountains, China. Sustain. Cities Soc. 2019, 51, 101729. [Google Scholar] [CrossRef]
  40. Gao, P.; Geissen, V.; Ritsema, C.J.; Mu, X.M.; Wang, F. Impact of climate change and anthropogenic activities on stream flow and sediment discharge in the Wei River basin, China. Hydrol. Earth Syst. Sci. 2013, 17, 961–972. [Google Scholar] [CrossRef]
  41. Song, J.X.; Xu, Z.X.; Liu, C.M.; Li, H.E. Ecological and environmental instream flow requirements for the Wei River—The largest tributary of the Yellow River. Hydrol. Process. 2007, 21, 1066–1073. [Google Scholar] [CrossRef]
  42. Zuo, D.; Xu, Z.; Yao, W.; Jin, S.; Xiao, P.; Ran, D. Assessing the effects of changes in land use and climate on runoff and sediment yields from a watershed in the Loess Plateau of China. Sci Total Environ. 2016, 544, 238–250. [Google Scholar] [CrossRef] [PubMed]
  43. Huang, S.Z.; Huang, Q.; Chang, J.X.; Leng, G.Y.; Xing, L. The response of agricultural drought to meteorological drought and the influencing factors: A case study in the Wei River Basin, China. Agric. Water Manag. 2015, 159, 45–54. [Google Scholar] [CrossRef]
  44. Li, C.; Zhang, H.B.; Gong, X.H.; Wei, X.W.; Yang, J.T. Precipitation Trends and Alteration in Wei River Basin: Implication for Water Resources Management in the Transitional Zone between Plain and Loess Plateau, China. Water 2019, 11, 2407. [Google Scholar] [CrossRef]
  45. Zhao, Q.; Wang, L.; Liu, H.; Zhang, Q. Runoff and sediment variation and attribution over 60 years in typical Loess Plateau basins. J. Soils Sediments 2019, 19, 3631–3647. [Google Scholar] [CrossRef]
  46. Gao, H.D.; Wu, Z.; Jia, L.L.; Pang, G.W. Vegetation change and its influence on runoff and sediment in different landform units, Wei River, China. Ecol. Eng. 2019, 141, 105609. [Google Scholar] [CrossRef]
  47. Kendall, M.G. Rank Correlation Measures; Charles Griffin: Oxford, UK, 1975. [Google Scholar]
  48. Mann, H.B. Nonparametric test against trend. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  49. Yue, S.; Wang, C.Y. The Mann-Kendall Test Modified by Effective Sample Size to Detect Trend in Serially Correlated Hydrological Series. Water Resour. Manag. 2004, 18, 201–218. [Google Scholar] [CrossRef]
  50. Pettitt, A.N. A Non-parametric Approach to the Change-point Problem. J. R. Stat. Soc. 1979, 28, 126–135. [Google Scholar] [CrossRef]
  51. Allen, D.G.; Pereira, L.O.; Raes, D.; Smith, M.; Allen, R.; Pereira, L.; Smith, M.A. Crop evapotranspiration: Guidelines for computing crop water requirements. Irrig. Drain 1998, 56, 300. [Google Scholar]
  52. Zhao, G.J.; Tian, P.; Mu, X.M.; Jiao, J.Y.; Wang, F.; Gao, P. Quantifying the impact of climate variability and human activities on streamflow in the middle reaches of the Yellow River basin, China. J. Hydrol. 2014, 519, 387–398. [Google Scholar] [CrossRef]
  53. Yang, H.; Yang, D. Derivation of climate elasticity of runoff to assess the effects of climate change on annual runoff. Water Resour. Res. 2011, 47, W07526. [Google Scholar] [CrossRef]
  54. Hu, W.; Si, B. Technical Note: Improved partial wavelet coherency for understanding scale-specific and localized bivariate relationships in geosciences. Hydrol. Earth Syst. Sci. 2021, 25, 321–331. [Google Scholar] [CrossRef]
  55. Zhao, G.J.; Mu, X.M.; Tian, P.; Wang, F.; Gao, P. Climate changes and their impacts on water resources in semiarid regions: A case study of the Wei River basin, China. Hydrol. Process. 2013, 27, 3852–3863. [Google Scholar] [CrossRef]
  56. Zuo, D.; Xu, Z.; Wu, W.; Zhao, J.; Zhao, 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]
  57. Deng, H.J.; Chen, Y.N. Influences of recent climate change and human activities on water storage variations in Central Asia. J. Hydrol. 2017, 544, 46–57. [Google Scholar] [CrossRef]
  58. Pang, J.T.; Duan, J.L.; Zhang, R.; Lu, J.C. Characteristics of Spatiotemporal evolution and climate response of vegetation cover in the Wei River Basin from 2000 to 2019. Res. Soil Water Conserv. 2021, 28, 230–237. [Google Scholar] [CrossRef]
  59. Gao, P.; Mu, X.; Wang, F.; Li, R. Changes in streamflow and sediment discharge and the response to human activities in the middle reaches of the Yellow River. Hydrol. Earth Syst. Sci. 2011, 15, 1–10. [Google Scholar] [CrossRef]
  60. Luo, K. Contribution of ecological conservation programs and climate change to hydrological regime change in the source region of the Yangtze River in China. Reg. Environ. Chang. 2022, 22, 1–24. [Google Scholar] [CrossRef]
  61. Sun, Q.; Miao, C.; Duan, Q.; Wang, Y. Temperature and precipitation changes over the Loess Plateau between 1961 and 2011, based on high-density gauge observations. Glob. Planet. Chang. 2015, 132, 1–10. [Google Scholar] [CrossRef]
  62. Chen, Y.; Wang, K.; Lin, Y.; Shi, W.; Song, Y.; He, X. Balancing green and grain trade. Nat. Geosci. 2015, 8, 739–741. [Google Scholar] [CrossRef]
  63. Wang, D.B.; Hejazi, M. Quantifying the relative contribution of the climate and direct human impacts on mean annual streamflow in the contiguous United States. Water Resour. Res. 2011, 47, W00J12. [Google Scholar] [CrossRef]
  64. Xu, G.C.; Zhang, J.X.; Li, P.; Li, Z.B.; Lu, K.X.; Wang, X.K.; Wang, F.C.; Cheng, Y.T.; Wang, B. Vegetation restoration projects and their influence on runoff and sediment in China. Ecol. Indic. 2018, 95, 233–241. [Google Scholar] [CrossRef]
  65. Qiu, D.; Xu, R.; Wu, C.; Mu, X.; Zhao, G.; Gao, P. Vegetation restoration improves soil hydrological properties by regulating soil physicochemical properties in the Loess Plateau, China. J. Hydrol. 2022, 609, 127730. [Google Scholar] [CrossRef]
  66. Zhou, J.; Fu, B.J.; Gao, G.Y.; Lu, Y.H.; Liu, Y.; Lu, N.; Wang, S. Effects of precipitation and restoration vegetation on soil erosion in a semi-arid environment in the Loess Plateau, China. Catena 2016, 137, 1–11. [Google Scholar] [CrossRef]
  67. Yang, X.; Sun, W.; Mu, X.; Gao, P.; Zhao, G. Run-off affected by climate and anthropogenic changes in a large semi-arid river basin. Hydrol. Process. 2020, 34, 1906–1919. [Google Scholar] [CrossRef]
  68. Liu, X.; Jin, S.; Liu, C.; Gao, Y.; Yang, S.; Gao, Y. Influences of shrubs-herbs-arbor vegetation coverage on the runoff based on the remote sensing data in Loess Plateau. Acta Geogr. Sin. 2014, 69, 1595–1603. [Google Scholar] [CrossRef]
  69. Mo, X.; Liu, S.; Hu, S. Co-evolution of climate-vegetation-hydrology and its mechanisms in the source region of Yellow River. Acta Geogr. Sin. 2022, 77, 1730–1744. [Google Scholar] [CrossRef]
  70. Zhang, M.F.; Liu, N.; Harper, R.; Li, Q.; Liu, K.; Wei, X.H.; Ning, D.Y.; Hou, Y.P.; Liu, S.R. A global review on hydrological responses to forest change across multiple spatial scales: Importance of scale, climate, forest type and hydrological regime. J. Hydrol. 2017, 546, 44–59. [Google Scholar] [CrossRef]
  71. Zheng, H.X.; Zhang, L.; Zhu, R.R.; Liu, C.M.; Sato, Y.; Fukushima, Y. Responses of streamflow to climate and land surface change in the headwaters of the Yellow River Basin. Water Resour. Res. 2009, 45, W00A19. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions, and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions, or products referred to in the content.
Figure 1. Location of the study region and hydrological and meteorological stations in the Wei River Basin.
Figure 1. Location of the study region and hydrological and meteorological stations in the Wei River Basin.
Water 15 02081 g001
Figure 2. Spatial distribution of multi-year average precipitation and intra-annual distribution of precipitation in the Wei River Basin.
Figure 2. Spatial distribution of multi-year average precipitation and intra-annual distribution of precipitation in the Wei River Basin.
Water 15 02081 g002
Figure 3. Observed annual runoff depth (R); precipitation (P); potential evapotranspiration (PET); temperature (T) and parameter n in the upstream of the Wei River (UWR), Jing River (JR), and Beiluo River (BLR). The solid red lines indicate the trend; the Z value was detected using the Mann–Kendall test (**: significant at p < 0.01; NS: no significance at p > 0.05). The 95% confidence interval is denoted by the pink strips.
Figure 3. Observed annual runoff depth (R); precipitation (P); potential evapotranspiration (PET); temperature (T) and parameter n in the upstream of the Wei River (UWR), Jing River (JR), and Beiluo River (BLR). The solid red lines indicate the trend; the Z value was detected using the Mann–Kendall test (**: significant at p < 0.01; NS: no significance at p > 0.05). The 95% confidence interval is denoted by the pink strips.
Water 15 02081 g003
Figure 4. Trend analysis of hydro-meteorological variables in 16 catchments of the Wei River Basin. Trends in annual runoff depth (R); precipitation (P); potential evapotranspiration (PET); temperature (T); and parameter n. △ and ▽ denote increasing and decreasing trends, respectively.
Figure 4. Trend analysis of hydro-meteorological variables in 16 catchments of the Wei River Basin. Trends in annual runoff depth (R); precipitation (P); potential evapotranspiration (PET); temperature (T); and parameter n. △ and ▽ denote increasing and decreasing trends, respectively.
Water 15 02081 g004
Figure 5. The elasticity of runoff to precipitation (a), elasticity of runoff to potential evaporation (b), elasticity of runoff to parameter catchment characteristics of parameters (c) in 16 catchments of the Wei River Basin. The elasticity of runoff means the percent change in runoff coming from the change in precipitation/potential evaporation/catchment characteristics of parameters by 1%.
Figure 5. The elasticity of runoff to precipitation (a), elasticity of runoff to potential evaporation (b), elasticity of runoff to parameter catchment characteristics of parameters (c) in 16 catchments of the Wei River Basin. The elasticity of runoff means the percent change in runoff coming from the change in precipitation/potential evaporation/catchment characteristics of parameters by 1%.
Water 15 02081 g005
Figure 6. Comparison of the observed and simulated runoff change. The solid red lines mean the regression line of the observed and simulated runoff, and the black dotted line represents the regression line with a slope of 1.
Figure 6. Comparison of the observed and simulated runoff change. The solid red lines mean the regression line of the observed and simulated runoff, and the black dotted line represents the regression line with a slope of 1.
Water 15 02081 g006
Figure 7. Runoff change induced by P (a), PET (b) and parameter n (c) in 16 catchments of the Wei River Basin. P: precipitation; PET: potential evapotranspiration; n: catchment characteristic parameter.
Figure 7. Runoff change induced by P (a), PET (b) and parameter n (c) in 16 catchments of the Wei River Basin. P: precipitation; PET: potential evapotranspiration; n: catchment characteristic parameter.
Water 15 02081 g007
Figure 8. Attribution of the annual runoff change in the 16 catchments in the Wei River Basin.
Figure 8. Attribution of the annual runoff change in the 16 catchments in the Wei River Basin.
Water 15 02081 g008
Figure 9. Spatial and temporal variability and net change of land cover in the Wei River Basin from 1995 to 2018. (a) Land cover change (1995–2018). (b) Net change of land cover.
Figure 9. Spatial and temporal variability and net change of land cover in the Wei River Basin from 1995 to 2018. (a) Land cover change (1995–2018). (b) Net change of land cover.
Water 15 02081 g009
Figure 10. The variations in the mean annual NDVI and maximum NDVI in the Wei River Basin from 1982 to 2015.
Figure 10. The variations in the mean annual NDVI and maximum NDVI in the Wei River Basin from 1982 to 2015.
Water 15 02081 g010
Table 1. The change in the hydro-meteorological characteristics and the elasticity coefficients of runoff to precipitation (P), potential evapotranspiration (PET), and parameter.
Table 1. The change in the hydro-meteorological characteristics and the elasticity coefficients of runoff to precipitation (P), potential evapotranspiration (PET), and parameter.
BasinAbrupt YearPeriodChange before and after the Abrupt YearElasticity Coefficient
R (mm)P (mm)PET (mm)nεPεPETεn
WRB1993 *1970–1993–12.3–3.228.70.43.15–2.15–2.12
1994–20183.28–2.28–2.38
UWR1993 **1970–1993–15.311.255.70.52.57–1.57–1.69
1994–20183.07–2.07–2.19
JR1996 **1970–1996–18.3–1.152.50.53.05–2.05–2.37
1997–20183.64–2.64–2.94
BLR1994 **1970–1994–24.0–19.044.80.33.31–2.31–2.41
1995–20183.34–2.34–2.46
Notes: abrupt year: **: significant at p < 0.01; *: significant at p < 0.05.
Table 2. The runoff change in hydro-meteorological characteristics and the contribution of climate change and human activities to runoff changes in the WRB and three sub-basins.
Table 2. The runoff change in hydro-meteorological characteristics and the contribution of climate change and human activities to runoff changes in the WRB and three sub-basins.
RegionRunoff Change Induced by P/PET/n (mm)Contribution Rate to Runoff Change (%)Contribution Rate of Climate Change (%)Contribution Rate of Human Activities (%)
∆RP∆RPET∆RnCPCPETCn
WRB–6.0–5.0–18.620.316.962.837.262.8
UWR–11.2–9.3–21.226.822.350.949.150.9
JR–0.2–3.4–15.31.117.980.920.879.2
BLR–3.4–3.2–9.720.719.559.840.259.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, R.; Gu, C.; Qiu, D.; Wu, C.; Mu, X.; Gao, P. Analysis of Runoff Changes in the Wei River Basin, China: Confronting Climate Change and Human Activities. Water 2023, 15, 2081. https://doi.org/10.3390/w15112081

AMA Style

Xu R, Gu C, Qiu D, Wu C, Mu X, Gao P. Analysis of Runoff Changes in the Wei River Basin, China: Confronting Climate Change and Human Activities. Water. 2023; 15(11):2081. https://doi.org/10.3390/w15112081

Chicago/Turabian Style

Xu, Ruirui, Chaojun Gu, Dexun Qiu, Changxue Wu, Xingmin Mu, and Peng Gao. 2023. "Analysis of Runoff Changes in the Wei River Basin, China: Confronting Climate Change and Human Activities" Water 15, no. 11: 2081. https://doi.org/10.3390/w15112081

APA Style

Xu, R., Gu, C., Qiu, D., Wu, C., Mu, X., & Gao, P. (2023). Analysis of Runoff Changes in the Wei River Basin, China: Confronting Climate Change and Human Activities. Water, 15(11), 2081. https://doi.org/10.3390/w15112081

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