Next Article in Journal
Water Security Assessment of Groundwater Quality in an Anthropized Rural Area from the Atlantic Forest Biome in Brazil
Next Article in Special Issue
Simulation of the Potential Impacts of Projected Climate Change on Streamflow in the Vakhsh River Basin in Central Asia under CMIP5 RCP Scenarios
Previous Article in Journal
Disentangling the Effects of Multiple Stressors on Large Rivers Using Benthic Invertebrates—A Study of Southeastern European Large Rivers with Implications for Management
Previous Article in Special Issue
Climate Change Impacts on Cold Season Runoff in the Headwaters of the Yellow River Considering Frozen Ground Degradation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Streamflow into Beijing and Its Response to Climate Change and Human Activities over the Period 1956–2016

1
China Institute of Water Resources and Hydropower Research, Beijing 100038, China
2
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, Beijing 100038, China
3
Beijing Institute of Water, Beijing 100048, China
4
MWR General Institute of Water Resources and Hydropower Planning and Design, Beijing 100120, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(3), 622; https://doi.org/10.3390/w12030622
Submission received: 15 January 2020 / Revised: 18 February 2020 / Accepted: 21 February 2020 / Published: 25 February 2020

Abstract

:
Streamflow is likely affected by climate change and human activities. In this study, hydro-meteorological data from six rivers upstream of Beijing, namely, the Yongdinghe, Baihe, Heihe, Chaohe, Juhe, and Jumahe Rivers, were analyzed to quantify the spatial and temporal variability of streamflow and their responses to climate change and human activities over the period of 1956–2016. The Mann–Kendall test and moving t-test were used to detect trends and changing points of the annual streamflow. Results showed that the streamflow into Beijing experienced a statistically significant downward trend (p < 0.05), abruptly changing after the early 1980s, owing to climate and human effects. The climate elasticities of the streamflow showed that a 10% decrease in precipitation would result in a 24.5% decrease in total streamflow, whereas a 10% decrease in potential evapotranspiration would induce a 37.7% increase in total streamflow. Human activities accounted for 87% of the reduction in total streamflow, whereas 13% was attributed to climate change. Lastly, recommendations are provided for adaptive management of water resources at different spatial scales.

1. Introduction

Streamflow is influenced by global climate change and human activities. Changes in streamflow, especially drastic increases and decreases, have increased the stress on ecological and socio-economic systems [1,2]. In recent decades, hydrological research has been focused on identifying the driving factors behind changes in streamflow [1,3,4,5,6,7]. The findings can support sustainable water resource use and are of practical significance. Climate change, such as changes in precipitation and potential evapotranspiration, results in streamflow changes [8,9]. Human activities can alter streamflow through land use/land cover change, soil and water conservation, reservoir operation, urbanization, water withdrawal, etc. [10,11,12,13]. However, the impacts of these activities on streamflow are interconnected and complex. Therefore, quantitative investigations of the effects of climate change and human activities on streamflow have been carried out [14,15,16,17,18]. As shown in the studies [15,16,17,18], investigating the characteristics of streamflow and quantifying its driving factors are imperative.
Two main approaches, model-based and statistics-based methods, are widely used for a quantitative evaluation of the effects of climate change and human activities on streamflow [14]. The model-based method uses a hydrological model to quantify the contribution of climate change to the streamflow change. Distributed hydrological models can accurately describe hydrological processes and provide high resolution spatial information that can be analyzed using a geographical information system and remote sensing technology [19]. Modelers can quantify the contribution of climate change to streamflow change by varying the meteorological inputs, given fixed land use/land cover conditions, and they can quantify the contribution of human activities to streamflow change by varying the land use/land cover, given fixed meteorological inputs [11]. This approach may be physically sound; however, collecting and handling data requires effort and can be time-consuming. Furthermore, several studies [2,14,20,21] have shown that a modeling approach can lead to different results between different models because of uncertainty in the model structure and parameter calibration. In addition, the applicability of the hydrological model could be affected in areas of intense human activity because it is difficult to fully describe soil and water conservation, reservoir operation, and water withdrawal status.
The statistics-based method mainly relies on the Budyko hypothesis and the concept of climate elasticity. The original Budyko equations [8] were based on the water and energy balance; thus, the change of catchment water storage is assumed to be zero over a long time period. However, Cavalcante et al. [22] considered water storage in the equation. The Budyko hypothesis provides only a rough approximation of the sensitivity of streamflow to climate change; it is most useful for basins with humidity ratios > 1 and for basins in which the moisture and energy balance are in phase [23]. Ma et al. [24] analyzed eight catchments in the Shiyang River basin in northwest China to investigate changes in annual streamflow over the past 50 years based on the Budyko hypothesis. They estimated that climate change accounted for over 64% of the reduction in the mean annual streamflow, owing to decreased precipitation, while the effect of human activities was limited. Tan and Gan [25] applied the Budyko framework to data from 96 Canadian watersheds and found that climate change increased the mean annual streamflow, while human activities decreased the mean annual streamflow.
The climate elasticity method proposed by Schaake [26] is a widely used approach for investigating the impacts of climate change on streamflow [20,23,27,28,29]. Sankarasubramanian et al. [23] proposed a nonparametric approach, median descriptive statistics, to estimate the climate elasticity directly from the observed data. They also found that the nonparametric estimator had a low bias and was as robust as, or more robust than, watershed model-based approaches for evaluating the sensitivity of long-term streamflow to climate change. Chiew [20] concluded that the nonparametric estimator was simple to use and estimated the elasticity directly from the historical data. Fu et al. [28] extended the nonparametric estimator to two-parameter climate elasticity to assess the climate change effect on annual streamflow. Their application indicated that the modified method could reflect the complicated non-linear relationship among streamflow, precipitation, and temperature and was appropriate to assess the climatic effect on annual streamflow of future climate change scenarios. The nonparametric approach does not require a model assumption or a calibration strategy and can be easily applied in comparative studies to data sets from numerous catchments across large regions [20,23]. This method requires a small amount of data, usually basic meteorological and hydrological data [30]. However, the climate elasticity method still demonstrates uncertainty. For instance, the framework used to quantify the impacts on streamflow are based on the assumption that climate change and human activities are independent [2,30]. In reality, however, the two factors interact, even during the baseline period. The detection of changing points involves subjective judgment, providing an additional source of uncertainty. Recently, many studies have suggested that climate elasticity can also be used to quantify and separate the contributions of climate change and human activities to streamflow, requiring less data than other methods and no parameter adjustments [9,14,30,31,32,33,34].
The attribution analysis could be an important reference for adaptive management of local water resources, especially in areas with decreasing streamflow. Beijing, located in the middle of the Haihe River basin, is the national political, cultural, international communication, and science and technology innovation center of China. The per capita water availability in Beijing is 119 m³, which is only 1/20 of the per capita water availability in China and 1/80 of the average per capita water availability in the world. Beijing is transitioning to a multi-source water supply system to ensure water security, and surface water is an important part of the system. Based on the results of the second water resources investigation and evaluation in Beijing [35], the mean annual amount of river water flowing into Beijing is 2.1 billion m3, which is 1.2 times the amount of water generated by local precipitation in Beijing. In the past 30 years, streamflow into Beijing has decreased, and the water storage and ecosystem health of Beijing have been affected. However, the spatial and temporal changes of streamflow have not been systematically analyzed, and their causes are not well understood. Consequently, the upstream area of Beijing is chosen as our study area. As mentioned above, each analysis method has merits and demerits. Considering the natural condition and available data of the study area, the climate elasticity method is used for the attribution analysis.
The objectives of this study are: (1) to detect statistically significant trends and changing points in the annual streamflow into Beijing, (2) to quantify the impacts of climate change and human activities on the streamflow changes into Beijing, by doing (1) and (2), (3) provide a reference for adaptive management of the water resources at different spatial scales.

