Next Article in Journal
Characterization of Regional Groundwater System Based on Aquifer Response to Recharge–Discharge Phenomenon and Hierarchical Clustering Analysis
Previous Article in Journal
Variation Characteristics of Summer Water Vapor Budget and Its Relationship with the Precipitation over the Sichuan Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contribution Analysis of the Streamflow Changes in Selected Catchments on the Loess Plateau, China, Using Multiple Budyko-Based Approaches

1
Department of Infrastructure Engineering, Faculty of Engineering and Information Technology, University of Melbourne, Parkville, VIC 3010, Australia
2
Department of Civil and Environmental Engineering, University of Tennessee, Knoxville, TN 37996, USA
3
Guangdong Open Laboratory of Geospatial Information Technology and Application, Guangzhou Institute of Geography, Guangzhou 510070, China
4
College of Geographical Science, Inner Mongolia Normal University, Hohhot 010022, China
*
Authors to whom correspondence should be addressed.
Water 2021, 13(18), 2534; https://doi.org/10.3390/w13182534
Submission received: 28 June 2021 / Revised: 5 September 2021 / Accepted: 11 September 2021 / Published: 15 September 2021
(This article belongs to the Section Hydrology)

Abstract

:
A better understanding of how streamflow interacts with climate change and human activities would contribute to the efficiency and effectiveness of water resources management. Specifically, quantifying the climate and human contributions has widely been used when attributing streamflow changes. However, only a few previous studies compared the results derived by different methods that are currently available, and even fewer studies have ever had a close look at the uncertainties induced by various estimations of evapotranspiration. This research first examined the streamflow changes for 12 catchments on the Loess Plateau in China during the period of 1961–2018 with Mann–Kendall test and relevant statistical measurements. Then, 8 Budyko-based climate elasticity methods, each with 13 estimations of evapotranspiration, were used to quantifying human and climate contributions to streamflow change in the study area (i.e., 104 pairs of values for human and climate contributions for one catchment). The results showed that significant declining trends could be found in 11 of the 12 catchments studied. In terms of contribution rates, human activity has been shown as the major contributor to the streamflow decrease (60–90%) compared to climate change (10–50%). By comparing the contribution results derived by possible combinations of attribution method and evapotranspiration estimation, the variability due to different Budyko-based methods being used seems to be related to geographical location and climate. Although the spatial pattern of variability due to different estimations of evapotranspiration is not obvious, it is necessary to consider the uncertainties induced when launching contribution analysis over specific regions.

1. Introduction

Water resources are closely related to agricultural production and human well-being. Meanwhile, both climatic changes and human interferences can influence the availability of water resources [1]. Climate change has been observed to have impacts on the global water cycle and alter the patterns of regional extremes (e.g., flood and drought) [2]. On the other hand, human activities also affect water dynamics through the water extraction and production activities that lead to changes in land use and land cover [3,4]. Climatic and anthropogenic factors sometimes come into play together, which leads to intensified impacts on natural flows and extreme conditions, such as drought [5]. In this manner, it is necessary to better understand the impacts of either internal or external factors on water resources and the underlying mechanisms to help decision makers manage water resources in a more sustainable way. Streamflow, also known as runoff, is one of the key components in hydrologic cycles. It is important to investigate the influence of environmental changes (due to either climatic or anthropogenic factors) on streamflow to assess how regional water resources would respond to the changes and to generate adaptation strategies for the changes [6]. Previous studies have investigated the changes in volume of runoff and its flow regime, aiming at identifying the characteristics of streamflow changes and the factors driving streamflow changes [7,8,9]. There are numerous approaches to describe the variations and characteristics of streamflow, such as indicators of hydrologic alteration (IHA), range of variability approach (RVA), and flow duration curve (FDC)- or hydrograph-based hydrologic alteration assessments [10,11,12]. Over a long period, both climate change and human activities contribute to streamflow changes [13]. Therefore, it is vital to separate impacts by different factors when analyzing the causal relationship between streamflow changes and influencing factors.
Previous studies have applied several different methods to partition climate and human contributions. Specifically, the partitioning methods that are currently available can be categorized into physically based hydrological models, statistical methods, and climate elasticity methods based on the Budyko hypothesis [13,14]. Li et all (2009) claimed a predominant impact of land use changes on hydrology within a catchment located at LP by applying SWAT model in their research [15]. Alike models are often affected by the lack of high-resolution datasets and the uncertainties raised by parameterization and model structures [16,17]. By using linear regression methods, Xu (2011) reported a higher contribution of climate change to the variability of streamflow [18], while other studies showed a relatively lower contribution of climate change, in which different study periods and research methods were applied [15,19]. Recently, Budyko-based climate elasticity methods have been increasingly used in the assessment of streamflow changes, because they are physically centered and only involve simple calculations [20,21].
In terms of the Budyko framework, previous studies have proposed several Budyko-based hypotheses that define the annual water balance (i.e., the ratio of actual evapotranspiration to precipitation) as a function of the ratio of potential evapotranspiration to precipitation (i.e., dryness index, φ) [22]. The climate elasticity methods based on those hypotheses can then be used to quantitatively separate the impacts of climate change and human activities on hydrology. However, few previous studies have investigated streamflow changes by quantifying human and climate contributions with multiple Budyko-based methods and assessing the uncertainties involved in the results of quantitative contribution [13]. One of the input parameters for the Budyko-based climate elasticity methods is potential evapotranspiration (E0) the value of which, unlike other parameters, can be obtained from gauging stations; the E0 being used here is often estimated using radiation-based methods, such as the Penman–Monteith (P–M for short) method recommended by Food and Agriculture Organization of the United Nations (FAO). Different estimations of E0 can lead to uncertainties in results but there is a lack of studies focusing on the comparison contribution analyses using different Budyko-based methods and varying inputs of E0.
The Loess Plateau (LP) in China, due to its long history of agricultural practice, low level of vegetative coverage, and the regional climate characterized by intense rainfall in summer, has high soil erosion rates and heavy sediment loads, which, in turn, make the ecosystem here extremely vulnerable to climate change [23]. Numerous studies have investigated the historical climate change using observational datasets [24,25]. For LP, Sun et all (2015) examined the mean and extreme values of temperature and precipitation from 1961 to 2011 and found that there was a general warming trend across the plateau while the magnitude varied spatially; on the other hand, more than one third of the plateau had witnessed a decrease in precipitation, which can be attributed to decreases in the frequency of the event or the intensity of the event [26]. Given the background of global warming, the temperature of LP increased at an annual rate of 0.02 °C since the 1950s, and along with the temporal inequality in precipitation, it contributed to the drying trends identified in almost half the plateau [27,28]. From the 1950s to the 1980s, efforts have been taken to increase the resilience of the LP ecosystem, including a series of water–soil conservation practices (e.g., afforestation and grassing to alter land surface conditions and micro-topography) and hydrological engineering events, such as dam constructions [29,30]. Then, in the 1990s, after experiencing severe drought and subsequent floods, two major conservation actions were launched, aiming at controlling soil erosion and sediment transport, known as the Natural Forest Conservation Program (NFCP) and the Grain-for-Green Program (GFGP) [31]. Those nationwide actions promoted the existing conservation measures, such as building check-dams, converting farmland to forest or grassland, and making investments on ecosystem services, termed as the Ecological restoration (ER) measures. Those ER activities have been proven to be effective in vegetation recovery and reducing sediment [32,33]. The rainfall interception and actual evapotranspiration in this region increased consequently. Meanwhile, an obvious decreasing trend in the streamflow within LP has been noticed, which affects the ecosystem in a negative way [34,35,36]. Thus, relevant studies are typically interested in the water resource changes within LP and the underlying mechanisms [21,37,38].
Considering the ongoing attempts, or those to be raised in the future, to protect the fragile ecosystem of LP, a better understanding of the impacts of human interference and climate change on the streamflow change would be essential. Therefore, this research employed eight Budyko-based climate elasticity methods, each with 13 estimations of E0, to analyze the contributions of climate change and human activities to the streamflow changes in 12 catchments on LP during the period 1961–2018. The objectives are to (1) have more insights into the variabilities of LP streamflow; (2) capture the characteristics of the interactions between climate, human activity, and streamflow change within the study area; (3) discuss the uncertainties in the contribution analysis results that associated with different methods and parameter values.