2. Study Area and Data

2.1. Study Area

Beijing is located between 39°26′–41°04′ N and 115°20′–117°30′ E in the middle reaches of the Haihe River basin, in the northwest part of the North China Plain. There are seven major rivers within Beijing: the Yongdinghe River (Yongdinghe River catchment), Baihe River (Chaobaihe River catchment), Heihe River (Chaobaihe River catchment), Chaohe River (Chaobaihe River catchment), Juhe River (Juhe River catchment), Beiyunhe River (Beiyunhe River catchment), and Jumahe River (Jumahe River catchment). Each is a tributary of the Haihe River. With the exception of the Beiyunhe River, which originates from the Jundu Mountain in Beijing, the other six rivers flow into the territory of Beijing. Therefore, these six rivers that belong to four catchments were chosen for analysis in this study. The Baihe, Heihe, Chaohe, Juhe, and Jumahe Rivers originate from the Hebei Province, and the Yongdinghe River originates from the Shanxi Province and Inner Mongolia Autonomous Region. The total study area is 58,194 km2. The Yongdinghe, Chaobaihe, Juhe, and Jumahe River catchments account for 69%, 20%, 8%, and 3% of the study area, respectively.
The study area is affected by human activities, and there are severe ecological problems. With economic development, urbanization, and population growth, demands for agricultural, industrial, and domestic water use have significantly increased since China’s land reform and opening-up policy. In addition, land use/land cover change alters the hydrological processes. Major land use types in the study area are cropland (38%), grassland (29%), forest land (28%), and urban and built-up land (3%).

2.2. Data Sources and Processing

Annual streamflow data at the Bahaoqiao, Xiabao, Sandaoying, Gubeikou, Sangyuan, and Dashadi stations were collected from the national hydrologic annals, published by the Hydrological Bureau of the Ministry of Water Resources. The total streamflow into Beijing consists of the streamflow of the Yongdinghe, Baihe, Heihe, Chaohe, Juhe, and Jumahe Rivers. Because the Xiabao (Baihe River), Sandaoying (Heihe River), Gubeikou (Chaohe River), and Sangyuan (Juhe River) stations are not close to the territory of Beijing, the streamflow into Beijing should be modified with a conversion coefficient considering the drainage area. The Beijing Water Authority [35] calculated the total streamflow into Beijing as
Q T = Q B H Q + 1.414 × ( Q X B + Q S D Y + Q G B K ) + 1.307 × Q S Y + Q D S D
where QT is the total streamflow into Beijing; QBHQ is the streamflow at the Bahaoqiao station on the Yongdinghe River; QXB is the streamflow at the Xiabao station on the Baihe River; QSDY is the streamflow at the Sandaoying station on the Heihe River; QGBK is the streamflow at the Gubeikou station on the Chaohe River; QSY is the streamflow at the Sangyuan station on the Juhe River; and QDSD is the streamflow at the Dashadi station on the Jumahe River.
Daily precipitation records and other meteorological data, including daily mean temperature, wind speed, sunshine duration, vapor pressure, and relative humidity, from 86 national meteorological stations within and around the study area were obtained from the National Meteorological Information Center of the China Meteorological Administration (CMA). The Thiessen polygon method was used to calculate the average precipitation and potential evapotranspiration for each catchment. Potential evapotranspiration represents the integrated effect of climate variables [15]. Based on the meteorological data, potential evapotranspiration was estimated using the Penman-Monteith equation recommended by the Food and Agriculture Organization of the United Nations (FAO) [36]. The Penman-Monteith equation calculates potential evapotranspiration with respect to a uniform reference surface, making it convenient to compare the potential evapotranspiration of different locations. Thus, potential evapotranspiration is calculated as
E T p = 0.408 ( R n G ) + γ 900 T m e a n + 273 U 2 ( e s e a ) + γ ( 1 + 0.34 U 2 )
where ETp is the daily potential evapotranspiration (mm day−1); Rn is the net radiation at the crop surface (MJ m−2day−1); G is the soil heat flux (MJ m−2day−1); Tmean is the average daily air temperature at a 2 m height (°C); U2 is the wind speed at a 2 m height (m s-1); es is the saturation vapor pressure (kPa); ea is the actual vapor pressure (kPa); (es- ea) is the saturation vapor pressure deficit (kPa); Δ is the slope of the saturated vapor pressure curve (kPa °C−1); and γ is the psychrometric constant (kPa °C−1).
Data over the period 1956–2016 were obtained. Figure 1 shows the locations of the catchments and hydro-meteorological stations.

3. Methodology

3.1. Detection of Trends and Changing Points

The nonparametric Mann–Kendall test developed by Mann [37] and Kendall [38] is used to detect the significance of the trends and changing points in the hydro-meteorological time series. This widely used method has the following advantages: (1) it is suitable for non-normally distributed data and censored data that are frequently encountered in a hydro-meteorological time series [39], and (2) it is robust to the presence of a small amount of abnormal data [40].
For the given time series X(x1, x2, . . . ,xn), the statistic S is defined as follows.
S = i = 1 n 1 j = i + 1 n s g n ( x j x i )
where
s g n ( x j x i ) { 1 x j > x i 0 x j = x i 1 x j < x i
Mann [37] and Kendall [38] showed that the statistic S is approximately normally distributed, with the expected value and variance given by the following.
E ( S ) = 0
V a r ( S ) = n ( n 1 ) ( 2 n + 5 ) 18
The standardized statistic is shown below.
Z = { ( S 1 ) / V a r ( S ) S > 0   0 S = 0 ( S + 1 ) / V a r ( S ) S < 0
In a two-sided test, the null hypothesis (no trend) is accepted when –Z1-α/2ZZ1-α/2. Here, Z < –Z1-α/2 indicates a decreasing trend, and Z > Z1-α/2 indicates an increasing trend, where α is the significance level [41]. In this study, α is set to 0.05.
To detect the changing points in the hydro-meteorological data Xt = (x1, x2,…, xn), the sequential Mann-Kendall test is used. To apply this method, the cumulative number ni of samples, where xi > xj (1 ≤ ji), should be calculated. The normally distributed statistic Sk can be calculated via the following formula.
S k = i = 1 k n i ( k = 2 ,   3 ,   ,   n )
The expected value and variance of the normally distributed statistic Sk are given by the following.
E ( S k ) = k ( k 1 ) / 4
V a r ( S k ) = k ( k 1 ) ( 2 k + 5 ) / 7
The normalized statistic UFk is estimated as follows:
U F k = S k E ( S k ) V a r ( S k ) ( k = 1 ,   2 ,   ,   n )
where UFk is the forward sequence, and the backward sequence, UBk, is calculated using the same equation but with a reversed series of data. A positive UFk denotes an upward trend, while a negative UFk denotes a downward trend. The null hypothesis (no step changing point) is rejected if any of the points in the forward sequence, UFk, are outside the confidence interval. An intersection point of UFk and UBk located within the confidence interval indicates the beginning of a step changing point [42,43].

3.2. Nonparametric Estimator of Climate Elasticity

The concept of streamflow elasticity was originally introduced by Schaake [26] to evaluate the sensitivity of streamflow to climate change [14]. The climate elasticity of streamflow (ε) is defined by the proportional change in streamflow (Q) divided by the proportional change in a climatic variable, such as precipitation or potential evapotranspiration (X), and is expressed as follows. [2]
ε = d Q / Q d X / X
Sankarasubramanian et al. [23] proposed and defined a nonparametric estimator ε at the mean value of the climatic variable. Similar to the study by Zheng [2], an estimator of climate elasticity (ε) of streamflow is expressed as follows.
Q i / Q ¯ = ε · X i / X ¯
Thus, the elasticity of streamflow ε can be regarded as the linear regression coefficient between ∆Xi/ X ¯ and ∆Qi/ Q ¯ . Then, the elasticity of streamflow can be estimated as [2]:
ε = X ¯ Q ¯ · ( X i X ¯ ) ( Q i Q ¯ ) ( X i X ¯ ) 2 = ρ X , Q · C Q / C X
where ρX,Q is the correlation coefficient of X and Q, and CX and CQ are the coefficients of variation of X and Q, respectively. These can be calculated as shown below.
ρ Q , X = ( Q i Q ¯ ) ( X i X ¯ ) ( Q i Q ¯ ) 2 · ( X i X ¯ ) 2
C Q = ( Q i Q ¯ ) 2 / n Q ¯
C X = ( X i X ¯ ) 2 / n X ¯

3.3. Attribution of Climate Change and Human Activities on the Streamflow

For a given catchment, the streamflow can be modeled as a function of climate change (C) and human activities (H) [30]:
Q = f ( C ,   H )
where Q represents streamflow; C represents the integrated effects of climate variables; and H represents the integrated effects of human activities. Therefore, changes in streamflow can be approximated as [14]
Q = Q C + Q H
where ΔQ is the total change in streamflow; ΔQC is the change in streamflow owing to climate change; and ΔQH is the change in streamflow caused by human activities.
The total streamflow change is determined as
Q = Q o b s , 2 Q o b s , 1
where Qobs,1 and Qobs,2 are the observed annual streamflow before and after the changing point, respectively.
Precipitation (P) and potential evapotranspiration (E0) are the dominant controls on the mean annual water balance [8,44]. Changes in the mean annual precipitation and potential evapotranspiration can lead to changes in the annual streamflow, and the relationship can be calculated as [2]
Q C = ( ε P · P / P o b s , 1 + ε E 0 · E 0 / E 0   o b s , 1 ) · Q o b s , 1
where εP and εE0 are the precipitation and potential evapotranspiration elasticities of the streamflow. Here, ΔP and ΔE0 are the changes in precipitation and potential evapotranspiration and are determined as
ε P = P o b s , 2 P o b s , 1
E 0 = E 0 o b s , 2 E 0 o b s , 1
where Pobs,1, Pobs,2, E0obs,1, and E0obs,2 are the annual precipitation and potential evapotranspiration before and after the changing point.

4. Results

4.1. Trends and Changing Points of the Annual Streamflow

The Mann-Kendall test was applied to detect the trends and changing points of the annual streamflow over the period 1956–2016. The UF value of the total streamflow into Beijing has been consistently lower than −1.96 since 1968. From then onwards, the streamflow began to decline. The value of the Mann-Kendall statistic ( Z = 7.52 ) indicated a statistically significant downward trend of the total streamflow into Beijing during the past six decades, as listed in Table 1. Figure 2a shows that the intersection point of the UF and UB curves was outside the confidence interval, based on the Mann-Kendall test. Thus, the confidence interval exceeded 95%. For more than one intersection point of the UF and UB curves, the detection of the changing point should be combined with other methods. In addition to the Mann-Kendall test, the moving t-test was also employed to further investigate the streamflow time series changing points (Figure 2h) in this study. An abrupt change in the total streamflow into Beijing was detected in 1982 by the integrated application of the Mann-Kendall test and moving t-test [45]. The detected changing points divided the hydro-meteorological time series into two periods. The first period represents the baseline period, when human impacts are assumed to be negligible, and the second period, referred to as the changed period, is associated with the influence of human activities [10,16]. The streamflow data for the baseline period were compared with data from the changed period. The streamflow into Beijing declined by 69% from the baseline period to the changed period.
The UF values for the annual streamflow at the six stations stayed below zero for most of the period. As listed in Table 1, the values of the Mann-Kendall statistic, Z, varied from −3.74 to −8.32; all of these values are lower than the critical value of −1.96. The largest absolute value of Z occurred in the Bahaoqiao station, highlighting that the annual streamflow in the Yongdinghe River has experienced the most significant downward trend during the past six decades. Figure 2b–g and Figure 2i–n show the UF, UB, and t-curves of the Bahaoqiao, Xiabao, Sandaoying, Gubeikou, Sangyuan, and Dashadi stations, respectively. Table 1 lists the detected changing points and change ratios of the mean annual streamflow. The change ratio is defined as the percent change in the mean annual streamflow before and after the changing point. The changing points at each of the stations, with the exception of the Gubeikou and Sangyuan stations, occurred in the early 1980s, consistent with the changing point of the total streamflow into Beijing. The Bahaoqiao station exhibited the largest change in streamflow, while the Yongdinghe River appeared to be the most influential with regard to the total streamflow into Beijing.

4.2. Climate Elasticity of the Streamflow

To assess the impacts of climate change on streamflow, the climate elasticity of the streamflow with respect to annual precipitation and potential evapotranspiration was calculated using Equation (14) and the annual precipitation and potential evapotranspiration data from 1956 to 2016. The elasticity values for each of the six rivers flowing into Beijing and the overall elasticity values are listed in Table 2. Figure 3 shows the streamflow time series at each station over the period 1956–2016.
The mean annual total streamflow into Beijing decreased from 2830 million m3 (48 mm) during the baseline period to 872 million m3 (15 mm) during the changed period. Generally, streamflow positively correlates with precipitation and negatively correlates with potential evapotranspiration. This was observed in the elasticity values, with precipitation elasticity, εP, and potential evapotranspiration elasticity, εE0, calculated as 2.45 and −3.77, respectively, over the long-term period 1956–2016. These elasticity values indicate that a 10% decrease in precipitation is associated with a 24.5% decrease in streamflow, while a 10% decrease in potential evapotranspiration is associated with a 37.7% increase in streamflow. Thus, the total streamflow into Beijing was more sensitive to potential evapotranspiration than to precipitation.
Evaluating each river separately, the precipitation elasticity of the streamflow varied from 0.75 to 3.05, and the potential evapotranspiration elasticity of the streamflow ranged from −0.47 to −4.37. The streamflow in the Yongdinghe, Chaohe, and Jumahe Rivers was more sensitive to potential evapotranspiration than to precipitation, whereas the streamflow of the Baihe, Heihe, and Juhe Rivers was more sensitive to precipitation than to potential evapotranspiration.

4.3. Impacts of Climate Change and Human Activities on the Streamflow

Based on the general framework described in Section 3.2 and the estimated streamflow elasticities (Table 2), the impacts of climate change on the streamflow were quantified. Table 3 lists the estimated contributions of climate change and human activities to the streamflow in the six rivers. Using Equation (21), with εP = 2.45 and given the observed 26 mm decrease in precipitation, the total streamflow depth should decrease by 6.1 mm. At the same time, because εE0 = −3.77, the observed 8 mm decrease in potential evapotranspiration should result in a 1.8 mm increase in the total streamflow depth. The cumulative effect of the changes in precipitation and potential evapotranspiration predicted by Equation (21) is a 4.3 mm decrease of the total streamflow depth or approximately 13% of the observed streamflow reduction. The remaining 28.7 mm decrease in the total streamflow is attributed to human activities based on Equation (19). Thus, human activities were estimated to have accounted for 87% of the total streamflow change.
As listed in Table 3, the combined effects of precipitation and potential evapotranspiration changes caused the streamflow to decrease in all rivers, except the Baihe River. Decreased precipitation was the major climatic factor associated with the decreasing streamflow trends in each of the rivers, except the Baihe River, where the increase in streamflow owing to decreased potential evapotranspiration outweighed the decrease in streamflow owing to decreased precipitation. The overall dominance of the precipitation change effects on the streamflow is a reasonable result because the change of the streamflow was mainly controlled by precipitation variability rather than potential evapotranspiration variability in the water-limited study region [6,46]. While precipitation had more influence than changes in potential evapotranspiration, as shown below, human activities had more influence on the streamflow than climatic factors.
As described in other studies [30,41,47,48], there is no remarkable climate change in the Haihe River basin, especially the Yongdinghe and Chaobaihe River catchments. In this study, the combined impacts of climate change and human activities caused a 56 to 978 million m³ streamflow decrease, with the Yongdinghe River experiencing the largest decline. The reduced precipitation caused decreases in the streamflow depth ranging from 0.1 mm to 68.8 mm, and the changes in potential evapotranspiration caused streamflow depth changes ranging from −3.4 mm to 7.9 mm. Therefore, climate change contributed −4% to 46% of the streamflow drop, and human activities were responsible for 54% to 104% of the streamflow decline. The contribution ratios of human activities to the streamflow change were larger than 80%, except for the Juhe and Jumahe River catchments. Human activities were the dominant cause of the streamflow decline across all catchments.

5. Discussion

5.1. Rationality of the Attribution Analysis

To reduce uncertainty of the climate elasticity method in this study, the Mann–Kendall test and moving t-test were jointly used to detect trends and changing points of the annual streamflow. An attribution analysis was performed. In this study, human activities accounted for the majority of streamflow reduction in the Yongdinghe River catchment, approximately 87%, in agreement with previous studies [41,49]. Table 4 lists the results of previous studies that attributed the streamflow changes to climate change and human activities in the Chaobaihe River catchment. Comparing the results of this study (Table 3) with those of the previous studies listed in Table 4, our estimate of the main driving force for the streamflow reduction in the Chaobaihe River catchment was consistent with previous estimates. However, there was a difference in the contribution of human activities to the streamflow changes. We attribute these differences to methodological choices; different hydro-meteorological stations, different study periods, and different methods of analysis were used in each study. The results are reasonable and could be a reference for the adaptive management of water resources. Nevertheless, previous studies [9,32,33,34,50] performed an attribution analysis using different methods in other regions. To compare and verify our results, other methods should be evaluated to further quantify the impact of climate change and human activities on streamflow.

5.2. Streamflow Change Response to Human Activities

Downward trends of the total streamflow into Beijing have been observed since the early 1980s, despite the non-significant decreasing trends in precipitation and potential evapotranspiration. The streamflow decline implies that streamflow could be affected by human activities in addition to climate change. In this study, we assume that the primary human activities affecting the streamflow into Beijing are water withdrawal and land use/land cover change. The early 1980s marked the beginning of China’s land reform. The reform motivated farmers to manage their lands more productively. This was accompanied by rapidly increasing agricultural water use [14]. Meanwhile, with the concurrent improvement of water infrastructure, the area of effective irrigation expanded, despite a slight decrease in the cropland area. The increase in agricultural irrigation triggered high evapotranspiration, hence decreased streamflow [52]. Figure 4 shows the land use categories for the study area in the summer of 1980 and 2015.
Urban and built-up land increased by 47% in the study area, and their distributions corresponded to the regions with a higher population density and GDP (Figure 5 and Figure 6). Dense population and socio-economic development may imply high industrial and domestic demand for water resources. Increasing demands for agricultural, industrial, and domestic water use have led to over-exploitation of surface and groundwater. Withdrawal from surface water would reduce the streamflow directly. Over-exploitation of groundwater leaves a gap to be filled by surface runoff, thus reducing streamflow [10]. A series of soil and water conservation practices have been carried out in the study area since the 1980s, reducing streamflow by increasing evapotranspiration and improving water conservation [17,52]. Although land use changes, such as an increase of forest and grassland, are not distinct in Figure 4, the degree of coverage has grown. In addition, the low efficiency of water use in the study area is considered to be responsible for the sharp decline in streamflow.
Because of the lack of long-term water consumption data, population and GDP maps (Figure 5 and Figure 6) in the study area were evaluated to determine how they might drive the human impact on streamflow. Sun et al. [50] showed that the population density correlated to a higher impact of the human activities on the streamflow. Figure 5 and Figure 6 show that the GDP and population density are concentrated in the Yongdinghe and Chaobaihe River catchments. The impacts of human activities were the most intensive in these two catchments, consistent with the results in Table 3. Massive water withdrawals could considerably reduce the streamflow of the Yongdinghe and Chaobaihe River catchments. In the Juhe and Jumahe River catchments, the population was small owing to the comparatively inhospitable natural conditions. Therefore, human activities were low. Human activities contributed to only 54% and 57% of the streamflow reduction in the Juhe and Jumahe River catchments, respectively.
Cropland is associated with water withdrawals for agriculture. Yang and Tian [48] found that when farmland exceeded 25% of the catchment area, runoff in the Haihe River basin declined significantly. In this study, the cropland in the Yongdinghe River catchment occupied an area of 18149 km2, accounting for 45% of the catchment area. The mean annual water withdrawal from the surface water for agricultural use was approximately 600 million m3, accounting for 80% of the total surface water withdrawal in the Yongdinghe River catchment. The cropland in the Chaobaihe River catchment occupied an area of 2724 km2, accounting for 25% of the catchment area. The mean annual water withdrawal from the surface water for agricultural use was approximately 33 million m3, accounting for 90% of the total surface water withdrawal in the Chaobaihe River catchment. The Cetian Reservoir, with a water storage capacity of 580 million m3 located in the Yongdinghe River, is the only large reservoir in the study area upstream of Beijing. Reservoir storage impounds water and allows more evaporation to occur, causing the streamflow to decrease [52]. The Cetian Reservoir regulates the streamflow of the Yongdinghe River to meet a range of social, economic, and environmental demands, mainly agricultural irrigation, flood control, and ecological water demand.

5.3. Adaptive Management of Water Resources

A quantitative analysis using the climate elasticity method can enhance our understanding of the streamflow changes into Beijing and the streamflow response to climate change and human activities over time. The adaptive management of water resources should be considered at different spatial scales to achieve harmony of resources, ecology, and socio-economics. The excess emphasis placed on human water use crowds out ecological water use that might cause ecological system degradation. The ecological water demand can be estimated as the minimum amount of water to prevent the river from running dry and allow aquatic biology to survive [14]. To maintain ecosystem health, we suggest that not only water for living and production, but also ecological water demand should be guaranteed in the future. This could be achieved by a reservoir regulation strategy in the Yongdinghe River catchment. Since the early 1980s, numerous soil and water conservation measures, such as an afforestation and planting grass, have been implemented. The proportion of forest and grassland area exceeds 50% of the total area in the Chaobaihe, Jehu, and Jumahe River catchments. Conservation measures prevented the loss of water and soil erosion; however, these may be responsible for the significant reduction of streamflow. Therefore, soil and water conservation measures should be practiced more reasonably and moderately. In addition, because of low water use efficiency in the study area, enforcement of water conservation is required. Cropland accounted for 40% of the study area; thus, water infrastructure should be increased to improve irrigation water use efficiency. For regions with a high population density and socio-economic development, especially the upstream area of the Yongdinghe River catchment, there are growing demands for industrial and domestic water use. The use of water saving devices and unconventional water resources, such as rainwater resources and recycled water, should be promoted.
The streamflow into Beijing has experienced a statistically significant decline since the 1980s, impacting the city’s water supply. In recent years, Beijing has had to rely on water diversion and overexploitation of groundwater to meet water demand. Enhanced cooperation between the upstream and downstream areas of the catchment is required to achieve integrated watershed management. Among the four catchments, the Yongdinghe and Chaobaihe River catchments provide the majority of Beijing’s surface water. The drainage area of these two catchments accounts for 90% of the total study area. The mean annual streamflow of these two catchments over the period 1956–2016 was 1.2 billion m3, 72% of the total streamflow into Beijing. Compared to the other catchments, the Yongdinghe and Chaobaihe Rivers have more influence on the total streamflow into Beijing, and their limited water resources have a significant impact on sustainable water resource use and environmental protection in Beijing. Therefore, it is important to continue to investigate the impacts of climate change and human activities in these two catchments.

6. Conclusions

With the Mann-Kendall test and t-test, the total streamflow into Beijing, as well as the six rivers, experienced a statistically significant decline over the period 1956–2016. Almost all the changing points occurred in the early 1980s. The mean annual total streamflow into Beijing declined by 69% from the baseline period to the changed period. The precipitation and potential evapotranspiration elasticities of the total streamflow into Beijing were 2.45 and −3.77, respectively, indicating that the response of the total streamflow was more sensitive to potential evapotranspiration than to precipitation. The precipitation elasticity of the streamflow of each of the six rivers varied from 0.75 to 3.05, while the potential evapotranspiration elasticity ranged from −0.47 to −4.37.
Climate change and human activities accounted for 13% and 87% of the reduction in the total streamflow into Beijing, respectively. The catchments were subjected to intensive human activities owing to the social and economic developments since the 1980s. Human activities were the dominant factor affecting the streamflow decline for all catchments. However, there were considerable differences in the contribution rates of the human activities among four catchments owing to various degrees of water withdrawal, land use/land cover change, and soil and water conservation practices. Human activities accounted for over 80% of the streamflow reduction in the Yongdinghe and Chaobaihe River catchments, with a higher density of the GDP and population and a larger area of cropland than those of the other catchments. Beijing is faced with a water resource shortage; therefore, enhancing the cooperation between upstream and downstream management to achieve integrated watershed management is required.

Author Contributions

H.W. and Y.Z. conceived and designed the study; X.M. and H.L. collected and analyzed the data; X.M. wrote the manuscript; G.H. and J.L. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2016YFC0401407.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, C.; Xia, J. Water problems and hydrological research in the Yellow River and the Huai and Hai River basins of China. Hydrol. Process. 2004, 18, 2197–2210. [Google Scholar] [CrossRef]
  2. Zheng, H.; Zhang, L.; Zhu, R.; Liu, C.; 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. [Google Scholar] [CrossRef]
  3. 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]
  4. De Fraiture, C. Integrated water and food analysis at the global and basin level. An application of WATERSIM. Water Resour. Manag. 2007, 21, 185–198. [Google Scholar] [CrossRef]
  5. Hejazi, M.I.; Moglen, G.E. Regression-based approach to low flow prediction in the Maryland Piedmont region under joint climate and land use change. Hydrol. Process. 2007, 21, 1793–1801. [Google Scholar] [CrossRef]
  6. Yang, D.; Shao, W.; Yeh, P.J.F.; Yang, H.; Kanae, S.; Oki, T. Impact of vegetation coverage on regional water balance in the nonhumid regions of China. Water Resour. Res. 2009, 45. [Google Scholar] [CrossRef] [Green Version]
  7. Gai, L.; Baartman, J.E.M.; Mendoza-Carranza, M.; Wang, F.; Ritsema, C.J.; Geissen, V. A framework approach for unravelling the impact of multiple factors influencing flooding. J. Flood Risk Manag. 2018, 11, 111–126. [Google Scholar] [CrossRef]
  8. Budyko, M.I. Climate and Life, English ed.; Academic Press: New York, NY, USA, 1974; Volume 508. [Google Scholar]
  9. Cavalcante, R.B.L.; Pontes, P.R.M.; Souza-Filho, P.W.M.; de Souza, E.B. Opposite effects of climate and land use changes on the annual water balance in the amazon arc of deforestation. Water Resour. Res. 2019, 55, 3092–3106. [Google Scholar] [CrossRef]
  10. Miao, C.; Ni, J.; Borthwick, A.G.; Yang, L. A preliminary estimate of human and natural contributions to the changes in water discharge and sediment load in the Yellow River. Glob. Planet. Chang. 2011, 76, 196–205. [Google Scholar] [CrossRef] [Green Version]
  11. 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]
  12. Choi, W.; Nauth, K.; Choi, J.; Becker, S. Urbanization and rainfall–runoff relationships in the Milwaukee River basin. Prof. Geogr. 2016, 68, 14–25. [Google Scholar] [CrossRef] [Green Version]
  13. Jonoski, A.; Popescu, I.; Zhe, S.; Mu, Y.; He, Y. Analysis of Flood Storage Area Operations in Huai River Using 1D and 2D River Simulation Models Coupled with Global Optimization Algorithms. Geosciences 2019, 9, 509. [Google Scholar] [CrossRef] [Green Version]
  14. Hu, S.; Liu, C.; Zheng, H.; Wang, Z.; Yu, J. Assessing the impacts of climate variability and human activities on streamflow in the water source area of Baiyangdian Lake. J. Geogr. Sci. 2012, 22, 895–905. [Google Scholar] [CrossRef]
  15. Ye, X.; Zhang, Q.; Liu, J.; Li, X.; Xu, C.Y. Distinguishing the relative impacts of climate change and human activities on variation of streamflow in the Poyang Lake catchment, China. J. Hydrol. 2013, 494, 83–95. [Google Scholar] [CrossRef]
  16. Zhao, G.; Tian, P.; Mu, X.; Jiao, J.; 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]
  17. Zhang, Q.; Liu, J.; Singh, V.P.; Gu, X.; Chen, X. Evaluation of impacts of climate change and human activities on streamflow in the Poyang Lake basin, China. Hydrol. Process. 2016, 30, 2562–2576. [Google Scholar] [CrossRef]
  18. Wang, M.; Du, L.; Ke, Y.; Huang, M.; Zhang, J.; Zhao, Y.; Gong, H. Impact of climate variabilities and human activities on surface water extents in reservoirs of Yongding River Basin, China, from 1985 to 2016 based on Landsat observations and time series analysis. Remote Sens. 2019, 11, 560. [Google Scholar] [CrossRef] [Green Version]
  19. Pontes, P.R.; Cavalcante, R.B.; Sahoo, P.K.; da Silva Júnior, R.O.; da Silva, M.S.; Dall’Agnol, R.; Siqueira, J.O. The role of protected and deforested areas in the hydrological processes of Itacaiúnas River Basin, eastern Amazonia. J. Environ. Manag. 2019, 235, 489–499. [Google Scholar] [CrossRef]
  20. Chiew, F.H.S. Estimation of rainfall elasticity of streamflow in Australia. Hydrol. Sci. J. 2006, 51, 613–625. [Google Scholar] [CrossRef]
  21. Meng, D.; Mo, X. Assessing the effect of climate change on mean annual runoff in the Songhua River basin, China. Hydrol. Process. 2011, 26, 1050–1061. [Google Scholar] [CrossRef]
  22. Cavalcante, R.B.; Pontes, P.R.; Tedeschi, R.G.; Costa, C.P.; Ferreira, D.B.; Souza-Filho, P.W.; de Souza, E.B. Terrestrial water storage and Pacific SST affect the monthly water balance of Itacaiúnas river basin (eastern Amazonia). Int. J. Climatol. 2019. [Google Scholar] [CrossRef]
  23. Sankarasubramanian, A.; Vogel, R.M.; Limbrunner, J.F. Climate elasticity of streamflow in the United States. Water Resour. Res. 2001, 37, 1771–1781. [Google Scholar] [CrossRef]
  24. Ma, Z.; Kang, S.; Zhang, L.; Tong, L.; Su, X. Analysis of impacts of climate variability and human activity on streamflow for a river basin in arid region of northwest China. J. Hydrol. 2008, 352, 239–249. [Google Scholar] [CrossRef]
  25. 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] [PubMed]
  26. Schaake, J.C. From climate to flow. In Climate Change and US Water Resources; Waggoner, P.E., Ed.; Wiley: New York, NY, USA, 1990; pp. 177–206. [Google Scholar]
  27. Vogel, R.M.; Wilson, I.; Daly, C. Regional regression models of annual streamflow for the United States. J. Irrig. Drain. Eng. 1999, 125, 148–157. [Google Scholar] [CrossRef]
  28. Fu, G.; Charles, S.P.; Chiew, F.H.S. A two-parameter climate elasticity of streamflow index to assess climate change effects on annual streamflow. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  29. 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]
  30. 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]
  31. Ma, H.; Yang, D.; Tan, S.K.; Gao, B.; Hu, Q. Impact of climate variability and human activity on streamflow decrease in the Miyun Reservoir catchment. J. Hydrol. 2010, 389, 317–324. [Google Scholar] [CrossRef]
  32. Li, H.; Zhang, Y.; Vaze, J.; Wang, B. Separating effects of vegetation change and climate variability using hydrological modelling and sensitivity-based approaches. J. Hydrol. 2011, 420–421, 403–418. [Google Scholar] [CrossRef]
  33. Li, Y.; He, D.; Li, X.; Zhang, Y.; Yang, L. Contributions of climate variability and human activities to runoff changes in the upper catchment of the Red River Basin, China. Water 2016, 8, 414. [Google Scholar] [CrossRef] [Green Version]
  34. Zhou, Y.; Lai, C.; Wang, Z.; Chen, X.; Zeng, Z.; Chen, J.; Bai, X. Quantitative evaluation of the impact of climate change and human activity on runoff change in the Dongjiang River Basin, China. Water 2018, 10, 571. [Google Scholar] [CrossRef] [Green Version]
  35. Beijing Water Authority (BWA). Report of the Second Water Resources Investigation and Evaluation in Beijing; Internal Report; BWA: Beijing, China, Unpublished work; 2003. [Google Scholar]
  36. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop. Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; Food and Agricultural Organization of the United Nations (FAO): Rome, Italy, 1998; p. D05109. [Google Scholar]
  37. Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  38. Kendall, M.G. Rank Correlation Measures; Charles Griffin and Company, Ltd.: London, UK, 1975; Volume 202, p. 15. [Google Scholar]
  39. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann–Kendall and Spearman’s rho tests for detecting monotonic trends in hydrological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  40. Burn, D.H.; Elnur, M.A.H. Detection of hydrologic trends and variability. J. Hydrol. 2002, 255, 107–122. [Google Scholar] [CrossRef]
  41. Xia, J.; Zeng, S.; Du, H.; Zhan, C. Quantifying the effects of climate change and human activities on runoff in the water source area of Beijing, China. Hydrol. Sci. J. 2014, 59, 1794–1807. [Google Scholar] [CrossRef]
  42. Moraes, J.M.; Pellegrino, G.Q.; Ballester, M.V.; Martinelli, L.A.; Victoria, R.L.; Krusche, A.V. Trends in hydrological parameters of a southern Brazilian watershed and its relation to human induced changes. Water Resour. Manag. 1998, 12, 295–311. [Google Scholar] [CrossRef]
  43. Zhang, Y.; Guan, D.; Jin, C.; Wang, A.; Wu, J.; Yuan, F. Analysis of impacts of climate variability and human activity on streamflow for a river basin in northeast China. J. Hydrol. 2011, 410, 239–247. [Google Scholar] [CrossRef]
  44. Zhang, L.; Dawes, W.R.; Walker, G.R. Response of mean annual evapotranspiration to vegetation changes at catchment scale. Water Resour. Res. 2001, 37, 701–708. [Google Scholar] [CrossRef]
  45. Buishand, T.A. Some methods for testing the homogeneity of rainfall records. J. Hydrol. 1982, 58, 11–27. [Google Scholar] [CrossRef]
  46. Yang, D.; Sun, F.; Liu, Z.; Cong, Z.; Ni, G.; Lei, Z. Analyzing spatial and temporal variability of annual water-energy balance in nonhumid regions of China using the Budyko hypothesis. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  47. Wang, G.; 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. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, Y.; Tian, F. Abrupt change of runoff and its major driving factors in Haihe River Catchment, China. J. Hydrol. 2009, 374, 373–383. [Google Scholar] [CrossRef]
  49. Mo, C.; Ruan, Y.; Mo, G. Studying Response of runoff to climate change and human activities based on elastic coefficient. J. Hydrol. 2018, 38, 41–45. (In Chinese) [Google Scholar]
  50. Sun, Y.; Tian, F.; Yang, L.; Hu, H. Exploring the spatial variability of contributions from climate variation and change in catchment properties to streamflow decrease in a mesoscale basin by three different methods. J. Hydrol. 2014, 508, 170–180. [Google Scholar] [CrossRef]
  51. Zhan, C.; Niu, C.; Song, X.; Xu, C. The impacts of climate variability and human activities on streamflow in Bai River basin, northern China. Hydrol. Res. 2012, 44, 875–885. [Google Scholar] [CrossRef]
  52. Liu, J.; Zhang, Q.; Singh, V.P.; Shi, P. Contribution of multiple climatic variables and human activities to streamflow changes across China. J. Hydrol. 2016, 545, 145–162. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Locations of the catchments and hydro-meteorological stations.