2. Materials and Methods

2.1. Study Area and Data

The Loess Plateau (approximately 35–41°N, 102–114°E) is located in the middle reaches of the Yellow River basin in North China. Approximately 8.5% of the Chinese population live on this plateau, which covers an area of 632,520 km2 (approximately 6.6% of the entire land area of China). The plateau is dominated by arid and semiarid climates with relatively higher evaporation rates, while it is affected by both East and South Asian monsoons that mostly bring rainfall during the period from June to September. The high-intensity rainstorms, along with the coverage of highly erodible loess, accelerate the soil erosion process, especially for the regions with slopes that exceed 10° [21,39,40]. This study will investigate the hydrology in 12 catchments (Figure 1) across the regions with varying precipitation, land cover types, and land surface gradients [13,23].
For the study area, monthly streamflow data from 1961 to 2018 were acquired from the Yellow River Conservancy Commission (YRCC). Daily precipitation data was provided by the China Meteorological Data Service Centre (CMDC). The potential evapotranspiration was calculated by P-M method and 6 other radiation-base methods, each of which was processed with two different parameterizations (more details are discussed in Section 3.2.) [41,42]. ArcGIS was applied to conduct linear interpolation to spatially average meteorological variables across the whole region of interest. In addition, the interannual variabilities in normalized difference vegetation index (NDVI) with respect to individual catchments were checked to examine the correlation between vegetation cover and streamflow changes. The data were obtained from MODIS with a spatial resolution of 1 km.

2.2. Trend Test and Streamflow Change Measurements

2.2.1. Mann–Kendall Test

The Mann–Kendall (M–K) test has been widely used to detect trends in the time series of hydrological variables [8,43,44] as well as vegetation indices [45,46]. This test is based on a test statistic S, defined as:
S = 1 n 1 j = i + 1 n s g n x i x j ,
where xi and xj are observations, n is the sample size, and sgn is the sign function. When n is larger than 8, S can be approximated as the normal distribution with a mean value of zero [47]. The variance of S can be calculated as:
V a r S = n n 1 2 n + 5 i = 1 n t i i i 1 2 i + 5 18 ,
where ti is the number of ties of extent i. Then, the standard test statistic Z can be calculated:
Z = S 1 V a r S   S > 0 0                 S = 0 S + 1 V a r S   S < 0 ,  
By conducting the M–K test, there will be an increasing trend when Z is greater than zero and a decreasing trend when Z is less than zero.

2.2.2. Indicators of Hydrological Alteration (IHA) and Range of Variability Approach (RVA)

IHA has been widely used to measure streamflow changes and, based on that RVA, was developed to have more insights into this hydrological variability. IHA includes several groups of statistics describing different properties of stream flow changes, such as the magnitude of monthly water conditions and the magnitude and duration of annual extreme water conditions [48]. Here, the first group of statistics was selected considering the temporal resolution of data used in this research. Specific parameters included in the selected group are listed in Table 1.
For the baseline period of a long-term time series of hydrological data, RVA would typically set the range of IHA variations as the interval between the 25th and 75th percentile, or the interval defined by one standard deviation from the mean IHA [10]. Then, the IHAs for the evaluation period will be examined to see if they fall into the pre-defined range. Those IHAs falling outside the range reflect the alteration, which can be measure by the alteration ratio as:
D i = N o , i N e N e × 100 % ,  
N e = N T 2 ,  
where No,i is the number of years when the stream flow is within the range, the subscript i denotes the corresponding IHA statistics, and the Ne is the expected number of years when the stream flow is within the range, the value of which is typically half the total number of years.
Furthermore, the alteration ratios for the high season, low season, and total year can be computed as:
D h i g h = i = 6 9 D i 2 4 ,  
D l o w = i = 1 5 D i 2 + i = 10 12 D i 2 8 ,  
D a n n u a l = i = 1 12 D i 2 12
The absolute value of D between 0–33% is deemed as mild alteration, between 33 – 67% as moderate alteration, and between 67–100% as severe alteration.

2.3. Contribution Analysis on the Impacts of Climate Change and Human Activities on the Streamflow Variability

2.3.1. Budyko Hypothesis

On a long-term average, the Budyko hypothesis ignores the change in water storage of a catchment and uses physical principles to show that the water balance (represented by the ratio E/P, where E is actual evapotranspiration and P is precipitation) is largely controlled by the water–energy balance (represented by the ratio of E0/P). Various equations, known as Budyko-type equations, have been proposed along the years as the solution to the Budyko hypothesis. Wu et all (2017) tabulated a total of 8 mathematical functions that represent the Budyko hypothesis, respectively (Table 2) [13]. Some of these equations are non-parametric, while others include parameters that are connected with catchment characteristics, such as soil type, topography, and vegetation.

2.3.2. Budyko-Based Climate Elasticity Methods

Climate elasticity method was first proposed to identify the relationship between streamflow (Q) changes and their influencing factors ( ε X = Q / Q X / X ) [49]. Based on the Budyko hypothesis, numerous studies have attempted to explain the variability of streamflow by assuming that the changes in E0, P, and human activities are independent throughout the catchment of interest [50,51]. In this manner, the streamflow changes induced by climate change (Qc) can be expressed as:
d Q c ¯ = f P d P + f E 0 d E 0
According to the concept of climate elasticity, Equation (9) can be rewritten as:
d Q c ¯ Q = ε P d P P + ε E 0 d E 0 E 0
Then, the P elasticity of Q, 0, and the E0 elasticity of Q, εE0, can be given as:
ε E 0 = 1 + ϕ f ϕ 1 f ϕ ,   ϕ = E 0 P ,  
ε E 0 + ε P = 1
Following this, the climate-induced changes in streamflow can be calculated as:
Δ Q c ¯ = ε p Δ P P ¯ + ε E 0 Δ E 0 E 0 ¯ Q ¯
It is worth noting that all the precipitation, potential evapotranspiration, and streamflow in the equation above are annually averaged over a long period.

2.3.3. Decomposition Method

As suggested by Wu et all (2017) and Liang et all (2015) the study period (1961–2018) can be split into two sub-periods: the baseline period that was less disturbed by human activities (1961–1970) and the evaluation period (1971–2018), within which human activities were greatly intensified [13,23]. In this research, the difference between the mean annual streamflow in the baseline period and that in the evaluation period is regarded as the total change in the streamflow ( Δ Q ¯ ). Thus, the relationship between the two sources that contribute to the streamflow changes can be expressed as:
Δ Q ¯ = Δ Q c ¯ + Δ Q h ¯ ,  
where Δ Q h ¯ is the streamflow changes caused by human activities. Furthermore, the contributions of climate change (Cc) and human activities (Ch) can be given, in percentage, as:
C c = Δ Q c ¯ Δ Q ¯ × 100 % ,   C h = Δ Q h ¯ Δ Q ¯ × 100 %

2.3.4. E0 Estimation Methods

Potential evapotranspiration (E0) is an important process in both water balance and energy balance. The assessment of a climate’s influence on water resources could be associated with E0 estimation. However, E0 is controlled by many factors, including humidity, wind, and radiation fluxes, which may lead to uncertainties in the results of estimation models. While there are several estimation methods available, most of them were developed focusing on different regions or climate types and were based on specific assumptions [52,53,54,55,56,57,58,59]. Therefore, it is necessary to compare the estimation efficiency in various circumstances to determine the feasibility of individual methods. In this study, we would like to examine the sensitivity of contribution analysis to the choice of E0 estimation method. More details about the methods involved in this study can be found in Table 3.

3. Results and Discussion

3.1. Changes in Hydrometeorological Variables

The M–K test was applied to examine the trends in the time series of three hydrometeorological variables (i.e., Q, P, E0) for the 12 catchments at both interannual and intra-annual scales. Specifically, the intra-annual analysis includes trends for high-flow (June to September) and low-flow seasons (October to May), respectively. As shown in Table 4, no significant trend can be captured for P over the study area; in all the catchments, E0 increased significantly in the low-flow seasons at rates ranging from 2.1 mm/yr (Kuye) to 2.5 mm/yr (Beiluo), while only the E0 in Xinshui shows a significant upward trend in high-flow season. In terms of annual E0, again, not all catchments exhibit significant upward trends. On the other hand, streamflow shows downward trends in both high-flow and low-flow seasons in most catchments, except that no such trend was captured for Dali and Qingjian. Furthermore, Dali turned out to be the only catchment without a significant downward trend found at the interannual scale.
RVA was then applied to the streamflow time series of all 12 catchments for further investigation. The results are also shown in Table 4. In general, the annual streamflow in the catchment Fenhe and Jialu were severely altered; those in the catchment Jinghe and Beiluo were mildly altered; the rest of the catchments have witnessed moderate alternations in the streamflow. At a finer temporal resolution, alterations are more noticeable in high-flow seasons than in low-flow seasons. Specifically, in high-flow seasons, three catchments in the study area have their streamflow severely altered, and the streamflow of the other catchments were also moderately altered. In comparison, more than half the catchments have their streamflow only mildly altered in low-flow seasons. The temporal patterns of streamflow changes in the catchments Fenhe and Jialu are so different from the others’ that most server alterations to the streamflow in those two catchments appeared in the low-flow seasons. The time series of streamflow more explicitly depicted the changes (Appendix A, Figure A1). In general, the streamflow within most catchments studied in this research have fluctuated along the years, but the overall trends are downwards. Compared to low-flow seasons, the streamflow fluctuated more obviously in high-flow seasons, especially for the catchment Gushan.

3.2. Contribution Analysis of Streamflow Changes

In this section, the quantitative contributions of climate change and human activity to streamflow changes in LP are analyzed by comparing results derived by different Budyko-based methods as well as various estimations of evapotranspiration (E0). This research applied eight different Budyko-based methods and 13 different estimation methods for E0 (except for the FAO method, all other methods have corresponding versions that have been tuned by adjusting specific parameters involved). In other words, there are a total of 104 pairs of values for the contribution of climate change and human activity to stream change in each of the 12 catchments being investigated in this research (quantitative attribution of stream flow change for each catchment is not shown here). Figure 2a,b separately illustrate the variability of climate and human contributions at the catchment scale by boxplot.
To control the variables involved in analysis, the results with different E0 estimation methods involved will be averaged (i.e., there is only one pair of human and climate contributions left) when comparing various Budyko-based methods, and similarly, the results derived by different Budyko-based methods will also be averaged when comparing various E0 estimation methods.

3.2.1. Quantitative Attribution Using Various Budyko-Based Methods

As shown by Figure 3, for a specific Budyko-based method, the results were averaged over calculations using different E0. Note that the sum of human and climate contributions equals 100% by design. A red–white–blue color ramp was applied to the results of human contribution, and the same color ramp was reversed for the results of climate contribution. Consequently, in both upper and lower panel, the color that is closer to the red side of the color ramp indicates that human activity contributes more to streamflow change compared to climate change. From Figure 3, we can tell that, for catchments such as Tuwei, Gushan, and Kuye, the results of streamflow change attribution are quite consistent. In fact, except for the Budyko-based methods denoted by Fu, Wang, and Yang, all other methods have shown a predominant role of human activity in the decline of streamflow in those catchments studied. More specifically, most methods revealed that the contribution of human activity to streamflow decline ranges from 71 to 79% with the corresponding contribution of climate ranging from 29 to 21%. By applying ER measures, there has been a greening trend across LP, especially in the last 20 years, and the slope gradient has been altered due to the water-conservation constructions [60,61]. The tradeoffs between conservation objectives and the negative impacts of ER measures on the streamflow remain to be addressed in the future works. Even though some methods (e.g., Fu, Wang, and Yang) are showing higher climate contributions compared to other methods (~50%), human activities still outweigh climate change in most catchment in terms of their impacts on streamflow, the results of which are consistent with previous contribution analyses in the same region [13,23]. The three methods (i.e., FU, Wang, and Yang), with which higher climate contributions would be derived, can be termed as parametric methods (Table 2). They will typically put extra weights on factors such as topography, soil and vegetation conditions, and hydrological features at the catchment scale. However, considering more parameters may sometimes add more uncertainties, let alone geographical factors (e.g., location) and large-scale events, which can also affect the streamflow change. To draw conclusions on the sensitivity of various methods to parameters, elaborate research on a much larger sample size of catchments and corresponding hydrometeorological data is needed, but that is out of the research scope of this study. When comparing the performance of different Budyko-based methods at the catchment scale, the eight methods agreed to the largest degree in Jialu catchment with a standard deviation of 0.81%. On the other hand, the standard deviation of results derived by different methods in the Beiluo catchment reached 11.24%, which also indicates that further investigation is needed to determine the most suitable method to use when assessing human or climate contribution to streamflow change.
The bar plots of standard deviation in Figure 3 reveal that the catchments can be divided into two categories based on the magnitude of the standard deviation. The first six catchments (i.e., Xinshui, Beiluo, Yanhe, Qingjian, Jinghe, and Dali) in Figure 3 have higher standard deviations, on average, compared to the rest of the catchments. By projecting the locations of catchments with the corresponding standard deviations of human/climate contribution (Figure 4), we can tell that those catchments are typically located in southern parts of LP and lower stream regions, compared to the rest of the catchments that are studied in this research. An exception is presented by methods using Hargreaves (HARG) E0 estimation (Figure A4), which corresponds to consistent low standard deviations across the catchments. Moreover, the results calculated by methods using the index-adjusted version of Hargreaves (HARGi) estimation show the highest standard deviations in the six austral catchments. More detailed discussions about the variability due to different E0 estimations will be given in later sections.