Figure 1. Locations of the catchments and hydro-meteorological stations.
Water 12 00622 g001
Figure 2. Mann-Kendall test statistics (ag) and t-test statistics (hn) for the annual streamflow over the period 1956–2016. The horizontal dashed lines show the critical value (p < 0.05).
Figure 2. Mann-Kendall test statistics (ag) and t-test statistics (hn) for the annual streamflow over the period 1956–2016. The horizontal dashed lines show the critical value (p < 0.05).
Water 12 00622 g002aWater 12 00622 g002b
Figure 3. Mean annual streamflow change of the different rivers (ag) over the period 1956–2016. The dashed lines indicate multi-year averages over the different periods.
Figure 3. Mean annual streamflow change of the different rivers (ag) over the period 1956–2016. The dashed lines indicate multi-year averages over the different periods.
Water 12 00622 g003aWater 12 00622 g003b
Figure 4. Land use categories for the study area in 1980 and 2015.
Figure 4. Land use categories for the study area in 1980 and 2015.
Water 12 00622 g004
Figure 5. Spatial distribution of population in the study area.
Figure 5. Spatial distribution of population in the study area.
Water 12 00622 g005
Figure 6. Spatial distribution of the GDP in the study area.
Figure 6. Spatial distribution of the GDP in the study area.
Water 12 00622 g006
Table 1. Trends and changing points of the annual streamflow.
Table 1. Trends and changing points of the annual streamflow.
StationArea
(km2)
M-K testt-test
Changing Point
Detected
Changing Point
Mean Annual Streamflow
(Million m3)
Change Ratio (%)
ZChanging PointBaseline PeriodChanged Period
Total streamflow −7.521982–1983198219822830872−69%
Bahaoqiao39,339−8.321982–1983198319831185207−83%
Xiabao4015−8.061980–198219831982236102−57%
Sandaoying1600−3.741983–19881985198513878−44%
Gubeikou4701−5.081994–199519941994362175−52%
Sangyuan375−4.291994–19951996/199819958327−68%
Dashadi4446−5.991971–198019791979700208−70%
Table 2. Streamflow depth, precipitation, potential evapotranspiration, and climate elasticity over the different periods.
Table 2. Streamflow depth, precipitation, potential evapotranspiration, and climate elasticity over the different periods.
RiverHydrological StationData PeriodLong-Term Mean Value (mm)Climate Elasticity
Q PE0εPεE0
Total streamflow 1956–1981485018631.66−3.73
1982–2016154758542.14−4.03
1956–2016294868582.45−3.77
Yongdinghe RiverBahaoqiao1956–1982304678901.18−3.65
1983–201654398862.30−4.75
1956–2016164518882.19−3.89
Baihe RiverXiabao1956–1981595168411.12−1.94
1982–2016255157990.24−1.81
1956–2016405158170.75−0.47
Heihe RiverSandaoying1956–1984875647892.49−2.36
1985–2016495437382.37−1.59
1956–2016655527602.54−1.27
Chaohe RiverGubeikou1956–1993775967362.67−2.89
1994–2016375617122.51−1.37
1956–2016545767222.94−3.11
Juhe RiverSangyuan1956–19942227517862.39−1.11
1995–2016726757874.45−8.39
1956–20161687247863.05−1.77
Jumahe RiverDashadi1956–19781576628112.07−3.71
1979–2016476008152.93−4.14
1956–2016896238143.00−4.37
Table 3. Quantitative assessment of the impacts of climate change and human activities on the streamflow.
Table 3. Quantitative assessment of the impacts of climate change and human activities on the streamflow.
RiverΔQ
(mm)
ΔP
(mm)
ΔE
(mm)
ΔQP
(mm)
ΔQE
(mm)
ΔQC
(mm)
ΔQC/ΔQΔQH
(mm)
ΔQH/ΔQ
Total streamflow
Yongdinghe River
−33
−25
−26
−28
−8
−5
−6.1
−3.9
1.8
0.6
−4.3
−3.4
13%
13%
−28.9
−21.5
87%
87%
Baihe River−33−1−42−0.11.41.3−4%−34.7104%
Heihe River−38−22−52−8.47.2−1.23%−36.597%
Chaohe River−40−35−24−13.37.9−5.414%−34.586%
Juhe River−150−762−68.8−0.8−69.646%−80.754%
Jumahe River−111−624−44.0−3.4−47.443%−63.257%
Table 4. Previous attribution studies on the streamflow change in the Chaobaihe River catchment.
Table 4. Previous attribution studies on the streamflow change in the Chaobaihe River catchment.
CatchmentHydrological StationStudy PeriodMethodContributionSource
Climate Change (%)Human Activities (%)
BaiheXiabao1956–2005Budyko hypothesis16.583.5Xu, et al. [11]
HeiheSandaoying1956–2000Distributed time-variant gain model3070Xu, et al. [11]
ChaoheDaiying1956–2005Distributed time-variant gain model22.177.9Xu, et al. [11]
Baihe-1961–2001Distributed time-variant gain model29.670.4Wang et al. [47]
Chaohe-1961–2001Distributed time-variant gain model31.468.6Wang et al. [47]
BaiheZhangjiafen1986–1998SIMHYD model37.562.5Zhan et al. [51]
ChaoheDaiying1957–2000Two-parameter model method4654Wang et al. [30]
ChaoheDaiying1957–2000Hydrological sensitivity analysis method3466Wang et al. [30]
ChaoheDaiying1957–2000Climate elasticity method3565Wang et al. [30]

Share and Cite

MDPI and ACS Style

Mu, X.; Wang, H.; Zhao, Y.; Liu, H.; He, G.; Li, J. Streamflow into Beijing and Its Response to Climate Change and Human Activities over the Period 1956–2016. Water 2020, 12, 622. https://doi.org/10.3390/w12030622

AMA Style

Mu X, Wang H, Zhao Y, Liu H, He G, Li J. Streamflow into Beijing and Its Response to Climate Change and Human Activities over the Period 1956–2016. Water. 2020; 12(3):622. https://doi.org/10.3390/w12030622

Chicago/Turabian Style

Mu, Xing, Hao Wang, Yong Zhao, Huan Liu, Guohua He, and Jinming Li. 2020. "Streamflow into Beijing and Its Response to Climate Change and Human Activities over the Period 1956–2016" Water 12, no. 3: 622. https://doi.org/10.3390/w12030622

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