3.2.2. Discrepancies Raised by Various E0 Estimations

In Figure 5, for each of the 13 studied E0 estimation methods, results were averaged over calculations using different Budyko-based methods. The same scenario of using color ramp to illustrate human and climate contributions to streamflow changes was applied as Figure 3 (i.e., the cells marked by dark red indicates that human factors predominantly contribute to streamflow change and the human contribution decreases as the cells become bluish). Figure 5 has confirmed the findings in Figure 3, that human activities are the primary driver for the streamflow decline in all the catchments studied in this research, although the magnitude of human contribution might vary due to the method applied and catchment characteristics. For the quantitative attributions using E0 estimation methods, such as MAK, PT, and their coefficient-adjusted versions, the difference identified by the catchment location showed that boreal catchments (e.g., Kuye and Gushan) have higher human contributions to streamflow decline (~81%) than do austral catchments (e.g., Dali) (~74%). By calculating the standard deviation of contribution results at the catchment scale, different E0 estimations tend to cause the highest variability in Xinshui, while it is the catchment Beiluo that has the highest variability of quantitative contribution in the case of comparing various Budyko-based methods (Section 3.2.1). Again, Jialu has witnessed the lowest variability in contribution results; the catchment Jialu may not be sensitive to either different Budyko-based methods or different E0 estimations involved in this calculation.
The quantitative attributions using E0, estimated by methods including FAO, JENSEN, and JENSENi, typically show reduced human contributions to the streamflow decline, while the contributions are still maintained at a level of 69%. The impacts of adjusting coefficients vary with different E0 estimation methods; they either increase or decrease the final contributions calculated. However, for any particular coefficient adjustments, the sign of change (i.e., positive or negative) in the magnitude of contribution is always consistent in all catchments. For example, after adjusting the coefficient for the E0 estimation method, ABTEW, the adjusted version (i.e., ABTEWi) leads to a decrease in human contribution to streamflow decline for all catchments. In other words, coefficient adjustments are independent from catchment characteristics such as geographical location. The bar plots of standard deviation in Figure 4 confirm the relatively lower sensitivity of the catchment Jialu to the method and E0 estimation used in decomposing human and climate contributions. At least for austral catchments (e.g., Xinshui and Beiluo), the variabilities of contributions calculated by methods such as FU, WANG, and YANG are remarkably higher than those calculated by the other methods (Figure A5). For catchments such as Xinshui, the standard deviation of quantitative contributions derived with different E0 can be over 14% when using FU and WANG methods for calculation. At the catchment scale, the level of sensitivity of Budyko-based method to E0 estimation also varies.
E0 has played a key role in calculating the sensitivity and contributions of climate change (which can also be partitioned into the contributions of P and E0 separately) and human activities on streamflow variability with Budyko-based methods [1,13,21,23]. In the real world, the robustness of E0 estimations by a particular method typically depends on the research domain (wet or dry climate) over which the method is applied. There is no doubt that a precise E0 estimation would help to increase the accuracy of contribution analysis. Although identifying the optimal method to calculate E0 for catchments on LP is not the focus of this study, the finding that, when using some estimation methods, the standard deviation is relatively lower for the contributions to the streamflow declines could provide implications for making further hypotheses regarding which method works the best on the selected catchments.

3.2.3. Uncertainties in the Quantitative Attribution

Typical sources of the uncertainty of quantitative contribution analysis include assumptions violating water balance assumption, the non-linearity of the interactions between climate change and human activities, as well as various ways to determine the value of E0. In terms of E0 estimation, FAO has widely been used, while other methods could also be helpful when climate data does not feed the FAO’s appetite. When examining the discrepancy shown by Figure 5, the HARG method formed the celling of human contribution of seven catchments, and after adjusting the coefficient, HARGi formed the floor of human contribution in the same seven catchments. Moreover, in the catchments Tuwei and Gushan, the situation reversed (i.e., HARG estimation led to the lowest human contribution, and HARGi estimation led to the highest estimation). Those results indicate that the different E0 estimations that were involved in this research extended the uncertainty interval of quantitative contribution, which could be attributed to the deviation between the actual value and the estimation of E0. This study indicates that the uncertainty induced by E0 estimation could weight a great share in total uncertainties. It is worth launching further sensitivity analyses upon the uncertainties in subprocesses of E0 estimation, such as pre-process raw data (in situ ore remotely sensed data), equations for calculating E0, and interpolating point/station data, if necessary.
Surface runoff is an important part of the hydrologic cycle. This research attributed streamflow changes to the two general terms, climate change and human activities. However, it is the particular processes, such as changes in CO2 concentration and land use and land cover (LULC), that directly or indirectly influence streamflow via biogeophysical and biogeochemical mechanisms. As discussed in Section 3.2.1, different Budyko-based methods greatly vary over the six catchments located at the austral part of LP, which may indicate that catchments with wet and warm climates are more easily subject to bias in contribution analysis of streamflow changes. However, a similar spatial pattern of contribution variability is not obvious for different E0 estimations. This can also be demonstrated by mapping the standard deviation of quantitative contributions derived by a selected method (Budyko) and various E0 estimation (Figure 6a) and the standard deviation of quantitative contributions derived by a selected E0 estimation (FAO) with various calculation methods (Figure 6b). To address those uncertainties, LULC variables, such as normalized difference vegetation index (NDVI) and leaf area index (LAI), as well as meteorological variables, including temperature and precipitation anomalies, need to be investigated carefully.

4. Conclusions

In this research, changes in streamflow, evapotranspiration (E0), and precipitation were first investigated for 12 catchments within Loess Plateau (LP) during the period of 1961–2018. With analysis methods such as the Mann–Kendall (M–K) test and the Range of Variability Approach (RVA), significant declining trends were found for streamflow in most catchments: low-flow seasons for the catchments studied were characterized by a decrease in E0; and no significant trend was found for the change of precipitation in the study area on LP. Then, eight Budyko-based climate elasticity methods with 13 E0 estimations were used to calculate the quantitative contributions of climate change and human activity to the decline of streamflow happening in those 12 catchments. With the whole period divided into baseline (the first 10 years) and evaluation (the remaining 48 years) periods, the results showed that human activity was the major contributor to the streamflow decline from baseline period to evaluation period for all 12 catchments, while the ratio of climate contribution to human contribution varied among catchments.
The discrepancies in results due to different methods and E0 estimations applied were also examined. In conclusion, the sensitivity of a specific catchment to attribution methods or the key parameters in applied calculations vary and may be affected by geographical locations, landscape characteristics, and local climates. We argue that it would be difficult to find a combination of climate elasticity and E0 estimation that is suitable for every case of contribution analysis, although it might be possible for catchments with similar climate types and catchment characteristics by more detailed sensitivity assessments. However, comparing common features and differences between methods as well as parameter values involved helps to capture the uncertainty in contribution analysis and could have implications for improving methods.

Author Contributions

Conceptualization, Z.Y., C.J.; methodology, Z.Y., C.J. and L.Z.; software, Z.Y., J.S.; validation, Z.Y., J.S. and K.W.; formal analysis, Z.Y. and J.S.; investigation, Z.Y., J.S. and R.H.; resources, Z.Y., L.Z. and C.J.; data curation, Z.Y. and J.S.; writing—original draft preparation, Z.Y.; writing—review and editing, Z.Y., J.S., C.J., K.W. and R.H.; visualization, Z.Y. and J.S.; supervision, C.J. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used and analyzed in this study are publicly available. Streamflow data are available via the YRCC website at http://www.yrcc.gov.cn/ (accessed on 7 July 2019). Meteorological data, such as precipitation and other variables needed for estimating evapotranspiration, are available via the CMDC website at http://www.data.cma.cn/ (accessed on 7 July 2019).

Acknowledgments

The authors wish to thank the anonymous referees for their comments and suggestions. All individuals included in this section have consented to the acknowledgement.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Time series of streamflow change for 12 catchments at intra- and inter-annual scales.
Figure A1. Time series of streamflow change for 12 catchments at intra- and inter-annual scales.
Water 13 02534 g0a1
Figure A2. Time series of precipitation change for 12 catchments at intra- and inter-annual scales.
Figure A2. Time series of precipitation change for 12 catchments at intra- and inter-annual scales.
Water 13 02534 g0a2
Figure A3. Time series of evapotranspiration change for 12 catchments at intra- and inter-annual scales.
Figure A3. Time series of evapotranspiration change for 12 catchments at intra- and inter-annual scales.
Water 13 02534 g0a3
For Figure A1, Figure A2 and Figure A3, the time series for each time period is marked by different colors (blue lines for baseline period and yellow lines for evaluation period) and is fitted by linear regression model separately (subscript b refers to baseline period and subscript e refers to evaluation period).
Figure A4. Standard deviation of quantitative contribution to streamflow change calculated by fixed E0 estimations and varied Budyko-based methods.
Figure A4. Standard deviation of quantitative contribution to streamflow change calculated by fixed E0 estimations and varied Budyko-based methods.
Water 13 02534 g0a4
Figure A5. Standard deviation of quantitative contribution to streamflow change calculated by fixed Budyko-based methods and varied E0 estimations.
Figure A5. Standard deviation of quantitative contribution to streamflow change calculated by fixed Budyko-based methods and varied E0 estimations.
Water 13 02534 g0a5

References

  1. Liang, W.; Bai, D.; Jin, Z.; You, Y.; Li, J.; Yang, Y. 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]
  2. Zhao, Y.; Zou, X.; Gao, J.; Xu, X.; Wang, C.; Tang, D.; Wang, T.; Wu, X. Quantifying the anthropogenic and climatic contributions to changes in water discharge and sediment load into the sea: A case study of the Yangtze River, China. Sci. Total Environ. 2015, 536, 803–812. [Google Scholar] [CrossRef]
  3. Gao, Z.; Zhang, L.; Zhang, X.; Cheng, L.; Potter, N.; Cowan, T.; Cai, W. Long-term streamflow trends in the middle reaches of the Yellow River Basin: Detecting drivers of change. Hydrol. Process. 2015, 30, 1315–1329. [Google Scholar] [CrossRef]
  4. Wang, C.; Wang, S.; Fu, B.; Lü, Y.; Liu, Y.; Wu, X. Integrating vegetation suitability in sustainable revegetation for the Loess Plateau, China. Sci. Total Environ. 2021, 759, 143572. [Google Scholar] [CrossRef]
  5. Lu, X.X.; Li, S.; Kummu, M.; Padawangi, R.; Wang, J. Observed changes in the water flow at Chiang Saen in the lower Mekong: Impacts of Chinese dams? Quat. Int. 2014, 336, 145–157. [Google Scholar] [CrossRef]
  6. Li, Y.; Liu, C.; Yu, W.; Tian, D.; Bai, P. Response of streamflow to environmental changes: A Budyko-type analysis based on 144 river basins over China. Sci. Total Environ. 2019, 664, 824–833. [Google Scholar] [CrossRef] [PubMed]
  7. Tan, X.; Gan, T.Y. Contribution of human and climate change impacts to changes in streamflow of Canada. Sci. Rep. 2015, 5, 17767. [Google Scholar] [CrossRef]
  8. Li, D.; Long, D.; Zhao, J.; Lu, H.; Hong, Y. Observed changes in flow regimes in the Mekong River basin. J. Hydrol. 2017, 551, 217–232. [Google Scholar] [CrossRef]
  9. Zhang, H.; Meng, C.; Wang, Y.; Wang, Y.; Li, M. Comprehensive evaluation of the effects of climate change and land use and land cover change variables on runoff and sediment discharge. Sci. Total Environ. 2020, 702, 134401. [Google Scholar] [CrossRef] [PubMed]
  10. Richter, B.D.; Baumgartner, J.V.; Wigington, R.; Braun, D.P. How much water does a river need? Freshw. Biol. 1997, 37, 231–249. [Google Scholar] [CrossRef] [Green Version]
  11. Vogel, R.M.; Sieber, J.; Archfield, S.A.; Smith, M.P.; Apse, C.D.; Huber-Lee, A. Relations among storage, yield, and instream flow. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  12. Ma, Z.; Wang, Z.; Xia, T.; Gippel, C.J.; Speed, R. Hydrograph-Based Hydrologic Alteration Assessment and Its Application to the Yellow River. J. Environ. Inform. 2014, 23, 1–13. [Google Scholar] [CrossRef]
  13. Wu, J.; Miao, C.; Wang, Y.; Duan, Q.; Zhang, X. Contribution analysis of the long-term changes in seasonal runoff on the Loess Plateau, China, using eight Budyko-based methods. J. Hydrol. 2017, 545, 263–275. [Google Scholar] [CrossRef]
  14. Budyko, M.I.; Miller, D.H. Climate and Life; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  15. Li, Z.; Liu, W.; Zhang, J.; Zheng, F.-L. Impacts of land use change and climate variability on hydrology in an agricultural catchment on the Loess Plateau of China. J. Hydrol. 2009, 377, 35–42. [Google Scholar] [CrossRef]
  16. 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. 2012, 27, 1158–1174. [Google Scholar] [CrossRef]
  17. Long, D.; Longuevergne, L.; Scanlon, B.R. Uncertainty in evapotranspiration from land surface modeling, remote sensing, and GRACE satellites. Water Resour. Res. 2014, 50, 1131–1151. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, J. Variation in annual runoff of the Wudinghe River as influenced by climate change and human activity. Quat. Int. 2011, 244, 230–237. [Google Scholar] [CrossRef]
  19. McVicar, T.R.; Li, L.; Van Niel, T.G.; Zhang, L.; Li, R.; Yang, Q.; Zhang, X.; Mu, X.; Wen, Z.; Liu, W.; et al. Developing a decision support tool for China’s revegetation program: Simulating regional impacts of afforestation on average annual streamflow in the Loess Plateau. For. Ecol. Manag. 2007, 251, 65–81. [Google Scholar] [CrossRef]
  20. Jiang, C.; Xiong, L.; Wang, D.; Liu, P.; Guo, S.; Xu, C. Separating the impacts of climate change and human activities on run-off using the Budyko-type equations with time-varying parameters. J. Hydrol. 2015, 522, 326–338. [Google Scholar] [CrossRef]
  21. Wang, F.; Duan, K.; Fu, S.; Gou, F.; Liang, W.; Yan, J.; Zhang, W. Partitioning climate and human contributions to changes in mean annual streamflow based on the Budyko complementary relationship in the Loess Plateau, China. Sci. Total Environ. 2019, 665, 579–590. [Google Scholar] [CrossRef]
  22. Gan, G.; Liu, Y.; Sun, G. Understanding interactions among climate, water, and vegetation with the Budyko framework. Earth-Sci. Rev. 2021, 212, 103451. [Google Scholar] [CrossRef]
  23. 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]
  24. Zhai, P.; Zhang, X.; Wan, H.; Pan, X. Trends in Total Precipitation and Frequency of Daily Precipitation Extremes over China. J. Clim. 2005, 18, 1096–1108. [Google Scholar] [CrossRef]
  25. Wei, K.; Chen, W. Climatology and trends of high temperature extremes across China in summer. Atmos. Ocean. Sci. Lett. 2009, 2, 153–158. [Google Scholar]
  26. 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]
  27. Sun, W.; Song, X.; Mu, X.; Gao, P.; Wang, F.; Zhao, G. Spatiotemporal vegetation cover variations associated with climate change and ecological restoration in the Loess Plateau. Agric. For. Meteorol. 2015, 209-210, 87–99. [Google Scholar] [CrossRef]
  28. Piao, S.; Ciais, P.; Huang, Y.; Shen, Z.; Peng, S.; Li, J.; Zhou, L.; Liu, H.; Ma, Y.; Ding, Y.; et al. The impacts of climate change on water resources and agriculture in China. Nature 2010, 467, 43–51. [Google Scholar] [CrossRef] [PubMed]
  29. Peng, J.; Chen, S.; Dong, P. Temporal variation of sediment load in the Yellow River basin, China, and its impacts on the lower reaches and the river delta. Catena 2010, 83, 135–147. [Google Scholar] [CrossRef]
  30. Zhao, G.; Mu, X.; Jiao, J.; An, Z.; Klik, A.; Wang, F.; Jiao, F.; Yue, X.; Gao, P.; Sun, W. Evidence and Causes of Spatiotemporal Changes in Runoff and Sediment Yield on the Chinese Loess Plateau. Land Degrad. Dev. 2016, 28, 579–590. [Google Scholar] [CrossRef]
  31. Liu, J.; Li, S.; Ouyang, Z.; Tam, C.; Chen, X. Ecological and socioeconomic effects of China’s policies for ecosystem services. Proc. Natl. Acad. Sci. USA 2008, 105, 9477–9782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Feng, X.; Fu, B.; Piao, S.; Wang, S.; Ciais, P.; Zeng, Z.; Lü, Y.; Zeng, Y.; Li, Y.; Jiang, X.; et al. Revegetation in China’s Loess Plateau is approaching sustainable water resource limits. Nat. Clim. Chang. 2016, 6, 1019–1022. [Google Scholar] [CrossRef]
  33. Wang, S.; Fu, B.; Piao, S.; Lu, Y.; Ciais, P.; Feng, X.; Wang, Y. Reduced sediment transport in the Yellow River due to anthropo-genic changes. Nat. Geosci. 2016, 9, 38–41. [Google Scholar] [CrossRef]
  34. Jiang, C.; Guo, H.; Wei, Y.; Yang, Z.; Wang, X.; Wen, M.; Yang, L.; Zhao, L.; Zhang, H.; Zhou, P. Ecological restoration is not sufficient for reconciling the trade-off between soil retention and water yield: A contrasting study from catchment governance perspective. Sci. Total Environ. 2021, 754, 142139. [Google Scholar] [CrossRef] [PubMed]
  35. Li, R.; Zheng, H.; O’Connor, P.; Xu, H.; Li, Y.; Lu, F.; Robinson, B.E.; Ouyang, Z.; Hai, Y.; Daily, G.C. Time and space catch up with restoration programs that ignore ecosystem service trade-offs. Sci. Adv. 2021, 7, eabf8650. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, H.; Li, X.; Xiao, J.; Ma, M. Evapotranspiration components and water use efficiency from desert to alpine ecosystems in drylands. Agric. For. Meteorol. 2021, 298-299, 108283. [Google Scholar] [CrossRef]
  37. Wang, H.; Sun, F. Variability of annual sediment load and runoff in the Yellow River for the last 100 years (1919–2018). Sci. Total Environ. 2021, 758, 143715. [Google Scholar] [CrossRef] [PubMed]
  38. Yin, Z.; Ottlé, C.; Ciais, P.; Zhou, F.; Wang, X.; Jan, P.; Dumas, P.; Peng, S.; Li, L.; Zhou, X.; et al. Irrigation, damming, and streamflow fluctuations of the Yellow River. Hydrol. Earth Syst. Sci. 2021, 25, 1133–1150. [Google Scholar] [CrossRef]
  39. Fu, B.; Liu, Y.; Lü, Y.; He, C.; Zeng, Y.; Wu, B. Assessing the soil erosion control service of ecosystems change in the Loess Plateau of China. Ecol. Complex. 2011, 8, 284–293. [Google Scholar] [CrossRef]
  40. Lü, Y.; Fu, B.; Feng, X.; Zeng, Y.; Liu, Y.; Chang, R.; Sun, G.; Wu, B. A Policy-Driven Large Scale Ecological Restoration: Quantifying Ecosystem Services Changes in the Loess Plateau of China. PLoS ONE 2012, 7, e31782. [Google Scholar] [CrossRef] [PubMed]
  41. Zhao, L.; (Guangdong Open Laboratory of Geospatial Information Technology and Application, Guangzhou Institute of Geography, Guangzhou, China). Personal communication, 2019.
  42. Wu, H.; Zhao, L.; Xia, J.; Song, J.; Li, T. Applicability of 12 potential evapotranspiration estimation methods under changing environments in different Arid and Humid regions in China. South-to-North Water Trans. Water Sci. Tech. 2021, 19, 01–11. (In Chinese) [Google Scholar]
  43. Li, Z.; Xu, X.; Xu, C.; Liu, M.; Wang, K. Dam construction impacts on multiscale characterization of sediment discharge in two typical karst watersheds of southwest China. J. Hydrol. 2018, 558, 42–54. [Google Scholar] [CrossRef]
  44. Patil, R.; Wei, Y.; Pullar, D.; Shulmeister, J. Evolution of streamflow patterns in Goulburn-Broken catchment during 1884–2018 and its implications for floodplain management. Ecol. Indic. 2020, 113, 106277. [Google Scholar] [CrossRef]
  45. Wu, D.H.; Wu, H.; Zhao, X.; Zhou, T.; Tang, B.J.; Zhao, W.Q.; Jia, K. Evaluation of Spatiotemporal Variations of Global Frac-tional Vegetation Cover Based on GIMMS NDVI Data from 1982 to 2011. Remote Sens. 2014, 6, 4217–4239. [Google Scholar] [CrossRef] [Green Version]
  46. Duan, H.; Xue, X.; Wang, T.; Kang, W.; Liao, J.; Liu, S. Spatial and Temporal Differences in Alpine Meadow, Alpine Steppe and All Vegetation of the Qinghai-Tibetan Plateau and Their Responses to Climate Change. Remote Sens. 2021, 13, 669. [Google Scholar] [CrossRef]
  47. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann-Kendall and Spearman’s rho tests for detecting monotonic trends in hy-drological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  48. Richter, B.D.; Baumgartner, J.V.; Powell, J.; Braun, D.P. A Method for Assessing Hydrologic Alteration within Ecosystems. Conserv. Biol. 1996, 10, 1163–1174. [Google Scholar] [CrossRef] [Green Version]
  49. Schaake, J.C. From climate to flow. In Climate Change and U.S. Water Resources; Waggoner, P.E., Ed.; John Wiley: New York, NY, USA, 1990; Volume 8, pp. 177–206. [Google Scholar]
  50. 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. [Google Scholar] [CrossRef]
  51. Xu, X.; Yang, D.; Yang, H.; Lei, H. 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]
  52. Penman, H.L. Natural evaporation from open water, bare soil and grass. Proc. R. Soc. Lond. A 1948, 193, 120–145. [Google Scholar]
  53. Makkink, G.F. Testing the Penman formula by means of lysimeters. J. Inst. Water Engng. 1957, 11, 277–288. [Google Scholar]
  54. Jensen, M.E.; Haise, H.R. Estimating Evapotranspiration from Solar Radiation. J. Irrig. Drain. Div. 1963, 89, 15–41. [Google Scholar] [CrossRef]
  55. Monteith, J.L. Evaporation and Environment. 19th Symposia of the Society for Experimental Biology; University Press: Cambridge, UK, 1965; Volume 19, pp. 205–234. [Google Scholar]
  56. Priestley, C.H.B.; Taylor, R.J. On the assessment of surface heat flux and evaporation using large scale parameters. Mon. Wealthy Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef]
  57. Hargreaves, G.H. Moisture Availability and Crop Production. Trans. ASAE 1975, 18, 0980–0984. [Google Scholar] [CrossRef]
  58. Doorenbos, J.; Pruitt, W.O. Crop Water Requirements. Irrigation and Drainage Paper No. 24, (Rev.); FAO: Rome, Italy, 1977; 144p. [Google Scholar]
  59. Abtew, W. Evapotranspiration Measurements and Modeling for Three Wetland Systems in South Florida. JAWRA J. Am. Water Resour. Assoc. 1996, 32, 465–473. [Google Scholar] [CrossRef]
  60. Chen, J.; Wang, Y.; Xu, J. The Influences of Soil and Water Conservation on the Water Resource and Sediment in the Loess Plateau; Yellow River Conservancy Press: Zhengzhou, China, 2004. (In Chinese) [Google Scholar]
  61. Gao, J.; Mu, X.; Sun, W. Spatial-temporal changes of vegetation restoration on Loess Plateau from 2000 to 2012. Yellow River 2015, 11, 85–91. (In Chinese) [Google Scholar]
Figure 1. Study Area: the locations of Loess Plateau (red line), catchments (blue line, with catchment names identified), and Yellow River (green line, only the parts within Loess Plateau are shown here) are plotted on a topography basemap, sourced from Esri, USGS, and NOAA.
Figure 1. Study Area: the locations of Loess Plateau (red line), catchments (blue line, with catchment names identified), and Yellow River (green line, only the parts within Loess Plateau are shown here) are plotted on a topography basemap, sourced from Esri, USGS, and NOAA.
Water 13 02534 g001
Figure 2. Boxplots made by all contribution data: (a) the variability of climatic contribution to streamflow change; (b) the variability of human contribution to streamflow change.
Figure 2. Boxplots made by all contribution data: (a) the variability of climatic contribution to streamflow change; (b) the variability of human contribution to streamflow change.
Water 13 02534 g002aWater 13 02534 g002b
Figure 3. Quantitative contribution to streamflow change by climate change (%) and human activities (%): for each type of eight Budyko-based methods, contribution values are averaged over results calculated with 13 E0 estimations; at the catchment scale, the variability (line plot) and standard deviation (bar plot) of contribution values derived by different methods are also shown here.
Figure 3. Quantitative contribution to streamflow change by climate change (%) and human activities (%): for each type of eight Budyko-based methods, contribution values are averaged over results calculated with 13 E0 estimations; at the catchment scale, the variability (line plot) and standard deviation (bar plot) of contribution values derived by different methods are also shown here.
Water 13 02534 g003
Figure 4. Map of the standard deviation (Std) of different Budyko-based methods for quantifying human and climate contributions to streamflow.
Figure 4. Map of the standard deviation (Std) of different Budyko-based methods for quantifying human and climate contributions to streamflow.
Water 13 02534 g004
Figure 5. Quantitative contribution to streamflow change by climate change (%) and human activities (%): for each type of 13 E0 estimations, contribution values are averaged over results calculated with eight Budyko-based methods; at the catchment scale, the variability (line plot) and standard deviation (bar plot) of contribution values derived by different methods are also shown here.
Figure 5. Quantitative contribution to streamflow change by climate change (%) and human activities (%): for each type of 13 E0 estimations, contribution values are averaged over results calculated with eight Budyko-based methods; at the catchment scale, the variability (line plot) and standard deviation (bar plot) of contribution values derived by different methods are also shown here.
Water 13 02534 g005
Figure 6. Map of the standard deviation (Std) of quantitative contributions calculated by (a) various methods with FAO E0 estimation and (b) method Budyko with various E0 estimations.
Figure 6. Map of the standard deviation (Std) of quantitative contributions calculated by (a) various methods with FAO E0 estimation and (b) method Budyko with various E0 estimations.
Water 13 02534 g006
Table 1. IHA statistics (Group 1).
Table 1. IHA statistics (Group 1).
IHA StatisticsHydrological Parameters
Group 1Mean flow in JanuaryMean flow in February
Mean flow in MarchMean flow in April
Mean flow in MayMean flow in June
Mean flow in JulyMean flow in August
Mean flow in SeptemberMean flow in October
Mean flow in NovemberMean flow in December
Table 2. Budyko-type equations.
Table 2. Budyko-type equations.
Function Identifier f Parameters
Budyko ( tanh ( 1 / ) 1 e ) 1 / 2 -
Budyko–Fu (FU) 1 + 1 + m 1 / m m, catchment characteristics
Budyko–Ol’dekop (OLDEKOP) tanh ( 1 / ) -
Bukyko–Pike (PIKE) 1 / 1 + 2 -
Budyko–Schreiber
(SCHREIBER)
1 e -
Budyko–Wang (WANG) 1 + 1 + 4 2 2 2 , vegetation-related catchment characteristics
Budyko–Yang
(YANG)
1 / ( 1 ) n + 1 1 / n n, climate seasonality and catchment characteristics
Budyko–Zhang
(ZHANG)
1 + ω / 1 + ω + 1 ω , plant-available water coefficient
Table 3. E0 estimation methods.
Table 3. E0 estimation methods.
MethodFunction (E0)Parameter
Penman–Monteith (FAO) 0.408 Δ R n G + γ 900 T a + 273 U 2 e s e a Δ + γ 1 + 0.34 U 2 -
Priestley–Taylor (PT) f Δ Δ + γ · R n G λ f 1.26
1.1623 *
Makkink (MAK) g Δ Δ + γ · R s λ g 0.61
0.7314 *
Abtew (ABTEW) K R s λ K 0.53
0.4352 *
Hargreaves (HARG) n T + 17.8 R s λ n 0.0023
0.0105 *
Doorenbos–Pruitt (DOOR) α Δ Δ + γ R s + β β −0.3
−0.4517 *
Jensen–Haise
(JENSEN)
C t T T x R s λ C t 0.025
0.0163 *
In this table, Δ is slope of vapor pressure curve ( kPa / ), Rn is net radiation flux at surface ( kw / m 2 ), G is water heat flux ( kw / m 2 ), γ is psychrometric constant ( kPa / ), Ta is daily average temperature at 2 m height ( ), U2 is wind speed at 2 m height (m/s), e s e a is vapor pressure deficit (kPa), λ is latent heat of water vaporization of water ( MJ / kg ), Rs is solar radiation ( MJ / m 2 / d ), α = 1.066 0.0013 R H + 0.045 U d 0.2 × 10 3 R H · U d 0.135 × 10 4 R H 2 0.0011 U d 2 (RH is mean relative humidity, %; Ud is mean daytime wind speed at 2 m height, m/s), Tx is temperature constant ( T x = 3 ); all methods listed below, except FAO (Penman–Monteith), include specific parameters, and this table lists their original values and adjusted values, denoted by *.
Table 4. Results of the Mann–Kendall test and RVA test.
Table 4. Results of the Mann–Kendall test and RVA test.
PE0QRVA-Q
HighLowTotalHighLowTotalHighLowTotalHighLowTotal
Beiluo----↑ ***-↓ ***↓ ***↓ ***MLL
Dali----↑ ***--↓ ***-MMM
Fenhe----↑ **-↓ ***↓ ***↓ ***MHH
Gushan----↑ ***-↓ ***↓ ***↓ ***HMM
Jialu----↑ ***↑*↓ ***↓ ***↓ ***MHH
Jinghe----↑ ***-↓ ***↓ ***↓ ***MLL
Kuye----↑**↑*↓ ***↓ ***↓ ***MLM
Qingjian----↑ ***--↓ **↓ *MMM
Tuwei----↑ ***↑*↓ ***↓ ***↓ ***MMM
Wuding----↑ **-↓ ***↓ ***↓ ***HMM
Xinshui---↑ *↑ ***↑ ***↓ ***↓ ***↓ ***MMM
Yanhe----↑ ***↑ *↓ **↓ **↓ **HMM
In this table, ↓ denotes that the Z-statistic of Mann–Kendall test is smaller than zero and indicates a downward trend; ↑ denotes that the Z-statistic of Mann–Kendall test is larger than zero and indicates an upward trend; - indicates that there is no significant trend detected; *** indicates that p value is less than 0.001, ** less than 0.01, and * less than 0.05; for the RVA results of streamflow, L denotes the mild alteration, M denotes the moderate alteration, and H denotes the severe alteration.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Z.; Song, J.; Jiang, C.; Wang, K.; Zhao, L.; Hao, R. Contribution Analysis of the Streamflow Changes in Selected Catchments on the Loess Plateau, China, Using Multiple Budyko-Based Approaches. Water 2021, 13, 2534. https://doi.org/10.3390/w13182534

AMA Style

Yang Z, Song J, Jiang C, Wang K, Zhao L, Hao R. Contribution Analysis of the Streamflow Changes in Selected Catchments on the Loess Plateau, China, Using Multiple Budyko-Based Approaches. Water. 2021; 13(18):2534. https://doi.org/10.3390/w13182534

Chicago/Turabian Style

Yang, Zhiyuan, Jian Song, Chong Jiang, Kao Wang, Lingling Zhao, and Runmei Hao. 2021. "Contribution Analysis of the Streamflow Changes in Selected Catchments on the Loess Plateau, China, Using Multiple Budyko-Based Approaches" Water 13, no. 18: 2534. https://doi.org/10.3390/w13182534

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